Spatial Distribution of Arsenic in the Aksu River Basin, Xinjiang, China: The Cumulative Frequency Curve and Geostatistical Analysis

: The quality of drinking water is crucial for human health and the sustainable development of societies. The Aksu River Basin, a typical inland river system, has areas where groundwater arsenic levels exceed safe drinking water standards (i.e., arsenic concentrations greater than 10 µ g/L). Identifying the causes of high arsenic levels in the basin’s groundwater requires further study. Analyzing the hydrogeochemical composition of the Aksu River basin helps us to understand the spatial distribution of groundwater environments and locate areas with dangerously high arsenic levels. In this research, we collected 196 groundwater samples from along the river. Out of these, 38 samples had arsenic levels above 10 µ g/L, which represents 19.4% of the total samples collected. By examining the slope changes in the cumulative frequency curves of major ion ratios and employing geostatistics (specifically, the Kriging interpolation), and taking into account the environmental characteristics of the entire basin, we divided the study area into five sub-regions (Zone I through Zone V). The geostatistical analysis showed a significant spatial variability in groundwater arsenic levels, with a clear spatial correlation. Our findings demonstrate that arsenic concentrations in the Aksu River basin’s groundwater vary widely, with Zones II and III—mainly located in the northeastern part of the basin and in Awat County—being hotspots for high-arsenic water. Factors such as a weak reducing environment, intense evaporation, strong cation exchange, and the low-permeability recharge of surface water contribute to the accumulation of arsenic in the basin’s groundwater. The results of this study are vital for assessing the risk of arsenic contamination in groundwater in similar basins and for identifying critical areas for further investigation and research.


Introduction
Access to high-quality drinking water is crucial for maintaining human health and is a cornerstone of sustainable development [1,2].Arsenic, a prevalent and harmful element found in groundwater, enters the human body through drinking water and food, posing significant health risks.Long-term consumption of arsenic-contaminated water and the consumption of fruits and vegetables irrigated with such water can lead to irreversible health effects [3,4].The most common condition resulting from prolonged exposure is chronic arsenic poisoning.Alarmingly, continuous intake of arsenic can increase the risk of various organ cancers, including those of the liver and kidneys [5].Data from The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019 Risk Factors Collaborators highlights that unsafe drinking water is a significant global risk factor, ranking 11th for females and 13th for males in causing deaths [6].This emphasizes the need for increased vigilance and protective measures regarding drinking water safety to mitigate the health risks associated with arsenic contamination [3,5,6].
According to statistics, to date, high arsenic groundwater has been discovered in more than 70 countries worldwide, affecting a population of 94 million to 220 million people [7].The majority of affected individuals, about 94%, are located in Asia, primarily in the river and lake basins or river delta regions of arid or semi-arid climates [7,8].Examples of such regions include the Ganges Basin in India and the Red River Delta in Vietnam [9][10][11][12][13], the Bengal Basin in Bangladesh [14,15], and the Hetao Basin in the Yellow River Basin in China [16][17][18].In recent years, with the reporting of research findings in this field, groundwater pollution has received significant attention from researchers as well as decision-makers in various governments.
Arsenic, as a non-metallic harmful element, primarily exists in geological sediments and natural water bodies in the environment.It accumulates in various organisms through the food chain [19].Existing research indicates that the arsenic content in the Earth's crust is relatively low, with an average concentration of 1.5 mg/kg, while the majority of natural water bodies have arsenic concentrations mostly below 1.0 µg/L [7,19].With the development of the Second and Third Industrial Revolutions, human exploitation of natural resources has increased significantly and become more refined.
Currently, over 200 naturally occurring minerals containing arsenic are known to exist [20,21].Due to the development of mining industries, unquantifiable harmful elements, including arsenic, have transitioned from stable states to active states, entering the natural cycling system [22].They undergo processes such as leaching, dissolution, evaporation concentration, and the cycling of the food chain, leading to iterative enrichment [23,24].With a development philosophy centered around people's health, both China and the World Health Organization (WHO) have established and adopted safe guideline values for arsenic content in drinking water (less than 10 µg/L) in 2006 and 2011, respectively [25,26].
While some researchers have speculated about the potential impacts of groundwater arsenic contamination in Xinjiang, China, the region's intricate terrain and geological makeup have historically rendered it a "black-box" for in-depth study [27].Specifically, research has been notably scarce in the Aksu River basin, a critical area for both agriculture and industry in southern Xinjiang, which also boasts a dense population [28].The inhabitants of this region depend heavily on the basin for their daily living and agricultural needs.Alongside China's rapid development, Xinjiang's industrial and agricultural sectors have seen significant growth and enhancement.Nevertheless, this progress has come at a cost, particularly in the form of relentless and widespread groundwater extraction, culminating in a marked deterioration of water quality, most notably in the downstream regions of the basin.
This study utilizes geostatistical Kriging interpolation, guided by the slope of cumulative frequency curves for arsenic content and major ion ratios in the study area's groundwater, to achieve three primary objectives: Firstly, to examine the hydrogeochemical processes and spatial distributions of arsenic across various groundwater types within the study area.Secondly, to introduce a methodology for pinpointing the spatial distributions of high arsenic contents by analyzing the cumulative frequency curves of ion ratios.And lastly, to explore the mechanisms behind arsenic enrichment in distinct regions, classified according to the slopes of various ion ratios.

Geographic Feature
The Aksu River is located in Aksu Prefecture, Xinjiang, China, and it belongs to the northern edge of the Tarim Basin.The overall terrain of the region is higher in the northwest and lower in the southeast.The highest elevation in the study area is 2759 m, located in the northwest, which is a branch of the Tianshan Mountains called the Latajiayi Mountain.This mountain range extends in a northwesterly direction.The lowest elevation is a fluvial plain at 1100 m, located in the southern part of the study area, with a slope ranging from 1 to 3‰.The maximum elevation difference in the study area is 1659 m.

Hydrology and Meteorology
The Aksu River basin is characterized by a typical arid to semi-arid warm temperate continental climate.Within the study area, the annual average temperature and evaporation increase with decreasing terrain elevation, while precipitation increases with higher terrain elevation.The Aksu River is the main tributary of the Tarim River, accounting for 75% of its flow.The Aksu River is formed by the confluence of the Toshkan River (western branch) and the Kumalak River (northern branch) within Wensu County.After the confluence, it flows through Aksu City and eventually joins the Tarim River to the south.Both rivers originate from the western segment of the Tianshan Mountains, with water supply mainly originating from mountain precipitation and snowmelt (Figure 1).

Hydrogeological Condition
The Aksu River basin is characterized by the presence of thick Quaternary unconsolidated sediments, with a deposition thickness ranging from 1000 to 1500 m in the alluvial plain.In the piedmont alluvial fan area of the basin, the particle size of the Quaternary sediments gradually decreases from the piedmont towards the front edge of the alluvial plain.The lithology of the Quaternary sediments in the upper part of the alluvial plain is composed of conglomerates and gravel, while the middle part consists of gravel and coarse sand.The lower part of the alluvial plain is composed of fine sand, silt, silty clay, and clayey silt.In the Tarim River alluvial plain, the Quaternary sediments are mainly composed of fine sand, silt, silty clay, and clayey silt.

Sample Collection and Analysis
A total of 196 groundwater samples were collected with the assistance of the First Regional Geological Survey Brigade of the Xinjiang Uygur Autonomous Region Geology and Mineral Exploration and Development Bureau.Each groundwater sample was taken from below 0.5 m below of the groundwater surface, with the location of each sample recorded using GPS technology.Groundwater was collected from each site in duplicate, stored in plastic bottles with airtight caps, and labelled appropriately.Membrane filtration (0.45 µm) of all samples was conducted on site.One of the duplicate samples was acidified on site using 2-3 drops of concentrated nitric acid (HNO 3 , 69%, Electron-grade) to ensure a pH value lower than 2 and to stabilize arsenic and metal ions and prevent their precipitation.This acidified water sample was subsequently used to analyze the total arsenic content and other elements.The second duplicate sample, which was not acidified, was utilized to analyze various cations and anions.All samples were transported to the laboratory at 4 • C, and analyzed within 7 days.
The water temperature, pH and total dissolved solids (TDS) were tested using a sensION+MM150 (MM156) portable multi-parameter water quality analyzer produced by HACH Company on the sampling day, and the turbidity was tested using a WGZ-200B portable turbidimeter (meter) produced by the Shanghai Xingrui Company.At the laboratory of the School of Resources and Geosciences, China University of Mining and Technology, the arsenic content in groundwater was measured via inductively coupled plasma-mass spectroscopy (ICP-MS) (7500C, Agilent, Santa Clara, CA, USA).K + , Na + , Ca 2+ , Mg 2+ , and other cations were measured via inductively coupled plasma-atomic emission spectroscopy (ICP-AES) and inductively coupled plasma-mass spectroscopy (ICP-MS).Anions (SO 4 2− and Cl − ) were analyzed using ICS1500 (Dionex, Sunnyvale, CA, USA), and HCO 3 was determined via titration.The concentration errors of both cations and anions in all samples were within a range of less than 5%, ensuring the accuracy of all the data.

The Partition Method of Cumulative Frequency Curve
Cumulative frequency profiles have been proved to be a valuable tool for characterizing distributions, providing a visual indication of the cumulative distribution of different hydrogeochemical component ratios [29,30].With the increases in both the ratio and the number of samples, the cumulative frequency gradually increased from 0% to 100%.Depending on the point of inflection, i.e., the point showing significant changes in the ratios, all data are represented by multiple straight lines with different slopes.These changes in slope reflect variations in the groundwater environment and hydrogeochemical processes.In this study, cumulative frequency curves of ion ratios were plotted and used in conjunction with geostatistical Kriging interpolation to clarify differences in the hydrogeochemical compositions of groundwater samples and reveal specific chemical processes in groundwater across different regions.

Hydrogeochemical Distribution of Samples
The concentrations of ions and TDS in the 196 samples range from two to four orders of magnitude, indicating strong spatial variability in the groundwater (Table 1).Specifically, the pH of the groundwater in the Aksu River basin was neutral to weakly alkaline, ranging from 7.05 to 8.88, with a mean value of 8.15.The highest concentrations of anions and cations were Cl − and Na + , respectively, with median values of 266 and 197 mg/L.The TDS content ranged between 200 and 27,024 mg/L, with a median value of 1234 mg/L; brine (TDS > 1000 mg/L) accounted for 53.6%.
The Piper diagram revealed that the hydrogeochemical type of groundwater in the study area consists of Na-HCO 3 •Cl and Na-HCO 3 (SO 4 ), while the chemical type in the pre-mountain area and the middle and lower reaches of the river alluvial plain in the plain area of Aksu region is Na-HCO 3 •SO 4 .The chemistry of high arsenic groundwater is mainly of the Na-SO 4 -Cl type.In high arsenic groundwater, sodium (Na + ) concentrations vary between 41.0 and 7532 mg/L (median: 280 mg/L), bicarbonate (HCO 3 − ) between 21.4 and 488 mg/L (median: 143 mg/L), calcium (Ca 2+ ) from 20.7 to 943 mg/L (median: 89.2 mg/L), and magnesium (Mg 2+ ) from 17.4 to 8337 mg/L (median: 197 mg/L).Total dissolved solids (TDS) range from 236 to 27,024 mg/L, with median and mean values of 1326 and 3672 mg/L, respectively.These values indicate a high level of salinization in arsenic-rich groundwater, particularly in groundwater with arsenic concentrations exceeding 30 µg/L, where the dominant hydrogeochemical type is Na-SO 4 -Cl (Figure 2).

Spatial Distribution of Arsenic
The total arsenic concentration in groundwater ranged from 0.05 to 210 µg/L, with a mean value of 8.70 µg/L and a median value of 2.05 µg/L (Table 1).High arsenic (As > 10 µg/L) groundwater accounted for 19.4% of the total.To further understand the distribution of arsenic in groundwater in different geomorphological units, the 196 groundwater samples were divided into three groups based on arsenic concentrations (i.e., 0-10 µg/L, 10-30 µg/L, and >30 µg/L), and raw data were analyzed using ArcGIS.⃝ and 3 ⃝ have the widest distribution area, with 2 ⃝ mainly being mainly the area where the Aksu River, Yarkant River, and Hotan River converge with the Tarim River, and Zone 3  ⃝ mainly being mainly distributed downstream of the Kumala River and Taushgan Darya, extending to Aksu City.Zones 4  ⃝ and 5 ⃝ are mainly distributed in parts of Aksu City and Onsu County near the Aksu River basin (Figure 4a).The majority of groundwater samples in the study area exhibited (HCO 3 < 2.83, and the arsenic concentration was found to be positively correlated with the ratio, so the groundwater in this area is a weakly reducing environment overall.The results of sample tests suggest that such an environment is influenced by the weakly reducing conditions.The majority of groundwater samples in the study area exhibited (HCO3 − + CO3 2− )/SO4 2− < 2.83, and the arsenic concentration was found to be positively correlated with the ratio, so the groundwater in this area is a weakly reducing environment overall.The results of sample tests suggest that such an environment is influenced by the weakly reducing conditions.

Ca 2+ /(HCO3
The concentration of Ca2+ in groundwater can serve as an indicator of the degree to which surface water is recharging the groundwater [17].Moreover, the increase in the Ca 2+ /(HCO3 − + CO3 2− ) ratio can be attributed to the dissolution of silicate minerals [32].The ratio is significantly influenced by the mixing of surface water with groundwater or rainfall, as well as cation exchange.A higher ion ratio indicates stronger cation exchange and lower mixing ratios of surface water or rainfall.The correlation between arsenic levels and the Ca 2+ /(HCO3 − + CO3 2− ) ratio was positive for the study area, with a greater correlation for high arsenic groundwater than for all samples combined.The study area was classified into five zones (①, ②, ③, ④, and ⑤) based on the slope variation in the cumulative frequency curve for Ca 2+ /(HCO3 − + CO3 2− ) ratio (Figure 3b).Zone ① (0-0.62) was predominantly found in the upper reaches of the Aksu River, as well as in parts of Aksu City and Onsu County near the Aksu River basin.Zone ② (0.62-1.50) had the widest distribution area and was found in areas outside of the other groups.Zone ③ (1.50-3.00)was mainly distributed in the areas of the Duolanghu Reservoir, the Xinjingzi Reservoir, and the confluence area of the Yarkant River and Tarim River.Zones ④ (3.00-4.50)and ⑤ (>4.50) were mainly distributed in the downstream area of the Tailan River, i.e., parts of the Duolanghu Reservoir (Figure 4b).
For most of the underground samples, Ca 2+ /(HCO3 − + CO3 2− ) > 0.62, suggesting weak surface water recharge and strong cation exchange.The enrichment of arsenic in the groundwater may also be related to these factors.The concentration of Ca2+ in groundwater can serve as an indicator of the degree to which surface water is recharging the groundwater [17].Moreover, the increase in the Ca 2+ /(HCO 3 − + CO 3 2− ) ratio can be attributed to the dissolution of silicate minerals [32].The ratio is significantly influenced by the mixing of surface water with groundwater or rainfall, as well as cation exchange.A higher ion ratio indicates stronger cation exchange and lower mixing ratios of surface water or rainfall.The correlation between arsenic levels and the Ca 2+ /(HCO 3 − + CO 3 2− ) ratio was positive for the study area, with a greater correlation for high arsenic groundwater than for all samples combined.The study area was classified into five zones ( 1⃝, 2 ⃝, 3 ⃝, 4 ⃝, and 5 ⃝) based on the slope variation in the cumulative frequency curve for Ca 2+ /(HCO 3 − + CO 3 2− ) ratio (Figure 3b).Zone 1 ⃝ (0-0.62) was predominantly found in the upper reaches of the Aksu River, as well as in parts of Aksu City and Onsu County near the Aksu River basin.Zone 2 ⃝ (0.62-1.50) had the widest distribution area and was found in areas outside of the other groups.Zone 3 ⃝ (1.50-3.00)was mainly distributed in the areas of the Duolanghu Reservoir, the Xinjingzi Reservoir, and the confluence area of the Yarkant River and Tarim River.Zones 4 ⃝ (3.00-4.50)and 5 ⃝ (>4.50) were mainly distributed in the downstream area of the Tailan River, i.e., parts of the Duolanghu Reservoir (Figure 4b).For most of the underground samples, Ca 2+ /(HCO 3 − + CO 3 2− ) > 0.62, suggesting weak surface water recharge and strong cation exchange.The enrichment of arsenic in the groundwater may also be related to these factors.

Ca 2+ /Mg 2+
In different lithological aquifers, dissolution and leaching can result in Ca 2+ /Mg 2+ ratio differences.As the ratio approaches 1 to 2, groundwater is primarily dissolved by carbonates such as calcite and dolomite, while a further increase in the ratio reflects the fact that the leaching of silicate minerals also played a role [33].Therefore, the Ca 2+ /Mg 2+ ratio reflects the intensity of lateral recharge, as an increase in the former can indicate an enhancement of the latter.The study area was divided into four zones, labeled 1  ⃝, 2 ⃝, 3 ⃝, and 4 ⃝, according to the slope change in the Ca 2+ /Mg 2+ cumulative frequency curve (Figure 3c).Zone 1  ⃝ (0-3.80) is mainly distributed in the upper reaches of the Aksu River and part of the lower reaches of the Kumala and Taushgan Darya rivers, specifically parts of Aksu City and Onsu County near the Aksu River basin.Zone 2 ⃝ (3.80-7.60)has a narrower distribution range and is mainly distributed in the periphery of Zone.Zones 3 ⃝ (7.60-11.4)and 4  ⃝ (>11.4) have the widest distribution area, including the upper reaches of the Kumala River basin, the lower reaches of the Taushgan Darya, the upper reaches of the Tailan River, and parts of the Yarkant and Tarim River basins in China (Figure 4c).
In the study area, there are 51 samples (26%) with Ca 2+ /Mg 2+ ratios between 1 and 2, and 97 groups (49.5%) with Ca 2+ /Mg 2+ ratios exceeding 2. This suggests that the dissolution of carbonates such as calcite and dolomite, as well as the leaching of silicate minerals, is affecting the region.

Na + /Ca 2+
Poor local runoff conditions can lead to cation exchange between groundwater and aquifer minerals.The Na + /Ca 2+ ratio can be used to evaluate the extent of cation exchange in groundwater [34].As cation exchange intensifies, the Na content and Na + /Ca 2+ ratio gradually increase.Based on the Na + /Ca 2+ cumulative probability curve, the study area can be classified into five zones and labeled as 1  ⃝, 2 ⃝, 3 ⃝, 4 ⃝, and 5 ⃝ (Figure 3d).Zone ⃝ (>4.80) is primarily distributed in the lower reaches of the Tailan River, including the Dolanghu reservoir area (Figure 4d).

Relationship between Different Ion Ratios and Spatial Distribution of Arsenic
Hydrogeochemical processes in the groundwater environment are the primary drivers of arsenic presence.The spatial distribution of different ion ratios is key to identifying these processes and the mechanisms of arsenic distribution in the study area.By analyzing the concentration characteristics and spatial distributions of these ion ratios, the final boundaries of each zone are determined using geostatistical techniques, including kriging interpolation.This approach delineates five distinct zones (I-V) (Figure 5).Zone I was characterized by 0.01 < (HCO 3 − + CO 3 2− )/SO 4 2− < 0.4, 0.43 < Ca 2+ /(HCO 3 − + CO 3 2− ) <6.52, 0.09 < Ca 2+ /Mg 2+ < 11.5, 0.60 < Na + /Ca 2+ < 6.35.Zone I is more sporadically distributed, mainly in Xinjingzi Reservoir and Awat County (Figure 5), where the concentrations of Na + (mean and median: 2097 and 550 mg/L), Ca 2+ (483 and 519 mg/L), and Mg 2+ (353 and 219 mg/L) are the highest among the five areas (Figure 6a-c).In zone I, The concentrations of HCO 3 − (mean and median: 308 and 275 mg/L), SO 4 2− (4157 and 5240 mg/L) and Cl − (7715 and 6541 mg/L) are the highest among the five zones (Figure 6d-f).The mean and median arsenic concentrations in the area are 6.49 and 2.34 µg/L, respectively, and the majority of samples in the area had arsenic concentrations of less than 10 µg/L (82.4%) (Figure 7).Overall, the area appears to be in a weakly reducing, strong lateral/surface water recharge, and moderate cation exchange environment (Table 2).Zone II is located in the lower reaches of the Aksu River basin and the alluvial flood plain in the confluence basin of Aksu River, Yarkant River and Tarim River (Figure 5), for which the concentrations of Na + (mean and median: 661 and 302 mg/L), Ca 2+ (238 and 165 mg/L), and Mg 2+ (171 and 128 mg/L) (Figure 6a-c).It shows that the cation exchange capacity in this area is medium, and the groundwater runoff condition is general.The mean and median arsenic concentrations in the area are 10.51 and 3.97 µg/L, and in most samples, the value is smaller than 10 µg/L (Figure 7).
In the study area, Zone I exhibits the lowest (HCO 3 − + CO 3 2− )/SO 4 2− ratio, indicating the highest degree of groundwater oxidation [35].This finding suggests the possible influence of neighboring surface water quality on groundwater.The mixing of groundwater with surface water containing higher oxygen levels may result in increased groundwater oxidation.Moreover, the smallest Na + /Ca 2+ ratio in Zone V implies weak cation exchange between Ca 2+ and Na + within the aquifer and favorable runoff conditions.
As the topographic gradient declines from Zone V to I, groundwater runoff slows, and sediment becomes finer.This slowdown, combined with enhanced evaporation and water mineralization, leads to groundwater discharge mainly through evaporation, lateral runoff, and human extraction, diminishing river recharge.Arsenic concentrations in groundwater are highest in Zone III and lowest in Zone V, with zones II and III showing similar levels.
The areas experiencing abnormal arsenic enrichment are predominantly located in the northeastern part of the basin and in Awat County.Specifically, in the northeast, unsafe arsenic levels are primarily found between the Tailan River and the Karayulgun River, influenced by the stratigraphy of upstream outflows.In the northern sector, the terrain is steep with coarse sandy gravel, facilitating good water runoff.However, moving southward, the geological composition transitions to finer materials such as clay, mud, fine sand, and silty sand, resulting in poorer runoff conditions.This region's history of frequent marsh formation has led to a high content of organic matter, including humus, which acts as a significant source for arsenic accumulation.Additionally, the combination of flat terrain and intense evaporation further exacerbates arsenic concentrations in the groundwater.
Awat County, on the other hand, is an alluvial-lacustrine deposit soil plain with flat topography, featuring a slope of less than 1‰.Slow groundwater runoff often results in the formation of wet-sinking depressions on the surface.A series of factors, including poor water flow, strong evaporation, saline development, and long-term sediment saturation, can promote the release of arsenic from the sediments into the water and leads to a groundwater anomaly that is high with arsenic concentrations.

Hydrogeochemical Processes
Worldwide, high arsenic levels in groundwater are often accompanied by high HCO 3 − concentrations [10,15,17].Hence, HCO 3 − /SO 4 2− can indicate the degree of reducing conditions in the study area, a key factor in groundwater arsenic enrichment.Groundwater samples with lower arsenic levels typically have HCO 3 − /SO 4 2 ratios greater than 1, whereas samples with higher arsenic levels often show ratios less than 1.The concentration of groundwater arsenic is positively associated with HCO 3 − /SO 4 2− , which suggests that a weakly reducing environment plays a major role in generating the high concentrations of arsenic (Figures 6d,e and 8a).
Aside from reducing conditions, the high arsenic groundwater in the study area may also be related to the evaporative concentration.Due to its high solubility, Cl − is difficult to be removed from the groundwater environment, and the only way increase its concentration is through the dissolution of rock minerals in the aquifer.SO 4 2− is often considered to be a sensitive component of redox conditions, and Cl − /SO 4 2− can indicate the degrees of both the environmental redox and the groundwater evolution [36].The higharsenic groundwater in the study area has Cl − /SO 4 2− > 1 (47.4%)(Figures 6e,f and 8b), showing a weak oxidation environment.
Ca 2+ /(HCO 3 − + CO 3 2− ) reflects the extent to which Ca-HCO 3 type surface water recharges groundwater [37].The higher the ionic ratio, the stronger the cation exchange and the lower the ratio of surface water or rainfall mixing.Strong surface water/rainfall recharge is typically present in zones IV and V of the study area (Figures 6a-d and 8c).In contrast, the weakest lateral/surface water recharge effect is present in zones I and II, along with strong cation exchange.Weak recharge and strong cation exchange may also be associated with arsenic enrichment in groundwater in this study area.No significant correlation between F − and As was found in the study area (Figure 8d).Highly fluorinated groundwater is mainly distributed in Zones III and IV, where the Na + content is high.High Ca 2+ is not conducive to fluoride enrichment, as Ca 2+ and F − in groundwater can form fluorite precipitates.The enrichment of arsenic in Zones III and IV is facilitated by the shallow groundwater levels and significant evaporative concentration within the study area.However, the release of arsenic into groundwater is hindered in Zone V due to weak runoff and oxidizing environments.The characteristics and ion ratios in Zone IV indicate that the conditions in this transition zone are conducive to the release of arsenic (Figure 9).
The analysis of arsenic concentration and its spatial distribution within the study area indicates that the depositional environment plays a significant role in influencing arsenic levels.Notably, high arsenic concentrations in groundwater appear sporadically, primarily in the eastern part of the region.The Aksu area, characterized by an arid to semi-arid inland climate, experiences low rainfall and high evaporation rates, where evaporation surpasses precipitation.In zones of arsenic enrichment, the groundwater table is notably high, leading to predominantly vertical rather than lateral water flow.These conditions, along with intense weathering processes, accelerate decomposition rates and biogeochemical cycles, creating conducive conditions for arsenic migration and accumulation.Arsenic is mobilized from arsenic-rich soil particles and organic matter into the groundwater, where it becomes concentrated due to ongoing evaporation.Significant arsenic levels were detected in groundwater at various locations, including Awat County, Kezile Township, Kezile Kumu Village, and Gulawati Township.Notably, water samples from the vicinity of the Dolanghu Reservoir exhibited extremely high arsenic levels, reaching up to 210 µg/L, far exceeding the standard levels.Conversely, in several other regions, arsenic concentrations generally remain below 10 µg/L.

Conclusions
This study explored the spatial distribution of arsenic within China's Aksu River basin by employing cumulative frequency curves and geostatistical analysis.It found a wide range in groundwater arsenic levels, with the most severe contamination in the northeast and Awat County, highlighting the strong spatial variability and significant correlation of arsenic concentration.
Geostatistical analysis enabled the classification of the study area into three zones: areas with high arsenic anomalies, areas with moderate arsenic levels, and areas with low arsenic presence.The high arsenic zones, mainly in zones II and III, are characterized by a weak reducing environment, gentle terrain slopes (less than 1‰), slow groundwater movement, strong cation exchange, and significant evaporation, coupled with limited surface water/rainwater recharge capacity.Zone IV corresponds to the moderate arsenic area, while zones I and V are identified as low arsenic zones.The process of arsenic migration is intricate, and its exact pathways remain partially understood.No single mechanism adequately accounts for the migration and accumulation of arsenic; a combination of processes provides a more plausible explanation for arsenic release into groundwater.This research significantly contributes to the assessment of arsenic contamination risks in similar hydrogeological settings and prioritizes areas for further investigation.

Figure 1 .
Figure 1.Spatial distribution of As and the sedimentary structure of the Aksu River basin.

Figure 2 .
Figure 2. Piper diagram of shallow groundwater in the study area.

1 ⃝
(0-0.40) is primarily located within Awat County and the region stretching from the confluence of the Hotan River to the Tarim River.Zone 2 ⃝ (0.40-1.60) has a scattered distribution, mainly in the vicinity of Aksu City, as well as the Dolanghu Reservoir and Xinjingzi Reservoir areas.Zone 3 ⃝ (1.60-3.20)has the broadest distribution and is present in areas other than those covered by the other groups.Zone 4 ⃝ (3.20-4.80) is situated in the upper basin of the Tailan River, part of Onsu County that is close to the Kumala River basin, and in portions of Awat County near the Yarkant River basin.Zone 5

Figure 5 .
Figure 5.Ion ratio distribution and regional recognition in groundwater.

Figure 5 .
Figure 5. Ion ratio distribution and regional recognition in groundwater.

Figure 6 .
Figure 6.Main ions-box plot of different intensity zons in aquifers.

Figure 7 .
Figure 7. Box plot of arsenic concentration in different zones.

Figure 9 .
Figure 9. Hydrogeochemical model of groundwater in the study area.

Table 1 .
Distribution of main ion concentrations in groundwater.

Table 2 .
The meaning of different ion ratios and the distribution of the arsenic concentration.