Geophysical Journal International

SummaryTianshan Mountains in Central Asia are one of the largest and most active orogenic belts in the world, characterized by complex structures and strong seismic activity. In this paper, we use recently updated GNSS data to self-consistently estimate the slip rates of major faults in Tianshan region via the elastic block model. Our results indicate that crustal deformation in Tianshan region is predominantly manifested as crustal shortening, regulated by foreland thrust belts and intermontane basin boundary faults. The shortening rate decreases from 14.2 ± 3.4 mm/yr in the west to approximately 3 mm/yr in the east. By estimating the seismic moment accumulation rates of the major seismic belts and comparing them with the historical earthquake catalog, we identify six seismic belts with significant seismic moment deficits. This indicates a potential risk for earthquakes exceeding magnitude 7, including the Kash fault, Keping fault, the Maidan fault zone, and the North Tianshan seismic belt. The Nalati seismic belt exhibits a relatively small seismic moment deficit, indicating the potential of earthquakes in the magnitude range of 6 to 7. In contrast, Qiulitage fault, the West Tianshan seismic belt and the Manas fold-and-thrust belt show a moment surplus, suggesting a low likelihood of strong earthquakes occurring in the near future. This study provides critical data and theoretical support for the prediction and risk assessment of seismic activity in Tianshan region.

SummaryWe used the multiple-frequency seismic tomography method to image the upper mantle beneath Northeast Brazil by processing P-wave broadband seismograms in six frequency passbands. The data comprised 87 896 relative traveltime residuals for P and PKIKP phases, simultaneously inverted to obtain 3D models of P-velocity anomalies. We conducted resolution tests using checkerboard patterns with horizontal dimensions of 312×312 km and 390×390 km. For the 312×312 km structures, we observed good horizontal recovery beneath the areas of the Borborema province and northern São Francisco craton at the depths of 136 and 226 km, with the areas of the São Luís craton and Parnaíba basin also presenting recovery, albeit with reduced amplitudes. For the 390×390 km structures we observed good horizontal and amplitude recovery throughout the entire study area. Our model was unable to recover the sharp vertical transitions between the anomalies. We present the preferred model for Parnaíba basin and Borborema province. The model shows a fragmented basement for the Parnaíba basin, with two strong high-velocity anomalies consistent with the Parnaíba and Granja blocks, and another slight high-velocity anomaly north of the basin, consistent with the São Luís craton. To the west of the basin, the Parnaíba block appears separate from another high-velocity anomaly associated with the Amazonian craton. A strong high-velocity anomaly south of the Borborema province is interpreted as the northern portion of the São Francisco craton. The São Francisco craton anomaly presents strong high-velocity anomalies, interpreted as a thickening of certain portions of the craton, separated by a weaker positive anomaly, interpreted as the Paramirim Aulacogen. The Borborema province is characterized by a low-velocity anomaly. The central and northeastern portions of this anomaly presented even lower velocities, which was interpreted as lateral flow in the asthenosphere, originating from the passage of a plume to the north of the province. A low-velocity anomaly located west of the Borborema province strikes roughly NE-SW and separates the São Francisco craton from the Parnaíba block and the Amazonian craton. It is interpreted as the Transbrasiliano Lineament. To test the capability of our data to resolve the limits of large-scale structures, we created four synthetic models simulating the presence of different cratonic nuclei. The models show good horizontal recovery, with the fourth model, based on our findings, presenting the best correlation between the real and recovered models. Seismicity in the study region is mainly correlated to low-velocity anomalies.

SummaryThe reliable classification of seismic events is crucial to precise seismic cataloging and robust hazard evaluation. Recent advances in deep learning have achieved great success in seismic event identification, leveraging their exceptional ability to automatically extract and recognize features. However, existing deep learning approaches to seismic classification rely exclusively on deterministic models, which cannot quantify epistemic uncertainty, preventing the estimation of prediction confidence that is critical for reliability evaluation. To address this issue, in this paper, we develop two uncertainty-aware deep learning models for earthquake (EQ) vs. explosion (EP) classification using the DiTing 2.0 artificial intelligence training dataset: a Bayesian convolutional neural network (BCNN) and a dropout-based CNN (DropCNN). We also implement a conventional deterministic CNN as a baseline model for comparative analysis. The experimental results demonstrate that both the BCNN and DropCNN can achieve a classification accuracy comparable to the conventional CNN, while providing additional uncertainty metrics for estimation confidence of prediction. Crucially, their uncertainty scores increase markedly in terms of encountering misclassifications or out-of-distribution samples compared to correct classifications, enabling automatic rejection of unreliable predictions based on the uncertainty threshold setting, triggering human verification or alternative discrimination methods. We then apply the trained models to analyze suspicious explosion events in the DiTing 2.0 dataset. The BCNN and DropCNN results exhibits strong agreement, consistently identifying 79 EP and two EQ events and flagging the remaining samples as uncertain classifications needing further verification. Our findings demonstrate that deep learning methods incorporating uncertainty estimation not only maintain a high accuracy in seismic event discrimination but also provide uncertainty estimation. This capability significantly enhances the model’s reliability and decision-making value in practical applications.

SummaryThe central portion of the 2019 ± 2 Ma Vredefort (South Africa) impact structure comprises a 40-50 km diameter central uplift of Archean basement rocks surrounded by a 15-20 km wide collar of late Archaean to early Proterozoic Witwatersrand Supergroup sedimentary and volcanic rocks. The collar is characterized by a ring of strongly negative (←5 500 nT) aeromagnetic anomalies surrounding much of the structure where the strata dip steeply to overturned. To better understand the origin of this magnetic feature, we undertook a ground survey along 20 transects (340 km) in the Vredefort structure using a three-axis fluxgate magnetometer mounted on a mountain bicycle. Upward continuation of our profiles to 150 m matches the aeromagnetic data in shape and amplitude. From the bicycle measurements, we pinpointed the rocks responsible for the extremely negative anomalies. Field observations and microfabric analyses of the rocks from six outcrops substantiated that the magnetic signal correlates with 10-100 m thick metamorphosed banded iron formations (BIFs) at the base of the supergroup as the main producer of the anomalies. Paleomagnetic samples collected from the rocks at the surface that produce the most intense anomalies (up to -22 000 nT) have extremely high natural remanent magnetization intensities (up to > 1000 A·m−1) likely arising from lightning strikes. Stepwise demagnetization and rock magnetic experiments establish a new protocol to distinguish samples that escaped remagnetization from lightning and possess the established 2.02 Ga paleodirection at Vredefort. From a suite of thermoremanent magnetization (TRM) experiments, the best estimate for the paleofield intensity at the time of impact was 52 μT, corresponding to an average remanence of 32.5 A·m−1. The results of the TRM experiments together with the paleodirection enabled us to successfully model the prominent negative anomalies in the metasediments only when accounting for the post-impact orientation of the BIFs. We interpret the strongly negative magnetic anomalies in the collar region as being formed directly after crater exhumation and uplift of the rocks. This interpretation implies that Bushveld-related metamorphism at 2.06 Ga created the up to mm-sized magnetite and garnet crystals in the BIFs, which resided at temperatures higher than the Curie temperature of magnetite (580°C) until the impact rapidly brought the BIFs close to the surface, where magnetite cooled to acquire a thermal remanence in the 2.02 Ga field.

SummaryMagnetic fabric analysis of dikes is a powerful technique when assessing magma transfer processes. This study presents an integrated analysis combining magnetic susceptibility and anisotropy of magnetic susceptibility (AMS), magnetic mineralogy, geochemistry, and new ⁴⁰Ar/³⁹Ar dating of dikes intruding formations ranging from the Lower Cretaceous to the Miocene on the island of Maio, in the Cabo Verde archipelago. We show that the dikes, dated at ≈ 9.2 Ma, intruding the younger Miocene Casas Velhas formation, display a Ti-rich titanomagnetite composition, higher whole-rock TiO2 content and very high magnetic anisotropy. They are clearly distinguished from the dikes, ranging in age from ≈ 9.3 to 11.3 Ma, intruding older formations, which show a predominantly Ti-poor titanomagnetite composition with multiple magnetic phases, lower whole-rock TiO2 concentration, higher range of magnetic susceptibilities and very low anisotropy. Magnetic fabric is predominantly normal with no significant imbrication relative to the dike margins. Numerical analysis of fabric shows a dominant coaxiality between the magnetic lineation and the preferred orientation of opaques and phenocrystals suggesting that magnetic lineation is, therefore, the proxy of the magmatic flow axis orientation. Based on the orientation of the magnetic fabric, we infer that magmatic flow within the studied dikes is predominantly vertical. The differences observed between the younger dikes and all other dikes may be related to magma sourced from distinct magma chambers. One, probably shallow, underneath the Casas Velhas fm in the southwest of the island, which would explain the very high values of magnetic anisotropy and the inferred vertical flow, and another located in a central position in the island, responsible for the dikes intruding the older formations. The location of such magma reservoirs and the dikes ages suggest a hypothetical migration with age of the magmatic sources that fed the dikes from the central part of the island to the southwest region. The magnetic and mineralogical heterogeneities of the dikes intruding older Lower Cretaceous formations may also be a result of a wider age range of the intrusions.

SummaryIn Distributed Acoustic Sensing (DAS), a fiber-optic cable is used as a distributed seismic sensor, with channels representing successive short sections of the fiber, spaced at defined intervals along the 1D fiber axis. Typically, the positions of these channels are assumed to be a line projection along the cable's position. In reality, a fiber-optic cable contains many fibers that are not perfectly straight and are thus longer than the cable itself. Consequently, the real channel positions may not correspond to a simple interpolation along the cable axis. Moreover, the precise cable coordinates are often sensitive information and may not be provided to the end user who uses the cable for sensing applications. On land, a tap test is usually carried out before the start of a DAS acquisition to determine the exact channel locations. DAS with marine horizontal cables has recently been used for various offshore applications, including seismic imaging. To avoid errors in the seismic image, a precise receiver location is required. In this paper, we propose a traveltime-based inversion workflow to determine a more accurate channel position on the seafloor. Moreover, we show that we can resolve an unknown time shift between the acquisition and the recording system, in addition to the fiber position.

SummaryTransform faults are one of the major tectonic plate boundaries offsetting the global mid-oceanic ridge system. The topographic features within these transform faults provide crucial evidence for tectono-magmatic processes and crustal accretion in transform fault zones. These interesting features include median ridges, which are major bathymetric anomalies found within both slow-slipping and fast-slipping transform faults, often associated with exposures of ultramafic rocks on the seafloor. To explain the origin of median ridges, previous studies have invoked multiple processes such as serpentinite diapirism, thermal uplift at ridge-transform intersections, or transpressive uplift induced by global plate reorganization, without any knowledge of the seismic structure. Here, we present results from 2D travel time tomography of downward-continued multi-channel seismic data along and across an ∼80 km long median ridge that lies within the eastern end of the slow-slipping (∼3.4 cm/yr) Chain transform fault in the equatorial Atlantic Ocean. The data were acquired during the 2018 ILAB-SPARC survey using a 6-km long streamer. Our high-resolution P-wave velocity model of the median ridge shows distinct high and low velocities ranging from 2.5 to 5 km/s within 500 m below the seafloor, on either side of the presently active strike-slip fault trace that cuts through the ridge. The low velocity on the eastern side of the ridge could be due to the presence of highly fractured basalt (with porosity in the range of 28 to 36 per cent) due to transform fault motion, whereas the high velocity on the western flank could be due to the presence of gabbro or highly serpentinised peridotite. The basaltic origin of the median ridge is supported by the observation of a seismic triplication event, which we call the T-event. The depth at which the T-event maps is shallow (200–500 m below seafloor) in high-velocity regions and deeper (600–1400 m) in low-velocity regions. We also find that the currently active strike-slip fault has been active since at least 0.26 Ma and has sliced the ridge. We image low-velocity pockets at the northern and southern limits of the median ridge that could represent the expression of the currently less active strike-slip faults.

SummaryThe northeastern Tibetan Plateau is bounded by the left-lateral Altyn Tagh and Haiyuan faults. How crustal motion along these fault systems transitions to crustal shortening and uplift is key for deciphering the geodynamic link between the escape tectonics and the growth of the Tibetan Plateau. Here, we use the PS-InSAR observations, combined with GNSS and leveling data, to obtain a high-resolution 3D model of the present-day crustal motion in the northeastern Tibetan Plateau. The resolved deformation field covers the entire northeastern Tibetan Plateau with a spatial resolution of approximately 0.01° × 0.01 °. Our analysis of slip rates and strain partitioning reveals that crustal motion along the Altyn Tagh fault gradually diminishes eastward and is absorbed by thrusting and uplift in the Qilianshan orogenic belt within the plateau. A similar tectonic transition occurs between the Haiyuan fault and the Liupanshan orogen on the eastern margin of the plateau. Some of the eastward crustal motion is accommodated by the younger Xiangshan-Tianjingshan fault system to the north of the Haiyuan fault, indicating the ongoing northward expansion of the Tibetan Plateau. Our results align with geological evidence of crustal deformation in the past few million years, highlighting the continuing tectonic transition from eastward crustal motion along the left-lateral strike-slip faults to the growth of the Tibetan Plateau.

SummaryThe continuous NE-SW compression due to the Indo-Asian collision creates active and complex deformation in the northeastern (NE) Tibetan plateau. How the lithosphere of the NE Tibetan plateau deforms both vertically and laterally in response to the ongoing collision is still a question. Further investigations with refined lithospheric structure are required. Here we present a high-resolution radially anisotropic model of the lithosphere beneath the NE Tibetan plateau and surrounding areas that was constrained by the joint analysis of Rayleigh and Love wave dispersion at periods from 6 to 100 s using the methods of ambient noise cross-correlation for short periods and earthquake two-plane-wave for long periods. Results show relatively small regions of significant slow shear wave velocity and positive radial anisotropy (Vsh > Vsv) in the middle crust beneath the Qilian orogenic belt, suggesting the existence of partial melting and probably limited channel flow. Considering the variable lateral strength of shear wave velocity and radial anisotropy in the middle crust with large parts mechanically strong enough to pass the strain, vertical coherent lithospheric deformation could still work in the Qilian orogenic belt. Extensive low shear wave velocity anomaly in the uppermost mantle extend from the Qilian orogenic belt northward to the Alxa block and eastward to the southwestern Ordos block, implying a hot and weak mantle lithosphere. The observed negative radial anisotropy (Vsh < Vsv) in such warm mantle lithosphere beneath the Qilian orogenic belt, Alxa block and the southwestern Ordos block is ascribed to vertical deformation fabrics arising from the convergence between Indian and Asian plates. These observations imply that the lithosphere of Qilian orogenic belt, Alxa block and southwestern Ordos block deform coherently and the NE Tibetan plateau is expanding towards Alxa block and the southwestern Ordos block.

SummaryWe develop a new method for estimating the autocorrelation function (ACF) from segmented data with the assumption of stochastic stationarity. The ACF of a signal is represented as the summation of the cross terms of sub-segments of arbitrary length. To successfully remove undesired transients in the data, this method introduces a correction for the amplitude bias associated with the removal of sub-segments, based on the comparison between the expected stationary signal and the measured signal. The method reconstructs and accesses later lag times, provides finer frequency resolution, obtains a better signal-to-noise ratio, which enables the extraction of detailed temporal or spectral structures from noisy data sets. As an application, we successfully retrieved a spectrum of the Earth’s seismic hum on the vertical component with fine frequency resolution and compared it to synthetic autocorrelation for spatially isotropic and homogeneous excitation by random shear traction at the ocean bottom and random pressure at the Earth’s surface. Although both models can explain the observed fundamental spheroidal modes, shear traction is better at explaining the observed overtones above 3 mHz. From 2 to 3 mHz, the pressure source also contributes to the excitation of the overtones, and the shear traction becomes dominant again below 2 mHz. This new method is anticipated to be effective in extracting valuable information from rare records within the context of extraterrestrial seismology.

SummaryNortheast China, with its complicated regional tectonic evolution, situated within the eastern Central Asian Orogenic Belt, is a key region for understanding lithospheric deformation and mantle dynamics. However, the ongoing debate surrounding its lithospheric structure and evolutionary processes remains, largely attributed to data limitations and methodological constraints. In this study, we integrate topography, geoid height, surface heat flow, and Rayleigh wave phase velocity dispersion curves to conduct a detailed imaging of the lithospheric thermal and compositional structure in Northeast China. We find a significant east-west gradient in lithospheric thickness, ranging from approximately 60 km in the east to 140 km in the west, and a compositional transition in the lithospheric mantle from fertile peridotite in the east to refractory peridotite in the west. By integrating analyses of upper mantle anisotropy and the spatiotemporal distribution of Mesozoic basalts, we argue that the lithospheric delamination and mantle upwelling may have combined to cause the lithospheric thinning in the region. This study highlights the significance of joint inversion of multiple datasets and integrated multidisciplinary analysis.

SummaryTraditional methods for measuring geopotential difference using optical fiber frequency transfer or satellite-based time-frequency transfer, based on general relativity, require the use of two clocks and the calibration of these clocks. Here we present a simplified clock transportation experiment using a single hydrogen clock to measure the geopotential difference between two time-frequency stations, separated by 129 km with a height difference of 1,245 m, by GNSS precise point positioning time-frequency transfer. Taking the reference clock of International GNSS Service (IGS) time as a ‘bridge’, we extract the gravity frequency shift between the two stations by comparing the fractional frequency differences between the hydrogen clock and the ‘bridge’ before and after clock transportation. The determined geopotential difference between the two stations is 12075.9$\pm $118.5 m²/s², which closely aligns with the value computed by the EIGEN-6C4 global gravity field model, with a difference of -78.7 m²/s². These results validate the feasibility of geopotential difference measurements with a single clock and highlight several advantages compared to the dual-clock method: elimination of inter-clock calibration, low operational complexity and equipment cost, high data utilization efficiency but similar precision of geopotential difference measurement. Furthermore, this method can be extended to other similar techniques to measure geopotential differences, provided that they enable users to connect to a stable time-frequency reference.

SummaryShale gas extraction could produce underground stress perturbation and local seismicity, which could put a threat human casualty. The Weiyuan area in the Sichuan province, China, underwent massive gas production and a significant increase of earthquake since 2015. In this study, we focus on human-induced subsurface hydrofracturing, calculate cumulative underground Coulomb-stress changes using a 3D numerical model, and probe the main cause of recent seismic activity in the Weiyuan area based on continuous regional stress/displacement loading. The simulation reveals the regional extent of positive Coulomb stress change with fractures matches the distribution of the moderate and micro seismicity in the past ten years. Background regional tectonic stress in the vicinity of the active fault likely resulted in earthquake preparation within and around the active faults; hydraulic fracturing changes mainly the displacement and stress pattern in the vicinity of the fracturing wells, and enhanced fracturing intensity (fracture volume-to-model volume ratio (θ), causes more obvious difference; faults may be locked prior to fracturing, and even small fracturing intensity may trigger the earthquakes near faults and fracturing wells; the seismic risk will be significantly increased near the two faults and fracturing wells in the next 50 years.

SummaryIn this study, we compare the usability of a simplified microtremor-based empirical method and a conventional microtremor method based on an inversion analysis of a subsurface velocity structure model for constructing a map of average S-wave velocity (AVS) values. In the simplified (empirical) method, the phase velocities of Rayleigh waves, which can be obtained by processing a microtremor array, at wavelengths of 13, 25, and 40 m are regarded as AVS values from the ground surface to depths of 10, 20, and 30 m (${\overline {Vs} }_{10},{\overline {Vs} }_{20},\ {\rm{and}}\ {\overline {Vs} }_{30}$), respectively. Microtremor array surveys were conducted at 173 observation points within a 15 km × 17 km area east of Aso caldera, Kyushu, Japan (target area). AVS values are obtained by applying the empirical method to the phase velocities obtained at each observation point. The AVS values at an observation point (located near the centre of the target area) with velocity logging data are verified by a comparison with those based on the velocity logging data (i.e. overestimations by 6 per cent at maximum). It is found that for the entire target area, the spatial distribution of the obtained AVS values is consistent with the geological distribution. The AVS values within areas of the Aso-3 ignimbrite are 30–40 per cent larger than those within areas of thick soil and tephra on the strongly consolidated Aso-4 ignimbrite. In addition, the AVS values of the Aso-3 deposits are more than 10 per cent larger than those of the Aso-4 deposits and about 10 per cent smaller than those of geological units older than the Aso-3 deposits. We also apply a conventional (i.e. inversion) method to the phase velocity data at each observation point to obtain a one-dimensional S-wave velocity (Vs) structure model from which we deduce AVS values. The deduced AVS values at the velocity logging point are underestimated by -8 per cent, with differences from the AVS values obtained using the empirical method reaching 13 per cent. The average systematic difference between the two methods is 15 per cent, as determined from a statistical analysis. None the less, a strong correlation is found between the methods, with an average correlation coefficient of 0.94, with no evidence showing that either method is more accurate. The empirical method can be used to construct an AVS map if overestimation is carefully considered. This analysis also reveals that the average maximum survey depths of the one-dimensional Vs structures based on the inversion method are only 23±10 m, making them often insufficient to map ${\overline {Vs} }_{20}$ and ${\overline {Vs} }_{30}$ (the ratios of the available to total numbers of data points are only 60 and 21 per cent, respectively). In contrast, the empirical method can determine ${\overline {Vs} }_{10},{\overline {Vs} }_{20},\ \ {\rm{and}}\ {\overline {Vs} }_{30}$ at more than 80 per cent of all sites. The construction of AVS maps using the empirical method is effective in terms of the simplicity and reliability of planning, observational efficiency, and simplicity of data processing, which support a practical and objective approach to seismic assessments.

SummaryIn this study we obtain 35 903 high-quality P-wave receiver functions from 1737 teleseismic events recorded at 120 dense broadband TanluArray temporary stations deployed in and around the Tanlu fault zone (TLFZ). After station azimuth and sediment correction are made, a detailed Moho depth distribution is obtained by CCP stacking. Our results show a sharp change in the Moho depth across the TLFZ from the west to east, which well corresponds to the surface geological structure. The deepest Moho (38.0 ∼ 40.0 km) occurs beneath the Dabie orogenic belt and the Sulu orogenic belt. The Moho beneath the Luxi uplift, Jiangnan orogenic belt and Jiaodong uplift is deeper (36.0 ∼ 37.0 km), whereas the Subei basin and the southern basin of the South Yellow Sea have a shallow Moho (28.0 ∼ 30.0 km). There is an obvious Moho uplift near Weifang, which corresponds to the Changle ancient volcano on the surface and may be a channel for upwelling of hot mantle material. The Moho is unclear under the fault zone near Tancheng, which is speculated to be a channel for upwelling of hot mantle material. It may be related to upwelling of hot and wet flows in the big mantle wedge above the subducted Pacific slab that is stagnant in the mantle transition zone beneath East Asia, which is a possible cause of the 1668 M8.5 Tancheng earthquake.

SummaryThe Limpopo transform margin offshore southern Mozambique results from the separation of Gondwana along the East Africa continental margin. Over the last three decades, more than thirty different reconstruction models have been proposed, sometimes contradicting each other. Here, we present results from the travel-time tomography of wide-angle seismic data acquired during the second China-Mozambique Joint Cruise, allowing the interpretation of the crustal structure and magmatism in the Limpopo Corridor and the Mozambique Basin. Using these results, we determine the extent of the Continent Ocean Transition and the location of the Continent Ocean Boundary on the southern Mozambique margin. The seismic profile is 442-km long, extending from the eastern part of the North Natal Valley in the west and crossing the Limpopo Corridor and the Mozambique Basin to the east. Based on the tomographic velocity model, we delineated three distinct domains from west to east along the profile: (1) a western transitional domain with anomalous or mixed crust, bounded by the Mozambique Fracture Zone to the east, where the crust gradually thins eastward from ∼14 km at distance 45 km to ∼10.8 km at distance 140 km; (2) a domain of thickened oceanic crust resulting from enhanced magmatism, where the crust thins eastward, from ∼10.8 km to ∼8.5 km over ∼100 km distance; and (3) an eastern domain of normal oceanic crust, where the average crustal thickness is ∼8 km. We suggest that (1) the western transitional domain roughly corresponds to the Limpopo Corridor and is of continental crustal origin but was affected and modified by strike-slip motion and magmatic activity, resulting in anomalous or mixed crust. The eastern Continent Ocean Boundary of the Limpopo Margin is close to the Mozambique Fracture Zone; (2) The thickened oceanic domain thins eastward, and the crustal velocity and thickness change dramatically compared to the oceanic domain. This domain seems to have strongly interacted and contaminated by the Limpopo Corridor during the opening of Mozambique Basin and seafloor spreading; (3) The eastern oceanic domain shows a relatively uniform oceanic crust of ∼8 km and high velocity up to 7.4 km/s in the lower crust, suggestive of a hotter mantle that produces more MgO-rich melts probably due to the influence of a thermal mantle anomaly.

SummaryThe downward continuation of gravity field can provide valuable information for 3-D gravity-field modeling, shallow-layer geological interpretation, source depth estimation, and so on. However, downward continuation is ill-posed, and traditional approaches often suffer from computational instability, poor noise resistance, and limited continuation depth, making it a longstanding challenge in gravity data processing. We present a new approach for fast, stable and large-depth downward continuation of gravity anomalies by using frequency-domain 3-D imaging. First, we utilize the frequency-domain 3-D imaging approach to invert the gravity anomalies at the original observational plane to quickly obtain the equivalent density model in the subsurface. Then, we apply the optimized strategy of frequency-domain 3-D forward calculation on the equivalent density model to rapidly obtain high-precision gravity anomalies at the downward-continuation plane. The synthetic data tests prove the effectiveness of our approach, and demonstrate that our approach enables fast, stable, robust noise resistance and large-depth downward continuation of large-scale gravity anomalies data, and has superior performance compared to the traditional regularized filtering approach and spatial-domain equivalent-source approach. The real data test of the free-air gravity anomalies data in the central South China Sea also verifies the fast, stable and reliable downward continuations of large depths by our approach. The 3-D gravity-field model built by our approach will provide significant support for the tectonic studies and resource exploration in this area.

SummaryModelling and inversion of controlled-source electromagnetic data requires elaborate numerical tools. The major challenge is the high computational cost of computing solutions to numerous forward problems (for the forward responses as well as the sensitivity matrix). Forward modelling is accomplished using either a direct or an iterative solver. Current modelling suites predominantly employ direct solution methods in the forward operator since multiple solutions are easily accessible using inexpensive and quick forward-backward substitution after an initial resource-demanding matrix factorisation step. Iterative techniques, on the other hand, require little resources for single forward solutions, and are yet very time consuming if solutions for many right-hand sides are to be computed. Evaluations of different solution techniques for modelling and inverse problems are only sparsely investigated. In light of this, we integrated an iterative solver as alternative in the forward and and inversion operators of the open-source software custEM and pyGIMLi. In particular, we implemented a two-level iterative scheme where the outer solver employs a generalised conjugate residual algorithm preconditioned with a highly efficient block-based preconditioner for square blocks. The inner-level solver is either of the same type as the outer solver, but preconditioned with the auxiliary-space Maxwell preconditioner, or may alternatively be a direct solver. In this paper, we evaluate the described iterative forward operator for forward modelling tasks for the Marlim R3D model for a single as well as numerous right-hand side vectors and compare the performance to the direct solver MUMPS. We further investigate the solver’s applicability on small and medium-sized computing platforms. We then examine the iterative solver for inversions of synthetic land-based and semi-airborne data in terms of computational requirements. Our results demonstrate that forward modelling tasks are best performed using an iterative approach for single source problems. Moreover, simulations of large and complex problems are accessible on even on small computing platforms such as laptops in very reasonable time. For inversions, the iterative forward operator, in particular the mixed iterative-direct-based one, performs equally well in terms of time as the direct one while reducing the memory demands for the computations of the forward responses and the data sensitivities.

SummaryIn a region of complex geology, we examine the influence of spatial resolution of conductivity models on Geomagnetically Induced Currents (GICs) estimations. We focus on the southern region of Portugal mainland, for which magnetotelluric (MT) sounding measurements have been obtained with lower noise from human activity. Using two conductivity models inverted from sets of MT soundings with different sampling distance, we look for an interpretation of the differences in GIC estimations at substation grounding resistances. We make use of two different proxies, the Local Effective Field (LEF) and the Regional Electromotive Source (RES), built from the electric induced field at each substation site and the sum of electromotive forces along all transmission lines connected to that substation, respectively. We compare different time signals associated to GICs using a parameter that combines Pearson correlation and linear regression slope, the Correlation Regression Coefficient (CRC). Our main conclusion is that spatially detailed information on lateral heterogeneities of the conductivity associated to complex geology is crucial for a rigorous assessment of GIC hazard, leading to relative differences in GIC standard deviation and in GIC peak values that can amount to more than 100% in certain cases. Additionally, using LEF and RES, we emphasise the non-locality of GIC drivers and bring new input concerning the choice of proxies used to monitor and forecast this kind of hazard.

SummaryThe Scandinavian Peninsula and its vicinity comprise highly tectonically diverse blocks, including the Baltic Shield, the continental margin, and the North Sea Basin. The crustal rheology is a critical constraint to understanding the tectonic evolution in this region. Based on 19 416 Lg waveforms from 233 earthquakes and 560 broadband digital stations, using an inversion method combining both single- and two-station ray paths, we constructed a broadband (0.05 and 10.0 Hz) Lg wave attenuation model in the study region, with the resolution approaches to 110 km (∼1°) or higher in areas with dense ray path coverages. The QLg distributions correlate well with regional geological features. The Baltic Shield exhibits the highest QLg, consistent with its thick Precambrian crust and high rheological rigidity developed through Archean Svecofennian orogeny. In contrast, passive margins with crustal thinning, magmatic modification, and thick sedimentary sequences exhibit strong attenuation, reflecting a reduction in rheological strength resulting from interactions with mantle plumes and extensional tectonics. The North Sea Basin exhibits the lowest QLg values and the presence of hydrocarbon-bearing sediments. The extremely high QLg distribution reveals the ancient cratonic core of the Baltic Shield, particularly in areas where the surface rock dating sample cannot be collected due to seawater coverage.