Identification and quantification of contributions to karst groundwater using a triple stable isotope labeling and mass balance model

Multiple stable isotopes and geochemical approaches to elucidate groundwater reactive transport paths in mining cities: A case from the northern Anhui, China

Mining cities are ecotone areas where human and natural components interact. Indeed, the negative effects of mining activities on drinking water quality have become a serious public concern worldwide. To elucidate groundwater genesis and reactive transport path controlling the water pollution, a multi-bodies system in the Sunan Mine area in China was considered in this study. The results of the mineral phase characterizations, hydrochemical analysis, and multiple stable isotopes (δ2H/δ18O, δ34S and 87Sr/86Sr) indicated that calcite, dolomite, gypsum, quartz, halite, organic carbon, and gases (O2, CO2 and H2O) were the primary reactants in the aquifer system, accompanied by dissolution and precipitation of minerals, cation exchange, desulfation, and evaporation. An inverse hydrogeochemical model was employed to identify three paths, Path 1 demonstrated that mine water mainly originated from the Quaternary loose aquifer water (QLA), Permian fractured sandstone aquifer water (PFA), and Carbonifer fractured limestone aquifer water (CFA), accompanied by high K++Na+ and HCO3- concentrations due to the carbonate dissolution, halite dissolution, and cation exchange processes. Path 2 showed that the recharge of the CFA and Ordovician fractured limestone aquifer (OFA) occurred from the shallow recharge zone to the deeper OFA water through faults and fractures, mainly involving halite dissolution, carbonate dissolution, and gypsum dissolution. Path 3 demonstrated the interaction between the Hui River, collapsed pond water, and QLA, accompanied by gypsum dissolution, calcite dissolution, and cation exchange. Although the shallow QLA quality met the WHO drinking water standards, the pollution risk from the surface collapse pit water cannot be ignored. Therefore, effective approaches need to be considered in the study area to reduce the connection between the collapse pit water and QLA. The study results can help decision makers to predict water quality of complex water systems in ecotone areas and other similar regions worldwide.

The main strategies for soil pollution apportionment: A review of the numerical methods

Nowadays, a large number of compounds with different physical and chemical properties have been determined in soil. Environmental behaviors and source identification of pollutants in soil are the foundation of soil pollution control. Identification and quantitative analysis of potential pollution sources are the prerequisites for its prevention and control. Many efforts have made to develop methods for identifying the sources of soil pollutants. These efforts have involved the measurement of source and receptor parameters and the analysis of their relationships via numerical statistics methods. We have comprehensively reviewed the progress made in the development of source apportionment methodologies to date and present our synthesis. The numerical methods, such as spatial geostatistics analysis, receptor models, and machine learning methods are addressed in depth. In most cases, however, the effectiveness of any single approach for source apportionment remains limited. Combining multiple methods to address soil quality problems can reduce uncertainty about the sources of soil pollution. This review also constructively highlights the key strategies of combining mathematical models with the assessment of chemical profiles to provide more accurate source attribution. This review intends to provide a comprehensive summary of source apportionment methodologies to help promote further development.

River Sediments Downstream of Villages in a Karstic Watershed Exhibited Increased Numbers and Higher Diversity of Nontuberculous Mycobacteria

The impact of residential villages on the nontuberculous mycobacteria (NTM) in streams flowing through them has not been studied in detail. Water and sediments of streams are highly susceptible to anthropogenic inputs such as surface water flows. This study investigated the impact of seven residential villages in a karst watershed on the prevalence and species spectrum of NTM in water and sediments. Higher NTM species diversity (i.e., 19 out of 28 detected) was recorded downstream of the villages and wastewater treatment plants (WWTPs) compared to sampling sites upstream (i.e., 5). Significantly, higher Zn and lower silicon concentrations were detected in sediments inside the village and downstream of the WWTP's effluents. Higher phosphorus concentration in sediment was downstream of WWTPs compared to other sampling sites. The effluent from the WWTPs had a substantial impact on water quality parameters with significant increases in total phosphorus, anions (Cl-and N-NH3-), and cations (Na+ and K+). The results provide insights into NTM numbers and species diversity distribution in a karst watershed and the impact of urban areas. Although in this report the focus is on the NTM, it is likely that other water and sediment microbes will be influenced as well.

A simple multilevel sampler for synchronous collection of stratified waters

Stratified water collection plays a crucial role in water quality monitoring, as most water bodies are not perfectly mixed systems. In order to precisely collect stratified waters, we developed an inexpensive, simple, and high-resolution sampler to simultaneously collect and measure physical and chemical parameters along vertical water profiles. The water sampler predominantly consists of two parts: (1) an apparatus for measuring sampling depth below the water and (2) water sampling units secured below the water. Proof of concept water sampling was performed in Caohai wetland (Southwest China) at 40 cm intervals, as sampling depth and interval are adjustable. Stratified waters in four sampling sites were characterized by markedly different levels of major and trace elements as well as physicochemical parameters. Results indicate this simple multilevel sampler to be a cheap, precise, and portable option for simultaneously collecting water samples at different depths in a wide array of water body types.

Behavior of Sb and As in the hydrogeochemistry of adjacent karst underground river systems and the responses of such systems to mining activities

Through the investigation of Qinglong mining area and adjacent karst underground river system, mining activities and water-rock interactions are found to control the hydrogeochemical evolution of karst underground water. Along the flow direction of the karst underground river, the hydro-chemical type is converted from HCO₃-Ca type to SO₄-Ca type. The concentrations of Sb and As also gradually decrease. Using PHREEQC to calculate the SI shows that: in the karst underground river system, both gypsum and fluorite are unsaturated, indicating a high degree of water-rock interaction. LogPCO₂ is negatively correlated with pH, indicating that the karst underground river systems are both open systems. The dissolution of carbonate minerals and the alternate adsorption of ions are the main water-rock interactions that lead to the rapid decline of Sb and As concentrations. This research also applies principal component analysis to identify the types of pollution in adjacent karst underground river systems. The results show that the LongBaiwei underground river was mainly affected by coal mining activities, and Fe was more prominent; the ShuiYa underground river was more significantly affected by the leachate from the antimony tailings yard. This study provides a scientific basis for the evolution of the water environment as well as strategies for pollution prevention and control in typical karst underground river systems owing to the influence of mining activities.

The effects of rainfall on groundwater hydrogeochemistry and chemical weathering

Due to the special hydrogeological conditions in karst areas, groundwater responds quickly to rainfall. The covariation of ion concentrations and spring discharge can help better understand the hydrogeochemical process of groundwater occurring in the heterogeneous karst aquifers. In this study, high-resolution monitoring of groundwater discharge, hydrochemistry, and stable isotopes was conducted at the Qingjiangyuan (QJY), a spring of the Qingjiang watershed in Hubei, China. The purpose is to investigate the changes in hydrogeochemical processes and chemical weathering under the influence of rainfall. The dynamics of spring discharge indicate the presence of pipelines and fissures of different sizes. According to the spring discharge attenuation curves, there are at least three medium types in the aquifer, which account for 45.7%, 34.2%, and 20.1% of the total groundwater. Pearson correlation analysis shows that the main sources of the solute in the QJY are carbonate minerals (mainly calcite and dolomite), evaporites (mainly gypsum and sylvite), celestite, and strontianite. Anthropogenic activities have less impact on groundwater solutes. Although carbonate minerals dominate the hydrochemistry, the changes in hydrogeochemical behavior caused by rainfall may come from gypsum, which is supported by the ion concentrations. At the early rainfall stage, Ca²⁺ concentration increased from 42.9 to 45.6 mg/L, followed by the SO₄^2- from 15.2 to 16.6 mg/L. When the discharge increased to the maximum (2320 L/s), Ca²⁺ and SO₄^2- showed opposite trends, decreased to 39.7 mg/L and 10.4 mg/L, respectively. The results also suggest that carbonate rocks and evaporites have important roles in hydrochemistry. The contributions of these three end-members were quantified based on the law of mass conservation. The proportions of carbonate weathering and evaporite weathering were 83.4% (85.2-80.3%) and 11.6% (6.9-18.0%), respectively, and rain was 5.0% (0.1-10.4%). These results were integrated into a hydrogeological conceptual model that explains the hydrogeochemical processes, including rock weathering, piston, and dilution effects caused by rainfall. The proposed conceptual model helps to improve the understanding of hydrogeochemical processes and chemical weathering in karst areas.

Identifying strontium sources of flowback fluid and groundwater pollution using 87Sr/86Sr and geochemical model in Sulige gasfield, China

Hydraulic fracturing technology has made unconventional oil and gas development economically viable; however, it can lead to potential environmental issues such as groundwater pollution. Strontium isotope (⁸⁷Sr/⁸⁶Sr) is considered as a sensitive tracer to indicate potential groundwater contamination. In this study, strontium (Sr) and ⁸⁷Sr/⁸⁶Sr sources of hydraulic fracturing flowback fluid are identified with 87 flowback fluid samples and 5 borehole core samples. High Sr concentrations and ⁸⁷Sr/⁸⁶Sr values were found in fracturing flowback fluid. The hydrogeochemistry evidence shows high Sr and ⁸⁷Sr/⁸⁶Sr in fracturing flowback fluid mainly comes from formation water with high ion concentrations, while Sr and ⁸⁷Sr/⁸⁶Sr of formation water develop in diagenesis and long term water-rock interaction (e.g., feldspar dissolution and clay mineral transformations) under the high temperature and pressure. A complete evaluation system was executed to assess the sensitivity of ⁸⁷Sr/⁸⁶Sr indicating potential pollution on groundwater. The mixing curves which ⁸⁷Sr/⁸⁶Sr combined with Sr and Cl were also established by mixing models to indicate groundwater pollution. The modeling results show mineral dissolution/precipitation and cation exchange have little impact on ⁸⁷Sr/⁸⁶Sr in the mixing process between fracturing flowback fluid and groundwater, which ⁸⁷Sr/⁸⁶Sr can identify contamination when only 0.89% of fracturing flowback fluid mixes with groundwater. Finally, the potential contamination pathways are discussed. It is highly unlikely fracturing flowback fluid contaminates groundwater and soil through upward migration, whereas leakage is a more prevalent pollution pathway.

Determining nitrate and sulfate pollution sources and transformations in a coastal aquifer impacted by seawater intrusion-A multi-isotopic approach combined with self-organizing maps and a Bayesian mixing model

Over the past few decades, the La Paz aquifer system in Baja California Sur, Mexico, has been under severe pressure due to overexploitation for urban water supply and agriculture; this has caused seawater intrusion and deterioration in groundwater quality. Previous studies on the La Paz aquifer have focused mainly on seawater intrusion, resulting in limited information on nitrate and sulfate pollution. Therefore, pollution sources have not yet been identified sufficiently. In this study, an approach combining hydrochemical tools, multi-isotopes (δ²H_H2O, δ¹⁸O_H2O, δ¹⁵N_NO3, δ¹⁸O_NO3, δ³⁴S_SO4, δ¹⁸O_SO4), and a Bayesian isotope mixing model was used to estimate the contribution of different nitrate and sulfate sources to groundwater. Results from the MixSIAR model revealed that seawater intrusion and soil-derived sulfates were the predominant sources of groundwater sulfate, with contributions of ~43.0% (UI₉₀ = 0.29) and ~42.0% (UI₉₀ = 0.38), respectively. Similarly, soil organic nitrogen (~81.5%, UI₉₀ = 0.41) and urban sewage (~12.1%, UI₉₀ = 0.25) were the primary contributors of nitrate pollution in groundwater. The dominant biogeochemical transformation for NO₃^- was nitrification. Denitrification and sulfate reduction were discarded due to the aerobic conditions in the study area. These results indicate that dual-isotope sulfate analysis combined with MixSIAR models is a powerful tool for estimating the contributions of sulfate sources (including seawater-derived sulfate) in the groundwater of coastal aquifer systems affected by seawater intrusion.

Evaluating the genesis and dominant processes of groundwater salinization by using hydrochemistry and multiple isotopes in a mining city

The increasing salinization of groundwater renders it challenging to maintain the water quality. Moreover, knowledge regarding the characteristics and mechanism of groundwater salinization in mining areas remains limited. This study represents the first attempt of combining the hydrochemical, isotope (δD, δ¹⁸O, δ³⁷Cl, and ⁸⁷Sr/⁸⁶Sr) and multivariate statistical analysis methods to explore the origin, control, and influence of fluoride enrichment in mining cities. The TDS content of groundwater ranged from 275.9 mg/L to 2452.0 mg/L, and 54% of the groundwater samples were classified as class IV water according to China's groundwater quality standards (GB/T 14848-2017), indicating a decline in the water quality of the study area. The results of the groundwater ion ratio and isotope discrimination analysis showed that dissolution and evaporation involving water-rock interactions and halite were the main driving processes for groundwater salinization in the study area. In addition to the hydrogeological and climatic conditions, mine drainage inputs exacerbated the increasing salinity of the groundwater in local areas. The mineral dissolution, cation exchange, and evaporation promoted the F^- enrichment, while excessive evaporation and salinity inhibited the F^- enrichment. Gangue accumulation and infiltration likely led to considerable F^- enrichment in individual groundwater regions. Extensive changes in the groundwater salinity indicated differences in the geochemical processes that controlled the groundwater salinization. Given the particularity of the study area, the enrichment of salinization and fluoride triggered by mining activities cannot be ignored.

Numerical Study of the Strength and Characteristics of Sandstone Samples with Combined Double Hole and Double Fissure Defects

To explore the failure mechanism of rock with holes and fissures, uniaxial compression tests of sandstone samples with combined double hole and double fissure defects were carried out using Particle Flow Code 2D (PFC2D) numerical simulation software. The failure behaviour and mechanical properties of the sandstone samples with combined double hole and double fissure defects at different angles were analysed, and the evolution results of the stress field and crack propagation were studied. The results show that with a decrease in fissure angle, the crack initiation stress, damage stress, elastic modulus and peak stress of the defective rock decrease, while the peak strain increases, and the brittleness of the rock is weakened. Rocks with combined double hole and double fissure defects at different angles lead to different failure modes, crack initiation positions and crack development directions. After uniaxial compression, both compressive stress and tensile stress concentration areas are produced in the defective rock, but the compressive stress concentration is of primary importance. The concentration area is mainly distributed around the holes and fissures and the defect connecting line, and the stress concentration area decreases with the decreasing fissure angle. This study can correctly predict the mechanical properties of rock with combined double hole and double fissure defects at different angles and provide a reference for actual rock engineering.