Effect of Large-Scale Mining Drainage on Groundwater Hydrogeochemical Evolution in Semi-Arid and Arid Regions

Groundwater evolution and circulation in coal mining areas will be significantly affected by natural processes and human activities. However, the impacts of large-scale mining drainage on groundwater hydrochemistry are unclear in the semi-arid and arid inland coalfields in northwest China. In particular, for shallow buried areas, the spontaneous combustion of coal seam outcrops forms burnt rock that is rich in fractures. Being a strong water-yield aquifer after receiving recharge, burnt rock has become a potential source of mine water inrush hazards. Therefore, drainage from this aquifer is necessary to eliminate potential hazards, which also leads to the deterioration of the local ecological diversity and sustainability. The hydrogeochemical method is an effective way to study the source and evolution of groundwater in mining areas and to determine the long-term cumulative effect of mining and drainage on the hydrogeochemical evolution of burnt rock areas. It is, however, still poorly understood. In this study, we employed the hydrochemical and isotope (2H, 3H, and 18O) to investigate the long-term effects of drainage on the hydrogeochemical evolution in the coalfield of southern Xinjiang, China. The results showed that the hydrochemical environment became more complex as the effect of recharge of different tritium-based groundwater ages in multiple-layered aquifer system and leads to the changes in the concentration of the chemical components. Before large-scale mining drainage, groundwater flowed from west to east and was finally discharged into the Kuqa River. The major water–rock interactions that occurred were the dissolution of halite, carbonate and gypsum dissolution, cation exchange, and dedolomitization. After large-scale and long-term mining drainage from the burnt rock aquifer, the groundwater flow field had changed, and the Kuqa river turned to recharge into the burn rock aquifer, the hydrochemical evolution also changed from water–rock interaction to mixture. This was caused by the significant decline in the groundwater level and changes in groundwater circulation in the mining area, which further led to the loss of valuable surface water resources in this arid area. Based on the characteristics and recharge conditions of burnt rock, we suggested that grouting can effectively cut off the hydraulic connection between the river and burnt rock and thus protect water resources.


INTRODUCTION
Spontaneous combustion is commonly observed in coal seam outcrops all around the world (Oren and Sensogut, 2010;Song et al., 2017;Kruszewski et al., 2018;Onifade and Genc, 2018;Zhang et al., 2020). During spontaneous combustion, heat spreads upwards, which causes the overlying rocks to heat up and transform into a suite of thermally metamorphosed rocks (named burnt rock) through a process known as combustion metamorphism . As a result, the color and structure of the burnt rock change significantly (Huang et al., 2008) and eventually turn reddish and rich in fractures. Previous studies have mostly focused on the petrology and mineralogy of burnt rocks to describe mineral and textural transformations during spontaneous combustion (Baboolal et al., 2018;Laita et al., 2019). However, the burnt rock becomes a strong water-yield aquifer after receiving recharge due to fractures created by the combustion metamorphism processes, affecting the safety of the underground coal mining process. In northwest China, burnt rocks are widely formed in the coal outcrop area as the shallow buried Jurassic strata (Ma et al., 2021). Yuan et al. (2019) have documented that the structure of the burnt rock is damaged with a large number of pores and fractures, which is one of the key factors of instability in the surrounding rock. Xue et al. (2021) employed the AHP and entropy-weight method to evaluate the water abundance of BRA and to avoid the recurrence of water inrush during coal mining. Because of the strong water yield, a water inrush from BRA occurred in the underground coal mine of Yongxin, causing a tunnel flooding incident that killed two people (Xue et al., 2021). In order to avoid such hazards, it is necessary to drain water from the overlying aquifer before mining (Dong et al., 2021). Secondary problems such as groundwater pollution (Wang et al., 2021b), ground subsidence (Zeng et al., 2020), and soil salinization (Zhang et al., 2007) may also occur. Furthermore, excessive drainage may change the flow path of the groundwater (Dimitrova-Petrova et al., 2020), which can complicate the hydrogeological condition of the mining area. The hydrogeochemical process and evolution of groundwater can reflect the flow process (Carrillo-Rivera et al., 2007;Xiao et al., 2022a), and the hydrogeochemical characteristics can change significantly in the mining area after excessive drainage. Therefore, it is necessary to study the hydrogeochemistry of BRA under the influence of large-scale mining drainage, which is poorly understood.
The hydrochemical and stable isotope analysis can be used to investigate the groundwater circulation, hydraulic connection, and water-rock interaction (Keita and Tang, 2017;Emvoutou et al., 2018;Meyzonnat et al., 2021). By combining the Piper diagram, Durov plot, and Gibbs plot, groundwater hydrogeochemical processes can be studied (Salameh et al., 2018). Stable isotopes and radioactive isotopes are used to identify the source and age of groundwater (Huang et al., 2017;Wang et al., 2021a). With the help of multivariate statistical analysis (MSA), such as hierarchical cluster analysis (HCA) (Barzegar et al., 2017), principal component analysis (PCA) (Dehghanzadeh et al., 2015;Xiao et al., 2022b), and correlation analysis (CA) (Qiao et al., 2019), a more reasonable explanation for the hydrogeochemical changes can be obtained (Qu et al., 2021).
In this study, the integration of hydrogeochemistry and isotope ( 2 H, 18 O, and 3 H) analysis, MSA, hydrogeological condition, and coal mine exploitation was employed to identify the origins of groundwater and hydrogeochemical evolution under the influence of mining drainage in a typical arid coalfield located in southern Xinjiang that is affected by BRA in China. This work aimed to 1) identify the key factors controlling hydrogeochemical origin and evolution, 2) reveal the effect of large-scale mining drainage on the hydrogeochemistry for BRA in the coalfield, and 3) propose reasonable suggestions for water resource protection and mine safety production.

STUDY AREA
The A-ai coalfield is located in southern Xinjiang, China (Figure 1), which has a total area of 32 km 2 . It is dominated by a continental arid climate, characterized by strong evaporation, low precipitation, and frequent droughts (Sabit et al., 2008), with an average annual precipitation of 90.09 mm and mean annual pan evaporation of 2,219 mm. Approximately, 85% of the annual precipitation occurs between June and September . The landform is influenced by wind erosion with surface elevations varying from +1840 m to +2005 m, and the topography of the A-ai coalfield is generally a basin with high surroundings and a low middle. In addition, four gullies are perpendicular to this basin from the north to the south. The Kuqa River is the only perennial river in the region that acts as the main source of drinking water. A syncline runs through this area and stretches 12 km wide from west to east and has a length of 2-4 km from south to north. No other large faults are found here. The coal seam spreads along the two wings of the syncline and is exposed on the surface. Due to the subsurface spontaneous combustion of the coal seam, sandstone and mudstone on the roof have changed the appearance and petrologic characters into burnt rock with a great number of fractures ( Figure 2). Therefore, the A-ai coalfield is bounded by gully Ⅰ in the west, the Kuqa river in the east, and the coal seam and burnt rock outcrop in the north and south (Figure 1).
The sedimentary sequence above coal-bearing strata in the study area is shown in Figure 3. From top to bottom, the strata are Quaternary and Jurassic. According to the microscopic features of typical samples from Jurassic, quartz, feldspar, mica, carbonate, and gypsum account for more than 90% (Liu et al., 2011). The quaternary aquifer (QA) is only deposited near gullies and the Kuqa river and has a thickness of 0-40 m. The Jurassic Ahe Group (J 1a ) is dominated by fine sandstone and mudstone but is poor in water. The Jurassic Taliqike Group (J 1t ) is the coal-bearing stratum, and above the coal seam is a sandstone confined aquifer (SCA). In the coal outcrop areas, the overlying rock of the coal seam gradually turns to melted or honeycombed forms, generating plenty of cracks and fractures during the combustion of coal. Burnt rocks can easily receive recharges from precipitation and become a strong water-yield aquifer. The drilling unit water inflow (q) is a generally accepted index in China to evaluate the water yield of an aquifer (Qiao et al., 2017). According to the pumping test, the average hydraulic conductivity (K) and q value (flow rate when the water level drops 10 m for a 91 mm borehole) of the SCA are 0.049 m/d and 0.017 L/(sm), respectively, indicating a weak permeability and water yield. Therefore, this aquifer has little impact on coal mining. In contrast, the burnt rock aquifer (BRA) is distributed at the outcrops of coal seams on both sides of the syncline with great permeability and a very strong wateryield (with an average K of 17.874 m/d and q of 11.434 L/(sm), indicating that BRA is the main water-filled aquifer and the primary threat to coal mining safety. Under natural conditions or before large-scale coal mining, the groundwater recharge sources are precipitation and surface water (SW) from gullies in the rainy season. The direction of groundwater flow is from the west to the east, which then gets discharged to the Kuqa River in the east in the form of spring and direct contact.
Four coal mines are distributed in this area, namely, Yushuquan (YSQ), Dapingtan (DPT), Yongxin (YX), and Qiuci (QC). From 2012 to 2014, due to the production of new mines and upgrades in the mining technology (from the pillar mining method to the comprehensive mechanized mining method), the amount of coal production increased sharply from 0.12 to 1.59 million tons per year and has since been maintained to date. Several water inrush disasters have occurred in this area, causing the deaths of 10 people and the flooding of several tunnels in the A-ai coalfield. Therefore, with an increase in mining output, the large-scale drainage of BRA has been carried out by drilling holes in an effort to avoid disasters such as that which occurred in 2015.

Groundwater Sampling and Analytical Techniques
In this study, 46 water samples from groundwater and surface water were obtained from the A-ai coalfield from the years 2007-2021. The data included 22 samples from the BRA (including 1 sample from the spring), 17 samples from the SCA, 3 samples from QA, and 4 samples from SW (2 samples from the Kuqa river and 2 samples from the gully Ⅲ). Sampling locations were recorded using a portable GPS device and displayed in Figure 1. All of the groundwater samples were taken from springs or boreholes. Water samples were mainly analyzed at the Shaanxi Key Laboratory of Prevention and Control Technology for Coal Mine Water Hazard. The concentrations of Na + +K + , Mg 2+ , Ca 2+ , SO 4 2− , and Cl − ions were measured by ion chromatography, while the concentrations of HCO 3 − and CO 3 2were measured by the traditional titrimetric method. The total dissolved solids (TDS) were calculated based on all the ion concentrations determined previously (Table 1). pH was measured in situ with a portable pH meter, and tritium ( 3 H) was determined on electrolytically enriched water samples by liquid scintillation spectrometry. For the reliability of the chemical data , the error of the charge balance error (CBE) was calculated as follows: where all ions are measured in meq/L. The calculated results of the CBE showed that the uncertainty was less than ±5%, which is acceptable for this study. The stable isotopes of δ 18 O and δ 2 H were measured by using a liquid water isotope analyzer (LGR LWIA-24d) and compared with the Vienna standard mean ocean water standard, which had precisions of 0.21‰ and 0.04‰, respectively.
Pumping tests were conducted to obtain the aquifer hydraulic conductivity. A total of 10 single-hole steady-flow pumping tests were carried out in the study area in the past 10 years.

Multivariate Statistical Analysis
Multivariate statistical analysis has been widely used to extract valuable information from large hydrogeochemical datasets (Khanoranga and Khalid, 2019;Liu et al., 2020). In this study, hierarchical cluster analysis (HCA), principal component analysis (PCA), and correlation analysis (CA) were employed to analyze the hydrogeochemical classification and the temporal and spatial variation. The Ward linkage method and the Euclidean distance were performed in the HCA of this study, and a dendrogram was used to graphically show clusters of the datasets and their proximity (Cloutier et al., 2008). PCA is used to explain the majority of variance by extracting the eigenvalues from the correlation between the chemical variables (Singh et al., 2017), while CA can reveal a statistical relationship between two or more variables (Nagaraju et al., 2016). All statistical calculations were performed in SPSS software version 20.0 for Windows.

General Hydrogeochemical Properties
The variation of hydrogeochemical composition in coalfields is affected by many factors, such as hydrogeological condition, lithology, mineral composition, precipitation, temperature, and also by mining and other human activities . In Table 1, all the water samples were classified by the sampling location and arranged by the sampling time. The pH of all samples ranged from 7.11 to 8.88, with a mean of 7.82, which is slightly alkaline. It was closely related to alkaline mineral sources in the aquifer system, especially samples 40 and 41 from the Kuqa River (8 and 8.19, respectively) indicating the close relationship with soil salinization (Xuemei et al., 2009). The TDS reflects the strength of water-rock interaction and the hydraulic connection between multi-aquifers (Shvartsev et al., 2016;Hao et al., 2020). The TDS values of BRA ranged from 505.53 mg/L to 2,124.90 mg/L with a mean of 1570.23 mg/L, while those of SCA ranged from 1394.93 mg/L to 7058.26 mg/L with a mean of 3331.89 mg/ L. The ion concentrations in both BRA and SCA were all obviously larger than those in SW and QA. Since the  spontaneous combustion of coal seams leads to a significant increase in the fractures in BRA, there will be a good hydraulic connection and thus lead to a small standard deviation in the ion concentrations, the TDS, and pH in the aquifer. Figure 4 shows that several main hydrogeochemical types could be identified in all samples. The water samples from BRA and SCA showed a similar hydrogeochemical type but were different in mineralization. SCA samples were presented as a highly saline Na-Cl-SO 4 type. Although the major hydrogeochemical type in SCA was dominated by Na-Cl-SO 4 , the average concentration of SO 4 was just slightly excess 25%, while the concentrations of both Na and Cl exceeded 55%. Meanwhile, the standard deviation of TDS in SCA was 1667.68 mg/L, which gives an indication of the heterogeneity. A total of 22 water samples from BRA are Na-Cl-SO 4 water chemical type with a lower TDS, including a sample from the spring near the Kuqa River (the spring is origin from the BRA at a flow rate of 0.69 L/s in 2007 but stopped flowing in 2015, when large-scale mining drainage began). In contrast, mixed water chemical types were observed in samples from the surface water and QA with low mineralization (ranging from 259.16 to 968.14 with a mean of 634.43). As a local drinking water source, two samples from the Kuqa River over 10 years met the WHO guidelines for both drinking and irrigation. The samples of QA were collected from different gullies with the weak hydraulic connections.

Relationships Among Major Ions by Multivariate Statistical Analysis
Six major ions were classified by the R-mode HCA process, while the Ward method and squared Euclidian distances were employed to measure the similarity between them. As shown in Figure 5, these ions could be classified into two clusters. Ca, Mg, SO 4, and HCO 3 were classified into cluster Ⅰ, suggesting the dissolution of carbonate. Cluster comprised Na and Cl, indicating the dissolution of halite.
The Pearson correlation coefficients are presented in Table 2, which shows that the TDS correlates strongly with Cl and Na (greater than 0.9), indicating that the dissolution of halite plays an important role in the increment of TDS. In addition, Ca correlated strongly with Mg (0.846), indicating the dissolution of dolomite.
We further used the principal component analysis to investigate the hydrochemical processes in this area. The principal components were extracted from the correlation matrix from the seven parameters (Ca, Mg, Na, HCO 3 , SO 4 , Cl, and TDS). The Kaiser criterion is such that only the components with eigenvalues greater than one are retained (Yilmaz and Buyukyildiz, 2015). Table 3 shows the principal component loadings and their separate variances explained. The result shows that the first two components extracted have eigenvalues greater than one and account for 80.61% of the total variance in the dataset. Component 1 is characterized by highly positive loadings in Na, Cl, and TDS (all higher than 0.9) with 52.28% of the explained variance. Component 2 is characterized by highly positive loadings in Ca, Mg, and SO 4 . The loadings of the chemical parameters can be defined by descriptive terms according to the position in the plane  ( Figure 6). Component 1 was defined as the salinity component because of the close relationship of Na, Cl, and TDS, while component 2 was defined as the hardness component because of the association of Ca and Mg, which was used to calculate hardness. The factors that cause these high loads will be discussed in a later section. The results of the multivariate statistical methods (HCA, PCA, and CA), Na, Cl, and TDS are dominant in cluster 1 and those of Ca, Mg, and SO 4 are dominant in cluster 2. This means that the natural processes dominate the evolution of groundwater hydrogeochemistry in this study.

Relationship Between the Hydraulic Conductivity and TDS
The hydraulic conductivity K can reflect the residence time by affecting the seepage velocity (Loaiciga, 2004), while the length of residence time can reflect the degree of water-rock interaction, which causes changes in the characteristics of the hydrogeochemistry, especially TDS (Jia et al., 2019). Ten pumping tests were conducted to assess the hydraulic conductivity of BRA and SCA in the study area ( Figure 7). The values of K and TDS in BRA ranged from 3.11 to 27.86 m/d and 505.53-2,121.90 mg/L, with averages of 17.874 m/d and 1570.23 mg/L, respectively. In contrast, the values of K and TDS in SCA ranged from 0.0001 to 0.14 m/d and 1394.93-7058.26 mg/L, with averages of 0.0049 m/d and 3331.89 mg/L, respectively. Similarly, the values of q in SCA ranged from 0.00001 to 0.047 L/(m s) with an average of 0.017 L/(m s), while the values of q in BRA ranged from 2.13 to 27.226 L/(m s) with an average of 11.434 L/(m s). Obviously, the K and q values of BRA are significantly larger than those of SCA, while the TDS is lower in BRA than in SCA. The pumping test in the three boreholes (including samples 2 and 26 from the same borehole but in different aquifers, and samples 3 and 28, 6 and 30), the difference are shown in K and TDS values of the two aquifers. The four water samples (samples 2 and 3 from BRA and samples 26 and 28 from SCA) showed the same water chemical type but with difference in mineralization, while two samples (samples 6 and 30) from these two aquifers showed different chemical types. In addition, sample 25 showed the highest TDS but the lowest K. Under the same hydraulic gradient, the smaller value of K provided more sufficient residence time and a wider range of water-rock    interaction. As a result, the high velocity of the groundwater flow results in low concentrations of ions in BRA.

Geogenic Sources of Chemical Compositions in Groundwater
The source of groundwater and the process of hydrogeochemistry can be identified by the ratio coefficient of the major ions (Li et al., 2016). Combined with MSA, more valuable information about the relationship between hydrogeochemical components can be found. The ideal ratios of Cl and Na would be one, if such ions are released by halite dissolution (Ettazarini, 2005). Meanwhile, because of the close relationship between Na and Cl calculated by MSA, scatter plots of Na and Cl are drawn in Figure 8A. Most samples are plotted along the 1:1 line, indicating that the dissolution of halite could be the main hydrogeochemical process in this area. This result is consistent with the principal component analysis, and both refer to the salinity component.
Similar to the dissolution line of halite, carbonate and gypsum dissolution have the decided ratio relationship (Eqs. 1-3).
As the main parameter for calculating hardness, the relationship between Ca and Mg also indicates the hardness component in PCA.
Considering the significant correlations between Ca, Mg, and SO 4 , most water samples were below the 1:1 line of gypsum dissolution ( Figure 8B) and above the 2:1 line of carbonate dissolution ( Figure 8C), which point toward cation exchange. Therefore, the Schoeller indices (CAI-I and CAI-II) and the relationship between (Na-Cl) and (Ca + Mg-HCO 3 -SO 4 ) ( Figures 8D,E) were employed to verify that the cation exchange is one of the hydrogeochemical processes. CAI-I and CAI-II were calculated according to the following Eqs 4, 5, where ions are expressed in meq/L.   If cation exchange has a significant impact on hydrogeochemical evolution, the water samples should fall in the third quadrant in Figure 8D with negative CAI values and linear with a slope of −1 in Figure 8E (Fisher and Iii, 1997). Thus, most samples exhibit cation exchange plays an important role in this area. Furthermore, the slightly negative correlation between HCO 3 and Ca (-0.004 in Table 2) indicates the possible dedolomitization process. With the increase in the Ca concentration caused by the dissolution of gypsum, the dissolution of dolomite will also lead to a decrease in HCO 3 with the precipitation of calcite. The process could be simplified as follows (Zang et al., 2015).
Therefore, halite, carbonate and gypsum dissolution, cation exchange, and dedolomitization significantly affect the variation of the major cations in the groundwater of the A-ai coalfield.
The saturation index (SI) of minerals can be calculated from the equation: SI = log (K IAP /K SP ), where K IAP is the ion activity   product for the mineral equilibrium reaction, and K SP is the solubility product. The results of the SI suggest a trend of mineral equilibrium and water-rock interaction, which was calculated by the PHREEQC program and is shown in Table 1. The relationship between SI and TDS for each aquifer is shown in Figure 9. The halite and gypsum concentrations for all groundwater samples have negative SI values, while a positive trend with TDS can be seen in Figures 8C,D. This indicates that the mineral will be continuously weathered by the groundwater and cause increments of TDS. However, it should be noted that water samples from low permeability in SCA have a higher saturation index than those from the high permeability BRA and surface water. In contrast, the SI values of dolomite and calcite from all samples were greater than zero, indicating that these samples were oversaturated with carbonate minerals. The precipitation of calcite and dolomite in groundwater may be related to low rainfall and strong evaporation (Kumar and Singh, 2015).

Origin of Groundwater
Using stable isotopes of δD and δ 18 O is an effective method to analyze groundwater circulation and origin (Craig, 1961).
In this study, the local meteoric water line (LMWL) was used as a baseline, where δD = 7.70δ 18 O+5.74 (Wang and Zhang, 2017). The global meteoric water line (GMWL, δD = 8δ 18 O+10, Craig, 1961) is also plotted in Figure 10. The difference in the slope and intercept of LMWL can reflect the vapor sources and climatic conditions in different regions. Through comparison of the two lines, the slope and intercept of the LMWL (7.70, 5.74) are significantly smaller than those of the GMWL (8, 10), indicating the dry climate and strong evaporation in this area (Wu et al., 2010).
All groundwater samples are distributed close to the GMWL and LMWL, indicating that meteoric origin is the major recharge source. In addition, the deviation of groundwater isotope values from the meteoric waterline may be caused by the secondary evaporation of precipitation before infiltration (Peng et al., 2004). Meanwhile, most SCA water samples deviate significantly from LWML compared to BRA water samples, indicating that the longer water-rock interaction in SCA than that of BRA. Due to the low flow rate and the long flow path, the surface water from Gully Ⅲ experienced more intense evaporation and deviated from LMWL more than that from the Kuqa River.
According to the identification method from USGS (Lindsey et al., 2019), the tritium-based age of groundwater could be classified into the modern (post-1950s) and premodern (pre-1950s) thresholds. Groundwater with a tritium content of less than 0.5 TU is classified as premodern groundwater, while larger than 0.5 TU is classified as a modern one. As shown in Table 1, the 3 H concentrations of water samples from SCA range from 0.28 to 0.39 TU with a mean value of 0.34 TU, indicating that the age of groundwater from SCA is earlier than in the 1950s. In contrast, the 3 H concentrations from BRA (ranging from 3.03 to 8.96 TU with a mean value of 5.84) and surface water (ranging from 7.11 to 10.37 TU with a mean value of 8.74) are significantly higher than the threshold. Combined with the results of stable isotopes, one can conclude that the water in BRA is recharged since the 1950s. Furthermore, although the origins of groundwater are all from precipitation through rock outcrops, the circulation condition of BRA is better than SCA, which is consistent with the results of stable isotopic of δD and δ 18 O.

Effect of Large-Scale Mining Drainage on Hydrogeochemistry
According to the previous discussion, BRA is superior to SCA in terms of the recharge conditions, K, and flow velocity. Therefore, BRA is the main aquifer in the A-ai coalfield and has been drained by boreholes during the mining process. With increases in coal production, the average annual mine water inflow also increases ( Figure 11), which mainly comes from BRA (it should be noted that because the DPT coal mine has not been produced, only three mines are included). In particular, in the YSQ coal mine, large-scale mining drainage was carried out in 2015, and the spring was cut off in this year. This year is considered to be the start of large-scale mining drainage. Due to the lack of long-term and continuous water level observation data, the water level of BRA boreholes along the syncline in 2011 and 2021 is shown in Figure 12. Under natural conditions or before large-scale mining drainage in 2011, groundwater in BRA flowed from the west to the east and discharged into the Kuqa River. In contrast, after largescale mining drainage for many years, the water level of BRA declined sharply by almost 50 m. The water level of the Kuqa River is significantly higher than that of BRA and may cause the discharge of the river to the BRA. In order to further study the influence of water level decline on hydrogeochemical characteristics, the water samples of BRA were divided into two clusters, according to the sampling time (before or after the start year of largescale mining drainage in 2015), as shown on the planar graph (Figure 13), where the diameter of the circle represents the value of TDS. First, sampling locations near the outcrop area are easier to receive recharge from meteoric water thus with lower TDS. Second, the water chemical types of these two clusters are different and the concentration of HCO 3 slightly increased in the second cluster, indicating that additional origin or water-rock interaction occurs in these samples. In the YX coal mine area, the TDS of the second cluster samples is significantly lower than the first cluster, and TDS gradually increases from the east to the west. More attention should be paid to samples 21 and 22 from BRA with approximately 300 m away from the Kuqa River. These two samples were collected from underground boreholes in YX coal mine, which have similar hydrogeochemical characteristics and quite a low TDS (630.71 and 505.53 mg/L, respectively). Obviously, the decline of the water level resulted in changes in the groundwater circulation condition. Low mineralization surface water flowed into the BRA through outcrop from the Kuqa River, and thus, this river has transferred from the recharge boundary to the discharge boundary. Therefore, the major hydrogeochemical process of BRA near the Kuqa river is mixing. As drainage continues, this variation will spread more widely.

Comprehensive Analysis and Water Resource Protection
Combined with the aforementioned discussions, the origin and evolution of hydrogeochemistry and hydrogeology in the BRA of the A-ai coalfield is summarized in Figure 14. Before large-scale mining drainage, seasonal flow in gullies of surface water will directly recharge to BRA; then, groundwater flows from the west to the east and recharges to the Kuqa River. During this process, mineral dissolution, such as halite, carbonate and gypsum, cation exchange, and dedolomitization play a pivotal role in the hydrogeochemical evolution. With large-scale mining drainage, the spring is disappearing, and the groundwater in the BRA start receiving recharge from the Kuqa River. At this point, the effect of the mixing action is markedly stronger than in other hydrogeochemical processes. Meanwhile, much groundwater flowed into the coal mine water sump and was discharged to the gullies or the Kuqa River after treatment and purification. Finally, a greater amount of low mineralization surface water will enter into the aquifer. As a result, the hydrogeochemistry of the aquifer will be changed in a wider range. On the other hand, clean and valuable surface water resources would be contaminated continuously in this area.
Drainage is inevitable during the mining process, so obstructing the source of recharge to protect water resources is a better approach. Considering the dry climate and less effect of precipitation, cutting off the flow path from the Kuqa River to the BRA is a preferred method. According to the hydrogeological characteristic of the burnt rock, the high permeability provides a good foundation for grouting. On the east side of the YX and QC coal mine, which is along the Kuqa River, cement could be injected into the BRA through boreholes to construct an underground curtain to effectively separate the Kuqa River and BRA. After that, only limited storage in the BRA and little recharge from precipitation is necessary for drainage. This not only decreases the cost of drainage but also protects the water resource. More importantly, coal mines in the A-ai coalfield will not be threatened by the water hazards of the BRA and would avoid casualties related to tunnel flooding.

CONCLUSION
In this study, we analyzed the effect of large-scale mining drainage on the regional characteristics of hydrogeochemistry in an arid coalfield by analyzing hydrogeochemical (major ions, 3 H, and stable isotopes of δD and δ 18 O) and hydrogeological data (pumping tests, annual average mine water inflow, and water level). The results showed that precipitation is the main recharge source of groundwater and surface water. In addition, the SCA is dominated by the Na-Cl-SO 4 type with higher TDS but has poor permeability and water yield. The significantly low concentrations of 3 H and high mineralization in SCA indicate the old water age and the weak hydraulic connection with other aquifers. In contrast, as the main risk factor for coal mine water hazards, the BRA has a young age and similar water chemical type with SCA but with a lower TDS. Spontaneous combustion of the coal seam caused a remarkable increase in permeability and the good recharge condition caused the strong water yield, which resulted in casualties and tunnel flooding events. Under natural or small-scale mining drainage conditions, groundwater flowed from the west to the east under the water-rock interaction of dissolution of halite, carbonate and gypsum, cation exchange, and dedolomitization, and they were finally discharged to the Kuqa River. After the perennial large-scale mining drainage, the water level of the BRA decreased sharply, causing changes to the hydrogeological condition. Furthermore, the main hydrogeochemical process in the BRA gradually changed from being a water-rock interaction into a mixing action. The influence range will gradually expand with continuous drainage. On the basis of the high permeability of the BRA, we proposed practical advice for the sustainable development and protection of surface water resources by grouting. By doing so, we could avoid the loss of surface water caused by prolonged drainage.
Similar to other coalfield areas, long-term mining activities have largely destroyed the original groundwater circulation pattern and produced a large amount of acid mine water. The results of this study provided a comprehensive analysis of groundwater evolution for the A-ai coalfield and also could serve as a template for other coalfields in northwest China to enhance the awareness of water resource protection (Huang et al., 2017).

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
AK, SZ, WG, and WH contributed to the conception and design of the study. AK, JZ, and XJ organized the database. AK and DS performed the statistical analysis. AK and SZ wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.