Spatial variability of groundwater depth and quality parameters in the National Capital Territory of Delhi

Land use dynamics and its influences on groundwater depth levels in South region of National Capital Territory (NCT) of Delhi, India

The present study envisions the influences of land use dynamics on the spatial trend of groundwater depth levels over a period of two decades in the south region of NCT Delhi, India. The findings have inferred that among five major LULC categories, built-up has shown a sprawling trend (+2.17 km2/year) from the north, northeastern, and central portions to the confined patches observed in the south and southwest fringes of South Delhi from 2001 to 2021. Likewise, vegetation class has also witnessed significant increment (+1.91 km2/year) to the peripheral boundary, i.e., southern, southwest/eastern fringes, and central portions under recent initiatives of city forests, plantation drives, and urban green spaces programs. On the contrary, a negative rate of change has been observed in fallow land (-2.78 km2/year), agriculture (-1.22 km2/year), and water bodies (-0.07 km2/year). LULC transition matrix has also showed prominent conversion of fallow land into vegetation and built-up class, and change of vegetation into built-up, and fallow land category. Subsequently, corresponding LULC maps have been superimposed with long-term trends of groundwater depth levels (as spatial contours). For entire South Delhi region, depth to water levels ranged from 2.02 to 66.45 meters below ground level (mbgl) where shallower depths followed a negative trend and remain persistent throughout the time period in north and northeastern fringes. Higher fluctuations in groundwater depletion with positive trends directly get influenced with allied land use transitions such as a steady increase in built-up area and steeper depth levels (> 40 mbgl) as observed in central, southwestern, and southern parts. Moreover, buffer peripheries in the proximity of groundwater monitoring stations viz., Hauz Khas, Pushp Vihar, Jamali, Gadaipur, and Bhatti Kalan have observed deeper groundwater levels allied to built-up expansion. Thus, groundwater depletion trends can be ascribed to the incessant conversion of recharging areas into impervious zones along with uneven distribution of groundwater usage and supply. Conversely, expanding vegetative land has also shown improved groundwater depth levels. Therefore, land use influences must be managed in the long run for ensuring sustainable management of groundwater resources.

Spatio-temporal geostatistical modelling of sulphate concentration in the area of the Reocín Mine (Spain) as an indicator of water quality

Water stored in open-pit lakes can be a water resource when the mine is closed. This study aimed to develop a reliable model to evaluate the water quality, based on the sulphate concentration, in the Reocín Mine area (Spain) by using geostatistical algorithms. To this end, water samples were taken from the beginning of the flooding period in November 2004 until August 2020. The model showed that the sulphate concentration was highest between February 2009 and February 2012 and decreased as the flooding process progressed. The area with the highest concentration (2000 mg L^-1) was the central part of the study area, where the mine is located, while in the northeast and southwest, the values from the beginning of the flooding period were lower, below 500 mg L^-1. In the last obtained model, the values decreased considerably to 1300 mg L^-1 in the central area and below 250 mg L^-1 in the northeast and southwest areas. The modelling conducted to assess the water quality in the area of influence of the mine determined that the flooding process has little influence on the water in the rivers and streams in the area, since the sulphate concentration measured in the adjacent rivers and streams was less than 250 mg L^-1, indicating that anomalous concentrations were only found in the open-pit area. It was shown that geostatistical algorithms are useful tools that can be used to model the intensity and extension of water pollutants in space over time.

Capturing spatial variability of factors affecting the water allocation plans-a geo-informatics approach for large irrigation schemes

The livelihoods of poor people living in rural areas of Indus Basin Irrigation System (IBIS) of Pakistan depend largely on irrigated agriculture. Water duties in IBIS are mainly calculated based on crop-specific evapotranspiration. Recent studies show that ignoring the spatial variability of factors affecting the crop water requirements can affect the crop production. The objective of the current study is thus to identify the factors which can affect the water duties in IBIS, map these factors by GIS, and then develop the irrigation response units (IRUs), an area representing the unique combinations of factors affecting the gross irrigation requirements (GIR). The Lower Chenab Canal (LCC) irrigation scheme, the largest irrigation scheme of the IBIS, is selected as a case. Groundwater quality, groundwater levels, soil salinity, soil texture, and crop types are identified as the main factors for IRUs. GIS along with gamma design software GS + was used to delineate the IRUs in the large irrigation scheme. This resulted in a total of 84 IRUs in the large irrigation scheme based on similar biophysical factors. This study provided the empathy of suitable tactics to increase water management and productivity in LCC. It will be conceivable to investigate a whole irrigation canal command in parts (considering the field-level variations) and to give definite tactics for management.

Geospatial Assessment of Groundwater Quality with the Distinctive Portrayal of Heavy Metals in the United Arab Emirates

Groundwater is a valuable resource, and its quality is critical to human survival. Optimal farming and urbanization degraded groundwater reserves. This research investigates and reports the spatial variability of selected heavy metals developed in the Liwa area of the United Arab Emirates. Forty water samples were collected from existing wells and analyzed for different elements. Principal components analysis was applied to a subgroup of the data set in terms of their usefulness for determining the variability of groundwater quality variables. Geographic information systems were used to produce contour maps to analyze the distribution of heavy metals. Ordinary kriging was used with Circular, Spherical, Tetraspherical, Pentaspherical-Bessel, K-Bessel, Hole effect, and Stable models for better representation. The water quality index was constructed using heavy metal concentrations and other variables. This yielded a value of 900 beyond the limit stated by WHO and US EPA. Nugget analysis showed that Cd (0), K (7.38%), and SO4 (1.81%) variables exhibited strong spatial dependence. Al (27%), Ba (40.87%), Cr (63%), Cu (34%), EC (27%), HCO3 (56%), NO3(36%), Pb (64%), and TDS (53%) represented moderate spatial dependence. As (76%), Mn (79%), Ni (100%), pH (100%), Temp (93%), and Zn (100%) exhibited weak spatial dependence.

Hydrogeochemical Assessment of Groundwater Quality of Mokopane Area, Limpopo, South Africa Using Statistical Approach

Despite being a finite resource, both the quality and quantity of groundwater are under tremendous pressure due to rapid global changes, viz. population growth, land-use/land-cover changes (LULC), and climate change. The 6th Sustainable Development Goal (SDG) aims to “Ensure availability and sustainable management of water and sanitation for all”. One of the most significant dimensions of the SDG agenda is the emphasis on data and governance. However, the lack of good governance coupled with good observed data cannot ensure the achievement of SDG6. Therefore, this study strives to evaluate water quality status and hydrochemical processes governing it in the data-scarce Mokopane area of South Africa. Groundwater is the main source of fresh water supply for domestic usage, intensive agriculture, and mining activities in Mokopane. In this study, hydrogeochemical analysis of groundwater samples was employed to calculate the water quality index (WQI) and evaluate factors governing water quality evolution in the study area. Statistical and spatial analysis techniques were carried out to divide sampling sites into clusters and delineate principal factors responsible for determining water quality of the sampled groundwater. Results suggest that most of the physico-chemical parameters are within permissible limits for drinking water set by the World Health Organization (WHO), except for high fluoride in some samples. Na-HCO3 is the most abundant water type followed by Mg-HCO3, which indicates dominance of Na+, Mg2+, and HCO3±. Rock-water interaction is the prime factor responsible for fluoride enrichment in water. The alkaline nature of groundwater favors the release of exchangeable F− from minerals like muscovite. The WQI suggests that 80% of water samples fall into the good and excellent categories. Poor management of untreated domestic sewage and agricultural runoff is a main factor for the bad/very bad categories of water samples. As the area lacks any credible scientific/government work to report water quality and its management aspects, the findings of this study will definitely help both scientific communities and policy makers to do what is needed for sustainable water resource management in a timely manner.

Rice yield estimation based on forecasting the future condition of groundwater salinity in the Caspian coastal strip of Guilan Province, Iran

Irrigation water salinity is one of the factors that reduces agricultural production. Guilan Province is one of the most important rice-producing areas in Iran where groundwater is used for irrigation. The temporal and spatial variations of groundwater salinity were studied in the coastal strip covering 4285 km² of the province using data from 73 wells, as well as its estimated effect on the rice yield. Data on mean electrical conductivity (EC) for each 6-month period of 12 consecutive years, from the second half of 2002 until the end of 2014, were analyzed and resulted in 25 mean ECs. EC maps and maps of the probability of higher salinity areas were obtained by using ordinary kriging (OK) and indicator kriging (IK) in ArcGIS 9.3 software, respectively. Thereby, areas belonging to different salinity classes were outlined and places with higher salinity reducing the rice yield were identified. In addition, the Mann-Kendall test and Sen's slope were used to project future changes. The results indicated that due to the salinity of groundwater in the coastal strip area, the minimum and the maximum rice yields were 80% and 100%, respectively. Using the IK method, higher probability of groundwater salinity reducing the yield was found from the central parts toward the east. The Mann-Kendal test result showed significant temporal trends of the size of areas below the 100% yield (EC < 1 dS/m) and 90-100% yield (1 < EC < 1.34 dS/m) thresholds. The equations given by Sen's slope estimator indicated that the groundwater salinity will not be a limiting factor for achieving 100% rice yields from the year of 2021 onward in all of the Guilan coastal area. The trend of increasing precipitation in the area may be an important cause.

Improvement of water table interpolation and groundwater storage volume using fuzzy computations

The water table is an important piece of data for hydrogeological studies, particularly as input data to groundwater simulation models. Since the accuracy of groundwater simulation models significantly depends on input data, this study highlights the application of fuzzy kriging to improve the accuracy of water table interpolation. The results of the fuzzy kriging approach are compared with common methods in water table interpolation like ordinary kriging, inverse distance weighting (IDW), and Thiessen polygon methods to justify the suitability of the fuzzy kriging. The Gilan and Zanjan plains, located in the northwest of Iran, are used as case study areas. The Gilan Plain is characterized by a dense and regular piezometric network and gentle hydraulic gradient. The longitudinal plain of Zanjan has a sparse and irregular piezometric network and steep hydraulic gradient. Since these plains have different piezometric network configurations, the sensitivity of the interpolation methods to the monitoring point configuration is analyzed. The cross-validation method is employed to validate the accuracy of interpolation methods in water table interpolation. In control points, the average of root-mean-square errors associated with groundwater water table values estimated using fuzzy kriging, ordinary kriging, IDW, and Thiessen polygon methods are obtained to be respectively 1.36, 1.93, 3.49, and 9.10 in the Gilan Plain and 13.60, 22.86, 32.30, and 59.81 in the Zanjan Plain. The results indicate that the fuzzy kriging technique has greater precision in comparison with other methods, especially under the conditions of the sparse piezometric network and steep hydraulic gradient. The results also demonstrate that the used methods generally have higher accuracy in the Gilan Plain with a regular piezometric network than in the Zanjan Plain. Furthermore, Thiessen polygon, IDW, and ordinary kriging methods overestimated water table in comparison with the fuzzy kriging method in our cases. This overestimation may cause large error values in subsequent calculations such as water budget and aquifer storage which play a major role in the appropriate management of water resources.

Seawater intrusion vulnerability in the coastal aquifers of southern India-an appraisal of the GALDIT model, parameters' sensitivity, and hydrochemical indicators

An appraisal of seawater intrusion into the coastal aquifers is one of the major issues for groundwater resource management. The GALDIT model applies to the analysis of multiple parameters using systematic GIS techniques for mapping and assessment of seawater intrusion vulnerability. It demarcates the mapping of potential vulnerability that shows a higher vulnerability to seawater intrusion in various parts of the coast and the estimated vulnerability index value of 7.50 and 9.64. An area of 33.0 km² spread in the low-lying coastal area comprising estuaries, salt marshes, and saltpans shows the high vulnerability condition with an estimated vulnerability value of 6.42-7.50. An area of 73.20 km² spread over coastal and alluvial plains experiences moderate vulnerability (temporal salinity in the groundwater sources) with an estimated vulnerability index value of 5.46-6.42. Aquifers underlying coastal uplands (hard rock formations) and some parts of accretionary beaches (2.05 km²) are relatively protected fresh groundwater sources, wherein the estimated vulnerability index is 4.55-5.46. The vulnerability mapping of the GALDIT model using hydrochemical analysis of primary groundwater parameters such as TDS, Cl^-, HCO₃, and Cl^-/HCO₃ ratio is validated. Higher concentration of TDS (2637-4162 mg/l) and Cl^- (1268-2347 mg/l) is taken for the areas falling under higher vulnerability to seawater intrusion, especially in the placer mining sites and coastal areas facing erosion. Similarly, the groundwater sources of the low-lying areas including estuaries, salt marshes, saltpans, and backwater were noted to have higher values of Cl^-/HCO₃ with a rationality of 9.87-12.18. Hydrological facies shows the highest concentration of NaCl in the groundwater sources within the proximity of eroded beaches, saltwater bodies, and sand mining areas. A hydrochemical facies evolution (HFE) diagram represents the hydrochemical facies of groundwater elements that shows an intrusion of seawater into the coastal aquifers underlying the very high vulnerable zones. Higher bicarbonate concentration (233-318 mg/l) is noticed in the upland areas and some parts of dunes and accreted beaches, sandy coasts, and uplands. Vulnerability analysis reveals that those areas near saltwater bodies and eroding coasts are prone to lateral and vertical diffusion of saltwater. The geodatabase developed through such modeling studies can help in planning and developing activities for sustainable groundwater resource management in coastal areas.

Groundwater fluoride contamination, probable release, and containment mechanisms: a review on Indian context

Fluoride contamination in the groundwater has got great attention in last few decades due to their toxicity, persistent capacity and accumulation in human bodies. There are several sources of fluoride in the environment and different pathways to enter in the drinking water resources, which is responsible for potential effect on human health. Presence of high concentration of fluoride ion in groundwater is a major issue and it makes the water unsuitable for drinking purpose. Availability of fluoride in groundwater indicates various geochemical processes and subsurface contamination of a particular area. Fluoride-bearing aquifers, geological factors, rate of weathering, ion-exchange reaction, residence time and leaching of subsurface contaminants are major responsible factors for availability of fluoride in groundwater. In India, several studies have reported that the groundwater of several states are contaminated with high fluoride. The undesirable level of fluoride in groundwater is one of the most natural groundwater quality problem, which affects large portion of arid and semiarid regions of India. Rajasthan, Andhra Pradesh, Telangana, Tamil Nadu, Gujarat, and West Bengal are the relatively high-fluoride-contaminated states in India. Chronic ingestion of high doses of fluoride-rich water leads to fluorosis on human and animal. Over 66 million Indian populations are at risk due to excess fluoride-contaminated water. Therefore, groundwater contamination subject to undesirable level of fluoride needs urgent attention to understand the role of geochemistry, hydrogeology and climatic factors along with anthropogenic inputs in fluoride pollution.

Spatial and temporal variability of fluoride concentrations in groundwater resources of Larestan and Gerash regions in Iran from 2003 to 2010

There is discrepancy about intervals of fluoride monitoring in groundwater resources by Iranian authorities. Spatial and temporal variability of fluoride in groundwater resources of Larestan and Gerash regions in Iran were analyzed from 2003 to 2010 using a geospatial information system and the Mann-Kendall trend test. The mean concentrations of fluoride for the 8-year period in the eight cities and 31 villages were 1.6 and 2.0 mg/l, respectively; the maximum values were 2.4 and 3.8 mg/l, respectively. Spatial, temporal, and spatiotemporal variability of fluoride in overall groundwater resources were relatively constant over the years. However, results of the Mann-Kendall trend test revealed a monotonic trend in the time series of one city and 11 villages for the 8-year period. Specifically, one city and three villages showed positive significant Kendall's Tau values, suggesting an upward trend in fluoride concentrations over the 8-year period. In contrast, seven villages displayed negative significant Kendall's Tau values, arguing for a downward trend in fluoride concentrations over the years. From 2003 to 2010, approximately 52 % of the Larestan and Gerash areas have had fluoride concentrations above the maximum permissible Iranian drinking water standard fluoride level (1.4 mg/l), and about 116,000 people were exposed to such excess amounts. Therefore, our study supports for a close monitoring of fluoride concentrations from health authorities in monthly intervals, especially in villages and cities that showed positive trend in fluoride concentrations. Moreover, we recommend simultaneous implementation of cost-effective protective measures or interventions until a standard fluoride level is achieved.

Does irrigation with reclaimed water significantly pollute shallow aquifer with nitrate and salinity? An assay in a perurban area in North Tunisia

In Tunisia, reclaimed water is increasingly used for irrigation in order to mitigate water shortage. However, few studies have addressed the effect of such practice on the environment. Thus, we attempted in this paper to assess the impact of irrigation with reclaimed water on the nitrate content and salinity in the Nabeul shallow aquifer on the basis of satellite images and data from 53 sampled wells. Ordinary and indicator kriging were used to map the spatial variability of these groundwater chemical parameters and to locate the areas where water is suitable for drinking and irrigation. The results of this study have shown that reclaimed water is not an influential factor on groundwater contamination by nitrate and salinity. Cropping density is the main factor contributing to nitrate groundwater pollution, whereas salinity pollution is affected by a conjunction of factors such as seawater interaction and lithology. The predictive maps show that nitrate content in the groundwater ranges from 9.2 to 206 mg/L while the electric conductivity ranges from 2.2 to 8.5 dS/m. The high-nitrate concentration areas underlie sites with high annual crop density, whereas salinity decreases gradually moving away from the coastline. The probability maps reveal that almost the entire study area is unsuitable for drinking with regard to nitrate and salinity levels. Appropriate measures, such as the elaboration of codes of good agricultural practices and action programs, should be undertaken in order to prevent and/or remediate the contamination of the Nabeul shallow aquifer.

Impacts of urbanization on groundwater hydrodynamics and hydrochemistry of the Toluca Valley aquifer (Mexico)

The Toluca Valley is located on the high plains of Mexico, where there are significant industrial zones and large populations. Water needs are almost exclusively met by groundwater, which has brought about intense exploitation of the aquifer and indication of some contamination. The present study investigates the effect of urbanization, related to industrialization of the region, on groundwater in the central portion of the Toluca Valley aquifer--a zone with high population density and where the largest industrial park is located. A general decline in the groundwater level has been found over the years, at a rate of as much as 2.5 m/year. The appearance of a large drawdown cone was identified, indicating changes in the direction of groundwater flow. Also identified was the presence of several ground fissures, the location of which coincided with the drawdown cone. In hydrochemical terms, the water type is sodium-magnesium bicarbonate and this characteristic has not changed over time, although it has been possible to detect the presence of larger quantities of sulfates (up to 117 mg/L) and nitrates (up to 47 mg/L) in recent years, likely associated with contamination from industrial and urban wastewater. Factor analysis made it possible to identify ions that would characterize natural processes involving the acquisition of salts (HCO3 (-), Na(+), Mg(2+), and Si), as well as anthropic activities (SO4 (2-), NO3 (-), Cl(-), Ca(2+), and K(+)).

Categorical Indicator Kriging for assessing the risk of groundwater nitrate pollution: the case of Vega de Granada aquifer (SE Spain)

Groundwater nitrate pollution associated with agricultural activity is an important environmental problem in the management of this natural resource, as acknowledged by the European Water Framework Directive. Therefore, specific measures aimed to control the risk of water pollution by nitrates must be implemented to minimise its impact on the environment and potential risk to human health. The spatial probability distribution of nitrate contents exceeding a threshold or limit value, established within the quality standard, will be helpful to managers and decision-makers. A methodology based on non-parametric and non-linear methods of Indicator Kriging was used in the elaboration of a nitrate pollution categorical map for the aquifer of Vega de Granada (SE Spain). The map has been obtained from the local estimation of the probability that a nitrate content in an unsampled location belongs to one of the three categories established by the European Water Framework Directive: CL. 1 good quality [Min - 37.5 ppm], CL. 2 intermediate quality [37.5-50 ppm] and CL. 3 poor quality [50 ppm - Max]. The obtained results show that the areas exceeding nitrate concentrations of 50 ppm, poor quality waters, occupy more than 50% of the aquifer area. A great proportion of the area's municipalities are located in these poor quality water areas. The intermediate quality and good quality areas correspond to 21% and 28%, respectively, but with the highest population density. These results are coherent with the experimental data, which show an average nitrate concentration value of 72 ppm, significantly higher than the quality standard limit of 50 ppm. Consequently, the results suggest the importance of planning actions in order to control and monitor aquifer nitrate pollution.

Association of sociodemographic and nutritional factors with risk of neural tube defects in the North Indian population: a case-control study

OBJECTIVE

To assess the role of sociodemographic and nutritional factors in the incidence of births affected by neural tube defects (NTD) in the North Indian population.

DESIGN

Case-control study.

SETTING

Government hospitals of Delhi, India.

SUBJECTS

Subjects comprised 284 mothers of NTD children (cases) and 568 mothers of healthy children (controls).

RESULTS

Significant differences were found between case and control mothers with respect to maternal age (P = 0·005), type of drinking water (P = 0·03) and consumption of milk (P = 0·01). Univariate and multivariate analysis suggested an association of unpasteurized milk use, low consumption of vegetables, low consumption of fruits and vegetarian dietary habits with NTD births. Further, variation in the risk factors for upper and lower NTD types was also observed, pointing towards phenotypic heterogeneity in the aetiology.

CONCLUSIONS

The results of the present study suggest an increased risk of NTD infants in mothers with low consumption of vegetables, fruits and milk and having vegetarian dietary habits. So, in order to reduce these devastating birth defects in future offspring, better nutritional care should be provided to mothers by suggesting dietary modifications and augmenting additional micronutrient supplementation during the periconceptional period.

Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins

In sparsely monitored basins, accurate mapping of the spatial variability of groundwater level requires the interpolation of scattered data. This paper presents a comparison of deterministic interpolation methods, i.e. inverse distance weight (IDW) and minimum curvature (MC), with stochastic methods, i.e. ordinary kriging (OK), universal kriging (UK) and kriging with Delaunay triangulation (DK). The study area is the Mires Basin of Mesara Valley in Crete (Greece). This sparsely sampled basin has limited groundwater resources which are vital for the island's economy; spatial variations of the groundwater level are important for developing management and monitoring strategies. We evaluate the performance of the interpolation methods with respect to different statistical measures. The Spartan variogram family is applied for the first time to hydrological data and is shown to be optimal with respect to stochastic interpolation of this dataset. The three stochastic methods (OK, DK and UK) perform overall better than the deterministic counterparts (IDW and MC). DK, which is herein for the first time applied to hydrological data, yields the most accurate cross-validation estimate for the lowest value in the dataset. OK and UK lead to smooth isolevel contours, whilst DK and IDW generate more edges. The stochastic methods deliver estimates of prediction uncertainty which becomes highest near the southeastern border of the basin.

Spatial assessment of soil salinity in the Harran Plain using multiple kriging techniques

The Harran Plain is located in the southeastern part of Turkey and has recently been developed for irrigation agriculture. It already faces soil salinity problems causing major yield losses. Management of the problem is hindered by the lack of information on the extent and geography of the salinization problem. A survey was carried out to delineate the spatial distribution of salt-affected areas by randomly selecting 140 locations that were sampled at two depths (0 to 30 and 30 to 60 cm) and analyzed for soil salinity variables: soil electrical conductivity (EC), soluble cations (Ca(2+,) Mg(2+), Na(+), and K(+)), soluble anions (SO (4) (2-) , Cl(-)), exchangeable Na(+) (me 100 g(-1)) and exchangeable sodium percentage. Terrain attributes (slope, topographical wetness index) were extracted from the digital elevation model of the study area. Variogram analyses after log transformation and ordinary kriging (OK) were applied to map spatial patterns of soil salinity variables. Multivariate geostatistical methods-regression kriging (RK) and kriging with external drift (KED)-were used using elevation and soil electrical conductivity data as covariates. Performances of the three estimation methods (OK, RK, and KED) were compared using independent validation samples randomly selected from the main dataset. Soils were categorized into salinity classes using disjunctive kriging (DK) and ArcGIS, and classification accuracy was tested using the kappa statistic. Results showed that soil salinity variables all have skewed distribution and are poorly correlated with terrain indices but have strong correlations among each other. Up to 65 % improvement was obtained in the estimations of soil salinity variables using hybrid methods over OK with the best estimations obtained with RK using EC(0-30) as covariate. DK-ArcGIS successfully classified soil samples into different salinity groups with overall accuracy of 75 % and kappa of 0.55 (p < 0.001).

Spatial variability of the shallow groundwater level and its chemistry characteristics in the low plain around the Bohai Sea, North China

To characterize the spatial distribution of groundwater level (GWL) and its chemistry characteristics in the low plain around the Bohai Sea, shallow groundwater depth of 130 wells were determined. Water soluble ions composition, total dissolved solid (TDS), electric conductivity (EC), total hardness (TH), total alkalinity (TA), and total salt content (TS) of 128 representative groundwater samples were also measured. Classical statistics, geostatistical method combined with GIS technique were then used to analyze the spatial variability and distribution of GWL and groundwater chemical properties. Results show that GWL, TDS, EC, TH, TA, and TS all presented a lognormal distribution and could be fitted by different semivariogram models (spherical, exponential, and Gaussian). Spatial structure of GWL, TDS, EC, TH, TA, and TS changed obviously. GWL decreased from west inland plain to the east coastal plain, however, TDS, EC, and TS increased from west to east, TH and TA were higher in the middle and coastal plain area. Groundwater chemical type in the coastal plain was SO (4) (2-) ·Cl(-)-Na(+) while chemical types in the inland plain were SO (4) (2-) ·Cl(-)-Ca(2+)·Mg(2+) and HCO (3) (-) -Ca(2+)·Mg(2+).

Field-scale spatial variation of saline-sodic soil and its relation with environmental factors in Western Songnen Plain of China

The objectives of this study were to investigate the degree of spatial variability and variance structure of salinization parameters using classical and geostatistical method in Songnen Plain of China, which is one of largest saline-sodic areas in the World, and to analyze the relationship between salinization parameters, including soil salinity content (SC), electrical conductivity (EC), sodium adsorption ratio (SAR), and pH, and seven environmental factors by Pearson and stepwise regression analysis. The environmental factors were ground elevation, surface ponding time, surface ponding depth, and soil moistures at four layers (0-10 cm, 10-30 cm, 30-60 cm, and 60-100 cm). The results indicated that SC, EC, and SAR showed great variations, whereas pH exhibited low variations. Four salinization parameters showed strongly spatial autocorrelation resulting from the compound impact of structural factors. The empirical semivariograms in the four parameters could be simulated by spherical and exponential models. The spatial distributions of SC, EC, SAR and pH showed similar patterns, with the coexistence of high salinity and sodicity in the areas with high ground elevation. By Pearson analysis, the soil salinization parameters showed a significant positive relationship with ground elevation, but a negative correlation with surface ponding time, surface ponding depth, and soil moistures. Both correlation and stepwise regression analysis showed that ground elevation is the most important environmental factor for spatial variation of soil sanilization. The results from this research can provide some useful information for explaining mechanism of salinization process and utilization of saline-sodic soils in the Western Songnen Plain.