Estimating Annual Groundwater Evapotranspiration from Phreatophytes in the Great Basin Using Landsat and Flux Tower Measurements

Escalating concerns about water supplies in the Great Basin have prompted numerous water budget studies focused on groundwater recharge and discharge. For many hydrographic areas (HAs) in the Great Basin, most of the recharge is discharged by bare soil evaporation and evapotranspiration (ET) from phreatophyte vegetation. Estimating recharge from precipitation in a given HA is difficult and often has significant uncertainty, therefore it is often quantified by estimating the natural discharge. As such, remote sensing applications for spatially distributing flux tower estimates of ET and groundwater ET (ET_g) across phreatophyte areas are becoming more common. We build on previous studies and develop a transferable empirical relationship with uncertainty bounds between flux tower estimates of ET and a remotely sensed vegetation index, Enhanced Vegetation Index (EVI). Energy balance‐corrected ET measured from 40 flux tower site‐year combinations in the Great Basin was statistically correlated with EVI derived from Landsat imagery (r² = 0.97). Application of the relationship to estimate mean‐annual ET_g from four HAs in western and eastern Nevada is highlighted and results are compared with previous estimates. Uncertainty bounds about the estimated mean ET_g allow investigators to evaluate if independent groundwater discharge estimates are “believable” and will ultimately assist local, state, and federal agencies to evaluate expert witness reports of ET_g, along with providing new first‐order estimates of ET_g.     

Evaluating Flow Diversion Impacts to Groundwater‐Dependent Riparian Vegetation with Flow Alteration and Groundwater Model Analysis

An approach for assessing the potential ecologic response of groundwater‐dependent riparian vegetation to flow alteration is developed, focusing on change to groundwater. Groundwater requirements for riparian vegetation are reviewed in conjunction with flow alteration statistics. Where flow alteration coincides with groundwater‐related vegetation sensitivities, scenarios are developed for groundwater simulation. Groundwater depths and recession rates in the riparian zone are simulated for baseline and altered stream hydrographs, with changes to river stage and width represented with a transient, flow‐dependent boundary condition. Potential flow diversion from the Upper Gila River in New Mexico is examined. Statistical flow alteration analysis, applying prospective diversions to a 76‐year record of daily flow, shows that flows in the winter‐spring months and within the high‐pulse to small flood range are subject to greatest potential change. Groundwater simulation scenarios are developed for these flow conditions in representative dry, near‐average, and wet years. Differences in groundwater elevations, generally less than 0.25 m during the flow alteration period, dissipate rapidly following cessation of diversion. Relating groundwater depth, recession rates and range of fluctuations to riparian vegetation needs, we find adverse ecological response is not expected from groundwater impacts for the flow alteration examined.     

Assessment of Groundwater Depletion and Implications for Management in the Denver Basin Aquifer System

The Denver Basin Aquifer System (DBAS) is a critical groundwater resource along the Colorado Front Range. Groundwater depletion has been documented over the past few decades due to the increased water use among users, presenting long‐term sustainability challenges. A spatiotemporal geostatistical analysis is used to estimate potentiometric surfaces and evaluate groundwater storage changes between 1990 and 2016 in each of the four DBAS aquifers. Several key depletion patterns and spatial water‐level changes emerge in this work. Hydraulic head changes are the largest in the west‐central side of the DBAS and have decreased in some areas by up to 180 m since 1990, while areas to the northwest show increases in hydraulic head by over 30.5 m. The Denver and Arapahoe aquifers show the largest groundwater storage losses, with the highest rates occurring in the 2000s. The results highlight uncertainty in the volumetric predictions under various storage coefficient calculations and emphasize the importance of representative aquifer characterization. The observed groundwater storage depletions are due to a combination of factors, which include population growth increasing the demand for water, variable precipitation, and drought influencing recharge, and increased groundwater pumping. The methods applied in this study are transferable to other groundwater systems and provide a framework that can help assess groundwater depletion and inform management decisions at other locations.     

Groundwater Nutrient Concentrations and Mass Loading Rates at Iowa Golf Courses

Golf courses are often considered by the public to be significant nitrogen (N) and phosphorus (P) nonpoint sources but only limited information exists on nutrient concentrations and loads in golf course groundwater. In this study, we measured N and P concentrations in groundwater and available surface water at six randomly selected Iowa golf courses to assess the loading risk posed by these facilities to groundwater and local rivers. At each course, three shallow monitoring wells were installed, one each on representative tee, fairway, and rough locations. Wells and available surface water were sampled on eight occasions during 2015 and 2016. NO₃‐N concentrations were not detected above 1 mg/L at three of the six courses monitored in this study and the overall mean NO₃‐N concentration in Iowa golf courses was 2.2 mg/L. The mass of NO₃‐N recharged to groundwater averaged 3.3 kg/ha at the six courses, which represents approximately one‐tenth of the NO₃‐N load exported by the watershed that contains the course and represented approximately 0.1 to 8% of the fertilizer N applied. Groundwater orthophosphorus concentrations averaged 0.13 mg/L and were similar to those measured in a variety of settings across Iowa. Study results should prove useful in evaluating nutrient contributions from golf courses in Midwestern states where nutrient reduction strategies are being pursued.     

Forecasting the Probability of Future Groundwater Levels Declining Below Specified Low Thresholds in the Conterminous U.S.

We present a logistic regression approach for forecasting the probability of future groundwater levels declining or maintaining below specific groundwater‐level thresholds. We tested our approach on 102 groundwater wells in different climatic regions and aquifers of the United States that are part of the U.S. Geological Survey Groundwater Climate Response Network. We evaluated the importance of current groundwater levels, precipitation, streamflow, seasonal variability, Palmer Drought Severity Index, and atmosphere/ocean indices for developing the logistic regression equations. Several diagnostics of model fit were used to evaluate the regression equations, including testing of autocorrelation of residuals, goodness‐of‐fit metrics, and bootstrap validation testing. The probabilistic predictions were most successful at wells with high persistence (low month‐to‐month variability) in their groundwater records and at wells where the groundwater level remained below the defined low threshold for sustained periods (generally three months or longer). The model fit was weakest at wells with strong seasonal variability in levels and with shorter duration low‐threshold events. We identified challenges in deriving probabilistic‐forecasting models and possible approaches for addressing those challenges.     

Estimating Current and Future Groundwater Resources of the Maldives

The water resources of the atolls of the Republic of Maldives are under continual threat from climatic and anthropogenic stresses, including land surface pollution, increasing population, drought, and sea‐level rise (SLR). These threats are particularly acute for groundwater resources due to the small land surface area and low elevation of each island. In this study, the groundwater resources, in terms of freshwater lens thickness, total volume of fresh groundwater, and safe yield are estimated for the 52 most populous islands of the Maldives for current conditions and for the year 2030, with the latter accounting for projected SLR and associated shoreline recession. An algebraic model, designed in previous studies to estimate the lens thickness of atoll islands, is expanded in this study to also estimate volume of groundwater. Results indicate that average current lens thickness, groundwater volume, and per capita safe yield are approximately 4.6 m, 1,300 million liters, and 300 l/day, and that these values will decrease by approximately 10, 11, and 34%, respectively, by the year 2030. Based on results, it is demonstrated that groundwater, in terms of quantity, is a viable source of water for the islands of the Maldives both now and in coming decades, particularly for islands with large surface area and low population. Study results can provide water resource managers and government officials with valuable data for consideration in water security measures.     

Impact of Urban Growth and High Residential Irrigation on Streamflow and Groundwater Levels in a Peri‐Urban Semiarid Catchment

The impact of urbanization on groundwater is not simple to understand, as it depends on a variety of factors such as climate, hydrogeology, water management practices, and infrastructure. In semiarid landscapes, the urbanization processes can involve high water consumptions and irrigation increases, which in turn may contribute to groundwater recharge. We assessed the hydrological impacts of urbanization and irrigation rates in an Andean peri‐urban catchment located in Chile, in a semiarid climate. For this purpose, we built and validated a coupled surface–groundwater model that allows the verification of a strong stream–aquifer interaction in areas with shallow groundwater, higher than some sewers and portions of the stream. Moreover, we also identified a significant local recharge associated with pipe leaks and inefficient urban irrigation. From the evaluation of different future scenarios, we found a sustainable water conservation scenario will decrease the current groundwater levels, while the median flow reduces from 408 to 389 L/s, and the low flow (Q_95%) from 43 to 22L/s. Overall, our results show the relevance of integrating the modeling of surface and subsurface water resources at different spatial and temporal scales, when assessing the effect of urban development and the suitability of urban water practices.     

Wind‐Powered Desalination: An Estimate of Saline Groundwater in the United States1

James Androwski, Abraham Springer, Thomas Acker, and Mark Manone, 2011. Wind‐Powered Desalination: An Estimate of Saline Groundwater in the United States. Journal of the American Water Resources Association (JAWRA) 47(1):93‐102. DOI: 10.1111/j.1752‐1688.2010.00493.x Abstract: Increasing scarcity of freshwater resources in many regions of the world is leading water resource managers to consider desalination as a potential alternative to traditional freshwater supplies. Desalination technologies are energy intensive and expensive to implement making desalination using renewable energy resources a potentially attractive option. Unfortunately, saline groundwater resources are not well characterized for many regions hindering consideration of such technologies. In this assessment, we estimate the saline groundwater resources of the principal aquifers of the United States using a geographic information system and correlate these resources to wind resources potentially sufficient to supply the energy demand of desalination equipment. We estimate that 3.1 × 10¹⁴ m³ saline groundwater, total volume, are contained in 28 of the country’s principal aquifers known to contain saline groundwater. Of this volume, 1.4 × 10¹⁴ m³ saline groundwater are co‐located with wind resources sufficient for electrical generation to desalinate groundwater.     

Hydrogeochemical and Isotopic Investigations on Groundwater Salinization in the Datong Basin,Northern China1

Abstract: Analyses of major elements, environmental isotope ratios (δ¹⁸O, δ²H), and PHREEQC inverse modeling investigations were conducted to understand the processes controlling the salinization of groundwater within the Datong Basin. The hydrochemical results showed that groundwater with high total dissolved solid (TDS) concentrations was dominated by sodium bicarbonate (Na‐HCO₃), sodium chlorite (Na‐Cl), and sodium sulfate (Na‐SO₄) type waters, whereas low‐TDS groundwater from near mountain areas was dominated by calcium bicarbonate (Ca‐HCO₃) and magnesium bicarbonate (Mg‐HCO₃) type waters. The characterization of the major components of groundwater and PHREEQC inverse modeling indicated that the aluminosilicate hydrolysis, cation exchange, and dissolution of evaporites (halite, mirabilite, and gypsum) governed the salinization of groundwater within the Datong Basin. The environmental isotope (δ¹⁸O, δ²H) and Cl^?/Br^? ratios revealed the impact of fast vertical recharge by irrigation returns and salt‐flushing water on the groundwater salinization. According to the analyses of major hydrochemical components and PHREEQC inverse modeling, evaporite dissolution associated with irrigation and salt‐flushing practice was probably the dominant controlling factor for the groundwater salinization, especially in the central part of the basin. Therefore, groundwater pumping for irrigation and salt‐flushing should be controlled to protect groundwater quality in this area.     

Climatic Controls on Watershed Reference Evapotranspiration Varied during 1961–2012 in Southern China

Reference evapotranspiration (ET_o) is an important hydrometeorological term widely used in understanding and projecting the hydrological effects of future climate and land use change. We conducted a case study in the Qinhuai River Basin that is dominated by a humid subtropical climate and mixed land uses in southern China. Long‐term (1961–2012) meteorological data were used to estimate ET_o by the FAO‐56 Penman–Monteith model. The individual contribution from each meteorological variable to the trend of ET_o was quantified. We found basin‐wide annual ET_o decreased significantly (p < 0.05) by 3.82 mm/yr during 1961–1987, due to decreased wind speed, solar radiation, vapor pressure deficit (VPD), and increased relative humidity (RH). However, due to the increased VPD and decreased RH, the ET_o increased significantly (p < 0.05) in spring, autumn, and annually at a rate of 2.55, 0.56, and 3.16 mm/yr during 1988–2012, respectively. The aerodynamic term was a dominant factor controlling ET_o variation in both two periods. We concluded the key climatic controls on ET_o have shifted as a result of global climate change during 1961–2012. The atmospheric demand, instead of air temperature alone, was a major control on ET_o. Models for accurately predicting ET_o and hydrological change under a changing climate must include VPD in the study region. The shifts of climatic control on the hydrological cycles should be considered in future water resource management in humid regions.     

Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas1

Mittelstet, Aaron R., Michael D. Smolen, Garey A. Fox, and Damian C. Adams, 2011. Comparison of Aquifer Sustainability Under Groundwater Administrations in Oklahoma and Texas. Journal of the American Water Resources Association (JAWRA) 1‐8. DOI: 10.1111/j.1752‐1688.2011.00524.x Abstract: We compared two approaches to administration of groundwater law on a hydrologic model of the North Canadian River, an alluvial aquifer in northwestern Oklahoma. Oklahoma limits pumping rates to retain 50% aquifer saturated thickness after 20 years of groundwater use. The Texas Panhandle Groundwater Conservation District’s (GCD) rules limit pumping to a rate that consumes no more than 50% of saturated thickness in 50 years, with reevaluation and readjustment of permits every 5 years. Using a hydrologic model (MODFLOW), we simulated river‐groundwater interaction and aquifer dynamics under increasing levels of “development” (i.e., increasing groundwater withdrawals). Oklahoma’s approach initially would limit groundwater extraction more than the GCD approach, but the GCD approach would be more protective in the long run. Under Oklahoma rules more than half of aquifer storage would be depleted when development reaches 65%. Reevaluation of permits under the Texas Panhandle GCD approach would severely limit pumping as the 50% level is approached. Both Oklahoma and Texas Panhandle GCD approaches would deplete alluvial base flow at approximately 10% development. Results suggest periodic review of permits could protect aquifer storage and river base flow. Modeling total aquifer storage is more sensitive to recharge rate and aquifer hydraulic conductivity than to specific yield, while river leakage is most sensitive to aquifer hydraulic conductivity followed by specific yield.     

Historical Groundwater Trends in Northern New England and Relations with Streamflow and Climatic Variables

Water‐level trends spanning 20, 30, 40, and 50 years were tested using month‐end groundwater levels in 26, 12, 10, and 3 wells in northern New England (Maine, New Hampshire, and Vermont), respectively. Groundwater levels for 77 wells were used in interannual correlations with meteorological and hydrologic variables related to groundwater. Trends in the contemporary groundwater record (20 and 30 years) indicate increases (rises) or no substantial change in groundwater levels in all months for most wells throughout northern New England. The highest percentage of increasing 20‐year trends was in February through March, May through August, and October through November. Forty‐year trend results were mixed, whereas 50‐year trends indicated increasing groundwater levels. Whereas most monthly groundwater levels correlate strongly with the previous month's level, monthly levels also correlate strongly with monthly streamflows in the same month; correlations of levels with monthly precipitation are less frequent and weaker than those with streamflow. Groundwater levels in May through August correlate strongly with annual (water year) streamflow. Correlations of groundwater levels with streamflow data and the relative richness of 50‐ to 100‐year historical streamflow data suggest useful proxies for quantifying historical groundwater levels in light of the relatively short and fragmented groundwater data records presently available.     

Causes of decreased reference evapotranspiration and pan evaporation in the Jinghe River catchment, northern China

This study focused on the changes of reference evapotranspiration (ET₀) and pan evaporation (ET_pan) to study the impacts of climate change on the hydrological cycle in the Jinghe River catchment. Based on the Penman–Monteith equation, the ET₀ was calculated. The temporal trend and spatial distribution of ET₀ and E_pan measured with a 20-cm pan were examined at the 14 stations during 1957–2005. The effects of meteorological factors on the variation of ET₀ were determined by analyzing the trends in themselves with comparison between original climate and detrended climate scenarios and then their sensitivity to ET₀. Both the ET₀ and E_pan showed remarkable decreasing trends from 1957 to 2005 and their decreasing rate was 40.9 and 17.7 mm per 10 years, respectively. Trend analysis of meteorological factors exhibited that the reduction in ET₀ and ET_pan was principally caused by both significant decreases in wind speed and sunshine hours. Furthermore, the decreasing trend of ET₀ was mainly dominated by the significant decrease in wind speed with high sensitivity, to a less extent, by the decrease in net radiation. Although relative humidity is one of the most sensitive variables, its effect on ET₀ was negligible because of its temporal constancy. The contribution of wind speed reduction to decreased ET₀ has increased from 50 to 76.1%, but net radiation, by contrast, decreased from 50 to 23.9%.     

Groundwater Level Trends and Drivers in Two Northern New England Glacial Aquifers

We evaluated long‐term trends and predictors of groundwater levels by month from two well‐studied northern New England forested headwater glacial aquifers: Sleepers River, Vermont, 44 wells, 1992‐2013; and Hubbard Brook, New Hampshire, 15 wells, 1979‐2004. Based on Kendall Tau tests with Sen slope determination, a surprising number of well‐month combinations had negative trends (decreasing water levels) over the respective periods. Sleepers River had slightly more positive than negative trends overall, but among the significant trends (p < 0.1), negative trends dominated 67 to 40. At Hubbard Brook, negative trends outnumbered positive trends by a nearly 2:1 margin and all seven of the significant trends were negative. The negative trends occurred despite generally increasing trends in monthly and annual precipitation. This counterintuitive pattern may be a result of increased precipitation intensity causing higher runoff at the expense of recharge, such that evapotranspiration demand draws down groundwater storage. We evaluated predictors of month‐end water levels by multiple regression of 18 variables related to climate, streamflow, snowpack, and prior month water level. Monthly flow and prior month water level were the two strongest predictors for most months at both sites. The predictive power and ready availability of streamflow data can be exploited as a proxy to extend limited groundwater level records over longer time periods.     

Further Development of a Specific Conductivity Approach to Measure Groundwater Discharge Area within Lakes

Groundwater exchanges with most lakes are rarely quantified because there are many technical challenges to quantification. We investigated a lakebed mapping approach to infer the relative areas of groundwater exchange in 12 prairie shallow lakes and five Laurentian mixed forest shallow lakes in Minnesota, USA in 2011. We used a relatively common approach (seepage meters) to provide baseline information on the magnitude and direction of flow at four locations in each lake. To expand from point measurements to the whole‐lake scale, we explored use of specific conductivity as a cheaper and more time efficient proxy for groundwater discharge to lakes. We validated the approach at near shore stations in each lake where seepage meter measurements and specific conductivity surveys overlapped. Specific conductivity surveys provided a similar assessment of groundwater discharge compared to seepage meters for 50% of the lake‐sampling period combinations. The lakebed mapping approach, when validated for a lake with a limited number of seepage meter (or alternative methods) measurements, offers the advantages of being more time and labor efficient over the use of a similar number of seepage meter monitoring locations; seepage meters (or piezometers, for example) are costlier in terms of equipment and labor, even for single‐lake studies. We show the combined approach could provide useful baselines for understanding and mapping groundwater exchange in shallow lakes.     

Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin,Nevada1

Carroll, Rosemary W.H., Greg Pohll, David McGraw, Chris Garner, Anna Knust, Doug Boyle, Tim Minor, Scott Bassett, and Karl Pohlmann, 2010. Mason Valley Groundwater Model: Linking Surface Water and Groundwater in the Walker River Basin, Nevada. Journal of the American Water Resources Association (JAWRA) 46(3):554-573. DOI: 10.1111/j.1752-1688.2010.00434.x Abstract: An integrated surface water and groundwater model of Mason Valley, Nevada is constructed to replicate the movement of water throughout the different components of the demand side of water resources in the Walker River system. The Mason Valley groundwater surface water model (MVGSM) couples the river/drain network with agricultural demand areas and the groundwater system using MODFLOW, MODFLOW’s streamflow routing package, as well as a surface water linking algorithm developed for the project. The MVGSM is capable of simulating complex feedback mechanisms between the groundwater and surface water system that is not dependent on linearity among the related variables. The spatial scale captures important hydrologic components while the monthly stress periods allow for seasonal evaluation. A simulation spanning an 11-year record shows the methodology is robust under diverse climatic conditions. The basin-wide modeling approach predicts a river system generally gaining during the summer irrigation period but losing during winter months and extended periods of drought. River losses to the groundwater system approach 25% of the river’s annual budget. Reducing diversions to hydrologic response units will increase river flows exiting the model domain, but also has the potential to increase losses from the river to groundwater storage.     

Influence of Riparian Seepage Zones on Nitrate Variability in Two Agricultural Headwater Streams

Riparian seeps have been recognized for their contributions to stream flow in headwater catchments, but there is limited data on how seeps affect stream water quality. The objective of this study was to examine the effect of seeps on the variability of stream NO₃‐N concentrations in FD36 and RS, two agricultural catchments in Pennsylvania. Stream samples were collected at 10‐m intervals over reaches of 550 (FD36) and 490 m (RS) on 21 occasions between April 2009 and January 2012. Semi‐variogram analysis was used to quantify longitudinal patterns in stream NO₃‐N concentration. Seep water was collected at 14 sites in FD36 and 7 in RS, but the number of flowing seeps depended on antecedent conditions. Seep NO₃‐N concentrations were variable (0.1‐29.5 mg/l) and were often greater downslope of cropped fields compared to other land uses. During base flow, longitudinal variability in stream NO₃‐N concentrations increased as the number of flowing seeps increased. The influence of seeps on the variability of stream NO₃‐N concentrations was less during storm flow compared to the variability of base flow NO₃‐N concentrations. However, 24 h after a storm in FD36, an increase in the number of flowing seeps and decreasing streamflow resulted in the greatest longitudinal variability in stream NO₃‐N concentrations recorded. Results indicate seeps are important areas of NO₃‐N delivery to streams where targeted adoption of mitigation measures may substantially improve stream water quality.     

Groundwater characterization in major industrial and residential locations of Lagos metropolis

This study investigated the quality of groundwater collected from two industrial and residential locations, each within the Lagos metropolis. Prescribed standard procedures of the American Public Health Association (APHA) were used to measure the physicochemical parameters of each of the groundwater samples, which include pH, electrical conductivity (EC), dissolved oxygen, total dissolved solids (TDS), biological oxygen demand, chemical oxygen demand; the anions chloride (Cl^?), nitrate (NO₃^?), sulfate (SO₄^?), and phosphate (PO₄^?); and heavy metals copper (Cu), zinc (Zn), lead (Pb), manganese (Mn), iron (Fe), cobalt (Co), cadmium (Cd), and chromium (Cr). Based on the laboratory analysis, the physicochemical parameters that were measured were within the permissible ranges specified by the World Health Organization and the Nigerian Standard for Drinking Water Quality Standards Organization of Nigeria (SON), except for pH, TDS, EC, Pb, Mn, and Fe for groundwater samples from the industrial locations and for pH, Pb, Mn, and Fe for residential locations. The elevated concentrations of TDS and EC reported for groundwater samples from industrial locations were attributed to the heavy discharge of effluents from nearby industrial treatment plants as well as the dissolution of ionic heavy metals from industrial activities involving the use of heavy machines. Statistical analysis using Pearson's correlation revealed the physicochemical parameters to be moderately and strongly correlated with one another at either p < .05 or < .01. In conclusion, groundwater samples from residential locations are more suitable for drinking than those from industrial locations.     

Accuracy Assessment of NOAA Gridded Daily Reference Evapotranspiration for the Texas High Plains

The National Oceanic and Atmospheric Administration (NOAA) provides daily reference evapotranspiration (ET_ref) maps for the contiguous United States using climatic data from North American Land Data Assimilation System (NLDAS). This data provides large‐scale spatial representation of ET_ref, which is essential for regional scale water resources management. Data used in the development of NOAA daily ET_ref maps are derived from observations over surfaces that are different from short (grass — ET_os) or tall (alfalfa — ET_rs) reference crops, often in nonagricultural settings, which carries an unknown discrepancy between assumed and actual conditions. In this study, NOAA daily ET_os and ET_rs maps were evaluated for accuracy, using observed data from the Texas High Plains Evapotranspiration (TXHPET) network. Daily ET_os, ET_rs and the climatic data (air temperature, wind speed, and solar radiation) used for calculating ET_ref were extracted from the NOAA maps for TXHPET locations and compared against ground measurements on reference grass surfaces. NOAA ET_ref maps generally overestimated the TXHPET observations (1.4 and 2.2 mm/day ET_os and ET_rs, respectively), which may be attributed to errors in the NLDAS modeled air temperature and wind speed, to which reference ET_ref is most sensitive. Therefore, a bias correction to NLDAS modeled air temperature and wind speed data, or adjustment to the resulting NOAA ET_ref, may be needed to improve the accuracy of NOAA ET_ref maps.     

Spatial Variability of Groundwater Depth and Quality Parameters in the National Capital Territory of Delhi

The groundwater quantity and quality scenario is of much concern in the National Capital Territory of Delhi, India, which necessitates an investigation to envisage the extent of spatial variability of groundwater depth and pollutant concentration levels in this region. Therefore, in this study, an effort was made to generate the spatial variability map of groundwater depth and quality parameters (viz. chloride, electrical conductivity, fluoride, magnesium, and nitrate). Ordinary kriging was used to analyze the spatial variability of groundwater depth and quality parameters, whereas indicator kriging was used to analyze groundwater quality parameters equal to or greater than the pollution threshold values. It was observed that the semivariogram parameters fitted well in the exponential model for water depth and in the spherical model for water quality parameters. The generated spatial variability maps indicated that in 43% of the study area, groundwater depth was within 20 m. The salinity level was higher than 2.5 dS m⁻¹ in 69% of the study area and the nitrate concentration exceeded 45 mg l⁻¹ in 36% of the area. The probability maps showed that about 24% of the area had the highest probability (0.8–1.0) of exceedence of the threshold electrical conductivity value and an area of 2% exhibited the highest probability of exceedence of the threshold value of nitrate concentration in the groundwater. The generated spatial variability and probability maps will assist water resource managers and policymakers in development of guidelines in judicious management of groundwater resources for agricultural and drinking purposes in the study area.