Geophysical Journal International

SummaryEarth’s magnetic field has exhibited erratic polarity reversals over much of its history; however, the processes that cause polarity transitions are still poorly understood. Dipole reversals have been found in many numerical dynamo simulations and often occur close to the transition between dipole-dominated and multipolar dynamo regimes. However, the physical conditions used in reversing simulations are necessarily far from those in Earth’s liquid iron core because of the long runtimes needed to capture polarity transitions and because a systematic exploration of parameter space is needed to find the dipole-multipole transition. Here, we use the theory of distinguished limits in an attempt to simplify the search for the dipole-multipole transition at increasingly realistic physical conditions. We consider three limits that are all built from the requirements of a constant magnetic Reynolds number $\mathit {Rm}$; one limit further attempts to impose balance between Magnetic, Coriolis, and Archimedean forces (a QG-MAC balance) while the other two seek to constrain solutions to an inertia-MAC, or QG-IMAC, balance. The presence of inertia, although not geophysically realistic, allows us to build limits that more closely follow the conditions where simulated reversals have been found to date. Numerical simulations along paths in parameter space defined by these limits show some consistencies with the assumed dynamical balances within the accessible parameter space, but also important discrepancies from predicted behaviour for certain diagnostic quantities, particularly the magnetic field strength and the magnetic/kinetic energy ratio. Furthermore, the paths do not follow the dipole-multipole transition; starting from reversing conditions, simulations move into the dipolar non-reversing regime as they are advanced along the path. By increasing the Rayleigh number, a measure of the buoyancy driving convection, above the values predicted by the distinguished limit, we are able to bound the dipole-multipole transition down to an Ekman number $\mathit {E}\sim 10^{-6}$, comparable to the most extreme conditions reported to date. Our results, therefore, demonstrate that using distinguished limits is an efficient method for seeking the dipole-multipole transition in rapidly rotating dynamos. However, the conditions under which we bound the dipole-multipole transition become increasingly hard to access numerically and also increasingly unrealistic because $\mathit {Rm}$ rises beyond plausible bounds inferred from geophysical observations. Future work combining the theory of distinguished limits with variations in the core buoyancy distribution, as suggested by recent studies, appears a promising approach to accessing the dipole-multiple transition at extreme physical conditions.

SummaryLouis Néel’s theory of magnetism, which describes a rock’s magnetization as being carried by “ideal” uniformly magnetized or single domain (SD) particles, has been a cornerstone of paleomagnetic studies for over seven decades and has enabled paleomagnetists to make plate tectonic reconstructions, date geological and archaeological materials, and recover the history of Earth’s magnetic field. Unfortunately, many geological samples produce experimental results that disagree with the predictions made by Néel theory. This “non-ideal” behavior is often associated with larger particles that are non-uniformly magnetized. In this paper, we use simulations based on Néel theory to demonstrate that non-ideal behavior is also expected in assemblages of SD particles with a range of sizes and shapes previously assumed to be ideal. This effect occurs over a relatively small range of shapes and sizes for magnetite, but a much wider range for titanomagnetites. Our results call into question the typical interpretation of SD particles as ideal magnetic carriers. Instead, we suggest that the geological stability of a rock’s magnetization is influenced not only by the domain state of the magnetic carriers, but also by their shape and composition. This has important implications for paleointensities (and cooling rate corrections thereof), paleodirections, and the dating of viscous magnetizations.

SummaryThis paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

SummaryThe Mendocino Triple Junction (MTJ), where the Gorda, North American, and Pacific plates meet, is one of the most seismically active regions in California. The tectonic movements along the Mendocino transform fault zone (MTFZ), Gorda slab (GS), and northern San Andreas Fault systems (NSAF) lead to high background seismicity rates but relatively low aftershock productivity. To improve the understanding of earthquake processes in the area, we analyze relations between background seismicity, aftershock productivity, and stress parameters. We apply the nearest-neighbor approach to investigate the spatial distributions and properties of background and clustered seismicity, and invert focal mechanisms of events in Voronoi cells for features of the deviatoric stress field. The results indicate that the intensity of background seismicity and aftershock productivity decrease with distance from the MTJ, defined here for simplicity as the hypocenter of the 1992 Mw7.2 mainshock. We also find that the stress regime is the most compressive in the area directly surrounding the MTJ. In the MTFZ and GS, the compressive stress decreases with increasing distance from the MTJ, correlating with the reduced aftershock productivity and background seismicity. In the NSAF, the observed relations between the stress, aftershock productivity, and background seismicity are not clear, possibly due to crustal extension related to the slab window and elevated heat flow. Compared to the MTFZ and GS, the NSAF has a higher foreshock proportion, lower aftershock proportion, and small-to-medium mainshock magnitudes, indicating more swarm-like clusters in this region. The inverted stress regimes in the MTFZ and NSAF are dominated by strike-slip faulting. The GS exhibits mostly strike-slip and normal mechanisms despite the subduction environment, which may reflect slab bending and reactivation of preexisting normal faults.

SummaryUnderstanding when and where strong earthquakes occur is crucial for assessing seismic hazard. Changes in the b-value, which describes how frequently earthquakes of different sizes happen, have been investigated as possible indicators in the space-time vicinity of upcoming large earthquakes. In this study, we investigate short-term b-value variations in the central-northern Apennines (Italy) by comparing an observed earthquake catalogue spanning the period 1987–2025 with a 10 000-year synthetic catalogue generated by a physics-based earthquake simulator. The synthetic seismicity is produced using a three-dimensional seismotectonic model derived from the DISS database and an elastic-rebound framework for earthquake nucleation. We apply a stacking procedure to compute average b-values within symmetrical time windows of ± 15 days and 30 km distance from selected pivot events of moderate to large magnitude. The same methodology is consistently applied to both observed and simulated datasets, enabling a direct comparison of their temporal behaviour. The observed catalogue shows a statistically significant decrease in the b-value in the days preceding earthquakes with Mw ≥ 4.0, followed by a post-event recovery. A comparable pattern is reproduced by the synthetic catalogue, where a pronounced b-value drop precedes pivot events of Mw ≥ 4.5 and is systematically followed by an increase after rupture. The persistence of this behaviour across different magnitude thresholds in the simulated data supports its robustness. These results indicate that physics-based simulations can reproduce short-term b-value variations associated with earthquake nucleation, supporting the relevance of this parameter for investigating the physical processes governing seismicity.

SummaryBasaltic flows and sills of the Central Atlantic Magmatic Province (CAMP) along the eastern North American seaboard have been proposed to be present in buried Mesozoic basins. Their offshore distribution is poorly constrained, yet the strong magnetic and gravity signature produced by basaltic rocks means it should be possible to map them using magnetic and gravity surveys. We conducted forward modeling using existing magnetic and gravity data to identify Mesozoic basins and basaltic units offshore. Onshore and offshore basins containing CAMP basalts in forward models generally predict the best fit with observed magnetic and gravity data. A positive magnetic anomaly over the New York Bight Basin suggests it may contain multiple basalt flows at depths > 2500 m, and scenario testing indicates the Long Island Basin possibly hosts at least one flow. The newly identified Central Bight Basin is unlikely to contain basaltic units, although the adjacent East Coast Magnetic Anomaly may be overwhelming potential basalt signatures within the basin. Deeper basement structures and/or possible interbasinal basalt likely influence existing data, therefore higher-resolution aeromagnetic and marine gravity surveys are needed to constrain CAMP basalt presence in offshore basins.

SummaryGeocenter motion, defined as the displacement of Earth’s center of mass relative to its center of figure, is crucial for maintaining the International Terrestrial Reference Frame origin and quantifying large-scale mass redistribution. However, whether observing geocenter motion by tracking satellite orbits or inferring it using geophysical models, accurately acquiring such subtle motions imposes stringent requirements on the consistency and precision of both tracking data and geophysical models. This study improves geocenter motion estimates derived from the combination of GRACE/GRACE-FO time-variable gravity (TVG) and Ocean Bottom Pressure (OBP) models (the GRACE-OBP method) in two ways. First, we apply a forward modelling technique to mitigate land–ocean leakage in GRACE/GRACE-FO TVG fields, which demonstrably outperforms empirical coastline buffer-zone corrections in controlled simulation experiments. Second, we introduce the Bayesian Three-Cornered Hat (BTCH) method to optimally combine geocenter series derived from multiple GRACE solutions and two independent OBP models (ECCO2 and MPIOM), producing an improved geocenter product without requiring a ground-truth reference. Uncertainty analysis shows that the noise level is governed primarily by the GRACE solution, and that BTCH provides a clearer advantage over equal-weighted averaging when the number of input series is limited, reducing the noise level by about 30 per cent. After restoring atmospheric and oceanic contributions, our improved geocenter series shows good agreement with the CSR SLR-derived geocenter product. Although uncertainty levels vary among individual solutions, the estimated annual and secular trend signals are broadly consistent and show limited sensitivity to the choice of GRACE TVG solution and OBP model. Using the improved geocenter series, we revisit the annual geocenter oscillation and its drivers; the results indicate that cryospheric mass variability and land-ocean mass exchange (i.e. sea-level fingerprints) provide non-negligible contributions to the annual geocenter cycle and improve consistency with observations. Finally, the improved geocenter series yields the lowest uncertainty in degree-1 mass variations, with a global RMS of 0.55 mm. Incorporating these degree-1 terms into mass budget assessments yields secular trends of 38.8 Gt/yr for the Antarctic Ice Sheet and 0.57 mm/yr for global mean ocean mass, highlighting the need for accurate geocenter corrections to support reliable long-term climate monitoring.

SummaryAccurate earthquake hypocentres are fundamental to a wide range of geophysical studies, yet source depth remains poorly constrained in teleseismic earthquake catalogues. Near source surface reflections such as pP, sP, and sS (known as depth phases) provide critical information for resolving hypocentral depth, particularly for intermediate-depth earthquakes. The number of depth phases reported by global earthquake monitoring agencies has declined significantly in recent decades, potentially reducing the precision of resolved earthquake depths. To address this, we automatically detect P, pP, sP, S and sS phase arrivals using teleseismic ad-hoc arrays. We detect these phases for earthquakes in the South American Subduction Zone (SASZ) at depths of 40–350 km and between mb 4.7 to 6.5. The identified phases are integrated with the phases reported to the ISC Bulletin, and used to relocate earthquakes with ISCloc. We assess the impact of incorporating automatically detected, ad-hoc array-derived depth phases on earthquake relocations across the SASZ, and find an improvement in depth resolution for 88.8% of earthquakes. Using this enhanced catalogue we investigate the structure of the Wadati-Benioff zone, focusing on two significant earthquakes: the 2005 Mw 7.7 Tarapacá and 2019 Mw 8.0 Peru events. Finally, we successfully apply our methodology to deep focus earthquakes (350-700 km), which further define the deepest portion of the seismogenic slab. Our results demonstrate the potential for automatically detected, ad-hoc array-derived depth phases to substantially improve the accuracy of teleseismic earthquake hypocentres, and offer further constraint upon slab geometry and seismogenic structure.

SummaryThis study aims to retrieve P waves from seismic ambient noise recorded by a dense local broadband network at the Chémery underground gas storage site, where anticline deformation was previously identified through wells and seismic reflection survey. To this end, we propose a new approach for reconstructing P waves from ambient noise. We process the passive seismic data to reconstruct the body wave component of the empirical Green’s functions (EGF). The retrieved P-wave arrivals were identified and analyzed, revealing that in this dataset, the picked arrival times likely correspond to non-physical head waves rather than direct or diving P-wave arrivals between virtual sources and receivers. Numerical simulations support this idea of non-unique interpretation of the passively reconstructed P-wave arrivals. The simulations suggest the potential for mapping lateral heterogeneities above the critical refractor as a cumulative time-delay, although for this dataset the anticline-induced time-delay is likely within the measurement uncertainties. It is found that the dominance of non-physical head waves over diving waves is primarily due to the large distance from the network to ambient noise sources.

SummaryOceanic subduction zone is the dominant pathway for transporting carbon into the interior of the Earth, and thus plays a critical role in deep carbon cycling. Despite being recognized as a key mechanism for carbon release in subduction zones, the metamorphic decarbonation outflux and efficiency remain subjects of ongoing debate. The thermal structure of subduction zone is widely recognized as a primary dynamic control on metamorphic decarbonation, however, the quantitative relationship between metamorphic carbon outflux and simplified thermal parameters of subduction zones (here defined as φ= slab age × subduction velocity/100 in kilometer) remains poorly constrained. On the other hand, previous studies on metamorphic decarbonation have been conducted within two distinct scenarios: the P-T-dependent decarbonation (PTD) system versus P-T-H2O-dependent decarbonation (PTHD) system, yet a quantitative comparison between these two scenarios remains lacking. In order to investigate the metamorphic decarbonation behavior of subducting slab in the PTD versus PTHD systems, we develop a coupled thermo-petrological model by integrating the thermodynamic dataset of temperature-pressure-(H2O)-dependent CO2 content into the thermal model of subduction zones. Systematic numerical models indicate that the metamorphic carbon outflux in the PTHD system is about 50 per cent lower than that predicted in the PTD system. Meanwhile, the quantitative functional relationship has been built between the metamorphic carbon outflux and φ, which reveals that the decarbonation outflux and efficiency decrease exponentially with increasing φ in both systems. Under present-day widespread subduction thermal conditions with the φ values of around 30 km, both PTD and PTHD system models yield low metamorphic decarbonation efficiency, suggesting that a substantial proportion of slab carbon is likely retained in the slab and transported into the deeper mantle.

SummaryThe Sichuan Basin in China has experienced a number of devasting earthquakes in the past 20 years, particularly on the Longmen Shan fault (LMSFZ) with the 2008 Wenchan and 2013 Luschan events. This study employs a hierarchical, four-dimensional (latitude, longitude, depth, time) clustering framework to characterize seismic activity in the Sichuan Basin. After the identification of spatial features of the region (e.g., faults), we then apply two cluster algorithms on the Longmen Shan fault data and compare the identification of the 2008 and 2013 events. In particular, we apply and compare Density-Based Spatial Clustering of Applications with Noise (DBSCAN) and Bayesian Gaussian Mixture Model (BGMM) on the identification of mainshock-aftershock sequence (as well as any foreshock events). By applying temporal clustering to the dataset and comparing DBSCAN and BGMM methods, we find distinct differences between the results. Specifically, we find that DBSCAN identifies a simple mainshock-aftershock sequence, while BGMM produces a more complex foreshock-mainshock-aftershock sequence. However, both scenarios have been identified within previous work on these events, highlighting that additional analysis is required and that single cluster algorithms should be applied with caution. The work here in comparing machine learning techniques within an integrated clustering framework is timely and will serve as a guide for more in-depth analysis on earthquake patterns and fault dynamics using these methods.

SummarySurface fibre-optic distributed acoustic sensing (S-DAS) typically requires trenching to achieve adequate coupling, which can be challenging and costly in hardrock environments. As an alternative, untrenched S-DAS offers significant time and cost savings, though at the expense of data quality. In this study, we systematically evaluated the impact of untrenched deployment on DAS data quality through a multi-method comparison with conventional sensors, including three-component geophones and vertical-component accelerometers. All recorded data were converted to particle velocity to enable direct amplitude comparisons across arrays. Phase fidelity was assessed using the surface-wave method, cross-power spectral analysis, and ambient noise interferometry. Untrenched S-DAS recorded amplitude levels that differed from those of conventional sensors, particularly above 50 Hz. However, within the surface-wave band (up to 35 Hz), the data quality was sufficient to derive reliable shear-wave velocity profiles, especially when the data were resorted to common-receiver gathers. Sensitivity and coherence decreased significantly at higher frequencies and larger offsets, and a systematic time delay was observed for DAS in the surface-wave band. Ambient noise interferometry was ineffective for the untrenched DAS array, largely due to variable channel coupling and system-specific noise. This study provides a systematic field comparison of untrenched S-DAS and conventional sensors in hardrock settings, outlining both the limitations and the practical potential of this cost-effective deployment method.

SummaryFaults and fractures heal and seal over time, decreasing along-fault permeability, and increasing reactivation stress. This presents a dilemma in geothermal reservoirs as maintaining permeability is crucial for reservoir longevity, but the reactivation of faults to increase permeability can also cause hazardous seismicity. The healing rate of faults is temperature-dependent and shows significant differences under wet and dry conditions. We investigate the healing behavior of a bare, water-saturated fault surface at temperatures up to 163°C through slide-hold-slide experiments. The gneiss sample from the UtahFORGE geothermal demonstration project, is continuously actively probed with P-waves and monitored for passive acoustic emissions radiating from the fault. Our data show that with increased temperatures, the fault surface friction decreases, healing rate increases and the fault becomes more prone to unstable slip. The decrease in friction and increase in healing rate we measure are larger and occur at lower temperatures than previously demonstrated. P-wave amplitudes and P-wave velocities increase during healing, with amplitudes sensitive to temperature but velocities conversely insensitive. We attribute this to a sensitivity of the P-wave amplitude to changes in contact area with P-wave velocity correlating with mechanical compaction, off-fault microcracks, and the formation of wear products during sliding. The sample continues to creep throughout holds during our hotter experiments, but the creep motion does not erase continuous healing. Acoustic emissions spike upon slip reactivation, where higher event rates and higher slip velocities occur as healing progresses—after longer hold times and at higher temperatures. The amplitude of the P-wave, as well as the acoustic emission rate, show precursory signs of spontaneous reactivation and therefore might have potential in forewarning slip.

SummaryWe present a new, variational, fully nonlinear, probabilistic ambient noise tomography method, which estimates subsurface structure and quantifies the corresponding uncertainties directly in three dimensions (3D) from inter-receiver seismic surface wave dispersion data. We use the method to invert for high resolution 3D seismic velocity models of the upper crust beneath Great Britain using seismic ambient noise data recorded around the region – a task that proved too high-dimensional and hence computationally demanding for Monte Carlo sampling to converge to a stable solution. We compare the inversion results from the new method to those obtained from two standard, indirect inversion methods, in which 2D (geographical) surface wave velocity maps and 1D (depth) shear velocity profiles are estimated in two separate, consecutive steps. The results show that the direct-3D scheme preserves better lateral continuity and produces better data fit than the two-step methods, and provides information about lateral correlations that is absent from the two-step solutions. The inversion results are consistent with large-scale geology of Great Britain, and for the first time provide seismologically-imaged evidence of the Great Glen Fault and other major tectonic faults. We therefore propose that direct-3D inversion schemes should be used where possible for surface wave inversion as they provide improved results at little additional computational cost.

SummaryAccurate three-dimensional coseismic deformation fields are critical for fault mechanics analysis and hazard assessment, but the sparse distribution of Global Navigation Satellite System (GNSS) stations often limits reconstruction accuracy. This study proposes Integrated Dislocation and Strain Models (IDSM) that seamlessly integrate GNSS and Interferometric Synthetic Aperture Radar (InSAR) data. This is achieved by combining a surface-constrained strain model and a subsurface-constrained dislocation model, which adaptively optimizes multi-source data weights through Variance Component Estimation (VCE), challenging the traditional reliance on uniformly distributed observations. Simulation experiments demonstrate that under insufficient GNSS coverage, this method improves deformation recovery accuracy by 10 per cent to 70 per cent in the vertical, north, and east components compared to the ESISTEM-VCE (Extended Simultaneous and Integrated Strain Tensor Estimation from Geodetic and Satellite Deformation Measurements-VCE) method, with particularly significant enhancement in the north component. Applied to the 2021 Yangbi MW6.4 and Maduo MW7.4 earthquakes, the study reveals distinct deformation patterns: the Yangbi event exhibits right-lateral strike-slip rupture with a maximum east-west extensional displacement of 87 mm and vertical subsidence of 59.8 mm, showing antisymmetric horizontal deformation around the epicenter. In contrast, the Maduo earthquake is dominated by left-lateral strike-slip motion, with east-west displacement reaching 1.4 m, while north-south and vertical deformations display patchy distributions along the fault. Error analysis confirms accuracy improvements over the ESISTEM‑VCE method. For the Yangbi earthquake, the Root Mean Square Error (RMSE) decreased by 50 per cent (east), 64 per cent (north), and 44 per cent (vertical) at GNSS validation points. Corresponding improvements of 6.1 per cent (east) and 53.5 per cent (north) were achieved for the Maduo earthquake.

SummaryOceanic transform faults (OTFs) have long been viewed exclusively as vertical, strike-slip structures offsetting mid-ocean ridges, yet their deep geometry and structural complexity remain poorly constrained. Thus, key questions persist, including whether OTFs are single-stranded and continuous, whether they maintain vertical dip angles, if they accommodate mixed-mode slip, and what factors control their geometry. This study addresses these questions through a global statistical analysis of teleseismic earthquake focal mechanisms from 150 OTFs across diverse tectonic settings. We introduce stack maps, a novel method that quantifies fault dip and rake, providing a graphical representation of average focal mechanisms. Our findings reveal that while OTFs tend to conform to the standard vertical, strike-slip model, nearly half exhibit deviations, either in dip or motion, challenging the classical view of these plate boundaries. We identify four distinct OTF categories: (1) those adhering to the standard model, (2) non-vertical faults with transtensive/transpressive components, (3) non-vertical faults accommodating strike-slip motion, and (4) vertical faults with a vertical component of motion. Tectonic regime shifts emerge as a primary driver of structural changes, with non-vertical geometries persisting even after the regime reverts to pure strike-slip motion. This structural memory suggests that fault geometry, once established, remains stable over geological timescales of several tens of Myr. By reconciling previously ’unusual’ focal mechanisms with fault structure and dynamics, this work demonstrates that global seismic catalogues, when analysed statistically, offer robust insights into OTF geometry and tectonic regimes.

SummaryAs a critical category of geophysical data, magnetic anomalies play vital roles in geological interpretation, resource exploration and target detection. For most applications involving magnetic anomaly data, the ideal dataset should have uniformly distributed data points, high resolution and completeness without gaps. However, because of the environmental constraints and measurement limitations, magnetic anomaly data obtained from real-world measurements often fail to meet these requirements. Thus, interpolation techniques present effective and cost-efficient technical approaches for processing measured magnetic anomaly data to meet the aforementioned criteria. To our knowledge, current research on magnetic anomaly data interpolation has primarily focused on gridding methods for interpolating irregularly sampled data into gridded data and super-resolution interpolation methods aimed at enhancing spatial resolution. Meanwhile, studies on interpolation methods specifically designed to fill large-area data gaps remain relatively scarce. To address the challenge of reconstructing large-area missing magnetic anomaly data, we propose a data-driven method for magnetic anomaly data gap filling. First, based on the analysis of the characteristics of magnetic anomaly data, we construct an open-source magnetic anomaly interpolation dataset (MAID) specifically designed for magnetic anomaly data interpolation tasks. Subsequently, we develop a magnetic anomaly data gap-filling generative adversarial network (MADGF-GAN) tailored for magnetic anomaly data gap filling. Upon sufficient training on the MAID training set, MADGF-GAN can directly fill gaps in given magnetic anomaly data. Finally, the effectiveness of MADGF-GAN is validated using four test samples from the MAID test set and Afghan aeromagnetic data. Compared with four existing interpolation methods, MADGF-GAN demonstrates considerable advantages in terms of interpolation accuracy, computational efficiency and practicality. This study demonstrates the potential of data-driven approaches in magnetic anomaly data processing, providing crucial technical support for related geoscientific applications.

SummaryBioleaching is a biologically facilitated process that helps to dissolve valuable metals in order to extract them from the mineral gangue. Applied in the field to heap ores, its efficiency mainly depends on solution flow inside the heterogeneous heaps, which is often tortuous and can remain stagnant in the pores and crevices between the particles. Methodologies that can help to monitor the bioleaching processes are therefore needed to improve operational efficiency. In this article, we present for the first time preliminary laboratory-scale investigations on spectral induced polarization (SIP) during the bioleaching of chalcopyrite (CuFeS2) containing ore material from a mine in Chile. Two column experiments representing different stages of the bioleaching process were monitored under un-saturated and highly acidic environment (pH ~2). Our objective was to explore the feasibility of SIP for detecting changes in electrical properties potentially associated with bioleaching-induced mineral dissolution and alteration. The results show a rapid decrease in SIP phase shift and imaginary conductivity during the early stage of bioleaching, while the real conductivity remains relatively stable. At a more advanced stage of bioleaching, the phase response is weaker and more stable. A relaxation time distribution (RTD) analysis was applied to further investigate changes in polarization mechanisms. Prior to bioleaching, the RTD exhibits a well-defined peak consistent with polarization controlled by sulfide mineral grains, whereas after one month of bioleaching the RTD broadens and shifts toward larger relaxation times, accompanied by a decrease in chargeability. This combined evolution suggests bioleaching-induced modifications of electrochemically active surfaces, potentially related to mineral dissolution and the formation of passivation layers. Estimated particle sizes derived from the RTD analysis are consistent with scanning electron microscopy observations. Although, the absence of a dedicated abiotic control column prevents us from attributing these changes unambiguously to bioleaching alone, these results highlight the potential of SIP as a non-invasive, real-time and integrative tool to monitor leaching processes and to identify zones that may remain weakly affected by leaching.

SummaryModeling crustal deformation induced by fault slip is a fundamental problem in structural geology and seismology. However, the challenges of data sparsity and spatial discontinuity impose significant limitations on conventional forward and inverse methods, often resulting in low computational efficiency and limited accuracy. Although AI-based approaches such as Physics-Informed Neural Networks (PINNs) and Physics-Encoded Finite Element Networks (PEFEN) offer new solutions for sparse-data problems governed by physical laws, their underlying assumption of spatial continuity conflicts with the inherent displacement discontinuities of fault-slip fields. To address this limitation, we propose a novel method—the Split-Node Physics-Encoded Finite Element Network (SN-PEFEN)—which integrates the node-splitting mechanism into the PEFEN framework. By explicitly encoding spatial discontinuities into the nodal topology during mesh preprocessing, SN-PEFEN not only overcomes the theoretical limitations of existing PEFEN models in handling discontinuous fields but also maintains physical consistency. We apply SN-PEFEN to perform forward and inverse modeling of deformation fields induced by complex fault slip in both 2D and 3D heterogeneous media. For a model with over one million degrees of freedom, the forward simulation achieves over 40× speedup compared to traditional FEM (∼1,800 s vs. 42 s), while maintaining comparable accuracy. In inverse modeling, the solution converges within only 100 iterations, with a total runtime of approximately 2,000 s, demonstrating high computational efficiency. This method establishes a new high-efficiency paradigm for analyzing complex discontinuous deformation in geomechanics, offering promising applications in multi-fault system analysis and fault-slip inversion. Furthermore, SN-PEFEN facilitates rapid, physics-based assessments for emergency seismic response and disaster management, while laying the groundwork for next-generation data-driven regional earthquake early warning systems.

SummaryUltralow velocity zones (ULVZs) at the Earth’s core-mantle boundary (CMB) are marked by substantial reductions in seismic velocities. They are often associated with significant increases in density, providing important insights into deep Earth composition and dynamics. In this study, we investigate ULVZs beneath eastern and southern Asia, regions associated with long-term subduction, by analyzing high-frequency (~1 Hz) ScP waveforms recorded at the small-aperture KZ array. After correcting for attenuation along the ScP path, we perform a grid search to match the observed waveform complexities with synthetics generated for a comprehensive suite of 1D ULVZ models. The best-fitting models for each event constrain ULVZ thickness, P- and S-wave velocity reductions, and density anomalies, revealing widespread but laterally variable ULVZ structures, although the influence of finite ULVZ geometry cannot be entirely excluded. The correlations among these parameters point to iron-rich chemical heterogeneity as the dominant origin of the imaged ULVZs, likely reflecting iron enrichment associated with long-term subduction processes.