Physics of the Earth and Planetary Interiors最新文献

{"title":"Water in the lower mantle: constraints from mineral physics and phase equilibrium data","authors":"Li Zhang ,&nbsp;Qun-ke Xia","doi":"10.1016/j.pepi.2025.107475","DOIUrl":"10.1016/j.pepi.2025.107475","url":null,"abstract":"<div><div>Water cycle in the Earth’s deep interior plays a critical role in mantle dynamics. Some subducting materials can reach the core-mantle boundary delivering water along the subduction path. High water content in basalts of ocean islands and large igneous provinces, which were formed by mantle plumes, suggests the presence of water-rich sources in the lower mantle. The seismic observation of low-velocity anomalies in the uppermost lower mantle may imply partial melting associated with water transport out of the transition zone and the released hydrous melts can migrate into the overlying crust layer. Recent experiments indicate that hydrated aluminous silica phases are the best candidate hosts for water transport in slab crust in the lower mantle down to the core mantle boundary owing to their high water storage capacity and exceptional thermal stability. The presence of H₂O significantly lowers the stishovite-to-poststishovite transition pressure to a much shallower depth &lt;1000 km and further stabilizes alumina-rich NiAs-type silica phase (Nt) at the depth &gt;1800 km. In peridotitic compositions, the δ-H solid solution AlOOH–MgSiO₄H₂ is a potential water carrier. The alumina and water partitioning under the lower mantle conditions significantly changes the transition thickness of bridgmanite-to-postperovskite, which allows reasonable interpretation of the seismic scatters and discontinuity in the lower mantle. At the core-mantle boundary, water released upon heating may induce partial melting in subducted oceanic crust, possibly accounting for the ultra-low-velocity zones (ULVZs) detected in high-velocity region. Future studies on multicomponent systems under high pressure and high temperature conditions will be needed to understand how water is exchanged between various mantle reservoirs and how the subducted water and dehydration melting affects the dynamics at the core-mantle boundary.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107475"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145579815","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Source characteristics of the June 21, 2022 (Mw 6.1) Khost earthquake, Afghanistan","authors":"Mohammad Raeesi ,&nbsp;Aram Fathian ,&nbsp;Ali Misaghi ,&nbsp;Cristiano Tolomei","doi":"10.1016/j.pepi.2025.107467","DOIUrl":"10.1016/j.pepi.2025.107467","url":null,"abstract":"<div><div>The June 21, 2022, Khost earthquake in southeastern Afghanistan impacted a rural area housing multiple inadequately constructed towns and villages. Reports indicated over 3600 injuries and a minimum of 1163 fatalities. Historical records did not document any significant seismic events in the area prior to this incident. To examine the rupture process of the earthquake, we conducted multiple analyses using teleseismic and InSAR data, employing separate and joint slip inversion. Given that the 2022 earthquake only slightly ruptured the near-surface crust, the results imply that significantly larger earthquakes might occur again in the same location. The main stress was released as an impulsive rupture that lasted only two to three seconds.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107467"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520960","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The Earth’s magma ocean: Processes and current interpretations from an interdisciplinary perspective","authors":"Charles-Édouard Boukaré ,&nbsp;Laura K. Schaefer ,&nbsp;Hanika Rizo","doi":"10.1016/j.pepi.2025.107463","DOIUrl":"10.1016/j.pepi.2025.107463","url":null,"abstract":"<div><div>To celebrate over 55 years of Physics of the Earth and Planetary Interiors providing a venue for communicating advancements in the chemical and dynamical processes that lead to planetary differentiation, we revisit Earth’s magma ocean in light of the seminal works of Ohtani (1983); Abe and Matsui (1986). Models of Earth’s formation suggest a hot initial state in which much of the planet’s interior was substantially, if not entirely, molten. The global-scale molten silicate mantle is referred to as a magma ocean. Because elements of the periodic table show different affinities for liquid, solid, and gaseous phases, the transition from a molten to a solid mantle provides a key window for early chemical differentiation, with profound implications for Earth’s long-term evolution. The magma ocean hypothesis has been extensively studied in the context of the Moon’s evolution. Major advances in our understanding of the lunar magma ocean have been enabled by experimental access to relevant petrological conditions. Pioneering studies by Ohtani (1983); Abe and Matsui (1986); Solomatov and Stevenson (1993c); Abe (1997) explored magma ocean processes in the context of the Earth, and although high-pressure data for Earth’s mantle were limited at the time, these studies correctly anticipated much of the physics now central to early Earth models. Recent developments, including analyses of short-lived isotopic systems, high pressure experiments using diamond anvil cells, and ab-initio calculations are now providing new constraints on models of Earth’s magma ocean. This review summarizes these recent advances and how they change our understanding of the Earth’s magma ocean evolution. We also discuss the current challenges in developing an interdisciplinary yet coherent picture of the Earth’s earliest evolutionary stages.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107463"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145468950","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Sound velocities of (Al,Fe)-phase H and the seismic visibility of dense oxyhydroxides in Earth’s lower mantle","authors":"Benjamin T. Strozewski ,&nbsp;Ojashvi Rautela ,&nbsp;Johannes Buchen ,&nbsp;Wolfgang Sturhahn ,&nbsp;Takayuki Ishii ,&nbsp;Itaru Ohira ,&nbsp;Olivia Pardo ,&nbsp;Jiyong Zhao ,&nbsp;Thomas S. Toellner ,&nbsp;Jennifer M. Jackson","doi":"10.1016/j.pepi.2025.107466","DOIUrl":"10.1016/j.pepi.2025.107466","url":null,"abstract":"<div><div>Oxyhydroxide phases in the (Al,Fe)OOH–MgSiO₂(OH)<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> system may form within oceanic lithosphere and transport hydrogen in their crystal structures into the lowermost mantle via cold, subducted slabs. In this work, we present new measurements of the seismic wavespeeds of the dense oxyhydroxide (Al,Fe)-phase H (Al_0.84Fe<span><math><msubsup><mrow></mrow><mrow><mn>0</mn><mo>.</mo><mn>07</mn></mrow><mrow><mn>3</mn><mo>+</mo></mrow></msubsup></math></span>Mg_0.02Si_0.06OOH) to 100 GPa constrained by nuclear resonant inelastic X-ray scattering, incorporating previous constraints on the equation of state of this phase. At 300 K and pressure greater than 70 GPa, (Al,Fe)-phase H exhibits high P-wave speeds (<span><math><mi>δ</mi></math></span>lnV<span><math><msub><mrow></mrow><mrow><mi>P</mi></mrow></msub></math></span> <span><math><mo>&gt;</mo></math></span>+12%) and low S-wave speeds (<span><math><mi>δ</mi></math></span>lnV<span><math><msub><mrow></mrow><mrow><mi>S</mi></mrow></msub></math></span> <span><math><mo>&lt;</mo></math></span>−7%) relative to the preliminary reference Earth model (PREM). Experimentally determined sound velocities are incorporated into a model of a hydrous metabasalt including (Al,Fe)-phase H and compared with the seismic wavespeeds of pyrolitic mantle along appropriate adiabats. Hydrous metabasalt may reproduce an anti-correlation of negative shear wave velocity and positive bulk sound velocity at the upper edges of large, low velocity provinces when compared to pyrolitic mantle but has similar wavespeeds to PREM in this region. Hydrous metabasalt with conceivable concentrations of (Al,Fe)-phase H can be distinguished from PREM in V<span><math><msub><mrow></mrow><mrow><mi>S</mi></mrow></msub></math></span> at mid-mantle depths (1100–1700 km) and in V<span><math><msub><mrow></mrow><mrow><mi>P</mi></mrow></msub></math></span> at shallower depths (750–1000 km). Subducted hydrous metabasalt could contribute to scattering of seismic waves across the depth interval of the post-stishovite transition, which may be affected by the formation of (Al,Fe)-phase H.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107466"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520961","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Effective gauss-newton inversion for 3D magnetotelluric data with arbitrarily anisotropic conductivities","authors":"Junjun Zhou ,&nbsp;Bo Han ,&nbsp;Andreas Junge ,&nbsp;Xiangyun Hu ,&nbsp;Tiaojie Xiao ,&nbsp;Zhidan Long ,&nbsp;Ningbo Bai","doi":"10.1016/j.pepi.2025.107465","DOIUrl":"10.1016/j.pepi.2025.107465","url":null,"abstract":"<div><div>A three-dimensional (3D) magnetotelluric (MT) inversion framework for arbitrarily anisotropic media has been developed. This method integrates efficient forward modeling with Gauss-Newton (GN) minimization, utilizing an implicit formulation for sensitivity calculations. The forward problem is discretized using the edge-based finite element (FE) method, where a real-valued linear system, derived from the original complex-valued system, is solved with a direct-iterative hybrid solver to enhance computational efficiency. To address the storage and computational challenges posed by large, dense matrices during inversion, the preconditioned conjugate gradient (PCG) algorithm is employed for model updates. A horizontal anisotropic model is used to evaluate various inversion scenarios, including isotropic, principal axial anisotropic, and more complex anisotropic cases. The influence of inversion parameters on the data is analyzed, and the recoverability of these parameters is thoroughly investigated. Furthermore, the anisotropic inversion algorithm is applied to a complex model that includes both anisotropic and isotropic blocks to evaluate the algorithm's reliability. Finally, the inversion method is tested on the USArray MT dataset, demonstrating its robustness and effectiveness.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107465"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145521038","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Temperature and pressure dependence on slip systems in MgO: Insights from large-strain deformation experiments using the rotational diamond anvil cell","authors":"Keiya Ishimori ,&nbsp;Shintaro Azuma ,&nbsp;Kentaro Uesugi ,&nbsp;Masahiro Yasutake ,&nbsp;Keishi Okazaki ,&nbsp;Bunrin Natsui ,&nbsp;Eranga Gyanath Jayawickrama ,&nbsp;Kenji Ohta","doi":"10.1016/j.pepi.2025.107461","DOIUrl":"10.1016/j.pepi.2025.107461","url":null,"abstract":"<div><div>MgO is the end-member of the primary constituents of the lower mantle, and identifying its crystal preferred orientation (CPO) developments or slip systems is key to understanding seismic observations and the dynamics of the lower mantle. To investigate the temperature dependence on CPO developments in MgO under high-pressure conditions corresponding to the lower mantle, we conducted large-strain deformation experiments using the rotational diamond anvil cell (rDAC) combined with synchrotron X-rays, achieving pressures up to 80 GPa and temperatures up to 973 K. Our results revealed that the CPO developments in MgO under large-strain deformation are temperature-dependent even at relatively low temperatures. The crystal plane parallel to the shear plane changed from the {110} plane to the {100} plane with increasing temperature and pressure. Based on our experimental results, we constructed a temperature-pressure map that shows the CPO variation of MgO. The temperature-pressure map obtained in this study provides essential foundational information for advancing our understanding of rheology in the lower mantle.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107461"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145435476","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Detecting rapid lateral changes of upper mantle discontinuities using azimuth-dependent P-wave receiver functions and multimode surface waves","authors":"Kotaro Tarumi,&nbsp;Kazunori Yoshizawa","doi":"10.1016/j.pepi.2025.107468","DOIUrl":"10.1016/j.pepi.2025.107468","url":null,"abstract":"<div><div>Lateral structural variations in the upper mantle generate azimuthal dependence in receiver functions (RFs) based on incoming directions of body waves. Although these azimuthal variations in RFs have not been considered in earlier studies of RF inversions, they provide a means to image localized changes in upper mantle interfaces. In this study, we incorporate azimuth-dependent RFs into a joint Bayesian inversion with multimode surface waves, applying this approach to major permanent broadband stations in Australia. The resulting models reveal dependence on event directions, and by identifying P-to-S conversion depths, we constructed a localized map of conversion points, reflecting local lateral variations of upper mantle discontinuities beneath each station. At the CTAO station in northeastern Australia, the lithosphere thickens rapidly northwestward, from 70 km to 120-130 km depth, corresponding to the tectonic boundary with the western cratonic region. At stations in western and central Australia, lithospheric thickness also varies laterally within the stable cratons, though these changes are more gradual than those in northeastern Australia. In addition to the lithosphere–asthenosphere boundary (LAB), both mid-lithospheric discontinuities (MLDs) and X-discontinuities (X-Ds) are observed in the local 1-D profiles. The X-Ds, characterized by seismic velocity jumps below the LAB, are found at multiple depths around 170, 220, 260, and 310 km, depending on location, accompanying the weakened radial anisotropy across these depths. The multiple MLDs are also identified in the cratonic regions, showing substantial variations in their seismological properties, including both positive or negative S-velocity jumps, which vary with location and depth. Our approach, incorporating azimuth-dependent RFs, enables the detection of localized changes in the upper mantle discontinuities and associated elastic properties, providing new insights into the complex layering of the upper mantle.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"369 ","pages":"Article 107468"},"PeriodicalIF":1.9,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145468949","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A simple linear regression model for calculating the degree of melting of the upper mantle using the major element composition of basalts","authors":"Kenta Ueki ,&nbsp;Satoru Haraguchi ,&nbsp;Atsushi Nakao ,&nbsp;Hikaru Iwamori","doi":"10.1016/j.pepi.2025.107464","DOIUrl":"10.1016/j.pepi.2025.107464","url":null,"abstract":"<div><div>We present a simple linear equation to calculate the degree of melting of the mantle using the major element composition of basalts. We constructed a model based on compiled results from high-pressure mantle melting experiments. We used a model selection approach to objectively select the optimal equation from many potential models, based on statistical criteria. We found that the degree of mantle melting (<span><math><mi>F</mi></math></span>) can be predicted with a simple equation that uses only the concentrations of three major elements, which are <span><math><msub><mrow><mi>SiO</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span>, FeO<span><math><msup><mrow></mrow><mrow><mo>∗</mo></mrow></msup></math></span> (total iron as FeO), and MgO (wt.%), as follows:</div><div><span><math><mrow><mi>F</mi><mrow><mo>(</mo><mi>wt</mi><mo>.</mo><mtext>%</mtext><mo>)</mo></mrow><mo>=</mo><mn>4</mn><mo>.</mo><mn>020</mn><mo>×</mo><msub><mrow><mi>SiO</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>+</mo><mn>5</mn><mo>.</mo><mn>109</mn><mo>×</mo><msup><mrow><mi>FeO</mi></mrow><mrow><mo>∗</mo></mrow></msup><mo>+</mo><mn>1</mn><mo>.</mo><mn>436</mn><mo>×</mo><mi>MgO</mi><mo>−</mo><mn>244</mn><mo>.</mo><mn>924</mn><mo>.</mo></mrow></math></span></div><div>The model allows us to calculate the degree of melting of the uppermost upper mantle (spinel lherzolite and harzburgite), under both anhydrous and hydrous conditions. The equation yields the equilibrium degree of melting in the case of batch melting, and the weighted-mean degree of melting of accumulated melt in the case of fractional melting. We also describe the petrological and thermodynamic implications of the equation. The degrees of melting of natural basalts are calculated as examples of the application of the equation.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"368 ","pages":"Article 107464"},"PeriodicalIF":1.9,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145362498","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Core field changes from eleven years of Swarm satellite observations","authors":"C.C. Finlay ,&nbsp;C. Kloss ,&nbsp;N. Gillet","doi":"10.1016/j.pepi.2025.107447","DOIUrl":"10.1016/j.pepi.2025.107447","url":null,"abstract":"<div><div>The <em>Swarm</em> satellite trio has provided global vector magnetic field measurements, with high precision and absolute accuracy, for the past eleven years. Based on this consistent, high quality, dataset we describe here how Earth’s main (core-generated) magnetic field has evolved between 2014.0 and 2025.0.</div><div>At the Earth’s surface, we find that the region in the South Atlantic where the field strength is weakest (below 26,000 nT), has expanded by 0.9% of Earth’s surface area and that the minimum intensity has decreased by 336 nT from 22,430 nT to 22,094 nT. In the northern polar region, we find that in Canada the area of strong field (above 57,000 nT) has diminished, decreasing in size by 0.65% of Earth’s surface area and with the maximum field strength decreasing by 801 nT from 58,832 nT to 58,031 nT. In contrast the corresponding strong field region in Siberia has grown in size, increasing in area by 0.42% of Earth’s surface area, with the maximum field intensity increasing by 260 nT from 61,359 nT to 61,619 nT.</div><div>At the core-mantle boundary, reversed flux features under southern Africa have moved westward, converging towards reversed flux features that have moved eastwards under the mid-Atlantic. In the northern polar region a strong flux feature under the Bering strait has moved westwards along the inner-core tangent cylinder. At low latitudes, under Indonesia and the western Pacific, field features have surprisingly moved eastwards. Field accelerations, including oscillations, are found to be most intense at low latitudes.</div><div>The <em>Swarm</em> mission has for the past decade been an essential source of global information on the changes taking place in Earth’s main magnetic field. Due to the long timescales of the underlying core processes, extending the mission lifetime for as long as possible, in particular for the higher satellite <em>Swarm</em> Bravo, is expected to yield further scientific insights. A long mission for <em>Swarm</em> Bravo would be an efficient means of ensuring that the present era of high quality geomagnetic observations from space continues as new missions come online.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"368 ","pages":"Article 107447"},"PeriodicalIF":1.9,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145221140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Remarkable variations in geomagnetic disturbances preceding the 2017 M 7.1 Jiuzhaigou earthquake","authors":"Yiliang Guan ,&nbsp;Xiaona Dong ,&nbsp;Lili Feng ,&nbsp;Manqiu He ,&nbsp;Yingfeng Ji","doi":"10.1016/j.pepi.2025.107438","DOIUrl":"10.1016/j.pepi.2025.107438","url":null,"abstract":"<div><div>With the rapid development of geomagnetic observation networks and the increased availability of geomagnetic data in China over the past decade, multistation geomagnetic data analyses have become widely used in studies of earthquake precursors with magnitudes &gt;6. Building on the single-station polarization algorithm, we propose a multistation daily sliding polarization correlation algorithm for detecting regional geomagnetic anomalies. This method can be applied to analyze the spatiotemporal evolution of geomagnetic field anomalies that are related to seismogenic processes before major earthquakes. In this study, we calculate polarization correlation time series data before and after the 2017 M 7.1 Jiuzhaigou earthquake. The results show that following the appearance of polarization anomalies, the regional correlations increase, then decrease, and subsequently increase again. Earthquakes occurred at the inflection points of these correlation trends. Our analysis suggests that changes in multistation polarization correlations directly reflect regional geomagnetic field anomalies that are caused by tectonic activity, thereby capturing the dynamic evolution during earthquake preparation. This method offers valuable insights into the mechanisms linking geomagnetic polarization anomalies with seismic events.</div></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":"368 ","pages":"Article 107438"},"PeriodicalIF":1.9,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049616","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}