Progress in Earth and Planetary Science最新文献

The Ou Backbone Range in the Tohoku Region of Japan is bounded on its both sides by two major active fault systems: the Western Fault Zone of Kitakami Lowland in the east and the Eastern Fault Zone of Yokote Basin in the west. Although these two systems are primary active fault systems in the region, information on their long-term activity characteristics is still quite limited. Therefore, we analyzed the normalized channel steepness indexes of river valleys trending perpendicular to the range front along both flanks of the Ou Backbone Range. Our results show that the eastern flank has gentler river valleys, whereas rivers along the northwestern flank are steeper. Our field investigation shows that knickpoints in this area are mostly related to local lithologic boundaries or are check dams along the valleys, thus the river systems are likely under steady-state conditions. Hence, the steeper river valleys in northwestern Ou Backbone Range indicate a higher uplift rate of the area. Because both fault systems are primarily dip-slip reverse faults and do not have significant variations in their subsurface geometry, the faster uplift suggests that the northern segment of the Eastern Fault Zone of Yokote Basin has a higher slip rate. This is consistent with results of previous studies, and the fact that the rupture of the 1896 Rikuu earthquake, the only historical surface-rupturing event in this region, was only limited along the northern segment of this fault system.

The microstructure of realistic fluid–rock systems evolves to minimize the overall interfacial energy, enabling local variations in fluid geometry beyond ideal models. Consequently, the permeability–porosity relationship and fluid distribution in these systems may deviate from theoretical expectations. Here, we aimed to better understand the permeability development and fluid retention in deep-seated rocks at low fluid fractions by using a combined approach of high-resolution synchrotron radiation X-ray computed microtomography imaging of synthesized rocks and numerical permeability computation. We first synthesized quartzite using a piston-cylinder apparatus at different fluid fractions and wetting properties (wetting and non-wetting systems with dihedral angles of 52° and 61°–71°, respectively) under conditions of efficient grain growth. Although all fluids should be connected along grain edges and tubules in the homogeneous isotropic wetting fluid–rock system enabling segregation by gravitational compaction in natural settings, the fluid connectivity rapidly decreased to ~ 0 when the total fluid fraction decreased to 0.030–0.037, as the non-ideality of quartzite, including the interfacial energy anisotropy (i.e., grain faceting), became critical. In non-wetting systems, where the minimum energy fluid fraction based solely on the dihedral angle is ~ 0.015–0.035, the isolated (disconnected) fractions was 0.048–0.062. A streamline computation in the non-wetting system revealed that with decreasing total porosity, flow focusing into fewer channels maintained permeability, allowing the effective segregation of the connected fluids. These results provide insight into how non-wetting fluids segregate from rocks and exemplify the fraction of retained fluids in non-wetting systems. Thus, the findings suggest a potential way for wetting system fluids to be transported into the deep Earth's interior, and the amount of fluids dragged down to the Earth’s interior could be higher than what was previously estimated.

This study aims to quantify the topography of the Moho boundary, the lower crust and uppermost mantle contact of Egypt, in order to estimate the crustal thickness variation and its link to the distribution of thermal anomalies under Egypt. This is accomplished by modeling satellite gravity, supported by the passive seismic constraints throughout Egypt. However, when estimating the thickness of the crust in Egypt using just seismic data, substantial uncertainty and deviation are produced due to the sparsely dispersed stations. Integrating seismic and gravity data minimizes uncertainty and improves estimate accuracy. The investigation is broken down into four stages, the first involving utilizing the Sentinel-3B satellite to create land surface temperature maps. The subsequent steps consist of gravity and seismic data adjustments, inversion and forward modeling. We used seismically restricted nonlinear inversion to look at Goco06s satellite gravity data to model the Moho’s topographic surface. The data gathered from deep seismic refraction and receiver functions adjusted the analyzed data. The inversion process relies on the adapted Bott's approach and Tikhonov regularization, using the assumption of the sphericity of the Earth planet. Reference values for depth of Moho and density contrast were set at 35 km and 500 kg/m³, respectively. The average statistical difference for Moho depth between gravity-based model and seismic data is − 0.10 km. Through forward gravity modeling, five gravity profiles were chosen and interpreted in 2.5D models. The results indicated that the Moho depth in the south varies from 35 to 39 km and decreases in the north and the Mediterranean. In upper Egypt, the highest Moho depth is 39 km. The depth varies beneath the Sinai Peninsula as it is about 35 km in its south, reaches 30 km in the northern portion, and ranges along the Red Sea’s Rift Margin from 29 to 32 km. Moreover, the final model shows the relation between Moho coincides with the surface temperature anomalies approved by satellite images and hot springs. The model reveals a correlation between Moho discontinuity and surface temperature anomalies, revealing the highest geothermal potential in a rectangular area in central Egypt, between latitudes 25°N and 30°N, based on satellite imagery and hot springs distribution.

Garnet peridotite xenoliths have been rarely reported from suboceanic mantle. Petrographic and geochemical characteristics of garnet-bearing oceanic peridotite xenoliths provide precious information on dynamics of the suboceanic lithosphere and asthenosphere interaction. We examined a lherzolite xenolith included in olivine nephelinite lava from Aitutaki Island, a member of the Cook-Austral volcanic chain. The lherzolite xenolith contains reddish fine-grained (< 5 µm in size) mineral aggregates (FMAs) with size range of 0.5–6 mm, consisting of olivine, calcic and sodic plagioclases, aluminous spinel, native iron, and nepheline. Microstructural observations and chemical data corroborate that the FMA is a decomposed pyrope-rich garnet including chromian spinel grains with an irregular highly indented morphology in the center. The FMA is surrounded by pyroxene-poor and olivine-rich aureole. The spatial and morphological relationships of FMA and chromian spinel with pyroxene-depleted margin suggest a reaction of aluminous spinel + pyroxenes → pyrope-rich garnet + olivine, which requires a compression before decomposition of the garnet to FMA. An orthopyroxene grain shows slight but clear chemical zoning characterized by increase in Al, Ca, and Cr from the grain center to the rim. The zoning patterns of Al and Ca in the orthopyroxene grain can be modeled by diffusion-controlled solid-state reactions induced by pressure and temperature changes, keeping surface concentrations in equilibrium with the other coexisting mineral phases. The results indicate that the mantle, from which the lherzolite xenolith was derived, underwent isothermal decompression followed by a weak heating on a time scale of a few tenths of million years before the xenolith extraction. From the deduced compression and decompression histories, we hypothesize that the mantle beneath Aitutaki Island was once dragged down to a garnet-stable deep mantle region and brought up later by small-scale sublithospheric convection.

Atmospheric transport and dispersion models (ATDMs) are widely used to study and forecast pollution events. In the frame of the “Effect of Megacities on the transport and transformation of pollutants on the regional to global scales” (EMeRGe) project, ATDM forecasts were carried out to identify potential airborne sampling areas of perfluorocarbons (PFCs) emanating from controlled PFC releases. The forecasts involved short-distance transport over small-scale topographic maxima (Manilla; Philippines), short-distance transport over large-scale topographic maxima (Taipei, Taiwan) and long-distance transport over mixed topography (Nanjing, China, sampled over Taiwan). In situ aircraft measurements of PFC mixing ratios down to a few parts per quadrillion (ppqv) provide us with a unique dataset to explore the added benefits of dynamical downscaling. Transport simulations were repeated using FLEXPART driven by ERA5 and IFS meteorological data and FLEXPART-WRF with dynamically downscaled IFS data down to 1.1 km and four PBL parametrisations. Of the three cases studied, dynamical downscaling led to significant differences for the Manilla and Taipei releases that can be interpreted through changes in the modelled orographic flow regimes. The choice of PBL scheme also significantly impacted accuracy, but there was no systematically better-performing option, highlighting the benefits of ensemble forecasting. Results show how convergence and divergence between ensemble members can be utilised to help decision-making during field campaigns. This study highlights the role that dynamical downscaling can play as an important component in campaign planning when dealing with observations over orographically complex areas.

Methane emission in South Asia is poorly understood due to a lack of observations, despite being a major contributor to methane emissions globally. We present the first results of atmospheric CH₄ inversions using air samples collected weekly at Nainital, India (NTL), and Comilla, Bangladesh (CLA), in addition to surface background flask measurements by NOAA, CSIRO and AGAGE using the MIROC4-ACTM. Our simulations span from 2000 to 2020 (considering the fixed “edge” effect), but the main analysis period is 2013–2020, when both the NTL and CLA datasets are available. An additional flux uncertainty reduction of up to 40% was obtained (mainly in the northern part of the Indian subcontinent), which enhanced our confidence in flux estimation and reaffirmed the significance of observations at the NTL and CLA sites. Our estimated regional flux was 64.0 ± 4.7 Tg-CH₄ yr⁻¹ in South Asia for the period 2013–2020. We considered two combinations of a priori fluxes that represented different approaches for CH₄ emission from rice fields and wetlands. By the inversion, the difference in emissions between these combinations was notably reduced due to the adjustment of the CH₄ emission from the agriculture, oil and gas, and waste sectors. At the same time, the discrepancy in wetland emissions, approximately 8 Tg-CH₄ yr⁻¹, remained unchanged. In addition to adjusting the annual totals, the inclusion of NTL/CLA observations in the inversion analysis modified the seasonal cycle of total fluxes, possibly due to the agricultural sector. While the a priori fluxes consisted of a single peak in August, the a posteriori values indicated double peaks in May and September. These peaks are highly likely associated with field preparation for summer crops and emissions from rice fields during the heading stage (panicle formation). The newly incorporated sites primarily exhibit sensitivity to the Indo-Gangetic Plain subregion, while coverage in southern India remains limited. Expanding the observation network is necessary, with careful analysis of potential locations using back-trajectory methods for footprint evaluation.

Carbonaceous material (CM) undergoes progressive changes that reflect its thermal history. These changes are in general irreversible and provide valuable information for understanding diagenetic and metamorphic processes of crustal rocks. Among various approaches to quantify these changes, the R2 ratio, area ratio of specific peaks in CM Raman spectra, is widely used to estimate the maximum temperature of intermediate- to moderately high-grade metamorphism. The calculation of the R2 ratio requires peak deconvolution of the original spectrum, and the results depend on the details of how this is carried out. However, a clear protocol for selecting appropriate initial conditions has not been established and obtaining a reliable temperature estimate depends at least in part on the experience and skill of the operator. In this study, we developed a Python code that automatically calculates the R2 ratio from CM Raman spectra. Our code produces R2 ratios that are generally in good agreement with those of Aoya et al. (J Metamorph Geol 28:895–914, 2010, https://doi.org/10.1111/j.1525-1314.2010.00896.x) for the same Raman data, with much less time and effort than was the case in the previous studies. We have confirmed that the code is also applicable to other previous datasets from both contact and regional metamorphic regions. The overall trend of the recalculated data indicates that samples with R2 greater than ~ 0.7 are not sensitive to the changes in CM maturity and thus should not be used for the calibration of an R2-based geothermometer. We propose a modified geothermometer for contact metamorphism that is strictly applicable to samples with R2 from 0.023 to 0.516, with the proviso that a laser with a wavelength of 532 nm should be used. A slight extrapolation of the newly proposed geothermometer up to R2 of 0.57 provides a temperature estimate that is consistent with the geothermometer of Kaneki and Kouketsu (Island Arc 31:e12467, 2022; https://doi.org/10.1111/iar.12467); the boundary between the two geothermometers corresponds to a temperature of 391 °C.

The 1771 Meiwa tsunami which struck the southern Ryukyu Islands (Sakishima Islands) had greater than 22 m run-up height, leaving about 12,000 casualties in its wake. At many places, the tsunami inundation or lack of inundation is well recorded in historical documents. Several tsunami source models have been proposed for this event using historical records as constraints of tsunami calculations. Nevertheless, the source model remains under discussion. This study re-evaluated the tsunami wave source model of the 1771 Meiwa tsunami using high-resolution (10 m mesh) bathymetric and topographical data for tsunami calculation, the latest historical record dataset, and seismological knowledge. Results demonstrated that a tsunami earthquake along the southern Ryukyu Trench was the likely cause of the 1771 event. However, it is noteworthy that assumption of a large slip with 30 m is necessary for a shallow and narrow region (fault depth = 5 km, fault width = 30 km, Mw = 8.49) of the plate boundary in the Ryukyu Trench, which is far larger than previously thought. This requirement of very large initial water level change at the source might involve not only the fault rupture along the plate boundary but also deformation by splay faults, inelastic deformation of unconsolidated sediments near the trench axis, and/or giant submarine landslides. Results also show that the effects of fault parameters on the run-up were quite different depending on the offshore coral reef width. This phenomenon strongly constrained the fault width to 30 km. Our tsunami ray tracing analysis further revealed the effects of bathymetry on tsunami propagation. It is noteworthy that meter-long huge tsunami boulders tend to be distributed along the specific coasts at which the tsunami was concentrated by bathymetric effects. This finding suggests that past tsunamis, including the 1771 event, might have affected the specific coral reefs on Sakishima Islands repeatedly, which is crucially important for understanding the heterogeneous distribution of tsunami boulders. This feature might also be useful to elucidate the effects of large tsunamis on the corals and reefs because a direct comparison of coral reefs that are damaged and not damaged by tsunami waves is testable in narrow areas in the case of the Sakishima Islands.

During their formation, terrestrial planets underwent a magma ocean phase during which their metallic cores segregated from their silicate mantles and their early atmospheres formed. These planetary formation processes resulted in a redistribution of the abundances of highly volatile elements (HVEs, such as H, C, N, and S) between the planets’ metallic cores, silicate mantles, and atmospheres. This review presents the numerous experimental techniques used to simulate the conditions and identify the parameters that influenced the behavior of HVEs during planetary formation. We also review the analytical techniques used to characterize the different types of experimental samples and quantify the distribution of HVEs between metallic and silicate phases, as well as their solubilities in silicate glasses. This exhaustive review targets students and young researchers beginning their work on the subject, or, more generally, scientists seeking a better understanding of this field of research.

From 1998 to 2017, neutron monitors located at an altitude of 4300 m on the Tibetan plateau detected 127 long-duration bursts of high-energy radiation in association with thunderclouds. These bursts typically lasted for 10–40 min, and 89% of them occurred between 10:00 and 24:00 local time. They were also found to be more likely to occur at night, especially during 18:00–06:00 local time period. The observed diurnal and seasonal variations in burst frequency were consistent with the frequencies of lightning and precipitation on the Tibetan plateau. Based on 19 years of data, the present study suggests that an annual variation in burst frequency has a periodicity of (sim) 16 years and a lag of (sim) 3 years relative to solar activity.