In arid and semi-arid areas, groundwater flow is a potent geological agent. The typical profile of Aksu river basin was chosen as the research object. First, the survey region’s geological background and hydrological conditions were systematically analysed. Combined with geophysical and remote sensing characteristics, the thickness change of the Quaternary loose layer, boundary of the subaqueous system, and basement of the water circulation depth were revealed. Hydrogeochemistry and environmental isotopes were applied to explain the recharge-runoff-discharge process, water-rock interaction, movement law, and residence time of surface-groundwater. Owing to the dissolution of carbonate and gypsum minerals and evaporates, cation exchange between Ca (Mg ) and Na (K ), and the evaporation-concentration effect, specific ions (SO , Cl , Na ) and [total dissolved solids](javascript:;)(TDS) from the surface to groundwater gradually increased with the flow direction and gradually decreased with depth. This difference was more evident in the downstream discharge area, which indicated that they belonged to different groundwater flow systems(GFSs).Because of the elevation effect of the recharge source, the δ O values in different sections showed different degrees of depletion in the horizontal and vertical directions. In [association with](javascript:;) the structural ‘one convex and two concaves’ frameworks and the particle size variation of loose sediments, it reflected the distribution pattern of multiple local GFSs. The vertical zonation between the H and C isotope concentrations and the recorded groundwater [residence](javascript:;) time (modern-24000years) further illustrated the existence of intermediate and regional GFSs. Three surface-conversion boundary key zone (GFS ) were identified, and the GFS conceptual model was established. Finally, the corresponding relationship between the GFSs and the environmental effects, such as the distribution and aggregation of phreatic water with high F and As and soil salinization, were analysed, which had important theoretical significance for protecting the ecological balance of Aksu River basin.

Knowledge of the spatial distribution characteristics of hydraulic parameters is essential for the management and protection of karst groundwater resources. In this study, we propose a workflow integrating the transition probability geostatistics (T‐PROGS) and the discrete cosine transform (DCT) with the ensemble smoother with multiple data assimilation (ES‐MDA) method to map the hydraulic conductivity of karst aquifers that usually follow a multimodal distribution. The priori parameter ensemble is constructed from the T‐PROGS and transformed into an approximately Gaussian distributed coefficient field containing critical spatial structural features about the aquifer using DCT. The ES‐MDA method is employed to update the DCT coefficients by assimilating the measurements. In addition, a postprocessing process based on cumulative distribution function (CDF) mapping is used to address the problem of parameters gradually tending to a Gaussian distribution after updating using the ES‐MDA and the inappropriate selection of initial parameter values. In practice, the limited amount of available data makes it difficult to fully capture the spatial distribution of the parameters in the initial ensemble of a single realization. Therefore, we suggest a new strategy for structuring the initial ensemble by mixing samples from multiple realizations. We then apply the proposed approach to four single realization models and a combined multiple realizations model in a field hydraulic tomography investigation of a real karst aquifer. The computed results show that this method can effectively identify the characteristics of the spatial distribution of hydraulic conductivity in karst aquifers. Compared with an individual realization model, the uncertainty of the hydraulic conductivity estimated by the combined multiple realizations model is significantly reduced. Therefore, including more uncertainty of the aquifer in the initial ensemble is beneficial to improve the accuracy of the parameter estimation.

Studying groundwater flow systems is important for water resources management, for pollution prevention and for maintaining the ecological balance in arid and semi-arid areas. Systematic geophysics and hydrogeological investigations allow us to define the thickness of the Quaternary sedimentary layer, the lateral boundary of the groundwater system, and the depth and basement of water circulation. Hydrogeochemistry and environmental isotopes are used to gain insights into the recharge process, water-rock interactions, hydraulic characteristics and groundwater retention time and to identify groundwater flow systems at all levels in the Aksu River Basin. Owing to the dissolution of carbonate and gypsum minerals and evaporites, cation exchange between Ca (Mg ) and Na (K ), and the evaporation-concentration effect, concentrations of specific ions (SO , Cl , Na ) and [total dissolved solids](javascript:;) (TDS) gradually increase along the flow direction and decrease with depth (indicating that they belong to different groundwater flow systems (GFSs)). Furthermore, interpretation of stable isotope concentrations such as δ O values suggests different degrees of depletion in the horizontal and vertical directions. Combined with the unique structural framework (namely the Wensu uplift, Wushi sag, and Awat sag), the particle size variation of loose sediments and the distribution and aggregation of phreatic water with high F and As and soil salinization show the existence of the surface-ground water interaction and the distribution pattern of multiple local GFSs. The vertical zonation of H and C isotope concentrations and estimates of groundwater [residence](javascript:;) time (modern to 24000 years) further illustrate the hydrodynamic cycle of the local and regional GFSs. The hydrodynamic and hydrochemical characteristics confirmed the distribution of GFSs and the complex mixing relationships between GFSs in the Aksu River Basin under the tectonic conditions since the Neogene in the South Tianshan Mountains.