Physics of the Earth and Planetary Interiors最新文献

Earth's subduction zone processes and surface environments are intricately governed by mass transfer phenomena at plate convergent boundaries. The determination of their rates and timings from high-pressure metamorphic rocks (e.g., eclogite), or remnants of ancient convergent boundaries, remains an ongoing challenge. Here, we proposed the potential and versatility of ordering transformation kinetics of omphacite, an essential mineral found in eclogite, as a dynamic recorder of the metamorphic history. Through macroscopic phase-field simulation, we explored the growth of antiphase domains (APDs) in metastable disordered omphacite, discussing the feasibility of constraining metamorphic reaction kinetics based on the size and morphology of omphacite APDs in eclogitized oceanic crust. Our simulation corroborated that omphacite nucleating later during the prograde metamorphism can exhibit an incompletely ordered state with sparsely distributed ordered domains, which suggests their usefulness in estimating the recrystallization timing of the omphacite. Additionally, we confirmed that the APD formation dynamics are significantly influenced by the initial cation configuration of the disordered matrix. This implies the APD morphology in natural omphacite under slab-surface conditions may reflect their precipitation kinetics. These findings provide valuable insights into the microtextural evolution of omphacite due to its ordering transformation, thereby enhancing our ability to interpret morphological features.

Synchrotron X-ray diffraction (XRD) and Raman spectroscopy in laser heated diamond anvil cells and first principles molecular dynamics (FPMD) calculations have been used to investigate the reactivity of calcite and molecular hydrogen (H₂) at high pressures up to 120 GPa. We find that hydrogen reacts with calcite starting below 0.5 GPa at room temperature forming chemical bonds with carbon and oxygen. This results in the unit cell volume expansion; the hydrogenation level is much higher for powdered samples. Single-crystal XRD measurements at 8–24 GPa reveal the presence of previously reported III, IIIb, and VI calcite phases; some crystallites show up to 4% expansion, which is consistent with the incorporation of ≤ 1 hydrogen atom per formula unit. At 40–102 GPa XRD patterns of hydrogenated calcite demonstrate broadened features consistent with the calcite VI structure with incorporated hydrogen atoms. Above 80 GPa, the CO stretching mode of calcite splits suggesting a change in the coordination of CO bonds. Laser heating at 110 GPa results in the formation of CC bonds manifested in the crystallization of diamond recorded by in situ XRD at 300 K and 110 GPa and by Raman spectroscopy on recovered samples commenced with C¹³ calcite. We explored several theoretical models, which show that incorporation of atomic hydrogen results in local distortions of CO₃ groups, formation of corner-shared CO polyhedra, and chemical bonding of H to C and O, which leads to the lattice expansion and vibrational features consistent with the experiments. The experimental and theoretical results support recent reports on tetrahedral C coordination in high-pressure carbonate glasses and suggest a possible source of the origin of ultradeep diamonds.

The electrical conductivity of subsurface rocks is generally anisotropic. The anisotropy of the subsurface electrical conductivity provides important information on the stress-strain state and geodynamics. To quantitatively interpret anisotropic conductivity structures revealed by electromagnetic surveys, it is essential to use a mixing model considering the anisotropy. Although there exists a mixing model for transversely isotropic rocks with crack-shaped pores, the previous model seems inappropriate in interpreting conductive anomalies revealed by electromagnetic exploration because cracks are assumed to be isolated in the model. Therefore, this study develops a theoretical mixing model for transversely isotropic rocks with mutually interconnected cracks by a statistical approach. The derived mixing model considers the macroscopic tortuosity of a collection of cracks as well as the tortuosity of each crack. The derived model can represent general transverse isotropy and includes the isotropic and parallel models as special cases. I compare the developed model to previously proposed mixing models, showing that the developed model can reproduce a much wider range of anisotropy than the already-existing anisotropic mixing model. By applying the developed model to an example of the anisotropic conductivity in the oceanic upper crust inferred by electromagnetic exploration, I demonstrate that the developed mixing model enables us to quantitatively infer the crack orientation and fluid volume fraction that reproduce significant anisotropic conductivity found by field observations. Furthermore, I compare the developed model to the anisotropic seismic velocity model for fluid-filled cracks.

Geomagnetic excursion events have been widely studied in recent years as a key process for understanding the evolution of the Earth's magnetic field. The Santa Rosa geomagnetic excursion (SRE) event during the Matuyama chron has been globally recorded in sediment sequences and lava flows. Galapagos lavas distributed in near-equatorial with an ⁴⁰Ar/³⁹Ar age of 925.7 ± 4.6 ka display absolute paleointensity values of about 14% of the modern magnetic field, which is a valuable record of the Earth's magnetic field strength during the SRE event. However, the above extremely low estimates of paleointensities during the SRE were fitting from higher temperature segments (400 °C–575 °C) from previous paleointensity experiments, which is biased by the thermal instability of Galapagos lava samples during high-temperature heating treatments. From our comprehensive rock magnetic experiments in this study, Galapagos lava samples exhibit thermal instability after heating treatments higher than 400 °C. The severe thermal alteration occurred after the heating temperature reached 500 °C, mainly manifested as an increase in remanence-carrying capacities, such as the enhanced ability of paleointensity specimens to record partial thermoremanent magnetization, resulting in underestimated paleointensities during the SRE. In-depth experiments on rock magnetism and hysteresis parameters analysis provide a powerful method to detect the thermal instability of lava samples, which can help us confirm the biased geomagnetic field strength during this short-lived excursion period and prevent misinterpretations of the Earth's magnetic field evolution through erroneous low paleointensity records.

Episodic Tremor and Slip (ETS) events showcase dynamic interactions of oscillatory slow slips and tremors deep within subduction zones and offer a window into Earth's internal dynamics. However, the exact mechanisms driving these events remain unresolved. This study proposes a novel approach that goes beyond traditional explanations focused on fluid pressure from mineral dehydration. Existing models often neglect the intricate interplay between fluid and rock pressures across various depths and potential fluid sources. This calls for a more comprehensive understanding of how fluid release from reactions interacts with rock deformation. The present formulation captures the interplay between fluid and solid pressures providing a more rigorous picture of ETS events. It employs a minimalistic and efficient approach based on integrating dehydration reactions. The model thereby develops a generic framework for mineral dehydration, offering an enhanced perspective of the underlying processes without the need to trace down to specific minerals. It allows a refined fit to GPS data by including high-frequency components from linear and nonlinear stability analyses, giving rise to improved correlation coefficients. Through the inclusion of the dynamic interplay between fluid and rock pressure diffusion within subduction zones, we propose a unified model of ETS events.

In this study, we investigate the finite deformation of a polycrystalline mixture of bridgmanite (Br) and ferropericlase (Fp) by diffusion creep at the lower mantle-like temperature and pressure by using the self-consistent approach. We explore the influence of volume fraction of Fp, viscosity contrast, and strain dependence (effect of shape change) under both axial (coaxial deformation) and simple shear (non-coaxial deformation). Our present study shows: i) the strength (viscosity) contrast between Fp and Br increases with strain since the viscosity of Fp significantly decreases as Fp grain elongates, and (ii) deformation starts from nearly homogeneous strain to finally nearly homogeneous stress under simple shear whereas deformation behavior remains nearly homogeneous strain under axial deformation. A more substantial creep rate partitioning occurs in simple shear than in axial deformation. These results imply that strain localization via diffusion creep might occur in the lower mantle, particularly in regions where the simple shear is dominated (i.e., in the boundary layers (e.g., the D″ layer)).

The Early Miocene lacustrine sediments of the Most Basin in the Czech Republic preserve a European continental paleoenvironmental archive. A number of paleoenvironmental and magnetostratigraphic studies have been carried out on sediment cores from boreholes due to ongoing coal mining in the basin. However, the magnetic carriers of the studied sediments have not been identified clearly. Here, we present a detailed paleo-rock magnetic study from the Burdigalian sediments near the Bilina mining area, Most Basin. The studied clay sediments cover the period of local lakes and a basin-wide lake above the main coal seam. Our results suggest that the magnetic carriers of the studied section in the Most Basin are mixtures of authigenic greigite and magnetite magnetofossils with overlapping magnetic signatures. Greigite is formed by migration of pore water through the sediment column, where iron from siderite grains reacts with these fluids with limited H₂S, which then favors greigite precipitation. The co-existence of greigite and magnetite indicates a partial dissolution of magnetofossils due to H₂S deficiency. Diagenetic greigite has been problematic in paleomagnetic studies due to an unknown time lag between the depositional remanence and the chemical remanent magnetization (CRM). A ghost polarity interval reveals that greigite acquired at least ∼45 kyr delayed CRM. The revealed timing of remanence acquisition brings a new perspective to the chronostratigraphic structure of the Most Basin.

The statistical features of earthquake clusters in North-Central Iran (Tehran Region) are investigated, with the aim of quantitatively characterizing the properties of earthquake triggering and allow exploring their possible relations with the tectonic setting of the study area.

The nearest-neighbor approach is used for the identification of the earthquake clusters in the space-time-energy domain. This approach permits for a data-driven identification of clusters so that, within multi-event clusters, the features of secondary and higher orders dependent events can be explored. The study is based on a revised dataset that is extracted from the catalog compiled by the Iranian Seismological Center (IRSC) for the period of 1996–2022. In order to exclude the effect of non-tectonic events, which turn out quite numerous within the study region, explosions within quarry-rich areas are removed; the identification of non-tectonic events is performed by considering the normalized ratios of daytime to nighttime events in an iterative removal procedure. According to preliminary analysis of the resulting catalog, an area is selected, within which a satisfactory completeness level is assessed for events with magnitude >2.0. Robust values of the scaling parameters, namely the b-value and the fractal dimension of epicenters, are also computed and are used to calculate the nearest-neighbor distances and to identify the earthquake clusters.

The nearest-neighbor method also permits to investigate the internal structure of earthquake sequences, and to differentiate the spatial properties of seismicity according to the different topological features of the clusters structure. The obtained results allow us identifying two macro-areas, approximately separated by the 52°E meridian, which are characterized by different clustering features, namely: high complexity indexes, indicating simple (burst-like) structure of clusters, to the East; low complexity index, corresponding to complex multi-level (swarm-like) structure of clusters, to the West. The complexity measures, borrowed from network theory (i.e. the Closeness and Outdegree Centralization indexes), consistently capture the complexity of the identified clusters, and confirm that the cluster structures have distinct preferred geographic locations. The territorial heterogeneity of the examined clustering properties can be related with the spatial variability of tectonic, structural and geophysical features of the Alborz region, in good agreement with findings from the Alps-Dinarides junction (Northeastern Italy), a region also characterized by a contractional structural setting, mainly including reverse and strike-slip faulting systems, and by moderate to high seismic activity.

Palaeomagnetic and modern geomagnetic measurements indicate that the South Atlantic Anomaly (SAA) has undergone rapid changes over the past few hundred years. Its minimum intensity decreased at an average rate of 26 not yr⁻¹, accompanied by a continuous westward drift and spatial expansion. Recently, a secondary minimum of SAA emerged near southern Africa, leading to speculation that expansion of the SAA could indicate an impending geomagnetic reversal. Here, we focus on the evolution and disappearance of the paleo-West Pacific Anomaly (WPA), as another SAA-like structure, which may have implications for the future of SAA evolution. We regard the WPA as SAA-like due to its feature and its association with a reversal flux patch on the core-mantle boundary. Consequently, we suggest that the observed evolutionary pattern in the WPA can serve as a reference for other negative anomalies, such as the SAA. By analysing models that combine datasets of archaeomagnetic and historical records, such as gufm1 and HistKalmag, it is found that the WPA occurred between 1600 and 1820 CE. Over its duration, the WPA experienced phases of rapid expansion, drift, and division. Eventually, its primary component faded away, giving rise to a new segment that continued to expand. The initial two evolutionary phases of the WPA are similar to the evolution of the SAA over the past century. According to the WPA's evolution, it suggests that the current state of the SAA may correspond to an early stage of splitting. Forecasts based on the evolution of the WPA indicate a rapid expansion of the anomalous region in the short term, followed by a gradual reduction in its primary component and continued expansion of a new local minimum. This study provides valuable insight into the evolution of the SAA and highlights the potential utility of the WPA as an evolutionary reference for such geomagnetic phenomena.

Due to the gradual and constant accumulation of seismic energy, Peninsular India (PI) is typically considered seismically stable with low to moderate seismicity. The seismic studies in Peninsular India always resorted to synthetic ground motion simulations, because of the limited instrumentation and hence lack of recorded data. In the absence of a well-defined medium model for PI, the usual practice is to use simple site proxies or one-dimensional velocity structures for ground motion simulations. However, the region consists of multi-scale geometric complexities, significant topography, and sedimentary basins and is surrounded by deep oceans. Thus, the radiated seismic wave field in the region is influenced by the medium properties and in the absence of a well-defined tomography model the reliable estimation of seismic hazard is a challenging problem in PI. Therefore, the seismic wave propagation in PI can be investigated using numerical simulation with reliable 3D computational model for PI, incorporating the knowledge of the underlying Earth structure. Hence, the present study attempts to develop a sophisticated three-dimensional (3D) medium model of Peninsular India for physics-based ground motion simulations for regional earthquakes. This is aided by the availability of one-dimensional (1D) velocity models and the crustal structure from the receiver function analysis which provides valuable insight into the variation of material properties in the region. In the present study, >100 s of 1D velocity profiles are collected from various literature, which is then grouped under 23 different geological regions identified in PI (as per GSI (2000)). The averaged material properties are assigned per each geological region and the information on sediment depths, basin geometry, topography, and bathymetry are incorporated. We use the spectral element method (SEM) to calibrate our 3D computational model by simulating synthetic seismograms and comparing them to recorded ground motions for two past earthquakes: the 2001 Mw 7.6 Bhuj earthquake and the 1997 Mw 5.8 Jabalpur earthquake. Further, the seismic waveforms at the near field of 2001 Mw 7.6 Bhuj event are simulated using a refined regional model. The spatial variability of associated seismic intensities and peak ground velocity (PGV) amplification are investigated. In addition, a study of the impact of model depth truncation and sphericity on ground motion is also conducted. The implemented medium model is the first of its kind for Peninsular India and can reliably be used in seismic wave propagation studies in the region. The simulated outcomes from the model are of engineering importance as these results can be used for seismic hazard assessment of the region.