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.     

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.     

Long‐Term Trends in Streamflow and Precipitation in Northwest California and Southwest Oregon, 1953‐2012

Using nonparametric Mann‐Kendall tests, we assessed long‐term (1953‐2012) trends in streamflow and precipitation in Northern California and Southern Oregon at 26 sites regulated by dams and 41 “unregulated” sites. Few (9%) sites had significant decreasing trends in annual precipitation, but September precipitation declined at 70% of sites. Site characteristics such as runoff type (groundwater, snow, or rain) and dam regulation influenced streamflow trends. Decreasing streamflow trends outnumbered increasing trends for most months except at regulated sites for May‐September. Summer (July‐September) streamflow declined at many sites, including 73% of unregulated sites in September. Applying a LOESS regression model of antecedent precipitation vs. average monthly streamflow, we evaluated the underlying streamflow trend caused by factors other than precipitation. Decreasing trends in precipitation‐adjusted streamflow substantially outnumbered increasing trends for most months. As with streamflow, groundwater‐dominated sites had a greater percent of declining trends in precipitation‐adjusted streamflow than other runoff types. The most pristine surface‐runoff‐dominated watersheds within the study area showed no decreases in precipitation‐adjusted streamflow during the summer months. These results suggest that streamflow decreases at other sites were likely due to more increased human withdrawals and vegetation changes than to climate factors other than precipitation quantity.     

Statistical Models to Predict and Assess Spatial and Temporal Low‐Flow Variability in New England Rivers and Streams

In the northern hemisphere, summer low flows are a key attribute defining both quantity and quality of aquatic habitat. I developed one set of models for New England streams/rivers predicting July/August median flows averaged across 1985–2015 as a function of weather, slope, % imperviousness, watershed storage, glacial geology, and soils. These models performed better than most United States Geological Survey models for summer flows developed at a statewide scale. I developed a second set of models predicting interannual differences in summer flows as a function of differences in air temperature, precipitation, the North Atlantic Oscillation (NAO) index, and lagged NAO. Use of difference equations eliminated the need for transformations and accounted for serial autocorrelations at lag 1. The models were used in sequence to estimate time series for monthly low flows and for two derived flow metrics (tenth percentile [Q10] and minimum 3‐in‐5 year average flows). The first metric is commonly used in assessing risk to low‐flow conditions over time, while the second has been correlated with increased probability of localized extinctions for brook trout. The flow metrics showed increasing trends across most of New England for 1985–2015. However, application of summer flow models with average and extreme climate projections to the Taunton River, Massachusetts, a sensitive watershed undergoing rapid development, projected that low‐flow metrics will decrease over the next 50 years.     

Seasonal and Flood‐Induced Variations in Groundwater–Surface Water Exchange in a Northern Coldwater Fishery

Groundwater upwelling is important to coldwater fisheries survival. This study used stable isotopes to identify upwelling zones within a watershed, then combined isotope analyses with reach‐scale monitoring to measure surface water–groundwater exchange over time. Research focused on Amity Creek, Minnesota, a basin that exemplifies conditions limiting coldwater species survival along Lake Superior's North Shore where shallow bedrock limits groundwater capacity, lowering baseflows and increasing temperatures. Groundwater‐fed reaches were identified through synoptic isotope sampling, with results highlighting the importance of isolated shallow surficial aquifers (glacially derived sands and gravels) for providing cold baseflow waters. In an alluvial reach, monitoring well results show groundwater was stored in two reservoirs: one that reacts quickly to changes in stream levels, and one that remained isotopically isolated under most flow conditions, but which helps sustain summer baseflows for weeks to months. A 500‐year flood demonstrated the capacity of high‐flow events to alter surface water–groundwater connectivity. The previously isolated reservoir was exchanged or mixed during the flood pulse, while incision lowered the water table for years. The results here provide insight for streams that lack substantial groundwater inputs yet maintain coldwater species at risk in a warming climate and an approach for managers seeking to protect cold baseflow sources.     

Increased Frequency of Low‐Magnitude Floods in New England1

Armstrong, William H., Mathias J. Collins, and Noah P. Snyder, 2012. Increased Frequency of Low‐Magnitude Floods in New England. Journal of the American Water Resources Association (JAWRA) 48(2): 306‐320. DOI: 10.1111/j.1752‐1688.2011.00613.x Abstract: Recent studies document increasing precipitation and streamflow in the northeastern United States throughout the 20th and early 21st Centuries. Annual peak discharges have increased over this period on many New England rivers with dominantly natural streamflow – especially for smaller, more frequent floods. To better investigate high‐frequency floods (<5‐year recurrence interval), we analyze the partial duration flood series for 23 New England rivers selected for minimal human impact. The study rivers have continuous records through 2006 and an average period of record of 71 years. Twenty‐two of the 23 rivers show increasing trends in peaks over threshold per water year (POT/WY) – a direct measure of flood frequency – using the Mann‐Kendall trend test. Ten of these trends had p < 0.1. Seventeen rivers show positive trends in flood magnitude, six of which had p < 0.1. We also investigate a potential hydroclimatic shift in the region around 1970. Twenty‐two of the 23 rivers show increased POT/WY in the post‐1970 period when comparing pre‐ and post‐1970 records using the Wilcoxon rank‐sum test. More than half of these increases have p < 0.1, indicating a shift in flow regime toward more frequent flooding. Region wide, we found a median increase of one flood per year for the post‐1970 period. Because frequent floods are important channel‐forming flows, these results have implications for channel and floodplain morphology, aquatic habitat, and restoration.     

Streamflow Changes in the South Atlantic,United States During the Mid‐ and Late 20th Century1

Abstract: Repeated severe droughts over the last decade in the South Atlantic have raised concern that streamflow may be systematically decreasing, possibly due to climate variability. We examined the monthly and annual trends of streamflow, precipitation, and temperature in the South Atlantic for the time periods: 1934‐2005, 1934‐1969, and 1970‐2005. Streamflow and climate (temperature and precipitation) trends transitioned ca. 1970. From 1934 to 1969, streamflow and precipitation increased in southern regions and decreased in northern regions; temperature decreased throughout the South Atlantic. From 1970 to 2005, streamflow decreased, precipitation decreased, and temperature increased throughout the South Atlantic. It is unclear whether these will be continuing trends or simply part of a long‐term climatic oscillation. Whether these streamflow trends have been driven by climatic or anthropogenic changes, water resources management faces challenging prospects to adapt to decadal‐scale persistently wet and dry hydrologic conditions.     

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000 Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.     

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.     

Analyzing the Effects of Groundwater Pumping on an Urban Stream‐Aquifer System

When the cone of influence of a pumping well reaches a nearby river, the resulting hydraulic gradient can induce enhanced seepage of streamflow into the aquifer. The rate of seepage is often modeled using analytical solutions that are simple to apply but may not reproduce field data due to mathematical assumptions not being met in the field. Furthermore, the appropriateness of such models has not been investigated in detail due to difficulty in measuring streamflow loss in the field. In this study, a field experiment was conducted on a reach of the South Platte River near Denver, Colorado to estimate pumping‐induced streamflow loss. A network of stream gauges, monitoring wells, and in situ measurements was used to observe streamflow rates, groundwater levels, and temperature to assess if pumping wells have a significant impact on streamflow, and to compare observed streamflow depletion against analytical solutions. Data collected suggest that pumping wells have a noticeable impact on streamflow. The analytical solutions proved accurate if streamflow was low and constant but performed poorly if streamflow was high and variable. Therefore, for this reach, the use of analytical solutions to predict streamflow may only be appropriate under low‐flow, constant‐flow conditions. Methods and results can be used to guide other streamflow depletion studies and to inform cases of pumping‐induced streamflow depletion, particularly in regard to water rights.     

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.     

STREAMFLOW OF THE OCCOQUAN RIVER IN VIRGINIA AS RECONSTRUCTED FROM TREE-RING SERIES1

ABSTRACT Tree-ring indices representing seven sites were used to reconstruct monthly summer streamflow in the Occoquan River basin of northern Virginia from 1841 to 1975. Attempts were made to reconstruct flow for each of the months, April through August. Reconstructions for June, July, and August were judged most reliable. Major mid-summer flow minima persisting for more than one year were reconstructed as having occurred in the early 1870's, the early 1930's, and the mid-1960's. Aside from these major dry periods, a greater frequency of extreme low flow during individual years is indicated for the entire record than for the most recent 50 years.     

Streamflow Depletion Modeling: Methods for an Adaptable and Conjunctive Water Management Decision Support Tool

Groundwater pumping depletes streamflow, which can have significant ecological impacts depending on the magnitude of depletion relative to environmental flow needs. Streamflow depletion estimates from groundwater pumping have been quantified using both analytical and numerical methods, but are not routinely compared to environmental flow needs or used in practical water management tools. Decision support tools that incorporate groundwater dynamics are becoming increasingly necessary for water managers as groundwater regulations become more important in environmental policy, particularly concerning the preservation of environmental flow needs. We develop and apply methods for a web‐based decision support tool for conjunctive groundwater and surface water management, demonstrated using a case study watershed in British Columbia, Canada. Open‐source data are analyzed with a combination of spatial algorithms and previously developed analytical models, such that the tool can be applied to other regions. Streamflow depletion estimates are calculated in four regions in the largely undeveloped Bulkley Valley, British Columbia. Our transparent methodology has geographic and data input flexibility which is a significant improvement on currently existing water management tool methods.     

Streamflow Response to Climate as Influenced by Geology and Elevation1

Mayer, Timothy D. and Seth W. Naman, 2011. Streamflow Response to Climate as Influenced by Geology and Elevation. Journal of the American Water Resources Association (JAWRA) 47(4):724‐738. DOI: 10.1111/j.1752‐1688.2011.00537.x Abstract: This study examines the regional streamflow response in 25 predominately unregulated basins to warmer winter temperatures and snowpack reductions over the last half century in the Klamath Basin of California and Oregon. Geologic controls of streamflow in the region result in two general stream types: surface‐dominated and groundwater‐dominated basins. Surface‐dominated basins were further differentiated into rain basins and snowmelt basins on the basis of elevation and timing of winter runoff. Streamflow characteristics and response to climate vary with stream type, as discussed in the study. Warmer winter temperatures and snowpack reductions have caused significantly earlier runoff peaks in both snowmelt and groundwater basins in the region. In the groundwater basins, the streamflow response to changes in snowpack is smoothed and delayed and the effects are extended longer in the summer. Our results indicate that absolute decreases in July‐September base flows are significantly greater, by an order of magnitude, in groundwater basins compared to surface‐dominated basins. The declines are important because groundwater basins sustain Upper Klamath Lake inflows and mainstem river flows during the typically dry summers of the area. Upper Klamath Lake April‐September net inflows have decreased an estimated 16% or 84 thousand acre‐feet (103.6 Mm³) since 1961, with the summer months showing proportionately more decline. These changes will exacerbate water supply problems for agriculture and natural resources in the region.     

Conjunctive Water Use in Confined Basalt Aquifers: An Evaluation Using Geochemistry,a Numerical Model,and Historical Water Level

As withdrawals from deep compartmentalized aquifers increasingly exceed recharge throughout the western United States, conjunctive water use management alternatives have become an applied research priority. This study highlights both details and limitations of the role of irrigation canal seepage as groundwater recharge, revealing the regional limitations of canal seepage as a dependable source of recharge in overdrawn aquifers. A suite of geochemical indicators were used together with a numerical model to evaluate current and future management scenarios focused on recharge derived from seepage from a region‐wide irrigation canal system. Twenty‐five years of static groundwater level data were used to relate spatial trends determined using geochemistry and groundwater modeling with “on‐the‐ground” management practices, which vary based on acreage, crop, and irrigation scheduling. Increasing groundwater age determined using isotope analysis, and declines in potentiometric heads, each correlate with increasing distance from the canal reaches. Predictive modeling indicates that if pumping is gradually reduced, as has been suggested by management agencies, that recharge from canal seepage will be negligible by 2035 due to regional groundwater through‐flow and the pattern of potentiometric head recovery. Unfortunately, historic hydrographs suggest that under current groundwater development conditions most wells are not sustainable, irrespective of proximity to the canal.     

SUMMER LOW FLOWS IN NEW ENGLAND DURING THE 20TH CENTURY1

ABSTRACT: High springtime river flows came earlier by one to two weeks in large parts of northern New England during the 20th Century. In this study it was hypothesized that late spring/early summer recessional flows and late summer/early fall low flows could also be occurring earlier. This could result in a longer period of low flow recession and a decrease in the magnitude of low flows. To test this hypothesis, variations over time in the magnitude and timing of low flows were analyzed. To help understand the relation between low flows and climatic variables in New England, low flows were correlated with air temperatures and precipitation. Analysis of data from 23 rural, unregulated rivers across New England indicated little evidence of consistent changes in the timing or magnitude of late summer/early fall low flows during the 20th Century. The interannual variability in the timing and magnitude of the low flows in northern New England was explained much more by the interannual variability in precipitation than by the interannual variability of air temperatures. The highest correlation between the magnitude of the low flows and air temperatures was with May through November temperatures (r =?0.37, p= 0.0017), while the highest correlation with precipitation was with July through August precipitation (r = 0.67, p > 0.0001).     

Controls on Temperature in Salmonid‐Bearing Headwater Streams in Two Common Hydrogeologic Settings,Kenai Peninsula,Alaska

Headwater streams are the most numerous in terms of both number and length in the conterminous United States and play important roles as spawning and rearing grounds for numerous species of anadromous fish. Stream temperature is a controlling variable for many physical, chemical, and biological processes and plays a critical role in the overall health and integrity of a stream. We investigated the controls on stream temperature in salmon‐bearing headwater streams in two common hydrogeologic settings on the Kenai Peninsula, Alaska: (1) drainage‐ways, which are low‐gradient streams that flow through broad valleys; and (2) discharge‐slopes, which are high gradient streams that flow through narrow valleys. We hypothesize local geomorphology strongly influences surface‐water and groundwater interactions, which control streamflow at the network scale and stream temperatures at the reach scale. The results of this study showed significant differences in stream temperatures between the two hydrogeologic settings. Observed stream temperatures were higher in drainage‐way sites than in discharge‐slope sites, and showed strong correlations as a continuous function with the calculated topographic metric flow‐weighted slope. Additionally, modeling results indicated the potential for groundwater discharge to moderate stream temperature is not equal between the two hydrogeologic settings, with groundwater having a greater moderating effect on stream temperature at the drainage‐way sites.     

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.     

Estimating Contributions of Nitrate and Herbicides From Groundwater to Headwater Streams,Northern Atlantic Coastal Plain,United States1

Abstract: Groundwater transport often complicates understanding of surface‐water contamination. We estimated the regional flux of nitrate and selected herbicides from groundwater to nontidal headwater streams of the Atlantic Coastal Plain (New Jersey through North Carolina) based on late‐winter or spring base‐flow samples from 174 streams. Sampled streams were selected randomly, and flux estimates are based on resulting population estimates rather than on empirical models, which have been used previously for similar estimates. Base‐flow flux in the estimated 8,834 headwater streams of the study area are an estimated 21,200 kg/day of nitrate (as N) and 5.83, 0.565, and 20.7 kg/day of alachlor, atrazine, and metolachlor (and selected degradates), respectively. Base‐flow flux of alachlor and metolachlor is <3% of the total base‐flow flux of those compounds plus degradates. Base‐flow flux of nitrate and herbicides as a percentage of applications is typically highest in well‐drained areas and lowest in areas with abundant poor drainage and anoxic conditions. In Coastal Plain watersheds of Albemarle and Pamlico Sounds, <2% of applied nitrogen reaches headwater streams as base flow. On the Delmarva Peninsula part of the Chesapeake Bay watershed, however, more than 10% of such applications are transported through groundwater to streams, and base‐flow nitrate flux represents 70% of total nitrogen flux in headwater streams.     

Are Investments in Groundwater Irrigation Profitable? A Case of Rice Farms from South India1

Varghese, Shalet Korattukudy, Jeroen Buysse, Aymen Frija, Stijn Speelman, and Guido Van Huylenbroeck, 2012. Are Investments in Groundwater Irrigation Profitable? A Case of Rice Farms from South India. Journal of the American Water Resources Association (JAWRA) 1‐15. DOI: 10.1111/j.1752‐1688.2012.00690.x Abstract: This article examines the profitability of cultivating double rice under bore well irrigation, given the cumulative interference of and reduced life span of wells, and thus increases the cost of groundwater extraction and use. The overexploitation of groundwater is a common stock problem and the cultivation of water intensive crops, such as rice, further exacerbates the overdraft of groundwater. Under these circumstances, we quantify the marginal benefit of irrigation investments in rice farming by estimating the probability of having a double rice crop as a function of the investment made in wells. Using this information, we explore profit maximization behavior of farms with a mathematical programming model to derive individual economic optima of irrigation costs. The results demonstrate that the ongoing overexploitation of groundwater, and its use to cultivate an economically inefficient crop, such as rice, has resulted in low profitability at farm level. A sensitivity analysis found that even when the investment in irrigation wells is reduced by 70%, small farms are still not economically efficient, thereby confirming the Tragedy of the Commons. Raising awareness amongst farmers with regard to the economics of irrigation would facilitate the participatory implementation of control mechanisms to regulate groundwater extraction.