JGR–Solid Earth

The identification of the Grenvillian-age ophiolite suites in the Yangtze block in recent years suggested that the northern Yangtze subblock (NYB) and the southern Yangtze subblock (SYB) were once separated by an ocean in late Mesoproterozoic. Although some paleogeographic models advocated the Pre-Grenvillian connections between the south China blocks and Laurentia, none of them has been paleomagnetically tested. Here we report the new paleomagnetic results obtained from the ∼1,270 Ma purplish-red muddy dolomite beds of the Luanshigou Formation, Shennongjia Group, NYB, providing new constraints for reconstructing the paleogeographic positions of the south China blocks in late Mesoproterozoic. A total of 447 samples underwent stepwise thermal demagnetization. Two components were identified. The low-temperature component is interpreted as the recent viscous remnant magnetization. The high-temperature component was obtained from 64 samples below 580°C and from 177 samples below 690°C, directed northeast-up or southwest-down, antipodally, positioning the paleomagnetic pole at 18.5°S, 74.4°E (dm/dp = 2.5/1.6°). Rock magnetic results demonstrate that the magnetic carriers in purplish-red dolomite and pale-pink dolomite are predominated by hematite and magnetite, respectively. The data quality is supported by an inverse baked contact test, a B-class reversal test, and the paleomagnetic pole is distinct from any younger poles of the region. Based on the paleomagnetic results, aided by geological evidence, we propose a reconstruction in which the NYB was juxtaposed to southwestern Laurentia in the late Mesoproterozoic and suggest that the late Mesoproterozoic Miaowan-Shimian ophiolite zone in the Yangtze block was likely an extension of the Grenville belt of Laurentia.

In the present study, we use broadband seismic data recorded by 190 stations of the EarthScope program's Transportable Array to construct a 3-D shear wave velocity model for the upper 180 km using a non-linear Bayesian Monte-Carlo joint inversion of receiver functions (RFs) and Rayleigh wave dispersion curves. Ambient noise and teleseismic data are used for obtaining Rayleigh wave phase velocity dispersion curves. A resonance removal filtering technique is applied to the RFs contaminated by reverberations from the thick sedimentary layers that cover most of the region. Our observations of the higher crustal shear velocities (∼3.40 km/s) beneath the Sabine Block (SB), along with the estimated relatively thicker crust (∼34.0 km) and lower crustal V p/V s estimates (∼1.80) in comparison with the rest of the Gulf Coastal Plain (GCP) (∼3.10 km/s for crustal shear velocities, ∼29.0 km for crustal thickness, and ∼1.90 for crustal V p/V s estimates), indicating that this crustal block has different crustal properties from the surrounding coastal plain regions. The southern Ouachita Mountains have a thin crust (∼30.0 km) and low mean crustal V p/V s value (∼1.73), suggesting that lower crustal delamination has occurred in this region. Low velocities in the upper mantle beneath most of the GCP are interpreted as a combined result of thin lithosphere, higher-than-normal temperatures, and possibly compositional variations.

The thermal conductivity of bridgmanite, the primary constituent of the Earth's lower mantle, has been investigated using diamond anvil cells at pressures up to 85 GPa and temperatures up to 3,100 K. We report the results of time-domain optical laser flash heating and X-ray Free Electron Laser heating experiments from a variety of bridgmanite samples with different Al and Fe contents. The results demonstrate that Fe or Fe,Al incorporation in bridgmanite reduces thermal conductivity by about 50% in comparison to end-member MgSiO3 at the pressure-temperature conditions of Earth's lower mantle. The effect of temperature on the thermal conductivity at 28–60 GPa is moderate, well described as k=k300(300/T)a ${k={k}_{300}(300/T)}^{a}$, where a is 0.2–0.5. The results yield thermal conductivity of 7.5–15 W/(m × K) in the thermal boundary layer of the lowermost mantle composed of Fe,Al-bearing bridgmanite.

Recent observations beneath Kanto, Japan have shown that seismic activity and seismic attenuation within the overlying continental plate change with time due to drainage caused by slow-slip events (SSEs) along the upper boundary of the Philippine Sea plate. However, associated changes in rock properties have not been investigated. In this study, we estimate frequency-dependent shear-wave anisotropy to provide a detailed insight into the structural change associated with drainage. We perform shear-wave splitting analysis in frequency ranges of 1–4, 2–6, and 4–8 Hz for 306 earthquakes that occur during September 2009–August 2021 and recorded at the Metropolitan Seismic Observation network. Obtained time differences between fast and slow S waves (delay time) range from almost zero to 0.16–0.18 s, exhibiting spatio-temporal variation and frequency dependence. The fast S-wave polarization directions are generally consistent with the direction of the maximum horizontal compressional axis in the study region, which suggests that the observed anisotropy is probably caused by the NE–SW-oriented fractures developed under the regional stress field. The temporal variation in delay times is correlated with SSEs activity with a lag time of 0.0–0.1 year. Furthermore, comparisons between observed frequency-dependent delay times and numerical calculation of fracture-induced anisotropy suggest that the average fracture radius is almost constant (0.30–0.35 m) over time but fracture density temporally varies from 0.025 to 0.035. We infer that the fracture density is probably enhanced by opening of the NE–SW-oriented fractures during the upward migration of fluids that are expelled from the plate interface.

To gain a deeper understanding of the extensive and varied lithospheric deformations beneath northern Oman, we examine seismic anisotropy in this region using splitting analysis of teleseismic shear wave data. Our study utilizes data from a dense network consisting of 13 permanent and 45 temporary seismic stations, which were operational for approximately 2.5 years starting from October 2013. By examining the azimuthal distribution of shear wave splitting (SWS) parameters, we were able to divide the study area into three sub-regions. The stations located to the west of the Hawasina window exhibit relatively azimuthally invariant SWS parameters suggesting a single anisotropic layer. On the other hand, most of the stations located in the central and eastern regions display variations versus back-azimuth, indicating the potential presence of depth-dependent anisotropy. The General NW-SE trend of the Fast Polarization Directions (FPDs) of the one-layer anisotropy in the west and FPDs of the upper layers in the east is concordant with the strike of the structures resulting from the collision between the continental and oceanic plates. A clear contrast in SWS parameters is observed in the Semail Gap Fault Zone (SGFZ), suggesting that the SGFZ can be a lithospheric-scale structure that hampers the intrusion of mafic magma from the southeast. Furthermore, the FPDs of the lower layer in the east exhibit an NE-SW trend, which may be indicative of the large-scale mantle flow resulting from the present-day plate motion.

We developed an enhanced Kalman-based approach to quantify abrupt changes and significant non-linearity in vertical land motion (VLM) along the coast of Chile and the Antarctic Peninsula using a combination of multi-mission satellite altimetry (ALT), tide gauge (TG), and GPS data starting from the early 1990s. The data reveal the spatial variability of co-seismic and post-seismic subsidence at TGs along the Chilean subduction zone in response to the Mw8.8 Maule 2010, Mw8.1 Iquique 2014, and Mw8.3 Illapel 2015 earthquakes that are not retrievable from the interpolation of sparse GPS observations across space and time. In the Antarctic Peninsula, where continuous GPS data do not commence until ∼1998, the approach provides new insight into the ∼2002 change in VLM at the TGs of +5.3 ± 2.2 mm/yr (Palmer) and +3.5 ± 2.8 mm/yr (Vernadsky) due to the onset of ice-mass loss following the Larsen-B Ice Shelf breakup. We used these data to constrain viscoelastic Earth model parameters for the northern Antarctic Peninsula, obtaining a preferred lithosphere thickness of 115 km and upper mantle viscosity of 0.9 × 1018 Pa s. Our estimates of regionally-correlated ALT systematic errors are small, typically between ∼±0.5–2.5 mm/yr over single-mission time scales. These are consistent with competing orbit differences and the relative errors apparent in ALT crossovers. This study demonstrates that, with careful tuning, the ALT-TG technique can provide improved temporal and spatial sampling of VLM, yielding new constraints on geodynamic models and assisting sea-level change studies in otherwise data sparse regions and periods.

Unmixing of remanent magnetization curves, either isothermal remanent magnetization (IRM) or backfield IRM, is widely used in rock magnetic and environmental magnetic studies to discriminate between magnetic coercivity components of different origins. However, the wide range of physical properties of natural magnetic particles gives rise to an ambiguous interpretation of these components. To reduce this ambiguity and provide a straightforward interpretation of coercivity components in terms of domain state, interactions, and constituent magnetic phases, we combined backfield IRM unmixing with unmixing of nonlinear Preisach maps for two typical mid-latitude northern hemisphere loess-paleosol sequences. Both backfield IRM and nonlinear Preisach maps unmixing are based on the same non-parametric algorithm, and provide similar endmembers (EMs) in the two sections studied. The first EM (EM1) has a low median coercivity (∼21 mT) and is a non-interacting single domain (SD) magnetite/maghemite of pedogenic origin. The second EM (EM2) has a moderate median coercivity (∼60 mT) and is a mixture of pseudo-single domain/multidomain, SD magnetite/maghemite and non-interacting SD hematite, all of eolian origin. The same EM1 found in both sections suggests that this component's grain size and coercivity are independent of pedogenesis intensity. The same EM2 indicates that a similar magnetic population is being transported and deposited, irrespective of the dust source area and loess granulometry. The approach outlined here provides strong evidence that non-parametric backfield IRM unmixing isolates physically realistic EMs. Unmixing nonlinear Preisach maps elucidates these EMs in terms of domain states and their constituent magnetic phases.

Tsunamis propagate over long distances and can cause widespread damage even after crossing ocean basins. Prediction of tsunamis in distant areas based on observations near their sources is critical to mitigating damage. In recent years, the accuracy of numerical models of trans-oceanic tsunami propagation has improved significantly due to the incorporation of effects such as the solid earth response to tsunami loading and wave dispersion. However, these models are computationally expensive and have not been fully utilized for real-time prediction. Here, we derive the adjoint operator for the linear set of equations describing deep-ocean tsunami propagation and show how a pre-computed database of adjoint states can achieve rapid synthesis of tsunami waveforms at target sites from nonpoint arbitrary tsunami sources. The adjoint synthesis method allows for an exhaustive parameter search for tsunami source estimation. A method for simultaneous inversion of fault geometry and slip distribution using adjoint synthesis with Sequential Monte Carlo method was proposed and applied to the 2012 Haida Gwaii earthquake tsunami. The influence of model accuracy and the amount of observed data on the estimation of tsunami sources and waveforms was examined. It was found that with a highly accurate propagation model, using only a limited amount of observed data produced source and waveform estimates very similar to the final models obtained with much larger data sets. The final inferred fault model involved megathrust slip distributed between the Haida Gwaii trench and the Queen Charlotte fault. The proposed method can also quantify the uncertainty of the waveform forecasts.

To help meet emission standards, hydrogen sulfide (H2S) from geothermal production may be injected back into the subsurface, where basalt offers, in theory, the capacity to mineralize H2S into pyrite. Ensuring the viability of this pollution mitigation technology requires information on how much H2S is mineralized, at what rate and where. To date, monitoring efforts of field-scale H2S reinjection have mostly occurred via mass balance calculations, typically capturing less than 5% of the injected fluid. While these studies, along with laboratory experiments and geochemical models, conclude effective H2S mineralization, their extrapolation to quantify mineralization and its persistence over time leads to considerable uncertainty. Here, a geophysical methodology, using time-domain induced polarization (TDIP) logging in two of the injection wells (NN3 and NN4), is developed as a complementary tool to follow the fate of H2S re-injected at Nesjavellir geothermal site (Iceland). Results show a strong chargeability increase at +40 days, interpreted as precipitation of up to 2 vol.% based on laboratory relationships. A uniform increase is observed along NN4, whereas it is localized below 450 m in NN3. Changes are more pronounced with larger electrode spacing, indicating that pyrite precipitation takes place away from the wells. Furthermore, a chargeability decrease is observed at later monitoring rounds in both wells, suggesting that pyrite is either passivated or re-dissolved after precipitating. These results highlight that a sequence of overlapping reactive processes (pyrite precipitation, passivation, pore clogging and possibly pyrite re-dissolution) results from H2S injection and that TDIP monitoring is sensitive to this sequence.

We investigated the effects of the propagation path and site amplification of shallow tremors along the Nankai Trough. Using far-field S-wave propagation from intraslab earthquake data, the amplification factors at the DONET stations were 5–40 times against an inland outcrop rock site. Thick (∼5 km) sedimentary layers with V S of 0.6–2 km/s beneath DONET stations have been confirmed by seismological studies. To investigate the effects of thick sedimentary layers, we synthesized seismograms of shallow tremors and intraslab earthquakes at seafloor stations. The ratios of the maximum amplitudes from the synthetic intraslab seismograms between models with and without thick sedimentary layers were 1–2. This means that thin lower-velocity (<0.6 km/s) sediments just below the stations primarily control the estimated large amplifications. Conversely, at near-source (≤20 km) distances, 1-order amplifications of seismic energies for a shallow tremor source can occur due to thick sedimentary layers. Multiple S-wave reflections between the seafloor and plate interface are contaminated in tremor envelopes; consequently, seismic energy and duration are overestimated. If a shallow tremor occurs within underthrust sediments, the overestimation becomes stronger because of the invalid rigidity assumptions around the source region. After 1-order corrections of seismic energies of shallow tremors along the Nankai Trough, the scaled energies of seismic slow earthquakes were 10−10–10−9 irrespective of the region and source depth. Hence, the physical mechanisms governing seismic slow earthquakes can be the same, irrespective of the region and source depth.

The mid-lithosphere discontinuity (MLD), identified by a sharp velocity drop at ∼70–100 km depths within the cratonic lithosphere is key to comprehending the chemical composition and thermal structure of the cratonic lithosphere. The MLD is widely accepted to be caused by composition anomalies, such as hydrous minerals, which show low velocities and high electrical conductivities. However, noticeable high-electrical conductivity anomalies have not been detected in the most cratonic lithosphere. Dolomite has an electrical conductivity similar to olivine and can be originated by carbonatitic melts trapped at ∼80–140 km depths. Here we investigated the elasticity of dolomite under mantle conditions using ab initio calculations and found dolomite exhibits significantly lower velocities than the primary minerals in the lithospheric mantle. Therefore, the dolomite enrichment might provide a good explanation for the observed velocity drop of the MLD in cratonic regions where no high-conductivity anomaly has been detected, such as the northern Slave craton.

We use seismic ambient noise recorded by dense ocean bottom nodes (OBNs) in the Gorgon gas field, Western Australia, to compute time-lapse seafloor models of shear-wave velocity. The extracted hourly cross-correlation (CC) functions in the frequency band 0.1–1 Hz contain mainly Scholte waves with very high signal-to-noise ratio. We observe temporal velocity variations (dv/v) at the order of 0.1% with a peak velocity change of 0.8% averaged from all station pairs, from the conventional time-lapse analysis with the assumption of a spatially homogeneous dv/v. With a high-resolution reference (baseline) model from full waveform inversion of Scholte waves, we present an elastic wave equation based double-difference inversion (EW-DD) method, using arrival time differences between the reference and time-lapsed Scholte waves, for mapping temporally varying dv/v in the heterogeneous subsurface. The time-lapse velocity models reveal increasing/decreasing patterns of shear-wave velocity in agreement with those from the conventional analysis. The velocity variation exhibits a ∼24-hr cycling pattern, which appears to be inversely correlated with the diurnal variations in sea level height, possibly associated with dilatant effects for porous, low-velocity shallow seafloor and rising pore pressure with higher sea level. This study demonstrates the feasibility of using dense passive seismic surveys and wave-equation time-lapse inversion for quantitative monitoring of subsurface property changes in the horizontal and depth domain.

Antigorite is the high-temperature member of the serpentine group minerals and is broadly considered a primary carrier of water in the subducting oceanic lithosphere. It has a wavy crystal structure along its a-axis and several polysomes with different m-values (m = 13–24) have been identified in nature. The m-value is defined as the number of tetrahedra in one wavelength and is controlled by the misfit between the octahedral and tetrahedral layers. The degree of misfit primarily depends on the volumes of the MgO6 octehedra and SiO4 tetrahedra within the layers, which vary as a function of pressure and temperature. However, it is not well understood which m-values of antigorite are stable at different pressure and temperature conditions. To investigate the pressure dependence of the stability of different m-values in antigorite, we performed first-principles calculations for several polysomes (m = 14–19) at high pressure from 0 to 14 GPa and compared their enthalpies at static 0 K. We found that although the energy differences between polysomes are small, polysomes with larger m-values are more stable at ambient pressure, while polysomes with smaller m-values are more stable at elevated pressures. This suggests that the structure of antigorite in the oceanic lithosphere subducting into the deep Earth may gradually evolve into a different polysome structure than the antigorite samples observed at ambient or near-surface pressure conditions. These changes in the m-values are accompanied by a minor dehydration reaction. By modulating the available amount of free water in the system, antigorite polysomatism may influence the distribution of intermediate-depth seismicity, such as the observance of double seismic zones.

Tectonic plate motions feed the earthquake cycle—a process whereby stress along crustal faults slowly increases over decade– or century–long periods, to then suddenly drop during earthquakes. Steadiness of plate motions during such cycles has long been a central tenet in models of earthquake genesis and of faults seismic potential, and can be tested against measurements of contemporary plate motions available from Global Navigation Satellite Systems (GNSS). Here we present analyses of GNSS data from Central and Northern Italy that illuminate the motion of the Adria microplate over a period of 6 years preceding the M W 6.3, 6 April 2009 L’Aquila (Italy) earthquake. We show that the motion of the whole Adria microplate changed before the 2009 earthquake, and slowed down by around 20%. We demonstrate with quantitative models that the torque required upon Adria in order to drive such a kinematic change is consistent with what is imparted to Adria by temporal stress variations occurring during the late interseismic phase of the 2009 L’Aquila earthquake cycle. The inference that plate motions can be influenced by, and thus sensitive to, earthquake cycles offers an additional perspective to assessing the seismic potential of tectonic margins.

Scientists have observed the surface expression of creep events along the San Andreas Fault since the 1960s. However, the evolution of slip at depth has been examined relatively little. So here we probe that deep slip by analyzing strain observations just before and during hours- to day-long creep events at the northern end of the creeping section of the San Andreas Fault. We identify 71 strain offsets that are likely produced by few-hour bursts of slip at depth. Then, we grid search to determine the location, depth, and magnitude of these slip bursts. We find that the slip bursts occur at a range of along-strike locations, from 0 to 7 km away from the surface slip observations. Slip occurs at depths from 0 to 10 km; 42%–55% of the bursts are likely below 4 km depth. The bursts typically have moments equivalent to M w 3.2–4.1 earthquakes. These findings suggest that creep events are not just small shallow events; they are relatively large events that nucleate at significant depths and could play a prominent role in the slip dynamics of the creeping section.

Fine-grained authigenic magnetite has been recognized increasingly in iron-rich marine environments affected by methane seepage and is a major sedimentary magnetization source. However, it is unknown whether this magnetite forms via microbial or abiotic processes. We report here abundant fine magnetite crystals, in close association with goethite, within coarse-grained sediments from two adjacent methane seepage sites in the South China Sea. The magnetite- and goethite-rich horizons have sharply increased Zr/Ti, Zr/Rb, Ti/Al, and Fe/Al ratios, probably reflecting deposition by turbidity currents. Deeper intervals have elevated pyrite content, positive δ34S excursions of chromium reducible sulfur, and low magnetic susceptibilities, which is consistent with past sulfate-driven anaerobic oxidation of methane in environments with dynamically variable seepage intensity. In magnetically extracted aggregates (>63 μm), magnetite particles are mainly clustered euhedral crystals with 0.2–0.8 μm sizes, which will likely impact sedimentary magnetic signals. The fine, euhedral crystalline nature of the magnetite suggests formation in sulfide-free, ferrous iron-rich sedimentary environments. Based on 16S rRNA gene sequences, anaerobic methanotrophic archaea coincide with pyrite rich horizons. In contrast, two co-occurring methanogenic archaea groups of the Methanomicrobia class (mainly Methanosarcina and Methanocella) are particularly abundant in turbidites but have low abundance in all other horizons. Increased Methanomicrobia abundances suggest that this class of archaea may be involved in microbial iron reduction in turbidites with abundant goethite as a reactive iron source, and that they apparently trigger magnetite formation. Our findings provide new clues to microbial magnetite formation in iron-rich marine sediments.

How induced seismicity in deep geothermal project (enhanced geothermal systems, EGS) is controlled by fluid injection is of central importance for monitoring the related seismic risk. Here we analyze the relationship between the radiated seismic energy and the hydraulic energy related to the fluid injection during several hydraulic stimulations and circulation tests at Soultz-sous-Forêts geothermal site. Based on a harmonized database, we show that the ratio between these energies is at first order constant during stimulations and of the same magnitude independently of the stimulation protocol and injection depth. Re-stimulations are characterized by a sharp evolution of this ratio during injection which ultimately converges to the characteristic value of the reservoir. This supports that the seismicity is caused by the relaxation of the pre-existing strain energy in the stimulated volume, rather than by the deformation generated from fluid injection. The ratio appears as an intrinsic large-scale property of the reservoir that can be assessed at the very beginning of the first stimulation. Based on this property, we suggest a way to predict the largest magnitude of the induced seismic events knowing the maximum targeted hydraulic energy of the injection.

The Gakkel Ridge in the Eurasian Basin has the slowest seafloor spreading worldwide. The western Gakkel Ridge (3°W–85°E; 14–11 mm/a) alternate between magmatic and sparsely magmatic zones, while the eastern Gakkel Ridge (85–126°E; 11–6 mm/a) appears to be dominated by magmatic zones despite ultraslow spreading. Little is known about the seafloor spreading conditions in the past along the entire ridge. Here, we exploit the residual bathymetry and basement roughness to assess the crustal accretion process of the Gakkel Ridge over time using 23 published regional multichannel seismic reflection profiles. Full seafloor spreading rates were faster (20–24 mm/a) up to ∼45 Ma, and residual bathymetry for the older crust is deeper than the world average in the entire Eurasian Basin. There is a sharp transition to 300–400 m shallower residual bathymetry for seafloor <45 Ma in the eastern Eurasian Basin. The crustal roughness versus spreading rate of the western Eurasian Basin is on the global trend, while that of the eastern is significantly below. Both low roughness and shallow residual bathymetry of the eastern Eurasian Basin is close to that of oceanic crust for spreading rates above 30 mm/a, demonstrating increased magmatic production of the eastern Gakkel Ridge since ∼45 Ma. A recent mantle tomography model predicts partial melting in the upper mantle based on the low Vs anomaly underneath. The sedimentary pattern toward the Lomonosov Ridge indicates that this hot mantle anomaly started to cause dynamic uplift of the area at ∼45 Ma.

In the Earth's upper crust, rocks deform mostly by means of brittle fracturing processes. At the micro-scale these processes involve the formation and growth of microcracks in the vicinity of defects such as open fissures, pores and other cavities. Large defects can produce strong enough perturbations of the stress field to activate and/or intensify brittle damage around them. Here we considered the ideal cases of smooth cylindrical and spherical pores inside an infinite solid body subjected to remote triaxial compressive stresses (i.e., the intermediate and minimum principal stresses are assumed equal). We first established the resulting local stress field around one of those large pores and then verified whether certain brittle damage processes could be activated in these conditions. We mainly considered the formation of tensile microcracks and micro shear bands. The former requires the presence of tensile stresses in some regions around the pore, while the latter needs sufficiently large shear stresses on pre-existing optimally inclined microcracks to overcome their frictional resistance to sliding. We find that shear driven deformation remains localized in the vicinity of the pore. On the other hand, dilatational, tensile cracks can propagate large distances away from the pore but cannot form at and above some threshold ratio of the least to the largest principal stresses. Faulting associated with interacting tensile cracks is therefore suppressed with increasing depth. Our analysis leads to conclusions generally consistent with published experimental observations and provides some clues to discuss the physical cause of the brittle-ductile transition in rocks.

In this study, we measure seismic velocity variations during two cycles of crustal inflation and deflation in 2020 on the Reykjanes peninsula (SW Iceland) by applying coda wave interferometry to ambient noise recorded by distributed dynamic strain sensing (also called DAS). We present a new workflow based on spatial stacking of raw data prior to cross-correlation which substantially improves the spatial coherency and the time resolution of measurements. Using this approach, a strong correlation between velocity changes and ground deformation (in the vertical and horizontal direction) is revealed. Our findings may be related to the infiltration of volcanic fluids at shallow depths, even though the concurrent presence of various processes complicates the reliable attribution of observations to specific geological phenomena. Our work demonstrates how the spatial resolution of DAS can be exploited to enhance existing methodologies and overcome limitations inherent in conventional seismological data sets.