Abstract

Hydrogeochemistry of 32 hot springs in the western Sichuan Province after the Wenchuan 8.0 earthquake was investigated by analyzing the concentrations of cation and anion and the isotopic compositions of hydrogen and oxygen. The water samples of the hot springs were collected four times from June 2008 to April 2010. Hydrogeochemical data indicated the water samples can be classified into 9 chemical types. Values of δD and δ18O indicated that the spring waters were mainly derived from meteoric precipitation and affected by water-rock interaction and mixture of deep fluids. Concentrations of K+ and of the samples from the Kangding district exhibited evident increases before the Wenchuan earthquake, indicating more supplement of deep fluids under the increase of tectonic stress. The chemical and isotopic variations of the water samples from the area closer to the epicenter area can be attributed to variation of regional stress field when the aftershock activities became weak.

1. Introduction

Great earthquakes usually associate with the physical-chemical variations of groundwater and hot springs. The short-term hydrogeochemical precursors for earthquake [14], co-seismic response of hydrochemistry [5, 6] and postseismic geochemical and isotopic changes of hot springs have been reported [712].

The observed geochemical anomalies related to earthquakes are usually attributed to the alteration of groundwater in the circulating system under the action of increasing crustal stress before and after earthquakes [6, 7, 13, 14]. For understanding the hydrogeochemical anomalies related to earthquakes, some genetic mechanism models have been proposed, such as increased solubility of rocks under increased pressure and release of ions from rocks into water [15], pore collapse with fluids expulsion [8], the water-rock interactions at the enhanced reactive surfaces [16], aquifers rupture with mixing of different fluids [6] and expulsion of deep fluids by tectonic pumping [1].

Hydrogeochemical survey of hot springs in western Sichuan before the 8.0 Wenchuan earthquake mainly focused on the origin of waters, the pollutions, the chemical classifications and the heat reservoirs in the sites surrounding Kangding county [1725].

The 8.0 earthquake occurred in Wenchuan county, Sichuan province, Southwest China on 12 May, 2008, following by hundreds of aftershocks with magnitudes higher than 3.0.

In order to investigate the hydrogeochemical characteristics of the springs related to the seismic activities, the hydrogeochemical survey of 32 hot springs in the western Sichuan were performed four times from June 2008 to April 2010.

2. Geological Setting

The investigated area is located at the eastern margin of the Tibetan Plateau, in which there are four major fault zones named as Minjiang fault (MJF), Longmenshan fault (LMSF), Xianshuihe fault (XSHF) and Anninghe fault (ANHF) zones (Figure 1) where earthquakes frequently occurred [26].Triassic littoral-neritic clastic rocks interbeded with carbonate and intrusive granite are exposed in the west of the LMSF and ANHF zones, but the strata from Late Paleozoic to Cenozoic are widespread exposed in the LMSF and ANHF zones and on the east of the fault zones, and granite is found in the intersection of the LMSF, XSHF and ANHF zones [2730]. The fault zones act as the important passage for upward migration of thermal fluids from the deep earth, which is indicated by distribution of many hot springs in the MJF, LMSF, XSHF and ANHF zones (Figure 1). Historically, a number of great earthquakes have occurred in the investigated area since 1800. For instance, the 1850 7.5 Xichang earthquake occurred in the ANHF zone, two 7.2 earthquakes of 16 and 23 August 1976 in the MJF zone, the 7.5 one of 25 August 1933 in the LMSF zone, and the 7.5 one of 1955 and the 7.9 one of 1973 in the XSHF zone [12]. The 8.0 Wenchuan earthquake resulted in a 285 km surface rupture zone along the pre-existing Yingxiu-Beichuan, Guanxian-Anxian and Qingchuan faults, with the maximum vertical surface offset of about 6.2 m [31].

3. Methods

Water samples were repeatedly collected four times (in June and October 2008, June 2009, and April 2010) at 32 sites of spas, wells, and springs (Figure 1) in the MJF, LMSF, XSHF, and ANHF zones. The samples were sealed and stored in 500 mL glass bottles. The values of temperature were measured with a digital thermometer with an error of ±1% in the field. Isotopic compositions of δD and δ18O were measured with a Picarro L1102 mass spectrometer in the Laboratory of Gas Geochemistry, Institute of Geology and Geophysics, Chinese Academy of Sciences in Lanzhou, and the errors are 0.5 for δD and 0.1 for δ18O, respectively. The concentrations of cations (K+, Na+, Mg2+, and Ca2+) and anions (F, Cl, Br, , and ) were determined with the Dionex ICS-900 ion chromatography (reproducibility within ±2%) that is installed in the Institute of Earthquake Science, China Earthquake Administration. The and concentrations were measured by the standard titration procedures with a ZDJ-100 potentiometric titrator (reproducibility within ±2%). The ion balance (ib) was calculated according to the following equation [11]:

4. Results

The physicochemical parameters of the 32 springs were listed in Table 1. No water samples were collected from the springs of nos. 2 and 24–26 in June 2008 because the springs were damaged by the Wenchuan earthquake. No samples were collected from the springs nos. 8 and 29 in April 2010 and the spring no. 28 in June 2009 and April 2010 because of the postearthquake reconstruction. The ion balance values of measured chemical data were calculated to be less than 5% (Table 1).

4.1. Physicochemical Parameters of the Springs in the LMSF Zone

The water samples from in the LMSF zone have δ18O and δD values between −16.34 and −6.18, −116.65 and −61.66, respectively. The temperatures and TDS ranged from 10.5 to 58.0°C and 151.69 to 1569.81 mg/L. The concentrations of Na+, Ca2+ and Mg2+ varied from 29.82 to 202.49 mg/L, 4.74 to 252.78 mg/L, and 0.80 to 128.73 mg/L, respectively. The concentrations of Cl, , and varied from 2.31 to 64.74 mg/L, 21.33 to 878.08 mg/L, and 14.63 to 310.70 mg/L respectively (Table 1).

4.2. Physicochemical Parameters of the Springs in the MJF Zone

Fifteen water samples were collected from the sites nos. 6–9 in the MJF zone, of which the δ18O and δD were between −15.34% and −13.20%, −112.23% and −102.20%, the values of temperatures and TDS ranged from 8.8 to 21.7°C and 505.34 to 1657.61 mg/L. The concentrations of Na+, Ca2+, and Mg2+ ranged from 8.75 to 25.23 mg/L, 49.54 to 183.57 mg/L, and 23.23 to 183.57 mg/L, and the concentrations of Cl, , and from 0.93 to 10.78 mg/L, 3.56 to 17.78 mg/L, and 313.61 to 1324.27 mg/L separately (Table 1).

4.3. Physicochemical Parameters of the Springs in the XSHF Zone

Sixty-two water samples were collected from the sites nos. 10–25 in the XSHF zone, of which δ18O and δD varied from −18.84 to −10.04 and −141.98 to −78.77. The values of temperatures and TDS ranged from 27.8 to 83.0°C and 427.53 to 2159.43 mg/L, respectively. The concentrations of Na+, Ca2+, and Mg2+ ranged from 35.14 to 679.23 mg/L, 4.38 to 50.10 mg/L, and 0.12 to 45.71 mg/L, respectively. The concentrations of Cl, , and ranged from 2.35 to 336.16 mg/L, 6.56 to 161.87 mg/L, and 62.48 to 1894.64 mg/L, respectively (Table 1).

4.4. Physicochemical Parameters of the Springs in the ANHF Zone

Twenty-four water samples were collected from the sites nos. 26–32 in the ANHF zone, of which the δ18O and δD were between −11.64 and −96.10, −110.47 and −102.20. The temperatures of the spring waters were between 21.1 and 56.9°C, the TDS were between 206.70 and 1212.86 mg/L. The concentrations of Na+, Ca2+, and Mg2+ ranged from 50.17 to 213.50 mg/L, 2.00 to 81.48 mg/L, and 0.00 to 29.43 mg/L, and the concentrations of Cl, , and from 1.91 to 136.74 mg/L, 10.94 to 257.26 mg/L, and 81.95 to 449.06 mg/L, respectively (Table 1).

5. Discussion

5.1. Chemical Types of the Waters and Their Origins
5.1.1. The Chemical Types

The triangular diagrams in Figure 2 illustrated the proportions of the anionic (Cl, , and ) and cationic (K+ + Na+, Mg2+, and Ca2+) concentrations of the water samples from the 32 sites. The anions and cations of the water samples were mainly distributed in blocks i, ii, and iii, indicating that Na+, Ca2+, , and were the main chemical composition for most of the samples. Based on the Shoka Lev’s classification method, the water samples were classified into 9 chemical types according to the main chemical composition, which are Na(Ca)-HCO3(SO4), Na(Mg)-HCO3(SO4), Ca(Na)-HCO3(SO4), Mg(Ca)-SO4, Ca(Mg)-SO4, Ca(Mg)-HCO3, Mg(Ca)-HCO3, Na-Cl(HCO3), and Ca(Na)-SO4(HCO3) (Table 1). The cations of the samples of the springs nos. 1–5 occur in the LMSF zone were mainly distributed in the blocks 1, 2, and 4, and the anions were distributed in the blocks i and ii, which together formed 4 chemical types of Na(Ca)-HCO3(SO4), Ca(Na)-HCO3(SO4), Mg(Ca)-SO4, and Ca(Mg)-SO4. The cations of the samples of the springs nos. 6–9 occur in the MJF zone were mainly distributed in the blocks 3 and 4, the anions were distributed in the block i, and the chemical types for the samples were Ca(Mg)-HCO3 and Mg(Ca)-HCO3. Both of the cations and anions of the samples of the springs nos. 26–32 occur in the ANHF zone were distributed in the blocks i-ii and 1-2, these samples have the similar chemical type of Na(Ca)-HCO3(SO4). Half of the samples including the samples of the springs nos. 10–25 occur in the XSHF zone, the cations of these samples were distributed in the blocks 1, 2, and 5, with anions in the blocks i, ii, iii, and vi, and the chemical types were Na(Ca)-HCO3(SO4), Na(Mg)-HCO3(SO4), Na-Cl(HCO3), and Ca(Na)-SO4(HCO3).

5.1.2. The Origins

(1) δD and δ18O. The values of δ18O and δD of the 32 springs were plotted along the local meteoric water line (LMWL: δD = 7.06 δ18O + 4.06) with different extent shifts from the LMWL (Figure 3), which indicated that the spring waters were mainly originated from meteoric water and with different extent alteration. The values of δ18O and δD were plotted in three regions, the region I includes the samples of the springs no.1–3, with values of δD and δ18O ranged from −82.13 to −61.66 and −11.75 to −6.18, the region III includes all the samples from the XSHF zone (springs nos. 10–25), with values of δD and δ18O between −141.98 and −78.77, −18.84 and −10.04, respectively, and the others in the region II, with the values of δD and δ18O ranging from −116.81 to −96.10 and −16.34 to −11.64 (Table 1). The values of δD and δ18O of the springs nos.10–25 from the higher mountain area were more negative, while those of the springs nos. 1–3 collected in the lower altitude region were less negative, which is constant with the previous results [32].

(2) Na(Ca)-HCO3(SO4), Na(Mg)-HCO3(SO4), and Na-Cl(HCO3) Waters. The water samples of the springs nos. 4-5, 10–23 and 25–32 were belong to these types (Table 1). The springs nos. 4-5, 13–23 and 25 occur in granite, while the springs nos. 10-11 in Late Triassic clastic rock interbeded with volcanic rocks. The spring no. 12 is found in Late Permian volcanic tuff, while the spring no. 26 in Sinian pyroclastic rock. The spring nos. 27–29 and 31 are found in Cretaceous clastic rock, the spring no. 30 in Early Triassic sandstone, and the spring no. 32 in Paleozoic mixed layer (Figure 1). The springs nos. 4-5, 10–23 and 25–32 occur in the rocks enriched in feldspar resulting in similar higher concentrations of Na+ and because of rock-water interaction. In addition, the samples from the springs nos. 5, 19, 23, and 32 were characterized by the similar concentrations of Ca2+, Mg2+, and (Table 1), which should be attributed to the water-rock interactions between groundwater and Devonian carbonate (Figure 1). Cl is known to be conservative and derive from the deep earth mainly [33, 34]. The chemical type for the samples from the spring no. 16 was Na-Cl (HCO3), with the Cl concentration as 328.54 mg/l (Table 1), which suggested the input of deep fluid. Meanwhile, the high 3He/4He ratio (between 1.43 and 3.73RA, RA = 1.39 × 10−6) [30] together with the high temperatures (between 80.0°C and 70.2°C) also are the evidences for the upwelling of the deep-earth fluids into the spring (Table 1).

(3) Ca(Na)-HCO3(SO4), Ca(Mg)-HCO3, Mg(Ca)-HCO3, and Ca(Na)-SO4(HCO3) Waters. These included the samples from the springs nos.1, 6–9, and 24 (Table 1). The spring no. 1 occurs in Jurassic carbonate interbedded with sandstone. The springs nos. 6 and 9 occur in Late Triassic clastic rock interbedded with carbonate. The springs nos. 7 and 8 occur in Middle Triassic clastic rock interbedded with carbonate, and the spring no. 24 is found in Carboniferous carbonate (Figure 1). The springs nos. 1, 6–9, and 24 occur in strata enriched in carbonate, so that the main components of the waters were Ca2+, Mg2+ and because of the interaction between the groundwater and carbonate. The samples of the springs nos. 8 and 9 have the same chemical type of Mg (Ca)-, and similar Ca/Mg ratios (0.54 and 0.86, Table 1), which can be attributed to the dissolution of calcite combined with the low temperatures of the springs nos. 8 and 9 (8.8°C and 9.0°C respectively) (Table 1).

(4) Mg(Ca)-SO4 and Ca(Mg)-SO4 Waters. The springs nos. 2 and 3 occur in Devonian to Silurian pyrite-bearing sale and marlite (Figure 1), from which the samples had higher concentration (878.08 mg/L and 812.50 mg/L, resp.) and TDS (1571.81 mg/L and 1195.08 mg/L, resp.), with the chemical types of Mg (Ca)-SO4 and Ca (Mg)-SO4 (Table 1), indicating interaction between carbonate and groundwater enhanced by the sulfate from the oxidation of pyrite [35].

5.2. Chemical Variations of the Spring Waters Associated with Seismic Activities
5.2.1. Temporal Variations

The hydrochemical data of the water samples showed evident trend of temporal variations.

The geochemical parameters of the samples nos. 1, 3–5 from the LMSF zone and the samples nos. 15 and 20 from the XSHF zone varriedby more than 20% one month after the Wenchuan earthquake, and similar variations happened after October 2008. The concentrations of Na+, Ca2+, Cl, , and TDS of the sample of the spring no. 1 decreased by 42%, 39%, 53%, 31%, and 35%, respectively; for the sample number 2, the decreasing amplitudes for the concentrations of Na+, Mg2+, Ca2+, Cl, , and TDS were 23%, 45%, 25%, 38%, 22%, and 33%, separately. The concentrations of Mg2+, Ca2+, , and TDS of the sample from the spring no. 3 decreased by 33%, 42%, 95%, and 35%, respectively. The concentrations of Na+, , and TDS of the sample of the spring no. 4 decreased by 67%, 33% and 57%, respectively. For the sample of the spring no. 5, the concentrations of Na+, and TDS decreased by 25%, 41%, and 26%, respectively. The concentrations of Na+ and of sample of the spring no. 20 decreased by 28% and 26%, respectively. For the sample of the spring no. 15, however the concentrations of Na+, and TDS increased in 17 April 2010, with the amplitudes of 53%, 38%, and 38%, respectively. Meanwhile, the δ18O-δD plots of the samples from the springs nos.1, 3, and 4 approached the LMWL (Local Meteoric Line) after sizeable shift in June 2008, and the similar variation happened for the spring no. 2 after October 2008 (Figure 3), which indicated gradual input decrease of deep-earth fluids that formed by water-rock reaction after the Wenchuan earthquake

The supplement of deep water characterized by higher mineralization and enriched in δ18O is considered as a main factor for controlling pre- and co-seismic geochemical variations of the groundwater and springs [2, 6, 7, 11, 13, 14, 16, 36, 37]. The Coulomb stress in the middle-north segment of the LMSF zone, the southeast segment of the XSHF zone, and the south segment of the MJF zone had been enhanced, which resulted in the Wenchuan 8.0 earthquake [38, 39], and the geochemical variations of the samples from the springs nos. 1–5 and 20. Meanwhile, the input decrease of deep-earth water during the sampling period in which seismic activity decreased gradually. 747 aftershocks with higher than 4.0 occurred in the three fault zones, and most of which (including all the 12 events with amplitude ranged from 6.0 to 6.4) occurred before September 10, 2008 (China Earthquake Network Center, Figure 1). As the aftershock activities became weak, the TDS values of the samples from the springs nos. 1–5 and 20 decreased gradually to the normal value (Figure 4), indicating the supplement decrease of the deep-earth fluid as the Coulomb stress being released after the events as indicated by the data of oxygen and hydrogen isotope compositions (Figure 3).

Notably, the concentrations of Na+, , and TDS of the sample of the spring no. 15 increased by 53%, 38%, and 38% on April 17, 2010 (Figure 4), which was closely related to the 28 April 2010 an 5.4 aftershock 15 kilometers epicenter distance to the spring no. 15 (Figure 1). Therefore, it can be infered that the hydrogeochemical variations of the spring no. 15 might be the short-term hydrogeochemical precursors for the 5.4 event.

The geochemical parameters of the samples from the springs nos. 16-17 and 19 about 200 kilometers showed evident increase trend before the Wenchuan earthquake [20, 24, 25]. The K+ concentrations of the samples of the springs nos. 16-17 and 19 increased by 19.3% to 54.4% before the main shock, and rapidly decreased to the normal values after the main shock. The concentrations of for the samples from the springs nos. 16 and 19 increased by 32.0% and 59.6% respectively before the earthquake, and then dropped to the normal value after the main shock (Figure 5). The springs nos. 16-17 and 19 occur in the intersection of the LMSF, XSHF and ANHF zones where the extrusion stress increased yearly from 2004 to 2007 before the Wenchuan earthquake [40, 41]. Therefore, the increase of K+ and concentrations for the samples of the springs nos.16-17 and 19 could be attributed to the more supplement of the deep fluids with higher concentrations of K+ and [42], The 3He/4He ratios of the samples from the springs nos. 16-17 and 19 increased before the Wenchuan earthquake, but showed a remarkable drop after the shock, and then became normal indicating the more supplement of mantle-derived fluids before the main shock [30].

5.2.2. Spatial Variations

The geochemical and isotopic compositions of the samples from the springs nos. 1–5 occur in the LMSF zone performed evident variation after the Wenchuan earthquake, with the highest amplitude as 95% for of the sample from the spring no. 5. However, for the other samples from the other fault zones, just the samples from the springs nos. 15 and 20 presented obvious geochemical variations, with the highest amplitude as 53% for Na+ of the sample from the spring no. 15 (Figure 4). The Wenchuan earthquake happened in the LMSF zone, which enhanced the Coulomb stress of the zone [38, 39], and with most of the aftershocks occurred there (Figure 1). For the springs nos. 1–5 occurred in the LMSF zone, the epicenter distances are between 50 and 110 kilometers, which are smaller than those of the hot springs nos. 6-32 (between 190 and 410 km). Therefore, the variation of the tectonic stress induced by the Wenchuan earthquake and the aftershocks may resulted in, the geochemical and isotopic variations for the springs nos. 1–5 after the Wenchuan earthquake.

6. Conclusions

120 water samples of the 32 hot springs in the western Sichuan have been analyzed after the Wenchuan earthquake. The following conclusions can be drawn.(1)The waters of the 32 springs were mainly derived from meteoric precipitation. Because of the isotopic elevation effect, the δD and δ18O of the samples from the higher mountain area were lower, but those of samples from the lower altitude region were higher. all of the analyzed water samples from the 32 springs can be classified into nine chemical types: Na(Ca)-HCO3(SO4), Na(Mg)-HCO3(SO4), Ca(Na)-HCO3(SO4), Mg(Ca)-SO4, Ca(Mg)-SO4, Ca(Mg)-HCO3, Mg(Ca)-HCO3, Na-Cl(HCO3), and Ca(Na)-SO4(HCO3), which were mainly controlled by water-rock interaction and the input of deep fluid.(2)The concentrations of K+ and of the samples of the springs nos.16-17 and 19 in the intersection of LMSF, XSHF, and ANHF varied evidently, with the amplitude ranging from 19.3% to 59.6% before and after the Wenchuan earthquake which may be attributed to the interfusion of the deep fluids with high K+ and induced by the increased tectonic stress.(3)The hydrogeochemical variations of the springs closer to the epicenter performed more obviously after the Wenchuan earthquake. As the aftershock activities became weak, the geochemical parameters of the samples from the springs nos. 1–5, 15, and 20 located in the regions where the tectonic stress was enhanced before the Wenchuan 8.0 earthquake decreased by more than 20%, and with the δ18O the LMWL, which may be related to the change of tectonic stress when the aftershock activities became weak.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to the Seismological Bureau of Sichuan province and Professor Fang Du for helping in the field work. The authors sincerely thank the editor and reviewers for their valuable comments and grammar edits in improving the paper. This research was supported by the National Science and Technology Support Program (the characteristics of gas geochemistry of Longmenshan faults, WFSD-10), the Basic Science Research Plan of the Institute of Earthquake Science, China Earthquake Administration (The fluid geochemical survey of the mud volcanoes in the North of the Tianshan Mountain, Xinjiang Uygur Autonomous Region 02132417), and the Basic Science Research Plan of the Institute of Earthquake Science, China Earthquake Administration (02132417).