Geophysical Journal International

SummaryDeep learning methods such as neural networks offer a rapid alternative to geophysical inversion and are increasingly applied to airborne electromagnetic (AEM) data for a range of earth science problems. Such networks rely on both the network architecture as well as suitable training datasets to produce realistic earth resistivity models. Here we develop a Convolutional Neural Network (CNN) to invert AEM data and pay particular attention to the characteristics of the training dataset. As well as covering a realistic range of subsurface property distributions, the training data must emulate the noise and data gap characteristics of the target survey data. We find neural networks trained on noise free data, while performing very well on validation data, perform poorly on real world data, often inducing geologically unrealistic artifacts in the resistivity model. Conversely, training data conditioned with noise and data gap characteristics appropriate for the survey of interest, fit the validation data less well, but produce fewer artifacts when applied to survey data. Choice of loss function metric and data normalisation method is important and may vary between surveys. Despite strong advances in neural network design, developing training data with broad enough characteristics to be generalised across many surveys and data acquisition platforms remains an important challenge.

SummaryAccurate separation of the Earth’s internal and external geomagnetic fields requires robust modeling, especially during the geomagnetic storms. Because the comprehensive CHAOS model is primarily designed to represent the Earth’s magnetic field under geomagnetically quiet periods, residuals of satellite observations relative to the CHAOS model during storms exhibit systematic dependencies on magnetic local time (MLT) and the disturbance levels, highlighting limitations in its representation of external field morphology. To address this limitation, we developed an empirical correction for the CHAOS model that specifically targets the magnetic signatures of the residual ring current at low and mid-latitudes. Using the high-quality scalar magnetic data from Swarm A/B and the Macau Science Satellite-1 (MSS-1) at nighttime hours, between November 2023 and September 2025, we modeled the systematic perturbations remaining after CHAOS background removal. Our separable parameterization explicitly resolves MLT structure while incorporating dependencies on geomagnetic disturbance levels. Event analyses demonstrate that our correction model can successfully capture the pronounced dawn-dusk MLT asymmetries of the external field during storms main and recovery phases. Statistically, the empirical model robustly captures residual ring current signals across all evaluated nighttime sectors, yielding the greatest improvements at dusk. Ultimately, the proposed correction effectively mitigates the activity-dependent biases in the CHAOS model’s external field, yielding a more reliable representation for the accurate interpretation of storm-time magnetic residuals.

SummaryGeophysical joint inversion offers a robust approach for integrating diverse geophysical datasets, leveraging their complementary information to mitigate the nonuniqueness inherent in individual data inversion. This approach thereby enhances model consistency and improves subsurface characterization. In this study, we investigate the effectiveness of three-dimensional (3D) joint inversion of gravity, magnetic and magnetotelluric (MT) data for deep geothermal exploration. We first implemented a reweighting optimization strategy for gravity and magnetic data to increase the vertical resolution of jointly inverted models of potential data. Subsequently, we integrated the optimized gravity and magnetic inversion scheme with MT data using a cross-gradient constraint for comprehensive joint inversion. Synthetic experiments validated the effectiveness of our developed reweighting optimization strategy and joint inversion scheme, demonstrating improved structural resolution and model reliability compared to separate inversions. Finally, we applied the methodology to the Yanggao geothermal field in northern Datong Basin, China, to assess its applicability for deep geothermal exploration. The joint inversion approach, offering an advantage over separate inversions, provides enhanced subsurface characterization by more accurately delineating hydrothermal pathways and potential geothermal reservoirs of the geothermal system within the investigated area.

SummaryLocal seismic topographic amplification of rock slopes is a primary cause of earthquake-induced damage in mountainous engineering. To resolve the scattering of plane SH waves by individual rock slopes, this study introduces an analytical model derived from the Uniform Theory of Diffraction. Conventional wave function expansion methods (WFEM) often struggle with non-closed free boundary conditions across varying elevations and tend to suffer from series divergence under high-frequency incidence. This approach precisely partitions the slope into two semi-infinite ideal wedge domains sharing a common face, reconstructing the generalized geometric wave field through rigorous global ray tracing. Second-order coupled diffraction terms with a slope-diffraction correction are introduced to describe the two-way interaction between the slope crest and toe, while the UTD transition functions ensure continuity of the wavefield across shadow boundaries. Systematic numerical evaluations in both frequency and time domains reveal a frequency-selective interference mechanism with distinct spatial variability. Under horizontal grazing incidence, energy is strongly focused on the leeward face, while the back-slope forms a high-frequency seismic shadow zone due to intense destructive interference between direct and secondary diffracted waves. At certain oblique incidence angles, repeated scattering of broadband waves between the slope crest and toe leads to strong local resonance in the slope region. Findings indicate that the location of maximum dynamic response shifts sensitively with the incident wavefront inclination. The results show that the maximum amplification does not always occur at the slope crest, highlighting the importance of incidence angle in seismic microzonation and site planning for mountainous areas.

SummaryThe statistics of earthquake populations are governed by stress changes and event interactions at various scales which promote non-linear event cascades and seismicity clustering. Earthquake clustering in space and time is observed in both tectonic and volcanically dominated regions; however, notable differences exist for instance during eruptive sequences. Here, we characterize seismicity statistics in California and Hawaii and quantify differences in background rates and clustering. We compare three statistical approaches, i.e, i) Reasenberg declustering, ii) Gamma distribution fit, and iii) Nearest Neighbor Clustering which are applied to seismicity in Hawaii, southern California, and ETAS catalogs with known clustering characteristics. Tests with ETAS-catalogs suggest that background rate variations are more consistently resolvable using the Nearest Neighbor and Gamma distribution algorithm without requiring parameter adjustments to the specific study area. The Nearest Neighbor method tends to overestimate background rates, whereas the Gamma distribution algorithm shows systematic underestimation, suggesting that an ensemble model can improve the performance over any single method. Uncertainties in background rate estimates based on small magnitude events (e.g., M≥2.5) are high, in particular for the Reasenberg method, which overpredicts true rates by up to a factor of 4, and extrapolations to larger magnitudes need to be treated cautiously. Our analysis shows that the combination of the three approaches generates new insights into the processes that govern seismicity clustering. The analysis of southern California seismicity produces consistent results between the different methods, and the distribution of background seismicity rates is indistinguishable from Poissonian rates at magnitudes above 3.5. Conversely, the seismicity in Hawaii exhibits notable non-Poissonian background rate changes. Results from the Nearest Neighbor clustering analysis suggest notable differences between tectonic and volcano seismicity, including background rate changes and more pronounced spatiotemporal clustering around active volcanoes. Hawaiian seismicity is characterized by more localized spatial-temporal clustering, in particular during the 2018 eruption. This difference may be a result of magmatic, hydrothermal and eruptive processes, which make volcanic areas susceptible to stress perturbations and earthquake triggering at small scales.

SummaryMapping the internal structure of sedimentary basins, including the depth to basement, is valuable for a variety of geological applications, particularly the identification of natural resources such as groundwater and minerals, or of the geological structures associated with them. The magnetotelluric (MT) method has proven to be a reliable technique for imaging complete sedimentary sequences, especially in areas where thick sedimentary cover makes imaging challenging. However, due to the high regularisation required by MT inversion procedures, it is limited in its ability to precisely locate geological interfaces, which is crucial for exploration. Building on previous work where we imaged the basement using interface probabilities derived from 1D probabilistic inversion of MT data, we present a new method which allows for the simultaneous classification of multiple interfaces across an entire survey, incorporating constraints on the spatial relationships between models. As a result, we obtain spatially consistent model ensembles and classified interfaces corresponding to transitions between layers of consistent electrical resistivity. This classification effectively reduces the size of the model ensemble, decreasing uncertainty in the estimated depth of the interfaces of interest. We apply this method to the Eucla sedimentary basin in Western Australia, analysing 550 MT sites distributed across 12 profiles. The ensemble clustering clearly identifies three interfaces which are consistent with the sedimentary succession expected from this basin. The results show great consistency across all the survey, providing valuable insights into the geological setting of the area. This research demonstrates the capability to reliably image the structure of a sedimentary basin using MT, within a workflow that integrates probabilistic inversion and model ensemble classification.

SummaryThe Tournemire Underground Research Laboratory (URL) operated by ASNR (Autorité de Sûreté Nucléaire et de Radioprotection) is a French test site specifically developed to study and comprehend the confinement properties of low-porosity claystone formations. More specifically, this site allows developing an improved knowledge of the transport properties of the Toarcian (Early Jurassic) shale formation. This claystone shares indeed similar characteristics to the Callovo-Oxfordian claystone formation encountered at the Meuse/Haute-Marne underground research laboratory operated by ANDRA (Agence Nationale pour la gestion des Déchets RAdioactifs) to study the feasibility of the long-tern storage of nuclear wastes. In the present study, 18 rock samples of the Toarcian formation are used to study (i) the induced polarization properties of this formation at full water saturation and (ii) the desiccation of this claystone by looking at the evolution of the complex conductivity spectra with the decrease of the water saturation in the frequency range 1 mHz-45 kHz. In this second set of experiments, the complex conductivity spectra are fitted with a double Cole Cole complex conductivity model. We then obtain the values of the first and second Archie’s exponents and estimates of both surface conductivity and normalized chargeability. We also test relationships between the surface conductivity, the normalized chargeability, the quadrature conductivity, and both the cation exchange capacity and water content of these clay-rocks. Furthermore, the laboratory data are favorably explained by the so-called Dynamic Stern Layer (DSL) model developed in the last decade. In order to test this new knowledge in field conditions, an in situ experiment is performed in a gallery of the Tournemire URL to test the ability of the method to image the damage zone around a 80-cm hole drilled normal to the gallery axis. For this purpose, a set of 5 wells equipped with a total of 240 non-polarizing electrodes is successfully used to get both electrical conductivity and normalized chargeability tomograms, which are in turn used to image the water content around the main hole. The results are analysed in terms of stress distribution including the effect of plasticity.

SummaryAcoustic wave propagation under pressure is inherently nonlinear and often exhibits hysteresis in elastic properties. Understanding this behaviour is important in rock physics and exploration seismology because pressure-dependent elastic responses contain information on evolving rock microstructure. We present a petrophysical model that describes the relationship between acoustic P- and S-wave velocities and effective rock pressure during both loading and unloading phases. Instead of prescribing a single microscopic mechanism, the model introduces an effective extensive state variable that represents the cumulative effect of pressure-induced microstructural evolution. Depending on the dominant physical process, this variable may be associated with compliant pore and microcrack volume, grain-contact evolution, or other pressure-sensitive structural descriptors. The model distinguishes between open and closed microstructural states and describes their stress-dependent evolution through different sensitivities during loading and unloading, thereby reproducing irreversible velocity-pressure behaviour within a unified framework. Model parameters are estimated from laboratory-measured acoustic velocities using linearized joint inversion. The inversion further enables the reconstruction of an acoustically weighted effective descriptor that characterises the pressure-dependent contribution of compliant microstructural evolution to the observed velocity response. In parallel, we apply the Preisach theory of hysteresis to pressure-dependent acoustic velocities. In this framework, compliant microstructural features are represented by effective hysteretic units (hysteron) characterized by distributed closing and opening pressures. The resulting Preisach model reproduces the measured hysteresis with high accuracy and enables a model-dependent reconstruction of the Preisach density function, which provides an effective, non-unique representation of hysteretic switching contributions across closing and opening pressures. By combining the petrophysical and Preisach approaches, we introduce two diagnostic quantities, the Ratio of Hysteron Jumps (RHJ) and the Cumulative Hysteron Jumps (CHJ), which quantify the relative and cumulative contributions of hysteretic elements throughout loading and unloading. While both approaches reproduce the observed velocity–pressure hysteresis, they provide complementary information. The petrophysical model describes the macroscopic evolution of elastic properties via an effective state variable and its acoustically weighted representation, whereas the Preisach model resolves the distribution and activation of hysteretic switching contributions in Preisach space via the reconstructed density function. The results demonstrate that acoustic hysteresis measurements contain not only information on effective elastic properties but also recoverable signatures of pressure-dependent microstructural evolution relevant to seismic interpretation.

SummaryThe São Francisco Paleocontinental block (SFPB) and the Congo Paleocontinental block (CPB) composed, prior to the Atlantic opening in the Cretaceous, the São Francisco-Congo Paleocontinent (SFCP). The SFCP was one of the landmasses involved in the Neoproterozoic Brasiliano-Pan African Orogeny and the assembly of Western Gondwana. Both the SFPB and CPB underwent reworking during the amalgamation of Western Gondwana, forming the orogenic belts that surround the preserved portions of the SFPB and CPB that are known today as the São Francisco and Congo cratons. Although the crustal portions in the outermost areas of these paleocontinents are no longer visible in the surface due to this reworking, the limits of these paleocontinents remain relatively intact at subcrustal depths. The objectives of this study are to delimit the SFPB and to image the structures beneath the Araçuaí belt at mantle lithosphere depths. We present the results of our study which were obtained by the multiple-frequency seismic tomography method for the SFPB and adjacent structures between depths of 68 and 587 km. Data were obtained by processing broadband seismograms recorded between 11-August-2023 and 31-July-2024 for P and PKIKP phases. The processed database comprised 9 254 residuals, which were added to the database of a previous regional study resulting in 97 150 traveltime delays. Checkerboard resolution tests using patterns with horizontal dimensions of 390×390 km reveal very good recovery between depths of 68 and 226 km, although the sharp vertical transition between the checkerboard layers is not fully recovered between the depths of 316 and 384 km. The checkerboard pattern with horizontal dimensions of 312×312 km reveals that the best recovery is observed in the southern and eastern portions of the São Francisco craton (SFCr) and in the Borborema province (BP) between depths of 68 and 226 km. Our results show a high-velocity anomaly with roughly NE-SW trend that extends from the Paranapanema block (PaBl), through the entire extension of the SFCr up to the southern BP. This anomaly contains two major segmentations, one in the southern Brasilia belt, marking the division of the PaBl with the SFPB, and the other dividing the eastern and western arms of the SFPB with a roughly NW-SE trend, interpreted as the Paramirim aulacogen. We also observe a high-velocity anomaly beneath the Araçuaí belt (ArBe) that dips with a SE direction. This anomaly starts at the northern end of the Abre Campo suture at a depth of 136 km, deepens to the south, and is interpreted as the foundered cratonic lithosphere that originated from the delamination of the SFCr due to the passage of the Trindade plume. To reinforce our interpretation, we tested several synthetic models with geometries that were based on our results and previous studies in the area. The synthetic model which presents the most resemblance to the real data model contains a high-velocity anomaly beneath the Araçuaí belt between depths of 226 and 384 km. A strong low-velocity anomaly beneath the Ribeira Belt is interpreted as material rising from the Nazca Slab dehydration.

SummaryNon-Newtonian pore-fluid behaviour is incorporated into a generalized Lagrangian–dissipation framework for wave propagation in fluid-saturated porous media. In contrast to approaches that rely on prescribed pore-scale flow patterns, pore geometries, or other microscopic structural assumptions, the present formulation introduces nonlinear rheological effects directly at the macroscopic constitutive level. The resulting model remains consistent with continuum mechanics and Biot-type poroelasticity while allowing non-Newtonian effects to enter three fluid-related dissipation mechanisms: Darcy-type fluid–solid friction, local-deformation (LD) relaxation, and bulk-viscous dissipation. Nonlinear dissipation potentials are constructed for these mechanisms and then reduced, through a first-harmonic approximation, to effective frequency-dependent coefficients for harmonic wave analysis. Numerical examples show that the non-Newtonian Darcy-friction and LD mechanisms primarily modify the characteristic frequencies of the corresponding relaxation processes, whereas the non-Newtonian bulk-viscosity mechanism affects both the frequency range and the strength of attenuation and dispersion, with a strong sensitivity to the assumed magnitude of the pore-fluid bulk viscosity. The framework also allows alternative nonlinear representations of Darcy-type fluid–solid friction to be compared within the same Lagrangian–dissipation structure; an alternative nonlinear dissipation operator produces wave-propagation predictions nearly indistinguishable from those obtained using the power-law Darcy-friction model, suggesting that different nonlinear constitutive assumptions may lead to similar effective attenuation structures at the macroscopic scale. Calibration against laboratory measurements on tight carbonate rocks shows that the observed attenuation and velocity-dispersion behaviour can be reproduced using physically interpretable non-Newtonian parameters. These results suggest that non-Newtonian pore-fluid rheology may contribute to frequency-dependent wave-propagation effects in fluid-saturated porous media and provide a flexible macroscopic framework for investigating fluid-related dissipation beyond conventional Newtonian poroelastic models.

SummaryGravity data is widely applied to determine Moho topography due to globally dense and homogeneous gravity observations. Compared with gravity data, gravity gradients can reveal more details of Moho structure, which are more competitive. However, gravity gradients must be transformed into gravity data for common Vening Meinesz-Moritz (VMM) methods, which causes loss of some high-frequency signals. Thus, an improved method by constructing direct function relation between Moho topography and gravity gradients is proposed in this paper. Subsequently, A case study in Tibetan Plateau is conducted using this improved method, and a new Moho topography with higher spatial resolution is determined. More refined tectonic implications can be also revealed from this new Moho topography: (1) the prevailing wavelengths of detected Moho folds are approximately 578 km in the east-west direction and approximately 515 km - 708 km in the north-south direction; (2) two possible lower crustal mass flow channels can be identified clearly according to the traces leaving on Moho topography in the southeastern Tibetan Plateau (3) clearer Moho subsidence, approximately 5 km, can be observed in the central Tarim Basin.

SummaryThe northeastern margin of the Tibetan Plateau exhibits pronounced spatial variations in present-day vertical crustal deformation, yet the mechanisms controlling these differences remain debated. In particular, contrasting uplift patterns are observed between the Liupanshan tectonic belt and the Maxianshan Fault, two key tectonic structures along the northeastern expansion front of the plateau. Here we develop two three-dimensional viscoelastic finite element models constrained by geological structures, seismic profiles, and GNSS observations to investigate the controls on vertical deformation. The modeling results show that when the mid–lower crustal viscosity is laterally uniform between the Longxi Basin and the Ordos Block, vertical deformation of ~2 mm/a is concentrated within a ~50 km-wide zone near the Liupanshan tectonic belt, consistent with leveling observations. This deformation pattern mainly reflects crustal shortening and thickening dominated by pure shear. In contrast, the observed ~4 mm/a uplift along the hanging wall of the Maxianshan Fault can only be reproduced when weak mid–lower crust and enhanced crustal flow are introduced beneath the Linxia Basin, indicating that northeastward mid–lower crustal flow is converted into vertical uplift along the fault. These results suggest that the Maxianshan Fault represents the northeastern termination of mid–lower crustal flow within the Qilian Block, whereas deformation in the Liupanshan tectonic belt is primarily accommodated by crustal shortening. Our findings highlight the fundamental role of lateral rheological heterogeneity in controlling vertical deformation and provide new insights into the mechanisms of present-day lateral growth of the northeastern Tibetan Plateau.

SummaryAmbient noise based monitoring of subsurface velocity changes is possible without the explicit retrieval of Green’s functions by correlation. Velocity variations can directly be observed from the fluctuations in the spectrograms of ambient noise time series or their cross-spectra. This approach is more resource efficient than the conventional Green’s function based monitoring and ideally suited for edge processing and the analysis of large volumes of data for example from fibre-optic records. The spectral fluctuations result from wave propagation in the heterogeneous subsurface and interference between different scattering paths. The fluctuations are directly related to the occurrence of coda waves in the Green’s function. In contrast to coda wave interferometry in the time domain, monitoring the evolution of spectral fluctuations allows for measurements of velocity changes that are unbiased by changes in the source spectrum. Recognizing that the imprint of the subsurface heterogeneity is directly encoded in the spectral fluctuations opens a new perspective on the investigation of subsurface scattering and attenuation.

SummaryWe develop a discrete event modeling framework that captures the progression of geophysical systems toward catastrophic failure through sequences of distinct damage events. By representing geophysical system evolution as a succession of temporally accelerating and amplitude-varying events, the framework reveals how finite-time singularities, both logarithmic and power law types, naturally emerge from the interplay between shrinking interevent intervals and growing event magnitudes. This event-based perspective, which can be viewed as a discrete representation of progressive damage accumulation in heterogeneous geomaterials, provides an intuitive physical understanding of rupture processes, highlighting how precursory signals such as accelerating strain rate, event frequency, and energy release can be traced back to simple underlying mechanisms. A mean-field damage-based formulation further links the observed power law exponents to the evolving stiffness of the geophysical system under constant or time-varying stress. Incorporating stochastic fluctuations, the model captures the inherent randomness of natural systems leading to the emergence of stochastic finite-time singular behavior. Together, these results establish a simple yet powerful framework for interpreting the dynamics of catastrophic events, providing a common event-based perspective on observations from landslides, glacier breakoffs, volcanic eruptions, and other related processes, and strengthening the physical foundations of early warning and hazard forecasting.

SummaryThe Accurately Controlled Routinely Operated Signal System (ACROSS) is an artificial seismic source that generates highly repeatable and stable seismic waves for high-resolution temporal monitoring. While direct P- and S-waves from ACROSS have been widely used to monitor earthquakes, volcanic activity, and environmental changes, the scattered waves that arrive later, known as coda waves, have received limited attention. In this study, we used coda waves generated by ACROSS to monitor subsurface seismic velocity changes (dv/v) over 10 months. The ACROSS signals were recorded by a seismic array of 14 seismometers deployed approximately 3 km from the source in Morimachi, central Japan. By deconvolving the records with a known source function, we obtained transfer functions corresponding to band-limited Green’s functions. Coda wave interferometry applied to the coda portions of these transfer functions detected temporal variations in dv/v, exhibiting both seasonal long-term variations and rainfall-induced short-term variations. Comparative analyses using direct P- and S-wave travel-time changes from the same ACROSS data and ambient-noise interferometry show that, under the present observational conditions, these conventional methods primarily detect long-term variations; however, they do not resolve the rapid, transient short-term variations. In contrast, the stable ACROSS source combined with the broader sampling of coda waves provides superior sensitivity to environmental changes across both timescales over kilometer-scale distances, highlighting its effectiveness for high-temporal-resolution monitoring. To investigate the mechanisms responsible for the observed dual-timescale variations, we used a poroelastic model with precipitation input and considered thermoelastic effects and groundwater-induced changes as possible contributors to the long-term variations.

SummaryThe Alpine-Himalayan orogen preserves geological remnants of subducted lithosphere from the Paleotethys and Neotethys oceans and intervening microcontinents. This orogenic belt displays distinct segments separated by discontinuities aligned with paleo-transform faults, reflecting laterally varying ocean opening and closure histories. We investigate how upper and lower mantle slab remnants imaged by seismic tomography may correlate with Paleo- and Neotethyan subduction zones, focusing on the Anatolian, Aegean, and Iranian segments. Using plate tectonic reconstructions in a mantle reference frame, we predict slab subduction timing, amount, and location and make a semi-quantitative dimensional comparison with seismic tomographic images beneath the Eastern Mediterranean and Middle East. We identify three major Neotethyan slabs: the Pontides and Egypt slabs detached in the Late Cretaceous and now occupy the upper lower mantle. The Cyprus slab remains mainly in the upper mantle in which it overturns. To account for subducted lithosphere that reconstructions predict, we interpret that the Cyprus slab lies overturned in the lower mantle down to ∼1000 km. We interpret a large lower mantle anomaly volume between 2200-1500 km (the Herodotus anomaly) as representing Paleotethyan lithosphere that subducted between ∼240-180 Ma. Our reconstruction-tomography comparison suggests that current slab positions likely reflect past detachment locations, while geometries indicate paleo-trench absolute motions, including a Late Cretaceous same-dip, double subduction configuration. Slabs associated with Aegean, Anatolian, Iranian, and Tibetan segments define stable mantle provinces with boundaries aligned with transform-related orogenic segmentation, implying minimal paleo-longitudinal mantle flow since the Early Mesozoic. Our findings indicate upper and lower mantle structure primarily results from near-vertical slab sinking after detachment since the Triassic, without evidence for deflection by lateral components of mantle convection.

SummaryAs the demand for European domestic mineral resources increases, exploration is increasingly focused on deeper and covered targets, often located in populated regions. In such environments, the application of electromagnetic (EM) techniques is severely challenged by anthropogenic infrastructure. Metal-bearing structures, including power lines, pipelines, railway tracks, and mine shafts, can strongly distort the EM fields, producing significant artefacts in both the measured data and the resulting inversion models, thereby biasing or hindering geological interpretation. These effects currently limit EM investigations in inhabited areas, despite their potential for non-invasive and efficient exploration of challenging subsurface targets. Focusing on semi-airborne electromagnetic (sAEM) data, we develop a data-driven approach to address infrastructure-related effects. By inverting affected data for the distribution of infrastructure currents, we account for the full coupling between the EM transmitter, the conducting Earth, and the metal infrastructure. In sAEM field experiments, involving a grounded dipole transmitter and an airborne receiver system, we study the distortion by self-built infrastructure and investigate both measured currents in infrastructure as well as its EM coupling. We validate the functionality of our current inversion approach by reproducing the measured currents in a simplified 1D scenario. With a synthetic 3D study, mimicking a sAEM campaign in the presence of infrastructure, we find that the current inversion approach is capable of separating and correcting for infrastructure effects even if the true resistivity distribution of the subsurface is unknown. In particular, other conductive subsurface structures can be resolved well, even directly below infrastructure. On a small sAEM field data example, affected by the impact of a metal-built conveyor belt and pipeline, we demonstrate the applicability of this framework for real-world scenarios.

SummaryPermeability is a critical parameter for reservoir characterization and hydrocarbon development, yet its accurate prediction remains a challenge. Pore structure, as the intrinsic factor governing both the elastic and hydraulic transport properties of rocks, serves as a bridge between these properties and facilitates permeability prediction from well logs and seismic data. To accurately describe the variation in the physical properties of tight sandstone reservoirs with pressure, our study aims to construct a physical model that relates rock elastic properties with permeability. We developed a dual-porosity rock physics model by coupling David & Zimmerman’s pore-structure inversion method with Dienes’s percolation theory. Our model divides the pore space into pressure-insensitive stiff pores and pressure-sensitive compliant microcracks. By inverting the microcrack density and aspect ratio distribution-which evolve with pressure-from elastic wave velocities, we quantitatively predict permeability variations using Dienes statistical percolation model. To validate this model, we measured porosity, permeability, and P- and S-wave velocities on four tight sandstone samples under effective pressures of 5 to 50 MPa. The results show that the proposed model accurately captures the evolution of both elastic parameters and permeability with effective pressure, demonstrating strong predictive capability for the experimental data. The significance of this study lies in achieving a quantitative relation between elastic and transport properties by explicitly characterizing the pore structure and integrating percolation theory, which is then calibrated with real experiment data, thereby providing solid physical basis for better prediction of reservoir permeability using acoustic logs and seismic data.

SummarySatellite formation configuration plays a crucial role in formulating the next generation gravity mission aiming at improving the accuracy and spatiotemporal resolution of Earth’s gravity field modeling. For the proposed GRACE-Pendulum (GRACE-P) formation, technical challenges and limited accuracy gains from the third observation link led to the adoption of only two baselines, without incorporating the third link perpendicular to the GRACE-type baseline, resulting in potential limitations in modeling accuracy. With future technological advancements overcoming these constraints, this study first proposes a complete triangular GRACE-Pendulum-full (GRACE-PF) formation, which includes the omitted third link, to explore the upper limit of accuracy achievable with this formation geometry. Building on this, we further introduce a novel satellite formation, namely the Three-Observation-Link (Three-OL) formation, inspired by the Chinese TianQin gravitational wave detection project. While geometrically similar to GRACE-PF, Three-OL does not include a GRACE-type baseline. To systematically evaluate the performance of these formations, a closed-loop simulation experiment was conducted, followed by a comparative analysis of GRACE-P, GRACE-PF, and Three-OL formations in the spectral-spatial domain. Specifically, spatial-domain results show that GRACE-PF outperforms GRACE-P by 28.7 per cent globally, confirming the accuracy contribution of the third link. Three-OL exhibits superior performance to GRACE-PF, with higher accuracy for mid-high degree signals in the spectral domain, and an 11.3 per cent accuracy improvement over land areas in the spatial domain. Furthermore, compared with the mainstream Bender and GRACE formations, Three-OL reduces the residual RMS over land areas by 28.4 per cent and 75.9 per cent, respectively, while also demonstrating improved capability for recovering sub-monthly gravity signals. In summary, this study reveals the accuracy potential of the complete GRACE-P configuration and proposes a novel formation capable of supporting high-precision temporal gravity field modeling. This work provides a feasible and efficient formation design for gravity satellite missions dedicated to Earth’s gravity field detection.

SummaryThe North Maladeta fault system, reactivated after the Alpine orogeny as a normal fault (Ortuño et al., 2008; Ortuño and Viaplana-Muzas, 2018) extends for almost 75 km across the border between Spain and France. This system is considered the most probable source of the historical M6.3 Ribagorça earthquake (1373). Within the system, the Coronas Fault is a secondary fault subparallel to the North Maladeta Fault and located in the southern slope of the Maladeta massif. Recent seismic activity in the massif has markedly increased since 2020, with previously unknown clusters of shallow earthquakes detected. Using a Template Matching approach, we verified that this apparent increase is not an artifact of network evolution, identifying additional clusters during 2014–2016, with far fewer events however. Relocation and focal mechanism analysis suggest that most of the seismicity is associated with the Coronas Fault, which appears to be less vertical and flatter than previously assumed. The seismicity is distributed on both sides of the fault, surrounding a 5–6 km wide aseismic segment. Another cluster, active since 2022 (Gerbosa cluster), likely corresponds to a previously unknown parallel fault.Temporal correlations indicate that cluster activation coincides with periods of low surface water levels. In particular, one cluster at the Western part of the Coronas fault, beneath the most rapidly shrinking part of the Aneto glacier, reactivated in late June 2020 and 2024, suggesting a possible influence of glacier melting. In the eastern Coronas fault, the two largest clusters were triggered within days of major rainfall events recorded in the past five years, consistent with a combination of short-term infiltration processes, and longer-term surface unloading from reduced snow and water storage since 2019. These observations raise the possibility of a climate–seismicity link in the Maladeta Massif, whereby rainfall and snowmelt infiltrations, and/or surface unloading may contribute to triggering shallow seismicity. Longer-term monitoring and detailed studies are required to confirm these processes and assess their potential implications for seismic hazard in the Pyrenees.