Stable isotopes show that earthquakes enhance permeability and release water from mountains

Inorganic Hydrogeochemistry in the 21st Century

Chemical and isotopic processes occur in every segment of the hydrological cycle. Hydrogeochemistry-the subdiscipline that studies these processes-has seen a transformation from "witch's brew" to credible science since 2000. Going forward, hydrogeochemical research and applications are critical to meeting urgent societal needs of climate change mitigation and clean energy, such as (1) removing CO2 from the atmosphere and storing gigatons of CO2 in soils and aquifers to achieve net-zero emissions, (2) securing critical minerals in support of the transition from fossil fuels to renewable energies, and (3) protecting water resources by adapting to a warming climate. In the last two decades, we have seen extensive activity and progress in four research areas of hydrogeochemistry related to water-rock interactions: arsenic contamination of groundwater; the use of isotopic and chemical tracers to quantify groundwater recharge and submarine groundwater discharge; the kinetics of chemical reactions and the mineral-water interface's control of contaminant fate and transport; and the transformation of geochemical modeling from an expert-only exercise to a widely accessible tool. In the future, embracing technological advances in machine learning, cyberinfrastructure, and isotope analytical tools will allow breakthrough research and expand the role of hydrogeochemistry in meeting society's needs for climate change mitigation and the transition from fossil fuels to renewable energies.

Detecting groundwater level changes related to the 2016 Kumamoto Earthquake

This study presented the first attempt to detect precursory changes in groundwater level before the 2016 Kumamoto Earthquake. This detection was achieved by accurately determining the relationship between long-term groundwater level fluctuation and crustal deformation over 16 years through analysis of groundwater level time-series data acquired at 17 sites within the study area. Here, we show that the observed groundwater levels were lower than the modelled levels in aquifers composed of porous strata (Togawa lava and part of the pre-Aso volcanic rocks), and that there were larger differences until 2014, which diminished until the occurrence of the Kumamoto Earthquake. The initial reduction in the modelled groundwater level and the latter recovery were most likely caused by crustal strain relaxation associated with the large 2011 earthquake off the Pacific coast of Tohoku (Mw 9.0) and the strain accumulation prior to the 2016 Kumamoto Earthquake.

The 222Rn and CO2 soil gas distribution at Lembang Fault Zone, West Java - Indonesia

Series of ²²² Rn measurements had been employed around the Lembang Fault Zone, West Java - Indonesia to investigate its relationship with the fault zone. Furthermore, the mechanism of ²²²Rn transport is also investigated by performing soil gas CO₂ and its δ¹³C measurements. The evaluation of ²²²Rn concentration shows that 34% data are within threshold values with an average of 14,117 Bq/m³, with anomalous concentration greater than 20,000 Bq/m³. The concentration of CO₂ in soil gas is varied from 72 ppm to 13,241 ppm and consisted of three populations, with 40% of the data above 655 ppm. The spatial distribution pattern shows that most of Lembang Fault Zone segment coincides with high ²²²Rn concentration indicating high permeability zone. Furthermore, the average ²²²Rn concentration in western part of the fault is higher than the eastern part and this may be correlated with higher seismic activities. In contrast to ²²²Rn, CO₂ concentration shows less correlation to the fault structure. Based on δ¹³C values, the source of soil CO₂ is dominated by atmospheric CO₂, with minor mixing of biogenic origin. Although Lembang Fault Zone is located in the south of an active volcano, there is no indication of magmatic CO₂ origin. The lack of correlation between ²²²Rn and CO₂ suggests that there is no indication of ²²²Rn transport by CO₂ as carrier gas.

Self-organizing map improves understanding on the hydrochemical processes in aquifer systems

The holistic understanding of hydrochemical features is a crucial task for management and protection of water resources. However, it is challenging for a complex region, where multiple factors can cause hydrochemical changes in studied catchment. We collected 208 groundwater samples from such region in Kumamoto, southern Japan to explicitly characterize these processes by applying machine learning technique. The analyzed groundwater chemistry data like major cations and anions were fed to the self-organizing map (SOM) and the results were compared with classical classification approaches like Stiff diagram, standalone cluster analysis and score plots of principal component analysis (PCA). The SOM with integrated application of clustering divided the data into 11 clusters in this complex region. We confirmed that the results provide much greater details for the associated hydrochemical and contamination processes than the traditional approaches, which show quite good correspondence with the recent high resolution hydrological simulation model and aspects from geochemical modeling. However, the careful application of the SOM is necessary for obtaining accurate results. This study tested different normalization approaches for selecting the best SOM map and found that the topographic error (TE) was more important over the quantization error (QE). For instance, the lower QE obtained from min-max and log normalizations showed problems after clustering the SOM map, since the QE did not confirm the topological preservation. In contrast, the lowest TE obtained from Z-transformation data showed better spatial matching of the clusters with relevant hydrochemical characteristics. The results from this study clearly demonstrated that the SOM is a helpful approach for explicit understanding of the hydrochemical processes on reginal scale that may capably facilitate better groundwater resource management.

Hydrogeochemical Characteristic of Geothermal Water and Precursory Anomalies along the Xianshuihe Fault Zone, Southwestern China

Hydrogeochemical changes in association with earthquakes are considered as a potential means of identifying earthquake precursors. The Xianshuihe fault zone (XSHF) is considered one of the most active seismic fault zones in China; 43 hot springs were sampled and analysed in the laboratory for major elements, silica, stable isotopes (δD and δ18O) and strontium isotopes were investigated from 2008 to 2021. The meteoric water acted as the primary water source of the hot spring in the XSHF, and recharged elevations ranged from 1.9 to 4.8 km. The geothermometers method was used to estimate the region of thermal storage temperature and its temperature between 8 and 142 °C. And the circulation depth ranged from 0.1 to 6.9 km. Most of the hot spring water was immature water with a weak degree of water-rock reaction. However, the degree of water-rock reaction and the depth of hot spring water circulation were high in part of the Kangding and Daufu segments, which also had the highest reservoir temperature and the most frequent strong earthquakes. Temporal variations of hydrogeochemical showed that Na+, Cl− and SO42− decreased obviously following the 12 May 2008 Wenchuan Ms8.0 and existed abnormal value fluctuations from the 20 April 2013 Lushan Ms7.0 to 22 November 2014 Kangding Ms6.3 occurred and after 20 July 2017 returned to the normal levels. And the ion concentrations in hot springs increased by 5% to 35% three months before 22 November 2014 Kangding Ms6.3 with the obvious precursor anomaly. Hydrogeochemical anomalies could be useful for predicting an earthquake in the study area.

Hydrochemical Characteristics of Earthquake-Related Thermal Springs along the Weixi–Qiaohou Fault, Southeast Tibet Plateau

The Weixi–Qiaohou Fault (WQF) is considered an important zone of the western boundary of the Sichuan–Yunnan block, and its seismicity has attracted much attention after a series of moderate–strong earthquakes, especially the Yangbi Ms6.4 earthquake that occurred on 21 May 2021. In the present research, we investigate major and trace elements, as well as hydrogen and oxygen isotopes, of 10 hot springs sites located along the WQF, which are recharged by infiltrated precipitation from 1.9 to 3.1 km. The hydrochemical types of most analyzed geothermal waters are HCO3SO4-Na, SO4Cl-NaCa, and SO4-Ca, proving that they are composed of immature water and thus are characterized by weak water–rock reactions. The heat storage temperature range was from 44.1 °C to 101.1 °C; the circulation depth was estimated to range between 1.4 and 4.3 km. The results of annual data analysis showed that Na+, Cl−, and SO42− in hot springs decreased by 11.20% to 23.80% north of the Yangbi Ms5.1 earthquake, which occurred on 27 March 2017, but increased by 5.0% to 28.45% to the south; this might be correlated with the difference in seismicity within the fault zone. The results of continuous measurements of NJ (H1) and EYXX (H2) showed irregular variation anomalies 20 days before the Yangbi Ms6.4 earthquake. In addition, Cl− concentration is more sensitive to near-field seismicity with respect to Na+ and SO42−. We finally obtained a conceptual model on the origin of groundwater and the hydrogeochemical cycling process in the WQF. The results suggest that anomalies in the water chemistry of hot spring water can be used as a valid indicator of earthquake precursors.

Estimation of Groundwater Recharge in Kumamoto Area, Japan in 2016 by Mapping Land Cover Using GIS Data and SPOT 6/7 Satellite Images

Agricultural fields, grasslands, and forests are very important areas for groundwater recharge. However, these types of land cover in the Kumamoto area, Japan, were damaged by the Kumamoto earthquake and heavy rains in 2016. In this region, where groundwater provides almost 100% of the domestic water supply for a population of about 1 million, quantitative evaluation of changes in groundwater recharge due to land cover changes induced by natural disasters is important for the sustainable use of groundwater in the future. The objective of this study was to create a land cover map and estimate the groundwater recharge in 2016. Geographic information system (GIS) data and SPOT 6/7 satellite images were used to classify the Kumamoto area into nine categories. The maximum likelihood classifier of supervised classification was applied in ENVI 5.6. Eventually, the map was cleaned up with a 21 × 21 kernel filter, which is larger than the common size of 3 × 3. The created land cover map showed good performance of the larger filter size and sufficient validity, with overall accuracy of 91.7% and a kappa coefficient of 0.88. The estimated total groundwater recharge amount reached 757.56 million m3. However, if areas of paddy field, grassland, and forest had not been reduced due to the natural disasters, it is estimated that the total groundwater recharge amount would have been 759.86 million m3, meaning a decrease of 2.30 million m3 in total. The decrease of 2.13 million m3 in the paddy fields is temporary, because the paddy fields and irrigation channels have been improved and the recharge amount will recover. On the other hand, since the topsoil on the landslide scars will not recover easily in natural conditions, it is expected to take at least 100 years for the groundwater recharge to return to its original state. The recharge amount was estimated to decrease by 0.17 million m3 due to landslides. This amount is quite small compared to the total recharge amount. However, since the reduced recharge amount accounts for the annual water consumption for 1362 people, and 12.1% of the recharge decrease of 1.41 million m3 each year to fiscal year 2024 is expected by municipalities, we conclude that efforts should be made to compensate for the reduced amount due to the disasters.

Changes in groundwater trace element concentrations before seismic and volcanic activities in Iceland during 2010-2018

We analysed temporal variations of trace element concentrations in groundwater from a 101 m-deep borehole (HA01) in northern Iceland during 2010-2018 and compared them with seismic and volcanic events that occurred in the same period to identify potential hydrogeochemical precursors. An increase of B, Al, V, Li and Mo concentrations started from eight months to one month before the 2014 Bárðarbunga eruption (~115 km from HA01), a major rifting event in central Iceland, while Ga and V concentrations began to increase one day and one month after the onset of the event, respectively. We also found that concentrations of some trace elements (Li, B, Ga, Mo, Sr, Rb and Fe) significantly increased before an M_w 5.0 earthquake that occurred ~80 km from the borehole in 2018. However, other notable hydrogeochemical changes were detected during the monitoring period without apparent correlation with the seismic and volcanic events in the region. This study shows that the systematic long-term hydrogeochemical monitoring in seismic and volcanic areas is critical to advance the science of seismic and eruptive precursors. Furthermore, the use of statistical tools, such as Principal Component Analysis (PCA) and Change Point (CP) detection can help identify the most useful chemical elements and validate the trend variability of those elements in the time series, reducing arbitrary choices of pre-seismic and pre-volcanic hydrogeochemical anomalies as potential precursors.

Acesulfame as a suitable sewer tracer on groundwater pollution: A case study before and after the 2016 Mw 7.0 Kumamoto earthquakes

On April 14th and 16th, 2016, two large-scale earthquakes (M_w 6.2 and 7.0) occurred in Kumamoto, Japan. The sewer system was seriously damaged and there were concerns about groundwater pollution by sewer exfiltration. In this study, artificial sweeteners including acesulfame (ACE) in groundwater were analyzed before and after the earthquakes to evaluate sewage pollution and its temporal variation. Before the earthquakes, ACE was detected in 31 of 49 groundwater samples analyzed, indicating that wastewater may have leaked into groundwater. Groundwater was sampled from the same locations 2, 7, 12, and 30 months after the earthquakes. The detection frequency and median concentration of ACE in groundwater increased significantly 7 months after the earthquakes, from several tens to maximumly 189 times greater than the pre-earthquake concentrations. This suggests the earthquakes caused serious damage to sewer pipes and groundwater may be polluted. However, ACE concentrations drastically decreased or remained low 30 months after the earthquakes, probably due to the recovery and restoration work of sewer infrastructure. This study shows that ACE is an excellent tracer for evaluating sewer exfiltration to groundwater. In addition, it is important to obtain data on sewage tracers under normal condition as part of preparations for large-scale earthquakes.

Effects of the Japanese 2016 Kumamoto Earthquake on Nitrate Content in Groundwater Supply

The 2016 Kumamoto earthquake had a significant impact on groundwater levels and quality. In some areas, the groundwater level increased significantly due to the release of groundwater from upstream mountainous regions. Conversely, the groundwater level in other areas greatly decreased due to the creation of new fracture networks by the earthquake. There were also significant changes in certain groundwater quality variables. In this study, we used clustering based SOM (self-organizing maps) analysis to improve the understanding of earthquake effects on groundwater quality. We were especially interested in effects on groundwater used for drinking purposes and in nitrate concentration. For this purpose, we studied groundwater nitrate (NO3− + NO2−–N) concentrations for the period 2012–2017. Nitrate concentration changes were classified into seven typical SOM clusters. The clusters were distributed in three representative geographical regions: a high concentration region (>4 mg/L), a low concentration region (<1.6 mg/L) with minimal anthropogenic loading area, and an intermediate concentration region (2–4 mg/L). Depending on these regions, the nitrate concentration changes just before and after the earthquake had both increasing and decreasing trends between 2015–2017. This points to complex physiographical relationships for release of stored upstream groundwater, promotion of infiltration of shallow soil water/groundwater, and nitrate concentration as affected by earthquakes. We present an analysis of these complex relationships and a discussion of causes of nitrate concentration changes due to earthquakes.