Spatial Analysis of the Chemical Regime of Groundwater in the Karatal Irrigation Massif in South-Eastern Kazakhstan

: Understanding the groundwater chemical regime on irrigated lands is relevant for assessing hydrogeological and reclamation processes. Water chemistry, natural and climatic conditions, and drainage characteristics deﬁne the complex hydrochemical connections which are formed between irrigation waters, soil, and water of the unsaturated zone and groundwater. This research paper presents the results of studies of seasonal changes in the chemical composition of groundwater in the Karatal irrigated massif. A detailed analysis of the salt regime and chemical composition of groundwater have been performed based on a large number of samples for 2019. Besides, data related to the Karatal River and drainage waters quality inside the basin have also been analyzed. Results demonstrate that the salinization processes in the massif are determined by hydrogeological conditions and drainage availability. In the north-eastern part of the massif, a high risk of soil salinization has been observed as a consequence of all the salinization processes which are determined by its hydrogeological conditions and drainage availability. Data analysis over soil samples and irrigation and drainage waters show that the hydrochemical regime of groundwater on Karatal irrigation massif is inﬂuenced by the presence of saline soil-forming rocks, the groundwater depth, the volume, and the quality of irrigation waters, the characteristics of the natural drainage, and evaporation processes.


Introduction
Increasing salt concentrations is one of the most common environmental and reclamation problems in irrigated agriculture [1][2][3]. Due to the variability of water management and economic conditions on irrigation systems, groundwater levels and geochemistry transformations occur [4]. These changes have an impact on the reclamation state, which determines the productivity of the irrigated lands. The knowledge of the groundwater chemical regime is critical for assessing the reclamation status of irrigated land. Groundwater chemical composition depends on: (i) the irrigation water quality; (ii) the water supply volume; (iii) the climatic and hydrogeological conditions; (iv) the lithological composition of the aeration zone and (v) the drainage characteristics of the territory [5][6][7][8][9].
Hydrogeological conditions play an important role in the formation of all processes of salt and water exchange on irrigated lands. Some recent studies describe that irrigation with superficial water induces a salinity increase in groundwater [10,11]. However, other studies show that this type of irrigation reduces the salinization impacts over the aquifer [12,13]. So, groundwater levels and soil geochemistry are closely interrelated with the soil cover characteristics and the irrigation conditions. Despite water resources being very limited in arid regions, groundwater quality is determined by its interaction with rocks and soil profiles [14].
Groundwater, depending on its depth and degree of mineralization, drainage conditions, lithological features, and other factors, represents the main source of salinization of the upper soil cover. Therefore, the study of the groundwater regime on irrigated lands, the identification of all factors influencing them, the determination of the conditions of formation by type of regime, and the spatial boundaries of the impact should be given important importance [15][16][17]. Based on these data, it is possible to predict changes in hydrogeological and meliorative conditions, develop projects to improve the condition of irrigated lands, and engineering and ecological reconstruction of irrigation systems.
Monitoring of irrigated lands is a system of stationary observations of the reclamation condition of irrigated lands, carried out in order to timely identify changes occurring on them, assess and make recommendations for improving and rational use of irrigated lands. Efficient and well-managed irrigation water practice including quantifying the groundwater table depth and monitoring water quality is of utmost importance for decisionmakers, growers, and experts [18]. Indicators of the reclamation state of irrigated lands are: (i) hydrogeological conditions (level, mineralization, and chemical composition of groundwater); (ii) salinity and degree of salinity of soils; (iii) engineering and geological processes (subsidence, overexposure, irrigation erosion).
Monitoring irrigated lands allow assessing and forecasting their reclamation conditions, identifying the causes of their deterioration as well as changes in the territories in the zone of direct influence of irrigation.
The reclamation state of irrigated lands determines the state of the system: soils, soils of the aeration zone, groundwater formed as a result of soil development, hydrogeological and engineering-geological processes under the influence of various reclamation measures and characterizing the degree of suitability of land for agricultural production [16]. The content of chemical elements such as Cl, Na, as well as pH and mineralization indicators do not affect the use of land at certain stages but can affect the productivity of land [19][20][21]. Lime can be used to increase soil pH, reducing trace element mobility and bioavailability (like Cd, Co, Cu, Ni, Pb, Sb, and Zn) due to enhanced adsorption and/or precipitation [22]. Liming enhances the physical, chemical, and biological properties of soil through its direct effect on the amelioration of soil acidity and its indirect effect on the mobilization of some of the major and trace element nutrients, immobilization of toxic metals, and improvements in soil structure and hydraulic conductivity [23].
The impact of agricultural practices on water quality has always been studied primarily with an emphasis on surface waters. Groundwater is more difficult to study. However, the impact on ground and surface waters is equally important. In non-irrigated areas, agriculture often leads to an increase in recharge, which sometimes leads to the leaching of salts from the unsaturated zone into groundwater. In irrigated areas, groundwater salinization may be the result of salty water irrigation, salt movement down in an unsaturated zone of salt dissolution, as well as salt accumulation due to water uptake by plants.
The processes of redistribution of elements occur in the soil aeration zone. Rising groundwater levels and salinity are considered to be the main cause of secondary salinization of irrigated lands [24,25]. To date, 210,400 hectares of the Almaty region in Kazakhstan are saline to varying degrees which is 36% of the available irrigated area. The deterioration of the reclamation status of rice irrigation systems as a result of changes in irrigation and economic conditions is a matter of great concern [26][27][28]. Methods for monitoring the groundwater quality parameters and chemical element immigration are effective if there is a high-quality source and regularly updated material in the database [19,29]. The use of spatial analysis techniques significantly simplified the assessment of the hydrogeological and reclamation state [30]. Groundwater studies are used to analyze site suitability, manage available data and resources, assess soil vulnerability to salinization, model groundwater runoff, solute transport, and leaching [31,32]. The results of this analysis are necessary for creating systems to support spatial solutions using GIS [33,34]. Thus, in order to identify the salinity distribution area, mapping of surface and groundwater quality was required. The quality of groundwater depends on the evolution over time of various parameters, such as pH, salinity, Ca, Mg, or Cl [35]. The map is often compiled using the inverse distance weighting (IDW) method. Much of the research and development related to Earth sciences has been done using spatial technologies. GIS software is an effective tool in mapping the chemical composition of groundwater and is important for monitoring and managing hydrogeological and reclamation processes [19,29,33,36].
Under this working scope, the objective of this research is to describe and understand the chemical regime of groundwater in the Karatal irrigation massif.

Description of the Study Area
The study area is located on the territory of the Karatal district of the Almaty region of the Republic of Kazakhstan (77 • 53'17 east longitude and 45 • 23'26 north latitude, 78 • 00'40 east longitude and 45 • 06'46 north latitude). The irrigation system is located in the Karatal River valley, in the eastern marginal part of the South Balkhash depression. The massif stretches along the Karatal River for 20 km and is about 12 km wide ( Figure 1). The climate in the study area is sharply continental, characterized by severe winters, hot summers, and short springs and autumns. In 2019, the climatic conditions of the region during the growing season were characterized by a dry and hot summer. The maximum average monthly air temperature was observed in July (26.9 • C). The coldest month was December (−9.3 • C). The sum of active temperatures (the sum of the average daily air temperatures above +15 • C) for the period May-August was 2755 • C.
The total area of irrigated land was 11,845 ha. The terrain is calm, the slope does not exceed 0.001 m/m and decreases in a north-westerly direction. The absolute elevation of the terrain surface varies from 435 to 417 m. Soils are characterized as light chestnut gray soils, composed mainly of sandy loam, slightly saline. Soda-sulfate by anions, calcium-sodium by cations type of salinity chemistry is common in the eastern and northeastern parts of the massif. Here, the salt content in the soil horizon is higher than 0.6%. In the area of the PR-35 system (well 246), soils with chloride-sulfate magnesium-sodium with a salt content of 0.6% were identified. Well 179 is located on lands with sulfate sodium-calcium chemistry, where the salt content reaches 0.6%. In the northern and southeastern parts of the massif, soils are salted with medium sulfate by anions and a type of chemistry mixed by cations. On these soils, salt content varies from 0.3 to 0.5%.
On the territory of the massif, underlying rocks are alluvial deposits of the Quaternary age, represented by alternating layers of loam, sand, and sandy loam. The thickness of the cover deposits varies from 0.5 to 6 m and filtration coefficients from 0.02 to 0.3 m/day. The water-bearing rocks of the upper water-saturated stratum include sandy loam, loam, powdery and clay fine-grained sands, and lower-gravel-pebbles with sandy or sandy loam aggregates.
Land use has significantly changed in recent years due to market demands, subsidies, etc. The irrigated area is used entirely for agricultural crops. The water management conditions of the district, due to the physical wear and tear of irrigation channels and structures and the insufficient volume of repair and restoration work, are characterized as unsatisfactory. For these reasons, the destruction of the main structures of the structures, siltation, and overgrowth of weeds of the channels occur. The existing on-farm irrigation and collector-drainage channels are unattended and repair work is practically not carried out on them. All this leads to violations of the operating conditions of the irrigation system.
The research work described below was carried out in 2019 and included (i) observations of the chemical regime of groundwater in fields with rice crops and dry crops; (ii) determining the location of areas with the spread of brackish groundwater and (iii) determining the water quality in the Karatal River at the site of water intake and after mixing with collector and drainage waters.

Data Collection
Observations of the mineralization and chemical composition of groundwater were carried out by sampling water from 10 m deep observation wells following the requirements established by the monitoring rules of evaluation of irrigated lands in force at the Republic oz Kazakhstan [37]. The data collection process on the irrigated lands of the Karatal massif was done by the Zonal Hydrogeological and Reclamation Center of the Ministry of Agriculture of the Republic of Kazakhstan.
This process was done two times, in spring and autumn of 2019. Groundwater samples were taken from the existing 62 observation wells, irrigation wells located at the entrance to the massif in its southern part, and at the discharge point of the idle K-4 collector in the northern part ( Figure 2).

Data Analysis
Groundwater quality has been assessed by analyzing its degree of mineralization and chemical composition. Chemical composition was determined by the content of anions in the samples: HCO 3 , SO 4 , Cl, and cations-Ca, Mg, Na, K, which are expressed in mg/dm 3 or mg/eq. The sum of anions and cations determines the salinity. According to the degree of salinity, groundwater is classified into fresh (salt concentration up to 1 g/dm 3 ), slightly salty (1-3), strongly salty (3)(4)(5), salty (5-50), and brines (>50). The chemical composition and mineralization of water are expressed by Kurlov's equations (see Supplementary Materials). According to this formula, the chemical type of water is determined by the predominant anions and cations. The type of water is set at the beginning according to the content of anions, and then cations in descending order. If the content of components is less than 10%, they do not participate in determining the type of water [37].
There are two main groups in surface and underground waters: macro-and microcomponents. Macro components are determined by the predominant cations and anions. Their content in water determines the amount of dry residue, specific gravity, and the nature of the physical properties of water. For the general characteristics of the composition and properties of water, an abbreviated chemical analysis of water has been performed.
When conducting water laboratory tests, the requirements for the types of analyses performed and the applicable regulatory documents are observed. The types of analyses and current standards for determining the chemical composition of natural waters in the Republic of Kazakhstan are shown in Table 1. The value of the Sodium adsorption ratio (SAR) of groundwater to evaluate the possible salinization of soil in contact with these waters was calculated according to: The assessment of water quality by the SAR value is given in Table 2. Due to the large number of variables, the data analysis was based on a geostatistical analysis method. Results of the water analysis were processed and interpreted through cartographic methods. Spatial maps of pH, SAR, and mineralization, as well as maps of the concentration distributions of each ion in groundwater, were obtained using the IDW interpolation method implemented inside the ArcGIS software. This interpolation method is based on the assumption that the values that are closer to each other are more similar than those that are located further away. The measured values surrounding the unmeasured location were used for prediction. As mentioned above, mapping helped to understand field information but it was not enough to take into account all the influencing factors. Once data were displayed on a map a quantitative study was performed to establish correlations between the parameters leading to a better understanding of the chemical regime of groundwater in the study area. The zonal statistical tool was used for calculating areas with different groundwater mineralization and calculations were made according to the accepted gradation. Previously, it was necessary to apply a layer of farm boundaries on the constructed map, within which the spatial calculation of the area should be made. The Spatial Analyst tools were the main ones used in creating the maps. Therefore, spatially integrated data and analytical tools were used to assess the hydrogeological and reclamation status of irrigated land. Results obtained from the analysis were compared with the actual Kazakh normative data of monitoring irrigated lands.

Results and Discussion
The generated spatial and attribute databases are integrated to create maps of spatial variations of the main chemical parameters. Groundwater quality maps are shown in the figures below to visualize the spatial distribution of maximum ion concentrations.
The generation of groundwater quality maps helps to understand the current state of groundwater in the area and to assess the hydrochemical regime as well as to visualize the trend of salt accumulation processes [38].

Surface Water Quality
Results of laboratory analyses show that the water used for irrigation of agricultural crops is fresh and suitable for irrigation. During the growing season, the salinity of irrigation water varied up to 424 mg/dm 3 . The chemical composition of these waters is bicarbonate-sulfate sodium-calcium. The SAR index did not exceed 3.4, which indicates a low risk of soil salinization.
The salinity of the Karatal River water after mixing with the K-1 collector-drainage waters varied from 293 mg/dm 3 to 424 mg/dm 3 during the entire vegetation period. It should be noted that the collector-drainage waters currently entering the river do not lead to a deterioration in the quality of the river water.
The salinity of reservoir and drainage waters in the K-1 reservoir at the beginning of the growing season varied from 711 mg/dm 3 to 725 mg/dm 3 . In August, the mineralization decreased to 677 mg/dm 3 . After the water supply was stopped in September, the salt content in the collector-drainage waters decreased to 429 mg/dm 3 , and in October, when the intensive removal of salts from the soils of the aeration zone occurred with a decrease in the groundwater levels, the salinity of the collector-drainage waters reached 1370 mg/dm 3 . The chemical composition of the water is bicarbonate-sulfate calcium-sodium. The SAR index did not exceed 8.0 during the entire observation period. Figure 3 shows the dynamics of salinity of the irrigation and collector-drainage waters and the average salinity of groundwater over a five-year period. The calculation did not take into account samples of wells 86, 174, 175, 179 which are located in areas with saline soils and insufficiently provided with drainage in the north-eastern and central parts of the massif. Salt content observations in the Karatal River showed that the salinity of the river water at the beginning of the growing season is low. This effect is related to the amount of precipitation during this period. In the K-1 collector, there is a decrease in the mineralization of collector-drainage waters by the end of the growing season due to the flushing mode during irrigation.

Groundwater Quality and Chemistry
According to the hydrogeological conditions on the massif, the groundwater regime completely depends on agricultural and water management factors. Groundwater lies at a depth higher than 1.5 m. With the water supply to the fields, there is a widespread rise of the groundwater table, being the highest rate in May and the first days of June. After the end of the irrigation season and the discharge of water from rice checks, there is a decline in the groundwater levels, which continues throughout the inter-vegetation period. In the spring period, the groundwater levels are below 3 m. This indicates favorable conditions for the passage of redox processes in the soils of the aeration zone.
The qualitative characteristics of groundwater in the study area are shown in Figures 4 and 5.  The pH values for the growing season of 2019 ranged from 5.7 to 9.2. In most of the irrigated lands, alkaline groundwater is common, only in the eastern part, slightly acidic water is concentrated (see Figures 4A and 5A). At the beginning of the irrigation season, the salinity of groundwater varied from 237 to 1539 mg/dm 3 ( Figure 4B). Table 3 presents statistical data of the results of the analysis of groundwater samples (minimum, maximum, average, median, and standard deviation values). The representation of the main components for the pre-and post-vegetation periods is shown in Figures 6 and 7.
The chemical composition of water with minimal salinity is mainly bicarbonate or bicarbonate-sulfate sodium-calcium or sodium, with a salt concentration of more than 500 mg/dm 3 , the proportion of sulfates increases (see Figures 4B and 6C). Groundwater with salinity higher than 1000 mg/dm 3 has a predominantly bicarbonate sodium chemical composition.
The map of the spatial distribution of the main elements on the pre-( Figure 6) and postvegetation ( Figure 7) periods are very similar, which demonstrates the mineralization of groundwater in the massif (Figures 4 and 5). By the end of the growing season, the salinity of groundwater varied from 167 to 1407 mg/dm 3 ( Figure 5B), with the predominance of bicarbonate, sodium, and calcium ions in the chemical composition. Under rice fields, a decrease in the mineralization of groundwater is characteristic due to dilution and subsequent desalination with irrigation water [17,39].
In fields with rice crops (wells 4, 104, 25, 27, 52, 92, 79) groundwater salinity decreased from 323-1051 mg/dm 3 in April to 242-670 mg/dm 3 in August due to the flushing regime. This hydrochemical situation persists almost until the end of the irrigation season. After the termination of irrigation, there is a gradual increase in the mineralization of groundwater, due to the fact that in conditions of still high groundwater standing (from 1 to 2 m), evaporative processes manifest themselves. Their mineralization increases to 1000 mg/dm 3 , and the chemical composition is dominated by sulfates and bicarbonates from anions, and sodium comes first from cations, then magnesium.  The enrichment of groundwater with sodium cations is explained by the processes of ion exchange with soil solutions of rocks of the aeration zone. According to researchers [17,39], when solutions interact with soils, cation exchange most often occurs, as a result of which the content of sodium bicarbonate increases in groundwater. In the fields occupied by dry crops, both the accumulation of salts in the groundwater and their decrease is observed, depending on the degree of salinity of the soil, the drainage of the territory, the amount of water supply as well as the state of the collector and drainage network.
During the growing season, with a sufficiently close occurrence to the earth's surface, the processes of evaporative concentration increase, which leads to a further increase in their mineralization. At the end of the irrigation season, dissolved salts under the influence of climatic and biological factors partially crystallize and move to higher horizons, up to the surface of the earth [15,39,40]. For example, in well 225, the salinity of groundwater during the season increased from 314 to 699 mg/dm 3 , and in well 37, decreased from 1259 to 204 mg/dm 3 .
High salinization in wells 86, 174, 175, 179 is caused by the accumulation of salts in groundwaters due to the dissolution of salts of the soil complex. Here, the mineralization in April varied in the range of 5334-16,354 mg/dm 3 and 3849-13,632 mg/dm 3 in August (see Figures 4B and 5B). The chemical composition is sodium sulfate-chloride or sodiummagnesium chloride-sulfate. According to the survey data, as well as the comparison of the land salinity map, salinization processes are manifested here. According to the salinity patterns in this area, soils with soda and sulfate-soda types are dominated by anions and sodium-by cations. The weighted average salt content in the upper meter thickness reaches 1% or more.
The state of drainage systems largely determines the nature of the distribution of the salt in the soil profile. Long-term observations in areas occupied by rice crops show unsatisfactory drainage where salts during the irrigation period are washed out to a depth of 50-100 cm. These salts are not removed outside the massif but accumulate in the lower part of the soil profile.
At the end of the irrigation period, an intensive rise of salts in the surface horizons begins. The maximum salt content during this period is located at 0-30 cm depth. The unsatisfactory condition of the collector and drainage networks does not allow for the removal of drainage water. This leads to an increase in the level of mineralized groundwater, which is aggravated by high evaporation and, subsequently, to an intensive influx of toxic salts into the root layer. High mineralization rates in these wells may also be associated with their location in low terrain areas, where the water exchange process is difficult. In some locations, abnormal results may also be due to technical difficulties and clogged wells (well 86).
Under stabilized groundwater regime conditions, the time to restore the initial state directly depends on the degree of salinity of the soil, its filtration properties, and, to no lesser extent, on the provision of artificial drainage. The return to the original values usually occurs in February-March of the following year [15,39].
Assessment of water quality by the SAR value is mandatory when implementing hydrogeological and reclamation monitoring of irrigated land. The groundwater quality is assessed by determining the salinization risk degree defined by the SAR index. During the growing season, the SAR index in groundwater ranged from 0.5-10.2. In wells with high groundwater salinity, the SAR index varied over the season from 41.3 to 11.9 and from 39.9 to 14.1. These values indicate a high degree of soil salinization (see Figures 4C and 5C).
There is a close relationship between groundwater mineralization and SAR (correlation coefficient 0.71), which is explained by soda salinization.

Correlation Analysis of Chemical Parameters
In this study, a correlation analysis between chemical parameters was carried out to establish the relationships between ions in groundwater. Available data are shown in Table 4. A high positive correlation has been observed between Cl − and SO 4 2− (R > 0.91), Cl − and Mg 2+ (0.96), Na+K and Cl − (>0.97), Mg 2+ and SO 4 2− (>0.94), SO 4 2− and Na+K (>0.96), Mg 2+ and Na+K (0.97). Salinity shows a good correlation with Cl − , SO 4 2− , Mg +2 and Na+K (0.96-0.99). The N+K on the Cl bivariate plot (Figure 8a,a') is used to demonstrate the origin and amount of dangerous salt. Unusual high values are noted in the northeastern part of the massif. High values of Na+K and Cl have been observed as a result of the dissolution of soil salts during the growing season [41,42]. However, the existing irrigation and economic problems in this area do not allow them to be taken outside the massif. Therefore, their ratios change slightly from spring to autumn. The bivariate plot of Ca+Mg versus SO 4 +HCO 3 (Figure 8b,b') was used to determine the ion exchange processes, since the carbonate-calcium and the sulfate-calcium system is of particular importance in irrigation conditions. This is due to the fact that the processes of dissolution and precipitation of calcium and gypsum carbonates occur just within the limits of mineralization observed in pore waters with natural fluctuations in soil moisture and rocks of the aeration zone and in irrigation groundwater when their reserves change as a result of infiltration of irrigation waters, evaporation, and transpiration [30]. Figure 8 shows the predominance of Ca+Mg over SO 4 +HCO 3 and confirms the process. After infiltration of irrigation waters and dilution of the upper layer of groundwater, greater dissolution of Ca occurs (Figure 8b'), which is the main process responsible for the release of Ca +2 . Figure 9 shows the scatterplots of the distribution of anions and cations by mineralization. Seasonal patterns of changes in the chemical composition of groundwater are observed. At high mineralization ratios, Cl − and SO 4 2− predominate in anions while Na+K and Mg +2 predominate in cations. However, the values of the four samples differ in all graphs. Further analysis using ArcGIS when building maps was done to explain these anomalies. The main components analysis indicates the existence of mainly natural processes under the influence of irrigation. The HCO 3 − and CO 2− anions are formed due to the dissolution of carbonate-containing rocks. The content of the bicarbonate ion often prevails in freshwater (up to 1000 mg/dm 3 ).
The amount of sulfate ion increases in weakly and strongly saline groundwater, where it often dominates over other anions. The sulfate ion is the leading one in the formation of the chemical composition of groundwater with mineralization ranging from 1500 to 10,000 mg/dm 3 . The migration ability of this ion may be limited by a biochemical barrier and the use of chemical ameliorants (sulfate ion, copper sulfate, zinc sulfate). The chlorine ion content usually does not exceed 20 mg-eq%. It is most often manifested in the salinity of groundwater from 5000 to 10,000 mg/dm 3 . Chlorine is the most mobile ion and usually has no biochemical barriers due to the high solubility of its salts.
The content of sodium and potassium in the groundwater composition in most cases was determined jointly. The sodium ion is the main element in the formation of groundwater mineralization on the irrigation array. Most often, it acts as the first element for fresh (up to 1000 mg/dm 3 ) and slightly mineralized (1000-3000 mg/dm 3 ) waters [8]. The sodium ion has practically no solubility barrier due to its good ability to dissolve its salts [41]. It is also retained by plants and in animal organisms.
Calcium is also found in groundwater (15 to 40 mg-eq%). The leading role of this cation is manifested in two ranges of mineralization changes from 250 to 750 mg/dm 3 and from 1400 to 3400 mg/dm 3 . The first range is explained by the calcium carbonate solubility barrier. Its further inhibition within the mineralization range from 700 to 1400 mg/dm 3 is associated with its absorption by plants. The second range of mineralization, when the calcium content increases again, is obviously due to the cation exchange between soils and groundwater.
The concentration of magnesium in groundwater ranges from 59 to 83 mg-eq%. Its predominant role is manifested within the limits of mineralization from 200 to 900 mg/dm 3 . The content of this cation in water is controlled by a biochemical barrier, since it is well absorbed by plants, participates in photosynthesis, and affects phosphorus and protein metabolism [8].
The areal change in the distribution of fresh and brackish waters from spring to autumn confirms the established chemical regime. The areal distribution of irrigated land by the groundwater salinity degree is shown in Table 5. Results show that the irrigated area is mainly covered by fresh groundwater with a salinity of up to 1000 mg/dm 3 . Slightly saline and salty water is distributed over 10-13% of the area. The area of distribution of slightly saline waters decreased by 414 ha since April, desalinated by diluting them with fresh irrigation water. The maps of the concentration distribution for each chemical element, for all the tested points, facilitate the visual perception of the existing patterns. Based on the analysis of the actual material of hydrochemical observations shown in Figures 4-7, the existing relationships are noted. When irrigation begins salt ions are redistributed in the groundwater. The bicarbonate ions contents (Figures 6A and 7A) and magnesium ( Figures 6F and 7F) increase while sodium (Figures 6G and 7G) and sulfate ion ( Figures 6C and 7C) decrease.

Conclusions
On cultivated lands, irrigation reclamation has an impact on the chemistry of groundwater. As a result of such an impact, complex hydrochemical connections are formed between irrigation waters, soils of the aeration zone, and groundwater. The main factors that control this process are the chemistry of irrigation water, the natural and climatic conditions, and the drainage conditions.
Based on a set of multiple observations of the hydrochemical regime of groundwater in the Karatal massif, a seasonal pattern of changes in the chemical composition of groundwater has been distinguished. For the actual irrigation type of regime in the massif, the change in the ion-salt regime of groundwater directly depends on the ratio of the irrigated and non-irrigated area, the volume of supplied irrigation water, and the drainage capacity of the drainage systems. Depending on the type of crop, the following seasonal patterns are observed in the groundwater regime.
Groundwater mineralization has increased in a relatively small area along the K-4 collector on the PR-35 system and in the lower areas between the K-1 and K-4 reservoirs. The presence of mineralized waters along the collector is facilitated by the fact that the lands of the PR-33-35 system have not been washed in recent years and a significant excess of evaporation over soil wettability contributes to a greater accumulation of salts in waters and soils. The unsatisfactory condition of the drainage network does not compensate for the increased groundwater supply on the P-27 system between the K-1 and K-4 collector and the underground inflow associated with the influence of irrigation channels and groundwater of the surrounding irrigation lands located at higher elevations. This hydrodynamic influence is transmitted to the groundwater that saturates the sand column, from which the water is "squeezed" in the form of ascending currents into the groundwater of the cover loam, which is consumed by evaporation. As a result, salts accumulate in soils and groundwater and mineralization ratios of up to 500-600 mg/dm 3 have been observed in the immediate vicinity of the channels. Groundwater is chemically similar to irrigation water (i.e., mainly bicarbonate-sulfate sodium-calcium).
Data obtained from soil samples and irrigation and drainage waters showed that the formation of the hydrochemical regime of groundwater on irrigated lands is influenced by the presence of saline soil-forming rocks, the groundwater levels depth, the volume and quality of irrigation water, and the degree of provision of artificial and natural drainage, as well as by evaporation processes.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/w14030285/s1, Table S1: The available data for each one of the 62 sampling sites.