Progress in Earth and Planetary Science最新文献

Water is transported into the deep mantle by subducting slabs, playing important roles in mantle dynamics and evolution. An aluminous hydrous mineral, phase δ with a main component of AlOOH, has been considered an important water carrier in the lower mantle. Recent studies reported that SiO₂ stishovite can accommodate weight percent levels of water, indicating another important water carrier in the lower mantle. However, which mineral can mainly carry water is not clear yet. Recent hydrous phase relation studies reported that stishovite is depleted in alumina when coexisting with hydrous phase δ, in which water content of stishovite was not investigated. In this study, we investigated hydrogen partitioning between stishovite and hydrous phase δ at 24–28 GPa and 1000–1200 °C by means of Kawai-type multi-anvil press in combination with Fourier-transform infrared spectroscopy at ambient conditions on recovered samples. Fourier-transform infrared spectra of recovered stishovites showed that water contents of stishovite coexisting with hydrous phase δ were limited to up to ~ 500 ppm. This indicates that coexisting hydrous phase δ causes not only depletion in alumina but also in hydrogen in stishovite and therefore mainly transports water in a cold subducting slab. Once hydrous phase δ becomes thermally unstable, alumina and water contents in silica minerals are increased by the chemical reaction between SiO₂ and AlOOH, and aluminous silica minerals such as stishovite and CaCl₂-type phase will be a main water carrier in the lower mantle. Presence of small-scale seismic scatterers observed around 1900 km depth, which was considered to be caused by a transition from almost pure SiO₂ stishovite to CaCl₂-type phase, might also be able to be explained by the phase transition of stishovite coexisting with hydrous phase δ.

As a result of climate change, the pan-Arctic region has seen greater temperature increases than other geographical regions on the Earth’s surface. This has led to substantial changes in terrestrial ecosystems and the hydrological cycle, which have affected the distribution of vegetation and the patterns of water flow and accumulation. Various remote sensing techniques, including optical and microwave satellite observations, are useful for monitoring these terrestrial water and vegetation dynamics. In the present study, satellite and reanalysis datasets were used to produce water and vegetation maps with a high temporal resolution (daily) and moderate spatial resolution (500 m) at a continental scale over Siberia in the period 2003–2017. The multiple data sources were integrated by pixel-based machine learning (random forest), which generated a normalized difference water index (NDWI), normalized difference vegetation index (NDVI), and water fraction without any gaps, even for areas where optical data were missing (e.g., cloud cover). For the convenience of users handling the data, an aggregated product is provided, formatted using a 0.1° grid in latitude/longitude projection. When validated using the original optical images, the NDWI and NDVI images showed small systematic biases, with a root mean squared error of approximately 0.1 over the study area. The product was used for both time-series trend analysis of the indices from 2003 to 2017 and phenological feature extraction based on seasonal NDVI patterns. The former analysis was used to identify areas where the NDVI is decreasing and the NDWI is increasing, and hotspots where the NDWI at lakesides and coastal regions is decreasing. The latter analysis, which employed double-sigmoid fitting to assess changes in five phenological parameters (i.e., start and end of spring and fall, and peak NDVI values) at two larch forest sites, highlighted a tendency for recent lengthening of the growing period. Further applications, including model integration and contribution to land cover mapping, will be developed in the future.

Abstract

The Osaka Bay is situated at a seismically active region north of the Median Tectonic Line and east of Awaji Island in western Japan, known as part of the Kinki Triangle and the Niigata–Kobe Tectonic Zone. Dense distribution of active faults and high geodetic strain rates characterize the region, posing a major seismic hazard potential to the coastal and metropolitan areas of the Kansai region. To investigate the shallow structure and recent deformation history of active faults in the Osaka Bay, we acquired 15 high-resolution seismic profiles using a Mini-GI airgun and a Boomer as active sources, together with multi-beam bathymetry data across the Osaka Bay Fault. Our seismic sections image a ~ 0.1 to 3.7 km-wide asymmetric anticline forelimb above the Osaka Bay Fault at shallow depths, coupled with a ~ 2.6 km-wide syncline to the west, and a broad, ~ 11 km-wide syncline in the footwall to the east. The synclinal axial surface at shallow depths measured in this study ranges 75°–89°. We observe the vertical displacement of the Osaka Bay Fault increasing northwards along strike. The sediment thickness on the hanging wall, however, is variable, modified by non-tectonic processes such as by tidal currents, affecting the geometry of growth strata. The most recent deformation by the Osaka Bay Fault reaches to near the seafloor by active folding, with large vertical offsets of 8–14 m over the last ~ 11 ka, and 5–11 m over the last ~ 5 ka. By combining with previously reported borehole age data, the average uplift rate on the Osaka Bay Fault is estimated to be ~ 1.0 to 1.7 m/ka during the Latest Pleistocene to Holocene. The inferred slip of the Osaka Bay Fault during the Holocene is likely to account for > 5% of the regional geodetic strain accumulation within the Kinki Triangle. Further studies to evaluate the Holocene slip rates of regional faults are necessary to assess the seismic hazards and the internal strain budgets within the Kinki Triangle and the Niigata–Kobe Tectonic Zone.

Abiotic synthesis of hydrocarbon-bearing fluids during geological processes has a significant impact on the evolution of both the Earth's biosphere and the solid Earth. Aqueous alteration of ultramafic rocks, i.e., serpentinization, which forms serpentinite, is one of the geological processes generating abiotic methane (CH₄). However, abiotic CH₄ generation is not limited to the serpentinization of mafic and ultramafic lithologies. Metasedimentary dolomitic marble from the Hida Belt, Japan, is characterized by the presence of forsterite-rich olivine (Fo_~89–93), and olivine crystals contain abundant fluid inclusions (<1 to 10 μm in size). Raman spectroscopic analyses of olivine-hosted fluid inclusions found that both primary and secondary fluid inclusions contain CH₄, lizardite/chrysotile, and brucite. This indicates that micro-scale interactions between COH fluid and host olivine produced CH₄ through the reduction of CO₂ by H₂ released during local serpentinization within inclusions. Our observation implies that the dolomitic marble has the potential to be a key lithology for the synthesis and storage of abiotic CH₄ in a shallower crustal portion of orogenic belts.

GRACE and GRACE Follow-On (FO) missions provide time-variable gravity field models of unprecedented quality that allow for the hydrological, oceanic, and ice mass change studies on a global scale. However, the very low-degree coefficients derived from GRACE and GRACE-FO are of inferior quality due to thermal effects acting on satellites and malfunctioning of the onboard accelerometers. Therefore, C₂₀ and C₃₀ coefficients describing the Earth’s oblateness and the pear shape of the Earth, respectively, are being replaced by values derived from satellite laser ranging (SLR) in the standard GRACE solutions. This study assesses the impact of the replacement of low-degree gravity field coefficients in GRACE/GRACE-FO solutions by SLR data on the trend and seasonal signals of ice mass changes in Greenland and Antarctica. We found that the replacement of the low-degree gravity field coefficients changes the estimates of trends by 4, 8, and 22 Gt/year in Greenland, West, and East Antarctica, respectively, depending on the source of SLR coefficients and period for which the coefficients are replaced. In SLR and GRACE solutions, all coefficients of the same order and the same parity of degrees are strongly correlated. Therefore, replacing only two selected coefficients may lead to a biased solution. Thus, we propose to combine GRACE with SLR solutions up to a degree and order 10 × 10 to properly consider the sensitivity of each of the techniques to gravity field coefficients, instead of replacing two coefficients from SLR in GRACE solutions. The combined solution reduces the residual trend of post-glacial rebound from 1.2 to 0.9 Gt/year and from − 57.8 to − 57.0 Gt/year in Scandinavia and South Canada, respectively, when compared to GRACE/GRACE-FO solutions with the replacement of coefficients. The SLR-GRACE combination reduces the noise in the GRACE/GRACE-FO solutions by 8%, from 38 to 35 Gt, in the Fennoscandia region. In the periods when GRACE is at the end of its mission and observations are disrupted, the weights adjust the contribution from SLR and GRACE based on relative ratio of variances from each techniques. Thus, the combined solutions are more consistent with independent geophysical models of glacial isostatic adjustment, and the combinations are affected by smaller noise than the standard GRACE solutions and properly account for different sensitivities of SLR and GRACE techniques to low-degree time-variable gravity field coefficients.

Thermokarst development is a topographic change in the landscape that is commonly associated with permafrost degradation in ice-rich permafrost regions. The Lena-Aldan interfluvial area in Central Yakutia in eastern Siberia has undergone extensive thermokarst development in the last three decades, particularly in the vicinity of settlements. Despite the negative effects of thermokarst development on the inhabitants of these settlements, no quantitative observation methods have been developed to investigate the surface displacement due to thermokarst development over entire towns. This study utilized interferometric synthetic aperture radar to reveal ground-surface displacement associated with thermokarst near the settlements of selected towns. The findings showed that significant subsidence was detected in disturbed areas (farming and abundant arable land) near the towns. The magnitude of subsidence in the Tyungyulyu and Mayya areas was less than that in Churapcha and Amga. Polygon density in a defined area in each town was examined using high-resolution optical images. The polygon density in Churapcha was considerably lower than that in Mayya, whereas polygonal texture in some areas in Tyungyulyu and Amga was unclear. Spatial frequency analysis using satellite optical images showed clear differences in averaged spectrum models between well-developed and less-developed polygons, which may reflect trough depths and density of vegetation between polygons. Satellite-based subsidence maps and statistics describing polygon development may be useful for evaluating both initial and subsequent stages of thermokarst development.

Seismic techniques using earthquakes are powerful tools for exploring the Earth’s internal structure. However, the earthquake distribution limits the spatial resolution. In recent years, ambient noise surface wave tomography using ambient seismic wave field has resolved these limitations. A typical ambient seismic wave field is microseisms excited by ocean swell activities. Ambient noise surface wave tomography is a technique in seismic interferometry that extracts seismic wave propagation between pairs of stations by cross-correlating the seismic records. The cross-correlation function can be interpreted as an impulsive response at a station with a virtual source at the other station. This technique became standard with the development of modern dense seismic networks. This paper reviews a theory of seismic interferometry for ambient noise surface wave tomography and procedures for practical data processing to calculate cross-correlation functions. The tomographic method typically consists of four steps: (1) the construction of reference 1-D models, (2) phase velocity measurements for each path, (3) 2-D phase velocity inversions, and (4) the construction of a 3-D S-wave tomographic model obtained from series of local 1-D inversions at all the grids. This paper presents the feasibility of multimode surface wave dispersion measurements for improving depth resolution.

Evaluating the mechanical properties of gas (primarily methane) hydrate-bearing sediments is essential for commercial production as a next-generation resource and understanding the global carbon cycle. Triaxial and uniaxial compression tests have been conducted on synthetic gas hydrate and natural core samples recovered from deep-sea beds using pressure coring techniques. The results show that four factors are vital in establishing the strength of hydrate-bearing sediments: hydrate saturation, effective confining stress, porosity, and strain rate. However, no study has evaluated these factors in a unified and quantitative manner, and even if the physical properties of the reservoir are known in detail from logging, predicting the strength has been challenging. In this study, pressure cores were drilled and recovered from the Eastern Nankai Trough in April 2018 after Japan’s second offshore production test, and triaxial or uniaxial compression tests were performed on 12 pressure core samples brought back to the laboratory. The mechanical properties of the hydrate-bearing sediments were classified with previous obtained results from 53 pressure cores and 223 synthetic cores, and empirical equations for triaxial compressive strength and deformation modulus were proposed as functions of gas hydrate saturation, effective confining pressure, porosity, and strain rate. The obtained equations were found to correlate well with the experimental data and can predict the strength and deformation modulus from logging data.

Abstract

The occurrence of earthquakes is now understood as brittle shear fracture releasing the elastic potential energy stored in the earth. Since the 1950s, many studies on the energy balance in earthquake faulting have been done, but there seems to be some incoherence among them. The essential reason is because various changes in conceptual framework happened during the last six decades, specifically the introduction of the new paradigm of plate tectonics in the 1960s, the concept of moment tensor as source representation in the 1970s, and the fault constitutive law governing rupture growth in the 1990s. Therefore, it will be worthwhile to reconsider the energetics of earthquake faulting from a current perspective. For this purpose, first of all, we summarize the basic concepts of elastic potential energy and moment tensor and review the general representation of earthquake sources and the origin of background crustal stress to confirm that the effect of earth’s self-gravitation is negligible in the energetics of shear faulting. Next, as a starting point for discussion, we directly derive a basic equation of mechanical energy balance in dynamic shear faulting from the equation of motion for an elastic body subjected to tectonic-origin deviatoric stress. Then, we review the widely accepted formula for indirectly evaluating radiated seismic energy from a simplified energy balance equation and compare with the direct evaluation based on the analytical solution of displacement fields for a point dislocation source in order to call attention to inconsistency between them. The inconsistency comes from the omission of the effects of rupture growth rate in the simplified energy balance equation. So, finally, we review the energy balance at the tips of a propagating shear crack, which naturally leads to the introduction of the slip-weakening fault constitutive law as a fundamental equation governing earthquake rupture. Then, we discuss the whole process of earthquake rupture, consisting of initiation, acceleration, steady propagation, deceleration, and termination from the viewpoint of energy balance.

The purpose of this study is to clarify the vertical prokaryotic distribution in groundwater in a terrestrial subsurface sedimentary environment with a geological complex. Six groundwater samples were collected from a coastal 1200-m-deep borehole in which digging strata deposited between 2.3 and 1.5 Ma in Horonobe, Hokkaido, Japan. The studied succession was divided into three vertical zones that were geochemically differentiated according to their chloride contents and water-stable isotopes. The upper zone (UZ; shallower than 500 m) primarily contained fresh water supplied by penetrating meteoric water, the connate water zone (CWZ; deeper than 790 m) contained paleo-seawater, and the diffusion zone (DZ; 500–790 m depth) located between UZ and CWZ. Fluctuations in the prokaryotic density and constituents were observed across these three zones. The prokaryotic density decreased from UZ toward DZ, and the density of DZ was two orders of magnitude lower than that of UZ and CWZ. High prokaryotic activity was observed in CWZ below DZ. The upward expansion of prokaryotic distribution from CWZ, where high prokaryotic potential expressed by biomass can be maintained almost equivalent to that in the marine environment, probably occurred on a geological timescale from 80 ka to 1.3 Ma, as shown by the groundwater age of DZ. The DZ is a zone where the geochemistry has changed drastically owing to the mixing of penetrating meteoric water and the diffusion of deep paleo-seawater, preserving a unique subsurface environment. This chemically mixed zone might be considered as a buffering zone for prokaryotes to prevent the expansion of prokaryotic density and activity provided by diffusion and their in situ growth from both above and below the zones, which is expected to be maintained over a geological timescale. This zone is considered important for using subsurface space in the deep subsurface environment of the island arc.