Fault Geometry and Slip Distribution of the 2021 Mw 7.4 Maduo, China, Earthquake Inferred from InSAR Measurements and Relocated Aftershocks

A nearly 70 yr hiatus of major seismic activity in the central eastern Bayan Har block (BKB) ended on 22 May 2021, when a multislip-peak sinistral strike-slip earthquake struck western Maduo County, Qinghai. This earthquake, which ruptured the nearly 170 km long Kunlun Pass–Jiangcuo fault, is a rather unique event and offers a rare opportunity to probe the mechanical properties of the intraplate lithosphere of the central eastern BKB. Here, we inferred the fault geometry associated with the Maduo earthquake using Interferometric Synthetic Aperture Radar (InSAR), and relocated aftershocks and inverted the slip distribution through InSAR radar phases and range offsets. Our analysis revealed that the geometry of the fault varies along the strike: the southeastern end of the fault dips steeply to the northeast, whereas the northwestern end dips southwestward. Using the combined datasets to constrain a coseismic slip, we found that the 2021 Maduo event was dominated by sinistral strike-slip movement, with a slight normal-slip component at a shallow depth, rupturing the steep-dipping fault for nearly 170 km in length. Five asperities were detected along the fault strike in the shallow crust (0–12 km) with a peak slip of ∼4.2 m corresponding mostly to simple structures, namely, continuous and straight rupture segments, suggesting that the rupture propagated across geometrical barriers in a multiasperity way. Based on an analysis of the strain field and the focal mechanisms of both the 2021 Maduo earthquake and historical earthquakes that have occurred in the BKB, we propose that the fault zones within the BKB can also generate large earthquakes and have the ability to accommodate the ongoing eastward and northeastward penetration of the Indian plate into the Eurasian plate.

Dating folding beyond folding, from layer-parallel shortening to fold tightening, using mesostructures: lessons from the Apennines, Pyrenees, and Rocky Mountains

Dating syntectonic sedimentary sequences is often seen as the unique way to constrain the initiation, duration, and rate of folding as well as the sequence of deformation in the shallow crust. Beyond fold growth, however, deformation mesostructures accommodate the internal strain of pre-folding strata before, during, and after strata tilting. Absolute dating of syn-folding mesostructures may help constrain the duration of fold growth in the absence of preserved growth strata. Absolute dating of mesostructures related to early-folding layer-parallel shortening and late fold tightening provides an access to the timing and duration of the entire folding event. We compile available ages from the literature and provide new U–Pb ages of calcite cements from veins and faults from four folds (Apennines, Pyrenees, Rocky Mountains). Our results not only better constrain the timing of fold growth but also reveal a contraction preceding and following folding, the duration of which might be a function of the tectonic style and regional sequence of deformation. This study paves the way for a better appraisal of folding lifetime and processes and stress evolution in folded domains.

Separation of Split Shear Waves in Two Non-orthogonal Sets of Vertical Fractures

The phenomenon of S-wave splitting indicates the development of fractures in the shallow crust. Therefore, methods based on S-wave splitting have been established to predict the development of one set of parallel fractures. However, for rocks containing two non-orthogonal sets of vertical fractures, the mechanism of S-wave splitting is more complex, and the available methods cannot be applied. To resolve this inadequacy, we propose a two-way rotation method to separate split S-waves with the aim of restoring the split S-wave polarizations and predicting the fracture azimuths. First, we calculate the stiffness matrix of fractured media based on the linear slip theory and derive the phase velocities and polarizations of split S-waves induced by fractures using the Christoffel equation. Second, we clarify the S-wave splitting mechanism in this media by employing velocity analysis and deconstruct the S-wave polarizations on the horizontal components. Third, we deduce a two-way rotation matrix obtained by the S-wave splitting modes to separate the split S-waves. To solve for the angle parameters related to the fracture azimuths in the two-way rotation matrix, we superpose the subspace polarizations in two dimensions to determine the polarization azimuths of the split S-waves. Numerical model tests demonstrate that the proposed method is stable under noisy conditions. Finally, we apply the proposed method to real near-offset and walkaround VSP data, and the predicted fracture results are verified by imaging logs and prior knowledge.

The trisulfur radical ion S3•− controls platinum transport by hydrothermal fluids

Platinum group elements (PGE) are considered to be very poorly soluble in aqueous fluids in most natural hydrothermal–magmatic contexts and industrial processes. Here, we combined in situ X-ray absorption spectroscopy and solubility experiments with atomistic and thermodynamic simulations to demonstrate that the trisulfur radical ion S3•− forms very stable and soluble complexes with both PtII and PtIV in sulfur-bearing aqueous solution at elevated temperatures (∼300 °C). These Pt-bearing species enable (re)mobilization, transfer, and focused precipitation of platinum up to 10,000 times more efficiently than any other common inorganic ligand, such as hydroxide, chloride, sulfate, or sulfide. Our results imply a far more important contribution of sulfur-bearing hydrothermal fluids to PGE transfer and accumulation in the Earth’s crust than believed previously. This discovery challenges traditional models of PGE economic concentration from silicate and sulfide melts and provides new possibilities for resource prospecting in hydrothermal shallow crust settings. The exceptionally high capacity of the S3•− ion to bind platinum may also offer new routes for PGE selective extraction from ore and hydrothermal synthesis of noble metal nanomaterials.

Supergiant porphyry Cu deposits are failed large eruptions

Porphyry copper deposits, the principal natural source of Cu and Mo, form at convergent margins. Copper is precipitated from fluids associated with cooling magmas that have formed in the mantle and evolved at variably deep crustal levels, before raising close to the surface where they exsolve fluids and copper. Despite significant advances in the understanding of their formation, there are still underexplored aspects of the genesis of porphyry copper deposits. Here, we address the role played by magma injection rates into the shallow crust on the formation of porphyry copper deposits with different copper endowments. Using a mass balance approach, we show that supergiant porphyry Cu deposits (>10 Mt Cu) require magma volumes and magma injection rates typical of large volcanic eruptions. Because such volcanic events would destroy magmatic-hydrothermal systems or prevent their formation, the largest porphyry Cu deposits can be considered as failed large eruptions and this may be one of the causes of their rarity.

Seismic Velocity Response to Atmospheric Pressure Using Time-Lapse Passive Seismic Interferometry

ABSTRACT Seismic velocities in the shallow crust down to a few kilometers depth show a remarkable sensitivity to stress perturbations due to the presence of compliant pores, cracks, fractures, and faults. Monitoring temporal changes of seismic velocities can thus provide key insights on dynamic processes affecting the shallow crust such as those related to the atmosphere (rainfall, barometric pressure, and temperature) and those with deeper tectonic and volcanic origins. In this work, we investigate the specific response of the near surface down to 300 m depth to atmospheric pressure variations. We conduct a four month passive seismic monitoring experiment in the desert of Oman using continuous noise recorded at geophones located within five wells. The results show a clear, direct correlation between seismic velocities and barometric pressure variations for monthly transients. At a longer, seasonal temporal scale, seismic velocities are stable, whereas atmospheric pressure shows a clear positive trend. We use the undrained coupled poroelastic theory to model these observations and find that the lack of seasonal velocity changes can be partly explained by the atmospheric pressure that diffuses into the pores with a strong hydraulic diffusivity likely higher than 100 m2/s consistent with the local geology referring to carbonates. Finally, the comparison between the modeled and observed velocity changes leads to estimate a velocity–stress sensitivity on the order of 6.3×10−7 Pa−1 which is consistent with previous studies. Using this result for calibration, we find that a sudden step-change drop of velocity of 0.015% occurring in the beginning of October 2019 and corresponding to a stress perturbation likely larger than 240 Pa affected the entire studied area. This small change could be related to a perturbation at greater depth associated with variations in the production rates within the underlying reservoir.

Temporal and Petrogenetic Links Between Mesoproterozoic Alkaline and Carbonatite Magmas at Mountain Pass, California

Mountain Pass is the site of the most economically important rare earth element (REE) deposits in the United States. Mesoproterozoic alkaline intrusions are spatiotemporally associated with a composite carbonatite stock that hosts REE ore. Understanding the genesis of the alkaline and carbonatite magmas is an essential scientific goal for a society in which critical minerals are in high demand and will continue to be so for the foreseeable future. We present an ion microprobe study of zircon crystals in shonkinite and syenite intrusions to establish geochronological and geochemical constraints on the igneous underpinnings of the Mountain Pass REE deposit. Silicate whole-rock compositions occupy a broad spectrum (50–72 wt % SiO2), are ultrapotassic (6–9 wt % K2O; K2O/Na2O = 2–9), and have highly elevated concentrations of REEs (La 500–1,100× chondritic). Zircon concordia 206Pb/238U-207Pb/235U ages determined for shonkinite and syenite units are 1409 ± 8, 1409 ± 12, 1410 ± 8, and 1415 ± 6 Ma (2σ). Most shonkinite dikes are dominated by inherited Paleoproterozoic xenocrysts, but there are sparse primary zircons with 207Pb/206Pb ages of 1390–1380 ± 15 Ma for the youngest grains. Our new zircon U-Pb ages for shonkinite and syenite units overlap published monazite Th-Pb ages for the carbonatite orebody and a smaller carbonatite dike. Inherited zircons in shonkinite and syenite units are ubiquitous and have a multimodal distribution of 207Pb/206Pb ages that cluster in the range of 1785–1600 ± 10–30 Ma. Primary zircons have generally lower Hf (<11,000 ppm) and higher Eu/Eu* (>0.6), Th (>300 ppm), Th/U (>1), and Ti-in-zircon temperatures (>800°C) than inherited zircons. Oxygen isotope data reveals a large range in δ18O values for primary zircons, from mantle (5–5.5‰) to crustal and supracrustal (7–9‰). A couple of low-δ18O outliers (2‰) point to a component of shallow crust altered by meteoric water. The δ18O range of inherited zircons (5–10‰) overlaps that of the primary zircons. Our study supports a model in which alkaline and carbonatite magmatism occurred over tens of millions of years, repeatedly tapping a metasomatized mantle source, which endowed magmas with elevated REEs and other diagnostic components (e.g., F, Ba). Though this metasomatized mantle region existed for the duration of Mountain Pass magmatism, it probably did not predate magmatism by substantial geologic time (>100 m.y.), based on the similarity of 1500 Ma zircons with the dominantly 1800–1600 Ma inherited zircons, as opposed to the 1450–1350 Ma primary zircons. Mountain Pass magmas had diverse crustal inputs from assimilation of Paleoproterozoic and Mesoproterozoic igneous, metaigneous, and metasedimentary rocks. Crustal assimilation is only apparent from high spatial resolution zircon analyses and underscores the need for mineral-scale approaches in understanding the genesis of the Mountain Pass system.

Dating folding beyond folding, from layer-parallel shortening to fold tightening, using mesostructures: Lessons from the Apennines, Pyrenees and Rocky Mountains

Dating syntectonic sedimentary sequences is often seen as the unique way to constrain the initiation, kinematics and rate of folding and the sequence of deformation in the shallow crust. Beyond fold growth however, deformation mesostructures accommodate the internal shortening of pre-folding strata before, during and after strata tilting. Absolute dating of mesostructures developed during extension at fold hinge may help constrain the duration of fold growth in the absence of preserved growth strata, while dating of mesostructures related to layer-parallel shortening and late fold tightening provide a valuable access to the timing and duration of the entire folding event. We compile existing ages in the literature and provide new U-Pb ages of calcite cements from veins and faults from four folds (Apennines, Pyrenees, Rocky Mountains). Our results not only better constrain the timing of fold growth but also reveal a contraction preceding and following folding, the duration of which might be function of the tectonic style and regional sequence of deformation. This study paves the way for a better appraisal of folding lifetime and processes and of stress evolution in folded domains.