Deep-Seated Fluids in Thermal Waters Before and After the 2016 Kumamoto Earthquakes

Noble gases, oxygen-hydrogen isotope ratios, and ion compositions were measured at three sampling points (KUM, OTN, and ASO) from December 2013 to July 2021. The 3He/4He values at the three sampling points remained stable in the range of 3-4 Ra throughout the observation period, suggesting that the supply of deep-seated gases to the aquifer was stable. The 4He/20Ne values of KUM and OTN indicate that the supply of surface-source fluids to the aquifer decreased relative to that of deep-seated fluids at KUM and OTN. In contrast, in the ASO site, both the surface- and deep-seated fluids supplied to the aquifer were stable. The δD-δ18O relationship indicated the supply of deep-seated water to the KUM and OTN aquifers but not to the ASO aquifer. Nevertheless, the δD-δ18O relationship remained stable throughout the observation period, suggesting that the supply of deep-seated water to the three stations was stable. The Li/Cl and 1/Cl relationships for the three sampling points were plotted within a narrow range throughout the observation period, suggesting that the groundwater recharge was stable. Neither spikes nor step changes owing to the 2016 Kumamoto earthquake were observed in any of the data. These results indicate that the KUM and OTN aquifers are constantly supplied with deep fluids from the fluid-rich zone beneath the Kumamoto region, and that only deep-seated gas was supplied to the ASO aquifer. We also confirmed that these supply conditions were unaffected by the 2016 Kumamoto earthquake or the subsequent aftershock activity.

New insights on the active degassing system of the Lipari–Vulcano complex (South Italy) inferred from Local Earthquake Tomography

Seismic tomography is a very powerful and effective approach to look at depths beneath volcanic systems thus helping to better understand their behaviour. The P-wave and S-wave velocity ratio, in particular, is a key parameter useful to discriminate the presence of gas, fluids and melts. We computed the first 3-D overall model of Vp, Vs and Vp/Vs for the Lipari–Vulcano complex, central sector of the Aeolian volcanic archipelago (southern Italy). The investigated area has been characterized in recent times by fumaroles, hydrothermal activity and active degassing. In particular, in the Vulcano Island, several episodes of anomalous increases of fumarole temperature and strong degassing have been recorded in the past decades and the last “crisis”, started in September 2021, is still ongoing. For tomographic inversion we collected ~ 4400 crustal earthquakes that occurred in the last thirty years and we used the LOcal TOmography Software LOTOS. The results clearly depicted two low Vp and Vp/Vs anomalies located up to ~ 8 km depths below Vulcano and the western offshore of Lipari, respectively. These anomalies can be associated to the large presence of gas and they furnish a first picture of the gas-filled volumes feeding the main degassing activity of the area.

Seismic evidence of pop-up tectonics beneath the Shillong Plateau area of Northeast India

Our detailed 3-D seismic tomographic assimilation using high-quality phase arrival time data recorded by the local seismographic network demonstrated that heterogeneities in the crustal faults have contributed significantly to the pop-up tectonics beneath the Shillong Plateau, characterized by high-V and low-σ. The major seismogenic faults, namely, the north-dipping Dapsi thrust in association with Dauki fault in the south and south dipping Brahmaputra fault in the north, located either side of the Shillong Plateau that acted as the causative factors for the pop-up, which attributed to the lithostatic (high-V, low-σ) and sedimentary (low-V, high-σ) load, respectively. Seismicity is found confined to a depth ≤ 60 km. Uneven distribution of structural heterogeneities in the upper crust is responsible for earthquake genesis of varying strengths. It is intriguing to note that high-velocity anomalies and low-ϭ in the uppermost crust, interpreted as the Shillong Plateau that acted as a geometric asperity and the juxtaposition of high-V and low-V became the source zone of the 1897 Shillong earthquake (Ms 8.7) as a novel observation for the region. Structural heterogeneities are distinctly distributed between low-V, high-σ and high-V, low-σ in the lower crust plays a major role for future intense seismogenesis due to differential strain accumulation.

Thermal Monitoring of the Lithosphere by the Interaction of Deep Low-Frequency and Ordinary High-Frequency Earthquakes in Northeastern Japan

Deep low-frequency earthquakes (LFEs) are known to occur in dehydration phenomena from the subducting hydrous slab and in magmatic phenomena beneath Quaternary volcanoes in Japan. To realize the spatial and temporal characteristics of the magmatic deep low-frequency earthquakes, their hypocenters along with those of ordinary overhead high-frequency earthquakes are analyzed beneath six volcanic fields in northeastern Japan. This trial clarifies the rising basaltic magma conduits and rheological profiles of the lithosphere. Deep low-frequency earthquakes tend to form three vertical clusters corresponding to the rheological strength peak of the peridotite upper mantle, gabbroic lower crust, and granitic upper crust. Interactive aseismic gaps between low- and high-frequency earthquakes reveal the brittle–plastic transition as an isothermal indicator in the lithosphere. This relationship provides a tool to monitor the thermal evolution of the lithosphere and to explore sustainable geothermal resources with basaltic magma replenishment systems.

Fault stress inversion reveals seismogenic asperity of the 2011 Mw 9.0 Tohoku-Oki earthquake

We predict, with a model (earthquake stress model) that inverts the displacements documented at 163 GNSS onshore stations of the GEONET, the change of shear and normal stresses on the megathrust near the Japan Trench over the seven years before the 2011 Mw 9.0 Tohoku-Oki earthquake. We find three areas on the megathrust with greater accumulations of shear and normal stresses before the earthquake, which match the ruptured areas of the mainshock and two largest aftershocks (M_w 7.8 and 7.4) that occurred within half an hour after the mainshock. We also find that the change of normal stress on the fault before the earthquake is not uniform but increases in the up-dip portion (shallower depth) of the fault from the hypocenter and decreases in the down-dip portion. We infer that the occurrence of the giant earthquake at the shallow portion of the megathrust may be attributed to the increase of the normal stress there, which leads to an increase of fault shear strength and allows more elastic strain energy to accumulate to prepare for the next big earthquake. Based on these results we propose a new concept of the seismogenic asperity as the area of greater accumulations of shear and normal stresses. The method presented here may be useful for predicting the rupture zone of future large earthquakes.

Tremor activity inhibited by well-drained conditions above a megathrust

Tremor occurs on megathrusts under conditions of near-lithostatic pore-fluid pressures and extremely weakened shear strengths. Although metamorphic reactions in the slab liberate large amounts of fluids, the mechanism for enhancing pore-fluid pressures along the megathrust to near-lithostatic values remains poorly understood. Here we show anti-correlation between low-frequency earthquake (LFE) activity and properties that are markers of the degree of metamorphism above the megathrust, whereby LFEs occur beneath the unmetamorphosed overlying plate but are rare or limited below portions that are metamorphosed. The extent of metamorphism in the overlying plate is likely controlled by along-strike contrasts in permeability. Undrained conditions are required for pore-fluid pressures to be enhanced to near-lithostatic values and for shear strength to reduce sufficiently for LFE generation, whereas well-drained conditions reduce pore-fluid pressures at the megathrust and LFEs no longer occur at the somewhat strengthened megathrust. Our observations suggest that undrained conditions are a key factor for the genesis of LFEs.

Seismic Evidence for Deep-Water Transportation in the Mantle

We report seismic evidence for the transportation of water into the deep mantle in the subduction zone beneath northeastern Japan. Our data indicate that water is released from the hydrated oceanic crust at shallow depths (< approximately 100 kilometers) and then forms a channel of hydrated mantle material on top of the subducting plate that is the pathway for water into the deep mantle. Our result provides direct evidence that shows how water is transported from the ocean to the deep mantle in a cold subduction zone environment.

The role of fluids in lower-crustal earthquakes near continental rifts

The occurrence of earthquakes in the lower crust near continental rifts has long been puzzling, as the lower crust is generally thought to be too hot for brittle failure to occur. Such anomalous events have usually been explained in terms of the lower crust being cooler than normal. But if the lower crust is indeed cold enough to produce earthquakes, then the uppermost mantle beneath it should also be cold enough, and yet uppermost mantle earthquakes are not observed. Numerous lower-crustal earthquakes occur near the southwestern termination of the Taupo Volcanic Zone (TVZ), an active continental rift in New Zealand. Here we present three-dimensional tomographic imaging of seismic velocities and seismic attenuation in this region using data from a dense seismograph deployment. We find that crustal earthquakes accurately relocated with our three-dimensional seismic velocity model form a continuous band along the rift, deepening from mostly less than 10 km in the central TVZ to depths of 30-40 km in the lower crust, 30 km southwest of the termination of the volcanic zone. These earthquakes often occur in swarms, suggesting fluid movement in critically loaded fault zones. Seismic velocities within the band are also consistent with the presence of fluids, and the deepening seismicity parallels the boundary between high seismic attenuation (interpreted as partial melt) within the central TVZ and low seismic attenuation in the crust to the southwest. This linking of upper and lower-crustal seismicity and crustal structure allows us to propose a common explanation for all the seismicity, involving the weakening of faults on the periphery of an otherwise dry, mafic crust by hot fluids, including those exsolved from underlying melt. Such fluids may generally be an important driver of lower-crustal seismicity near continental rifts.

Regional crustal thickness and precipitation in young mountain chains

Crustal thickness is related to climate through precipitation-induced erosion. Along the Andes, the highest mountains and thickest crust (approximately 70 km) occur at 25 degrees south, a region of low precipitation. Westerly winds warm passing over the Atacama Desert; precipitation is modest in the High Andes and eastward over the Altiplano. Severe aridity, hence low erosion rates, helps to account for the elevated volcanogenic contractional arc and high, internally draining plateau in its rain shadow. Weak erosion along the north-central arc provides scant amounts of sediment to the Chile-Peru Trench, starving the subduction channel. Subcrustal removal might be expected to reduce the crustal thickness, but is not a factor at 25 degrees south. The thickness of the gravitationally compensated continental crust cannot reflect underplating and/or partial fusion of sediments, but must be caused chiefly by volcanism-plutonism and contraction. Contrasting climate typifies the terrain at 45 degrees south where moisture-laden westerly winds encounter a cool margin, bringing abundant precipitation. The alpine landscape is of lower average elevation compared with the north-central Andes and is supported by thinner continental crust (approximately 35 km). Intense erosion supplies voluminous clastic debris to the offshore trench, and vast quantities are subducted. However, the southern Andean crust is only about half as thick as that at 25 degrees south, suggesting that erosion, not subcrustal sediment accretion or anatexis, is partly responsible for the thickness of the mountain belt. The Himalayas plus Tibetan Plateau, the Sierra Nevada plus Colorado Plateau, and the Japanese Islands exhibit analogous relationships between crustal thickness and climate.