Study on the Impact of Closed Coal Mines on Groundwater in the Panlong River Basin (Shangdong Province, China) Based on Sulfur and Oxygen Isotopes

: To deeply investigate the impact of closed coal mines on groundwater sulfate contamination in the Panlong River basin, Zaozhuang, Shandong Province, China, the hydrochemical characteristics and controlling factors of groundwater were analyzed based on 64 surface water and groundwater samples. The sources of sulfate contamination were identified by sulfur and oxygen isotope analysis. The results indicate that the sulfate content in the pore groundwater exceeds the Standard for Groundwater Quality (China) Category III (250 mg/L), with a maximum content of 666.2 mg/L. Specifically, the exceedance rate of sulfate in pore water near the western boundary of the Taozao coal field is directly related to mine water. The exceedance rate of sulfate in karst groundwater is 28%, peaking at 1131 mg/L, and is mainly distributed outside the western boundary of the Taozao coalfield and near the Dingzhuang water source in the southeast of the coalfield, indicating a significant influence from high-sulfate mine water in the coal-bearing strata. The sulfur and oxygen isotope differences are significant between surface water, mine water, karst water, and pore water. Through Bayesian end-member analysis, it was determined that 40–83% of sulfate in the downstream pore water of the Taozao coal field originated from mine water, while 48–86% of the sulfate in the karst water originated from mine water, which greatly affects the groundwater in the western and southeastern parts of the Taozao coalfield. Therefore, measures should be taken as soon as possible to control the risk of sulfate pollution of neighboring groundwater, especially karst groundwater, by mine water from closed coal mines.


Introduction
Hydrochemical composition is a key element in evaluating the quality of surface water and groundwater of a region, as well as its development and utilization.The hydrochemical composition of surface water and groundwater in a region is influenced by multiple factors such as rock weathering, runoff, source conditions, human activities, and climate [1][2][3].
Meanwhile, based on the chemical parameters of various water bodies and statistical graphical interpretation, the sources and transformation patterns of hydrochemical ions in surface water and groundwater can be obtained [4].Using the hydrochemical statistical graphical interpretation method, Zhou et al. investigated the main ion characteristics of surface water in the upper reaches of the Shule River basin and the possible controls [5].Li et al. used a similar approach and analyzed the hydrochemical characteristics of karst groundwater during different hydrological periods.They also explored the natural factors and human factors affecting the major chemical ions [6].In addition to hydrochemical composition analysis, due to the exchange of substances and energy along with the transformation processes of different water bodies, comprehensive studies of hydrochemical and isotopic compositions of water bodies can also reveal the sources and cyclic process of surface water and groundwater [7,8].Gu et al. investigated the relationship between surface water and groundwater in the Liujiang Basin using hydrochemistry and environmental isotope analysis [9].Liu et al. applied hydrochemistry and hydrogen and oxygen stable isotope analysis to reveal that the ion content in different water bodies of the Donggong River Basin is affected by rainfall, temperature, and geological background [10].
Stable isotopes can be employed to identify the sources of water pollution.Sulfur isotopes in nature exhibit a unique isotopic composition that can be attributed to their source [11].Furthermore, apart from sulfate-reducing bacteria, no significant sulfur isotope fractionation occurs during the biogeochemical sulfur cycle [12].These characteristics render sulfur isotopes an invaluable tool for tracing the origins of sulfate in water bodies [13,14].Ma et al. [15] studied the source of sulfate in groundwater in Zaozhuang City (Shandong province, China) using the relationship between δ 34 SSO 4 -SO 4 , δ 34 SSO 4 -SI (gypsum), and δ 34 SSO 4 -SO 4 2− /Cl − .Based on the distribution characteristics of δ 34 SSO 4 -SO 4 , δ 34 SSO 4 -SI (gypsum), and δ 34 SSO 4 -SO 4 2− /Cl − , it is believed that sulfate pollution mainly originates from pyrite oxidation, gypsum dissolution and industrial and agricultural wastewater infiltration pollution, and in addition to the primary geological genesis, the infiltration of wastewater discharged by industrial and mining enterprises into the subsurface is also the main cause of sulfate pollution in groundwater.Zhang et al. [16] selected three groundwater potential sulfate source end elements, (A) sulfide oxidation SO 4 2− , (B) evaporite dissolution SO 4 2− , and (C) atmospheric precipitation SO 4 2− , and determined the degree of karst groundwater pollution by mine water using the relationship between δ 34 S (SO 4 ) and 1/SO 4 2− of water samples in the study area combined with the calculation of mixing ratio.Zhao et al. [17]  2− contamination may originate from the mine water of the coal system.The Taozao coalfield, situated in the central-eastern region of Zaozhuang City, Shandong Province, China, is a syncline extending in an east-west direction.It is part of the Carboniferous-Permian coalfield, which is bounded by faults of Beishan on its north side and seam outcrops on its east, west, and south sides.The coalfield extends for a length of 30 km in an east-west direction, with a width of 6 to 7 km in the north-south direction in the western part and 3 to 5 km in the eastern part.It covers an area of 90.10 km 2 .The coalfield has a long history of exploitation, with abundant coal resources.Following the closure of coal mines within the coalfield, there has been a significant rebound in formation water, which has changed the hydrodynamic conditions of the Taozao Basin.This may have an impact on the water quality of the surrounding karst water.The Panlong River is an important surface water body in Zaozhuang City.Therefore, this paper aims to clarify the impact of closed coal mine drainage on the Panlong River and regional groundwater.Firstly, surface water and groundwater samples were collected and tested.Next, the chemical characteristics of groundwater ions were analyzed, and analysis of the sulfur and oxygen Water 2024, 16, 1634 3 of 13 isotopes was applied to identify the pollution sources of groundwater sulfates.The findings of the study could serve as a scientific foundation for further assessment and prevention of groundwater sulfate pollution.

Study Area
The study area is situated within a warm and semi-humid continental monsoon climate, characterized by a moderate climate and distinct seasons.The precipitation is concentrated primarily from June to September, accounting for 80.5% of the total annual precipitation.The average precipitation over the period from 1980 to 2021 was 796.89 mm.The sunshine time is around 2380 h, the frost-free period is more than 200 d, and the average annual temperature is 13.9 • C. The study area is part of the Nansi Lake system in the Huai River Basin along the Grand Canal, with a well-developed surface water system.Major rivers include the Panlong River and its major tributaries, which flow from the north to the south, as well as the Dasha River in the Yicheng district.The Panlong River traverses the central and western regions, with its upper reaches divided into two branches: the south branch, which originates from Xingyu-Hengshankou in Xuecheng District, and the north branch, which originates from Baishan in Shizhong District.The drainage area of the north branch is 296 km 2 .As a rain-fed river, it receives surface runoff from the basin and a certain amount of industrial and mining enterprises' production and domestic sewage, which results in severe water pollution during the dry season.
The Taozao Basin exhibits pronounced topographical undulations.It is an irregular elliptical basin extending from east to west, displaying higher elevations in the north, south, and east and lower elevations in the west, with a relatively flat central section.The area comprises Quaternary unconsolidated porous aquifer groups, clastic rocks porous aquifer groups in the overlying strata of coal measure strata, the Carboniferous-Permian rocks fractured aquifer groups, old pits and water-filling roadway aquifer groups in coal seams (aquifer groups of coal measure strata) and karst aquifer groups of Cambrian-Ordovician carbonate rocks (Figure 1b).Among these aquifer groups, the water abundance of water in Quaternary unconsolidated porous aquifer groups is relatively low.The fractures in the overlying strata of the coal measure strata are poorly developed, resulting in a low water inflow per unit.The medium-bedded sandstone within the Carboniferous-Permian rocks fractured aquifer groups has a weak water abundance, and the shale has an extremely poor water abundance.While the shale layers between the coal seam floor and the Ordovician limestone act as relatively impermeable formations.The fractured rock karst aquifers of the Ordovician Majiagou group have undergone strong karst development, forming cool springs and sources of drinking water on the surface.In the Cambrian carbonate rock fractured karst aquifer groups, the karst features are well developed and the aquifers have abundant water, forming Shili Spring and Dingzhuang-Dongwangzhuang headwaters.Following the closure of mines, mine water infiltrates into the deep confined water through skylights, abandoned fractured wells, and faults, resulting in the pollution of the karst underground water.

Sample Collections
After a comprehensive examination of the characteristics of groundwater flow and the impact range of abandoned mines, a total of 64 sampling points were established in this survey (Figure 1a), comprising 7 surface water sampling points, 23 points of Quaternary pore water (referred to as pore water), 9 points of water in coal measure strata fissures (referred to as mine water), and 25 points of karst groundwater (referred to as karst water).Prior to sample collection, the wells were washed using portable water quality test kits, which were used to perform water quality tests every five minutes until at least three consecutive measurements of at least three testing indexes reached stable standards for groundwater sample collection.If these standards were not achieved after four hours of washing, Bailer's method was applied for sampling.Once the water quality stabilized, sample collection began and was completed within two hours.The samples were stored at 4 °C and sent to the laboratory within 24 h for analysis.

Sample Collections
After a comprehensive examination of the characteristics of groundwater flow and the impact range of abandoned mines, a total of 64 sampling points were established in this survey (Figure 1a), comprising 7 surface water sampling points, 23 points of Quaternary pore water (referred to as pore water), 9 points of water in coal measure strata fissures (referred to as mine water), and 25 points of karst groundwater (referred to as karst water).Prior to sample collection, the wells were washed using portable water quality test kits, which were used to perform water quality tests every five minutes until at least three consecutive measurements of at least three testing indexes reached stable standards for groundwater sample collection.If these standards were not achieved after four hours of washing, Bailer's method was applied for sampling.Once the water quality stabilized, sample collection began and was completed within two hours.The samples were stored at 4 • C and sent to the laboratory within 24 h for analysis.

Physical-Chemical Testing
A handheld GPS locator was employed to ascertain the coordinates of each sampling point.On-site pH and temperature measurements were conducted using a portable water Water 2024, 16, 1634 5 of 13 quality testing kit (Hach DS5, Loveland, CO, USA) and total dissolved solids (TDS) were quantified using a portable conductivity meter (SX-650).Anions, namely chloride (Cl − ), sulfate (SO 4 2− ), nitrate (NO 3 − ), and fluoride (F − ), were analyzed using a Thermo Scien-tific™ Dionex™ Aquion™ Ion Chromatography (IC) instrument (Waltham, MA, USA).The concentration of HCO 3 − was determined by titration with methyl orange as the indicator.K + and Na + were quantified using an atomic absorption spectrophotometer (GFA-6880 by Shimadzu Corporation, Kyoto, Japan).Ca 2+ and Mg 2+ were determined by titration with ethylenediaminetetraacetic acid (EDTA).Total nitrogen (TN) was quantified by alkaline potassium persulfate digestion UV spectrophotometry in accordance with the current national standard.The total phosphorus (TP) concentration was quantified by the ammonium molybdate spectrophotometric method.The total iron (Fe) and manganese (Mn) concentrations were determined by inductively coupled plasma spectroscopy (ICP).All samples were analyzed by the Water and Soil Testing Center of the Shandong Provincial Lunan Geology and Exploration Institute.The national standard [18] was used to conduct quality assurance and quality control (QA/QC) of the ion determinations.During the analysis, three replicate water samples were used to assess the precision of the analysis.The analytical precision of the ion charge balance was within 5% for all water samples.
Each sample for sulfate and oxygen isotope was collected at 1.5 L on the site and filtered through cellulose acetate membranes within 24 h.Thereafter, 1 mL of high-purity concentrated hydrochloric acid was added, mixed well, and allowed to stand for 0.5 h.Next, 10 mL of saturated BaCl 2 solution was added to the sample, which was then left to stand overnight.The precipitate of BaSO 4 was filtered by a PES membrane with a pore size of 0.22 µm (or a cellulose acetate membrane with a pore size of 0.45 µm) in the laboratory and then dried at 850 • C for 2 h.Subsequently, purification of BaSO 4 samples was carried out.The dried BaSO 4 solids were transferred and dissolved in a DTPA solution.Once all the precipitate had dissolved, high-purity concentrated hydrochloric acid was added to remove the carbonate.Then, a saturated barium chloride solution was added to re-precipitate the sulfate ions.After washing with ultrapure water to remove BaCl 2 and DTPA, the samples were dried at 60 • C for more than 24 h [19].After purification, the dried BaSO 4 precipitate underwent testing for sulfur and oxygen isotopes in sulfate at Nanjing University.Sulfur and oxygen isotopes were tested using the EA Isolink-Conflo IV-Delta V plus and EA-HT Isolink-Conflo IV-Delta V plus instruments, respectively.The international standard (NBS127, δ 18 OVMSOW = 8.6‰; δ 34 SVCDT = 20.3‰) was chosen to calibrate the samples to be tested, all with an analytical accuracy of ±0.3‰.

Analysis of Stable Isotopes
A Bayesian isotope mixing model was employed to determine the proportions of sulfates in water [20].The fundamental premise is elucidated by the following formula: In this context, k represents the number of potential sources; i represents the number of samples of the mixture, and i = 1, 2, 3, . .., N; j represents the type of stable isotope, and j = 1, 2, 3, . .., J; X ij represents the composition of isotope j in the ith mixture; P k represents the mixing ratios of each potential source, which are the final result calculated by the model.For multi-element isotopes, the model requires the input of the jth isotopic composition range S jk in the kth potential source, represented by mean µ jk and standard deviation ω jk .For defining the isotope fractionation during the process from source to mixture, the fractionation factor C jk of the jth isotope in the kth potential source is input to the model, represented by mean λ jk and variance τ jk .In the absence of isotope fractionation, the fractionation factor C jk may be disregarded.ε ij refers to additional unquantifiable differences between mixture i, typically represented by a mean of 0 and a variance of σ j .

Hydrochemical Characteristics of Water Samples
The statistical analysis of conventional hydrochemical parameters for the surface water and groundwater samples is presented in Table 1.The range of total dissolved solids (TDS) in the surface water samples is 385.00 to 1382.00 mg/L, with an average of 1049.80 mg/L.The most abundant cations are Ca 2+ and Na + , with mean mass concentrations of 183.10 mg/L and 89.00 mg/L, respectively.The main anions are SO 4 2− and HCO 3 − , with mean mass concentrations of 535.88 and 209.08 mg/L respectively.TDS in the pore groundwater samples ranges from 389.00 to 2156.00 mg/L, with an average of 985.78 mg/L, predominantly freshwater.The order of cation mass concentrations is as follows Ca 2+ > Na + > Mg 2+ > K + , with averages of 233.81 mg/L, 49.16 mg/L, 32.87 mg/L, and 1.12 mg/L, respectively.The order of anion mass concentrations is as follows: HCO 3 − > SO 4 2− > Cl − > NO 3 − , with mean values of 373.75 mg/L, 272.01 mg/L, 125.15 mg/L, and 18.48 mg/L, respectively, showing significant differences.NO 3 − concentrations in more than 50 percent of the wells exceeded the Standard for Groundwater Quality (China) [21] Category III, followed by SO 4 2− , exceeding 39%.The range of TDS is 670.00~3620.00mg/L in mine water, with an average of 2064.3 mg/L, and brackish water accounts for 34.7%.The order of cation mass concentrations is as follows: Ca 2+ > Na + > Mg 2+ > K + , with average of 279.08 mg/L, 128.50 mg/L, 51.84 mg/L, and 7.61 mg/L, respectively.The order of anion mass concentrations is as follows: , with mean values of 657.16 mg/L, 424.08 mg/L, 75.00 mg/L, and 14.99 mg/L, respectively, showing great differences.The anion with the highest rate of exceedance of the standard is SO 4 2− , with an exceedance of 78%, followed by NO 3 − , with an exceedance of 44%.The iron content in the mine water samples ranges from 0.044 to 22.032 mg/L, with an average of 2.731 mg/L, while manganese content ranges from 0.0007 to 5.478 mg/L, with an average of 1.898 mg/L.TDS in the karst groundwater samples ranges from 459.00 to 2139.00 mg/L, with an average of 868.96 mg/L.Brackish water accounts for 24%.The order of cation mass concentrations is as follows: Ca 2+ > Na + > Mg 2+ > K + , with the average of 192.70 mg/L, 47.99 mg/L, 32.49 mg/L, and 1.99 mg/L, respectively.The order of anion mass concentrations is as follows HCO 3 − > SO 4 2− > Cl − > NO 3 − (as N), with the main values of 332.18 mg/L, 296.53 mg/L, 66.06 mg/L, and 13.13 mg/L, respectively, showing significant differences.The anion with the exceedance rate is SO 4 2− with 24%, followed by NO 3 − , with 16%.The analysis shows that the dominant cations in groundwater are Ca 2+ and Na + , while the dominant anions in pore and karst groundwater are HCO 3 − , followed by SO 4 2− .In mine water, the dominant anion is SO 4 2− , followed by HCO 3 − .From the coefficient of variation, the coefficients of variation of total phosphorus in surface water, total nitrogen, K + and Na + in mine water, and K + and Na + in karst water were slightly higher than 1.The coefficients of variation of total iron and manganese were larger in pore water, mine water, and karst water (1.86-1.96,1.59-2.43),and the coefficients of variation of the rest of the indexes were all less than 1, which indicated that the spatial distribution of total iron and manganese was extremely unequal; in addition to being affected by mine mining activities, these indexes are also affected by stratigraphic factors and redox environment, etc.

Hydrochemical Types Analysis
A Piper trilinear diagram is a commonly used tool for presenting water chemistry data.Its purpose is to directly represent the relative abundance and distribution characteristics of major ions in water chemistry, thereby assisting in the identification of end-members of hydrochemical formation and evolution [22].As illustrated in Figure 2, the distribution of hydrochemical types is relatively concentrated, with cations predominantly distributed near the Ca 2+ and Mg 2+ ends, while anions are mainly distributed near the HCO 3 − and SO 4 2− ends.In particular, pore groundwater is primarily of the HCO 3 •SO 4 -Ca and HCO 3 -Ca types, accounting for 56.5% and 13.0%, respectively.Karst groundwater and mine water are mainly of the HCO 3 •SO 4 -Ca types, accounting for 80.0% and 77.7%, respectively.This indicates that each water body is mainly influenced by carbonate rock weathering and dissolution, in the form of sulfide mineral oxidation and sulfate mineral dissolution.

Analysis of Water-Rock Interaction of Water Samples
The primary sources of Na + and K + in natural water bodies are the weathering of silicate minerals and evaporite rocks.Cl − remains relatively stable in groundwater and is mainly derived from the weathering and dissolution of evaporite rocks.When Na + , K + , and Cl − in groundwater are derived from the weathering and dissolution of rock salts in

Analysis of Water-Rock Interaction of Water Samples
The primary sources of Na + and K + in natural water bodies are the weathering of silicate minerals and evaporite rocks.Cl − remains relatively stable in groundwater and is mainly derived from the weathering and dissolution of evaporite rocks.When Na + , K + , and Cl − in groundwater are derived from the weathering and dissolution of rock salts in the aquifer, the value of γ(Na + + K + )/Cl − should be equal to 1. Figure 3a illustrates the distribution of sampling points across the study area, which is generally scattered.In pore groundwater, the ratio of (Na + + K + ) to Cl − is less than 1 for most ions, indicating that except for the dissolution of rock salts, Na + and K + can more easily be absorbed by mineral surfaces and Cl − may be influenced by anthropogenic pollution.In contrast, for some karst groundwater, the ratio of (Na + + K + ) to Cl − is greater than 1, indicating that the Na + , K + , and Cl − mainly originate from the weathering and dissolution of evaporite rocks.Nevertheless, γ(Na + + K + )/Cl − is less than 1 in certain locations, suggesting the potential influence of other sources.For mine water and surface water, γ(Na + + K + )/Cl − is more than 1, indicating that Na + and K + in groundwater are not balanced by Cl − .This could be attributed to ion exchange among Na + , Ca 2+ , and Mg 2+ [23], or potential influences from other sources of Na + , such as anthropogenic inputs, etc.
Water 2024, 16, x FOR PEER REVIEW 9 of 14 drochemical characteristics of end-members of silicate rocks, evaporite rocks, and carbonate rocks.The influence of rock weathering on water solutes can be distinguished through relationships between Ca 2+ /Na + and HCO3 − /Na + and Ca 2+ /Na + and Mg 2+ /Na + .Figure 3c,d demonstrate that the sampling points are concentrated between the end-members of evaporite rocks and carbonate rocks.This indicates that the hydrochemical ions in this area are primarily derived from the weathering and dissolution of evaporite and carbonate rocks.

Sulfur and Oxygen Isotope Characteristics
The use of sulfur and oxygen isotope analysis allows for the effective identification of the sources of sulfate ions.The sulfur isotopic composition of sulfate ions, as indicated by δ 34 SSO4, can provide insights into the sources and chemical transformation processes.In contrast, the oxygen isotopic composition of sulfate ions, as indicated by δ 18 OSO4, can provide insights into the redox environment.Table 2 indicates that the average δ 34 SSO4 of pore groundwater, mine water, and karst water are 3.16‰, 0.23‰, and 2.44‰, respectively, while the average δ 18 OSO4 are +3.90‰,2.96‰, and 3.88‰, respectively.Except for mine water wells M04, M05, and M07, which have relatively shallow depths (wells less than 10 m deep) and low water table burial depths, resulting in larger δ 18 OSO4, all other mine water wells have relatively lower δ 34 SSO4 and δ 18 OSO4.The potential sources of SO4 2− include atmospheric precipitation, rock weathering (gypsum dissolution and pyrite oxi- In general, agricultural activities exhibit relatively high NO 3 − /Ca 2+ values and low SO 4 2− /Ca 2+ values, whereas mining activities exhibit the opposite trend [24].Figure 3b illustrates this trend, with surface water and mine water showing low NO 3 − /Ca 2+ values and high SO 4 2− /Ca 2+ values in the study area.The majority of points are distributed in the upper-left area of the diagram.This indicates that surface water and mine water are significantly impacted by mining activities.For pore water and karst water, some points exhibit similar distributions to surface water and mine water, being also located in the upper-left area of the diagram.This suggests that pore water and karst water are also influenced by mining activities in the study area.Based on the hydrochemical characteristics of single lithological small watersheds, Gaillardet et al. [25] summarized the hydrochemical characteristics of end-members of silicate rocks, evaporite rocks, and carbonate rocks.The influence of rock weathering on water solutes can be distinguished through relationships between Ca 2+ /Na + and HCO 3 − /Na + and Ca 2+ /Na + and Mg 2+ /Na + .Figure 3c,d demonstrate that the sampling points are concentrated between the end-members of evaporite rocks and carbonate rocks.This indicates that the hydrochemical ions in this area are primarily derived from the weathering and dissolution of evaporite and carbonate rocks.

Sulfur and Oxygen Isotope Characteristics
The use of sulfur and oxygen isotope analysis allows for the effective identification of the sources of sulfate ions.The sulfur isotopic composition of sulfate ions, as indicated by δ 34 S SO4 , can provide insights into the sources and chemical transformation processes.In contrast, the oxygen isotopic composition of sulfate ions, as indicated by δ 18 O SO4 , can provide insights into the redox environment.Table 2 indicates that the average δ 34 S SO4 of pore groundwater, mine water, and karst water are 3.16‰, 0.23‰, and 2.44‰, respectively, while the average δ 18 O SO4 are +3.90‰,2.96‰, and 3.88‰, respectively.Except for mine water wells M04, M05, and M07, which have relatively shallow depths (wells less than 10 m deep) and low water table burial depths, resulting in larger δ 18 O SO4 , all other mine water wells have relatively lower δ 34 S SO4 and δ 18 O SO4 .The potential sources of SO 4 2− include atmospheric precipitation, rock weathering (gypsum dissolution and pyrite oxidation), and anthropogenic factors [24], among others.Different sources of SO 4 2− have distinct δ 34 S SO4 values.Gypsum in marine sediments typically exhibits higher δ 34 S SO4 values, as exemplified by the Handan Middle Ordovician gypsum, with δ 34 S SO4 values ranging from 20‰ to 24‰ [26].In contrast, pyrite in sedimentary rocks typically exhibits lower δ 34 S SO4 values.The values range from 3.4‰ to 10.2‰ in coal seams in different regions of China [27], while those in magmatic rocks range from −5‰ to +5‰ [26,28].There is no significant isotopic fractionation of sulfur during the oxidation of pyrite, so the δ 34 S SO4 values of SO 4 2− formed from oxidation are similar to those of pyrite.Some studies indicate that the δ 34 S SO4 values of pyrite oxidation products range from −15‰ to 4‰ [29].The minimum δ 34 S SO4 and δ 18 O SO4 values of the samples are negative in the study area, indicating that some of the SO 4 2− may originate from the oxidation of sulfide minerals. 2− in various water bodies within the study area.It is evident that the values of the majority of sample points fall within the range of sulfide mineral oxidation, suggesting that nearby coal mining activities have a significant impact on the sulfate concentration in the regional water bodies.Furthermore, the values of δ 34 S-SO 4 of mine water are lower than those of the other two types of groundwater, indicating a greater influence of mineral oxidation on coal mine water.Based on the relationships between δ 18 O-SO 4 and SO 4 2− , the majority of mine water and some pore water and karst water samples are located in the lower-right area, indicating a significant influence of sulfide mineral oxidation.Some pore water and karst water samples are situated in the upper-right area of gypsum dissolution, indicating that the sulfate concentration in the groundwater is generally higher due to the combined effects of several factors.Moreover, it has been found that the isotopes undergo variations before and after rainfall even at the same sampling point, with an overall increase in atmospheric precipitation after rainfall, and a trend of increasing sulfur and oxygen isotopes.
Figure 4a,b illustrate the relationships between δ 34 S-SO4 and δ 18 O-SO4 and δ 18 S-SO4 and SO4 2− in various water bodies within the study area.It is evident that the values of the majority of sample points fall within the range of sulfide mineral oxidation, suggesting that nearby coal mining activities have a significant impact on the sulfate concentration in the regional water bodies.Furthermore, the values of δ 34 S-SO4 of mine water are lower than those of the other two types of groundwater, indicating a greater influence of mineral oxidation on coal mine water.Based on the relationships between δ 18 O-SO4 and SO4 2− , the majority of mine water and some pore water and karst water samples are located in the lower-right area, indicating a significant influence of sulfide mineral oxidation.Some pore water and karst water samples are situated in the upper-right area of gypsum dissolution, indicating that the sulfate concentration in the groundwater is generally higher due to the combined effects of several factors.Moreover, it has been found that the isotopes undergo variations before and after rainfall even at the same sampling point, with an overall increase in atmospheric precipitation after rainfall, and a trend of increasing sulfur and oxygen isotopes.

Source Identification of Sulfate
The previous analysis indicates that sulfate in surface water and groundwater in the study area primarily originates from sulfide oxidation, atmospheric precipitation, gypsum dissolution, and sewage discharge.The end-member values for sulfate sulfur and oxygen isotopes in atmospheric precipitation were derived from testing in North China, with δ 34 S SO4 = +3‰ and δ 18 O SO4 = +8‰ [30,31].The end-member values of the δ 34 S SO 4 and δ 18 O SO 4 in sulfide minerals were also calculated using the aforementioned average values, resulting in the following current test values for mine water, δ 34 S SO4 = −5‰ and δ 18 O SO4 = −0.5‰[32,33].For gypsum in carbonate-evaporite rock, δ 34 S SO4 is approximately +25‰, while δ 18 O SO4 typically ranges from +15‰ to +20‰.The δ 34 S SO4 value was found to be +25‰, while the δ 18 O SO4 value was 18‰ [34].Previous studies have indicated that for sewage discharge, the δ 34 S SO4 value is 8‰, while the δ 18 O SO4 value is 15‰ [35].A Bayesian stable isotope mixing model based on R programming language was employed to calculate the mixing ratios of different sources (https://cran.r-project.org/web/packages/simmr/vignettes/simmr.htmlaccessed on 20 January 2024).This model represents an enhanced version of the stable isotope mixing model installation package [36].Consequently, the mixing ratios of sulfate sources could be calculated within the study area.
Figure 5a,b illustrate the distributions and sources of karst water and pore water in the study area.The distribution contour of sulfate indicates that the concentration of sulfate in the groundwater is increasing in and surrounding the Taozao coalfield, exceeding the groundwater quality standard of 250 mg/L.In particular, sulfate concentrations at some points in the central part of the coalfield exceed 350 mg/L.The analysis of sulfate sources indicates that the average contribution percentages from sulfide minerals, gypsum dissolution, atmospheric precipitation, and sewage sources in karst water are 64.35%,16.00%, 12.62%, and 7.04%, respectively.This suggests that sulfate in karst water is primarily derived from the oxidation of sulfide minerals, with some locations exhibiting sulfide mineral contributions as high as 86%.The average contribution percentages from sulfide minerals, gypsum dissolution, atmospheric precipitation, and sewage sources in pore water are 65.13%, 18.30%, 8.55%, and 8.00%, respectively.This indicates that sulfate in pore water also primarily originates from the oxidation of sulfide minerals, with some locations showing sulfide mineral contributions as high as 83%.From the distribution of sampling points, it can be observed that in the vicinity of the western boundary of the Taozao coalfield, the ratio of sulfate in Quaternary pore water originating from the oxidation of sulfide minerals is relatively high, while the ratios in pore water are low in most areas.The values of the other three sources of sulfate pollution are generally similar, indicating that these sources of pollution also contribute to varying degrees of contamination in different water bodies.

Conclusions
The study area is the Panlong River Basin in Zaozhuang City, Shandong Province, China.The status of sulfate pollution and its sources are analyzed based on chemical and sulfur-oxygen isotope data.The results indicate that the sulfate in Quaternary ground-

Conclusions
The study area is the Panlong River Basin in Zaozhuang City, Shandong Province, China.The status of sulfate pollution and its sources are analyzed based on chemical and sulfur-oxygen isotope data.The results indicate that the sulfate in Quaternary groundwater exceeds the regulatory limits by 39%, reaching a peak of 666.2 mg/L.Specifically, the pore water near the western boundary of the Taozao coalfield is directly related to mine water recharge.In the downstream area of the Panlong River, Changcun Town, the sulfate exceedance is associated with emissions from paper mills, domestic sewage effluents, and closed coal preparation plants.The sulfate in karst groundwater exceeds regulatory limits by 28%, peaking at 1131 mg/L.It is distributed mainly outside the western boundary of the coalfield and near the Dingzhuang water source in the southeastern part of the coalfield.Its concentration is significantly influenced by mine water with high sulfate levels in coal-bearing strata.Isotope testing results show significant differences in sulfur-oxygen isotopes among surface water, mine water, karst water, and pore water.Bayesian endmember analysis was employed to ascertain the proportions of various sulfate sources.The results indicated that 40-83% of the sulfate in pore water downstream of the Taozao coalfield originated from coal seam water, while 48-86% in karst water was derived from coal seam water, exerting a notable influence on groundwater in the vicinity of the western boundary and southeastern part of the Dingzhuang area within the Taozao coalfield.

Figure 1 .
Figure 1.Distribution of sampling points and formations in the study area.(a) Distribution of sampling points (b) Hydrogeological profile of study area).

Figure 1 .
Figure 1.Distribution of sampling points and formations in the study area.(a) Distribution of sampling points (b) Hydrogeological profile of study area).

Figure 4 .
Figure 4. Distribution characteristics of sulfur and oxygen isotopes in water samples in the study area (The gray areas indicate values quoted from Refs.[26-29]).(a)the relationship of δ 34 S-SO4 and δ 18 O-SO4 (b) the relationship of δ 34 S-SO4 and SO4 2− .

Figure 4 .
Figure 4. Distribution characteristics of sulfur and oxygen isotopes in water samples in the study area (The gray areas indicate values quoted from Refs.[26-29]).(a)the relationship of δ 34 S-SO 4 and δ 18 O-SO 4 (b) the relationship of δ 34 S-SO 4 and SO 4 2− .

Water 2024 , 14 Figure 5 .
Figure 5. Distributions and sources of sulfate in karst water and pore water.(a) Karst water (b) Pore water.

Figure 5 .
Figure 5. Distributions and sources of sulfate in karst water and pore water.(a) Karst water (b) Pore water.
used δ 34 SSO 4 -SO 4 correlation to study the

Table 1 .
Statistical results of major hydrochemical indicators for different water bodies.

Table 2 .
Statistical results of δ 34 S SO4 and δ 18 O SO4 in different water bodies.Figure 4a,b illustrate the relationships between δ 34 S-SO 4 and δ 18 O-SO 4 and δ 18 S-SO 4 and SO 4