Calibrating a Basin‐Scale Groundwater Model to Remotely Sensed Estimates of Groundwater Evapotranspiration

Remotely sensed vegetation indices correspond to canopy vigor and cover and have been successfully used to estimate groundwater evapotranspiration (ET_g) over large spatial and temporal scales. However, these data do not provide information on depth to groundwater (dtgw) necessary for groundwater models (GWM) to calculate ET_g. An iterative approach is provided that calibrates GWM to ET_g derived from Landsat estimates of the Enhanced Vegetation Index (EVI). The approach is applied to different vegetation groups in Mason Valley, Nevada over an 11‐year time span. An uncertainty analysis is done to estimate the resulting mean and 90% confidence intervals in ET_g to dtgw relationships to quantify errors associated with plant physiologic complexity, species variability, and parameter smoothing to the 100 m GWM‐grid, temporal variability in soil moisture and nonuniqueness in the solution. Additionally, a first‐order second moment analysis shows ET_g to dtgw relationships are almost exclusively sensitive to estimated land surface, or maximum, ET_g despite relatively large uncertainty in extinction depths and hydraulic conductivity. The EVI method of estimating ET_g appears to bias ET_g during years with exceptionally wet spring/summer conditions. Excluding these years improves model performance significantly but highlights the need to develop a methodology that accounts not only on quantity but timing of annual precipitation on phreatophyte greenness.     

Annual Estimates of Recharge,Quick‐Flow Runoff,and Evapotranspiration for the Contiguous U.S. Using Empirical Regression Equations

This study presents new data‐driven, annual estimates of the division of precipitation into the recharge, quick‐flow runoff, and evapotranspiration (ET) water budget components for 2000‐2013 for the contiguous United States (CONUS). The algorithms used to produce these maps ensure water budget consistency over this broad spatial scale, with contributions from precipitation influx attributed to each component at 800 m resolution. The quick‐flow runoff estimates for the contribution to the rapidly varying portion of the hydrograph are produced using data from 1,434 gaged watersheds, and depend on precipitation, soil saturated hydraulic conductivity, and surficial geology type. Evapotranspiration estimates are produced from a regression using water balance data from 679 gaged watersheds and depend on land cover, temperature, and precipitation. The quick‐flow and ET estimates are combined to calculate recharge as the remainder of precipitation. The ET and recharge estimates are checked against independent field data, and the results show good agreement. Comparisons of recharge estimates with groundwater extraction data show that in 15% of the country, groundwater is being extracted at rates higher than the local recharge. These maps of the internally consistent water budget components of recharge, quick‐flow runoff, and ET, being derived from and tested against data, are expected to provide reliable first‐order estimates of these quantities across the CONUS, even where field measurements are sparse.     

Assessing Calibration Uncertainty and Automation for Estimating Evapotranspiration from Agricultural Areas Using METRIC

Agricultural irrigation accounts for a large fraction of the total water use in the western United States. The Mapping Evapotranspiration at high Resolution with Internalized Calibration (METRIC) remote sensing energy balance model is being used to estimate historical agricultural water use in western Nevada to evaluate basin‐wide water budgets. Each METRIC evapotranspiration (ET) estimate must be calibrated by a trained user, which requires some iterative time investment and results in variation in ET estimates between users. An automated calibration algorithm for the METRIC model was designed to generate ET estimates comparable to those from trained users by mimicking the manual calibration process. Automated calibration allows for rapid generation of METRIC ET estimates with minimal manual intervention, as well as uncertainty and sensitivity analysis of the model. The variation in ET estimates generated by the automated calibration algorithm was found to be similar to the variation in manual ET estimates. Results indicate that uncertainty was highest for fields with low ET levels and lowest for fields with high ET levels, with a seasonal mean uncertainty of approximately 5% for all fields. In addition, in a blind comparison, automated daily and seasonal ET estimates compared well with flux tower measurement ET data at multiple sites. Automated methods can generate first‐order ET estimates that are similar to time intensive manual efforts with less time investment.     

Assessing Satellite and Ground‐Based Potential Evapotranspiration for Hydrologic Applications in the Colorado River Basin

This study evaluates a remotely sensed and two ground‐based potential evapotranspiration (PET) products for hydrologic application in the Upper Colorado River Basin (UCRB). The remotely sensed Moderate Resolution Imaging Spectroradiometer product (MODIS‐PET) is a continuous, daily time series with 250 m resolution derived using the Priestley‐Taylor (P‐T) equation. The MODIS‐PET is evaluated against regional flux tower data as well as a synthetic pan product (Epan; 0.125°, daily) derived from the North American Land Data Assimilation System (NLDAS) and a Hargreaves PET derived from DAYMET variables (DAYMET‐PET; 1 km, daily). Compared to point‐scale PET computed using regional flux tower data, the MODIS‐PET had lower errors, with RMSE values ranging from 2.24 to 2.85 mm/day. Epan RMSE values ranged from 3.70 to 3.76 mm/day and DAYMET‐PET RMSE values ranged from 3.55 to 4.58 mm/day. Further investigation showed biases in temperature and radiation data contribute to uncertainty in the MODIS‐PET values, while bias in NLDAS temperature, downward shortwave (SW↓), and downward longwave (LW↓) propagate in the Epan estimates. Larger discrepancies between methods were observed in the warmer, drier regions of the UCRB, however, the MODIS‐PET was more responsive to landcover transitions and better captured basin heterogeneity. Results indicate the satellite‐based MODIS product can serve as a viable option for obtaining spatial PET values across the UCRB.     

Sensitivity and Uncertainty of Ground‐Water Discharge Estimates for Semiarid Shrublands1

Abstract: We proposed a step‐by‐step approach to quantify the sensitivity of ground‐water discharge by evapotranspiration (ET) to three categories of independent input variables. To illustrate the approach, we adopt a basic ground‐water discharge estimation model, in which the volume of ground water lost to ET was computed as the product of the ground‐water discharge rate and the associated area. The ground‐water discharge rate was assumed to equal the ET rate minus local precipitation. The objective of this study is to outline a step‐by‐step procedure to quantify the contributions from individual independent variable uncertainties to the uncertainty of total ground‐water discharge estimates; the independent variables include ET rates of individual ET units, areas associated with the ET units, and precipitation in each subbasin. The specific goal is to guide future characterization efforts by better targeting data collection for those variables most responsible for uncertainty in ground‐water discharge estimates. The influential independent variables to be included in the sensitivity analysis are first selected based on the physical characteristics and model structure. Both regression coefficients and standardized regression coefficients for the selected independent variables are calculated using the results from sampling‐based Monte Carlo simulations. Results illustrate that, while as many as 630 independent variables potentially contribute to the calculation of the total annual ground‐water discharge for the case study area, a selection of seven independent variables could be used to develop an accurate regression model, accounting for more than 96% of the total variance in ground‐water discharge. Results indicate that the variability of ET rate for moderately dense desert shrubland contributes to about 75% of the variance in the total ground‐water discharge estimates. These results point to a need to better quantify ET rates for moderately dense shrubland to reduce overall uncertainty in estimates of ground‐water discharge. While the approach proposed here uses a basic ground‐water discharge model taken from an earlier study, the procedure of quantifying uncertainty and sensitivity can be generalized to handle other types of environmental models involving large numbers of independent variables.     

A REEVALUATION OF THE GROUND WATER BUDGET FOR LAS VEGAS VALLEY,NEVADA, WITH EMPHASIS ON GROUND WATER DISCHARGE1

ABSTRACT: An essential component to the ground water budget for the Las Vegas Valley (LVV) in southern Nevada is discharge from the ground water system. Discharge for the LW has been based on estimates made more than 50 years ago of 35,524,224 m3 per year as evapotranspiration (ET) and 0 m³ per year as subsurface outflow. Newly published values for recharge based on a more robust data set (70,308,360 m³) indicate a large imbalance associated with the earlier discharge estimates, providing the basis for the reevaluation conducted in this study. ET estimates in this study, as opposed to previous studies, were assigned a range in values that included an approach that assigned higher weight to the unique soil, plant, water, and climatic conditions that existed in predevelopment (1905) LW. The earlier discharge estimates also assumed that the basin was hydrologically closed; however, based on our evaluation, a range in yearly discharge by subsurface outflow from 1,480,176 m³ to 19,735,680 m³ could be assigned. Likewise, a range in yearly ET from 20,475,768 m³ to 78,819,372 m³ could be assigned. Based on newly published recharge values, closure can only occur if higher values are assigned to both the subsurface outflow and/or ET components of ground water discharge. We cannot provide a complete water balance closure with our ground water discharge estimate of 64,140,960 m³. However our reevaluation gives support to the higher recharge estimates and provides the rationale for future studies to be conducted based on a more rigorous scientific assessment.     

Benchmarking Optical/Thermal Satellite Imagery for Estimating Evapotranspiration and Soil Moisture in Decision Support Tools

Generally, one expects evapotranspiration (ET) maps derived from optical/thermal Landsat and MODIS satellite imagery to improve decision support tools and lead to superior decisions regarding water resources management. However, there is lack of supportive evidence to accept or reject this expectation. We “benchmark” three existing hydrologic decision support tools with the following benchmarks: annual ET for the ET Toolbox developed by the United States Bureau of Reclamation, predicted rainfall‐runoff hydrographs for the Gridded Surface/Subsurface Hydrologic Analysis model developed by the U.S. Army Corps of Engineers, and the average annual groundwater recharge for the Distributed Parameter Watershed Model used by Daniel B. Stephens & Associates. The conclusion of this benchmark study is that the use of NASA/USGS optical/thermal satellite imagery can considerably improve hydrologic decision support tools compared to their traditional implementations. The benefits of improved decision making, resulting from more accurate results of hydrologic support systems using optical/thermal satellite imagery, should substantially exceed the costs for acquiring such imagery and implementing the remote sensing algorithms. In fact, the value of reduced error in estimating average annual groundwater recharge in the San Gabriel Mountains, California alone, in terms of value of water, may be as large as $1 billion, more than sufficient to pay for one new Landsat satellite.     

Evapotranspiration in the Meso‐Scale Forested Watersheds in Minjiang Valley,West China1

Abstract: Information on evapotranspiration (ET) can help us understand water balance, particularly in forested watersheds. Previous studies in China show that ET was relatively low (30‐40% of total precipitation) in the Minjiang Valley located in the upper reach of the Yangtze River Basin. However, this conclusion was derived from research on small‐scale watersheds (<100 km²). The objective of this paper was to present ET information on meso‐scale watersheds in the Minjiang Valley. Four meso‐scale watersheds (1,700‐5,600 km²) located in the Minjiang Valley were used to estimate ET using the water balance approach. We first generated forest vegetation variables (coniferous forest percentage, forest cover percentage, and derived forest vegetation index) using remote sensing data. Landsat 5 TM satellite images, acquired on June 26, 1994, were selected for the vegetation classification. Actual annual ET was calculated based on 11‐year estimated precipitation and measured streamflow data (1992‐2002). We also calculated potential ET (PET) using an improved Thornthwaite model for all four watersheds for the period of 1992‐1998. PET can provide additional information about potential capacity of water flux to atmosphere in the region. Seasonal (dry and rainy) PET and ET for all studied watersheds were also estimated for comparison purposes as the water balance approach, at shorter than annual scales, would likely provide inaccurate estimates of ET. The dominant vegetations in the Minjiang Valley were grasslands, conifer forests, and shrub‐lands. Our results confirmed that both ET and PET for three studied meso‐scale watersheds in the Minjiang Valley is relatively low (39.5‐43.8 and 28.2‐47.7% for ET and PET, respectively), with an exception of ET in the Yuzixi watershed being 71.1%. This result is generally consistent with previous research at small watershed scales. Furthermore, the low ET across various scales in the Minjiang Valley may be related to the unique deeply cut valley environment.     

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.     

GROUND-WATER RECHARGE THROUGH ACTWE SAND DUNES IN NORTHWESTERN NEVADA1

ABSTRACT: Most water-resouree investigations in semiarid basins of the Great Basin in western North America conclude that ground-water recharge from direct precipitation on the valley floor is negligible. However, many of these basins contain large areas covered by unvegetated, active sand dunes that may act as conduits for ground-water recharge. The potential for this previously undocumented recharge was investigated in an area covered by sand dunes in Desert Valley, northwestern Nevada, using a deep percolation model. The model uses daily measurements of precipitation and temperature th determine energy and moisture balances, from which estimates of long-term mean annual recharge are made. For the study area, the model calculated a mean annual recharge rate of as much as 1.3 inches per year, or 17 percent of the long-term mean precipitation. Model simulations also indicate that recharge would be virtually zero if the study area were covered by vegetation rather than dunes.     

HYDROLOGIC IMPLICATIONS OF GREATER GROUND‐WATER RECHARGE TO LAS VEGAS VALLEY,NEVADA1

ABSTRACT: Published estimates of natural recharge in Las Vegas Valley range between 21,000 and 35,000 acre‐feet per year. This study examined the underlying assumptions of previous investigations and evaluated the altitude‐precipitation relationships. Period‐of‐record averages from high altitude precipitation gages established in the 1940s through the 1990s, were used to determine strong local altitude‐precipitation relationships that indicate new total precipitation and natural recharge amounts and a new spatial distribution of that recharge. This investigation calculated about 51,000 acre‐feet per year of natural recharge in the Las Vegas Hydrographic Basin, with an additional 6,000 acre‐feet per year from areas tributary to Las Vegas Valley, for a total of 57,000 acre‐feet per year. The total amount of natural recharge is greater than estimates from earlier investigations and is consistent with a companion study of natural discharge, which estimated 53,000 acre‐feet per year of outflow. The hydrologic implications of greater recharge in Las Vegas Valley infer a more accurate ground‐water budget and a better understanding of ground‐water recharge that will be represented in a ground‐water model. Thus model based ground‐water management scenarios will more realistically access impacts to the ground‐water system.     

Isolated Spring Wetlands in the Great Basin and Mojave Deserts,USA: Potential Response of Vegetation to Groundwater Withdrawal

Desert springs, often the sole sources of water for wildlife and cattle, support wetland and wetland/upland transition ecosystems including rare and endemic species. In the basin and range province in Nevada, USA, springs in the Great Basin and Mojave deserts are sustained by interconnected deep carbonate and shallow basin-fill aquifers which are threatened by proposed groundwater withdrawal to sustain rapidly expanding urban areas, a common problem in arid regions worldwide. This paper draws on historic groundwater data, groundwater modeling, and studies of environmental controls of spring ecosystems to speculate on the potential effects of groundwater withdrawal and water table decline on spring-supported vegetation. The focus is on springs in the Great Basin and Mojave deserts representative of those that may be affected by future, planned groundwater withdrawal. Groundwater withdrawal is expected to reduce spring discharge directly through reduced flows from the shallow basin-fill aquifer or through reduction of the hydraulic head of the deep carbonate aquifer. This flow reduction will truncate the outflow stream, reducing the areal cover of wetland and wetland/upland transition vegetation. Lowering the local water table may also reduce the amount of upland phreatophytic vegetation by causing water levels to drop below plant rooting depths. Percolation of salts to surface soils may be reduced, eventually altering desert shrub cover from halophytes to nonhalophytes. The extent of these effects will vary among springs, based on their distance from extraction sites and location relative to regional groundwater flow paths. On-site monitoring of biotic variables (including cover of selected hygrophytes and phreatophytes) should be a necessary complement to the planned monitoring of local hydrologic conditions.     

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.     

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.     

Hydrological Variability and Uncertainty of Lower Missouri River Basin Under Changing Climate

The lower Missouri River Basin has experienced increasing streamflow and flooding events, with higher risk of extreme hydrologic impacts under changing climate. The newly available North American Regional Climate Change Assessment Program (NARCCAP) climate projections were used as atmospheric forcing for Soil and Water Assessment Tool (SWAT) model which runs with varying potential evapotranspiration (PET) methods to assess the hydrological change and uncertainty of 2040‐2069 over 1968‐1997. The NARCCAP temperature and precipitation predictions were refined using a bias correction method. The results show that, following the seasonal variability of precipitation, various water fluxes would increase in most seasons except the summer. Expected precipitation tends to increase in intensity with little change in frequency, triggering faster surface water concentration to form floods. The greatest streamflow increase would occur from November to February, increasing by around 10% on average. An increase of 3% occurs in the other months except for July and August in which river discharge decreases by around 2%. The climate predictions contribute more uncertainty annually, but PET algorithms gain more influence in winter or when other weather factors such as wind play a relatively more important role on evapotranspiration flux. This study predicts an even wetter environment compared to the historically very wet period, with the possibility of more flooding.     

Uncertainty Models for Estimates of Physical Characteristics of River Segments Over Large Areas

Spatially comprehensive estimates of the physical characteristics of river segments over large areas are required in many large‐scale analyses of river systems and for the management of multiple basins. Remote sensing and modeling are often used to estimate river characteristics over large areas, but the uncertainties associated with these estimates and their dependence on the physical characteristics of the segments and their catchments are seldom quantified. Using test data with varying degrees of independence, we derived analytical models of the uncertainty associated with estimates of upstream catchment area (CA), segment slope, and mean annual discharge for all river segments of a digital representation of the hydrographic network of France. Although there were strong relationships between our test data and estimates at the scale of France, there were also large relative local uncertainties, which varied with the physical characteristics of the segments and their catchments. Discharge and CA were relatively uncertain where discharge was low and catchments were small. Discharge uncertainty also increased in catchments with large rainfall events and low minimum temperature. The uncertainty of segment slope was strongly related to segment length. Our uncertainty models were consistent across large regions of France, suggesting some degree of generality. Their analytical formulation should facilitate their use in large‐scale ecological studies and simulation models.     

SIMULATION AND MAPPING OF SOIL-WATER CONDITIONS IN THE GREAT PLAINS1

ABSTRACT: Soil-water conditions provide valuable insight into the hydrologic system in an area. A soil-water balance quantitatively summarizes soil-water conditions and is based on climatic, soil, and vegetation characteristics that vary spatially and temporally. Soil-water balances in the Great Plains of the central United States were simulated for 1951–1980. Results of the simulations were mean annual estimates of infiltration, runoff, actual evapotranspiration, potential recharge, and consumptive water and irrigation requirements at 152 climatic data stations. A method was developed using a geographic information system to integrate and map the simulation results on the basis of spatially variable climatic, soil, and vegetation characteristics. As an example, simulated mean annual potential recharge was mapped. Mean annual potential-recharge rates ranged from less than 0.5 inch in much of the north-central and southwestern Great Plains to more than 10 inches in parts of eastern Texas and southwestern Arkansas.     

Comparison of Evapotranspiration Simulation Performance by APEX Model in Dryland and Irrigated Cropping Systems

Accurate estimation of evapotranspiration (ET) is essential to improve water use efficiency of crop production systems managed under different water regimes. The Agricultural Policy/Environmental eXtender (APEX) model was used to simulate ET using four potential ET (ET_p) methods. The objectives were to determine sensitive ET parameters in dryland and irrigated cropping systems and compare ET simulation in the two systems using multiple performance criteria. Measured ET and crop yield data from lysimeter fields located in the United States Department of Agriculture‐Agricultural Research Service Bushland, Texas were used for evaluation. The number of sensitive parameters was higher for dryland (11–14) than irrigated cropping systems (6–8). Only four input parameters: soil evaporation plant cover factor, root growth soil strength, maximum rain intercept, and rain intercept coefficient were sensitive in both cropping systems. Overall, it is possible to find a set of robust parameter values to simulate ET accurately in APEX in both cropping systems using any ET_p method. However, more computation time is required for dryland than irrigated cropping system due to a relatively larger number of sensitive input parameters. When all inputs are available, the Penman–Monteith method takes the shortest computation time to obtain one model run with robust parameter values in both cropping systems. However, in areas with limited datasets, one can still obtain reasonable ET simulations using either Priestley–Taylor or Hargreaves. Editor's note : This paper is part of the featured series on Optimizing Ogallala Aquifer Water Use to Sustain Food Systems. See the February 2019 issue for the introduction and background to the series.     

ADEQUACY OF SELECTED EVAPOTRANSPIRATION APPROXIMATIONS FOR HYDROLOGIC SIMULATION1

ABSTRACT: Evapotranspiration (ET) approximations, usually based on computed potential ET (PET) and diverse PET‐to‐ET conceptualizations, are routinely used in hydrologic analyses. This study presents an approach to incorporate measured (actual) ET data, increasingly available using micrometeorological methods, to define the adequacy of ET approximations for hydrologic simulation. The approach is demonstrated at a site where eddy correlation‐measured ET values were available. A baseline hydrologic model incorporating measured ET values was used to evaluate the sensitivity of simulated water levels, subsurface recharge, and surface runoff to error in four ET approximations. An annually invariant pattern of mean monthly vegetation coefficients was shown to be most effective, despite the substantial year‐to‐year variation in measured vegetation coefficients. The temporal variability of available water (precipitation minus ET) at the humid, subtropical site was largely controlled by the relatively high temporal variability of precipitation, benefiting the effectiveness of coarse ET approximations, a result that is likely to prevail at other humid sites.     

A System Dynamics Model for Conjunctive Management of Water Resources in the Snake River Basin

The Pacific Northwest is expected to witness changes in temperature and precipitation due to climate change. In this study, we enhance the Snake River Planning Model (SRPM) by modeling the feedback loop between incidental recharge and surface water supply resulting from surface water and groundwater extraction for irrigation and provide a case study involving climate change impacts and management scenarios. The new System Dynamics‐Snake River Planning Model (SD‐SRPM) is calibrated to flow at Box Canyon Springs located along a major outlet of the East Snake Plain Aquifer. A calibration of the model to flow at Box Canyon Springs, based on historic diversions (1950‐1995) resulted in an r² value of 0.74 and a validation (1996‐2005) r² value of 0.60. After adding irrigation entities to the model an r² value of 0.91, 0.88, and 0.87 were maintained for modeled vs. observed (1991‐2005) end‐of‐month reservoir content in Jackson Lake, Palisades, and American Falls, the three largest irrigation reservoirs in the system. The scenarios that compared the impacts of climate change were based on ensemble mean precipitation change scenarios and estimated changes to crop evapotranspiration (ET). Increased ET, despite increased precipitation, generally increased surface water shortages and discharge of springs. This study highlights the need to develop and implement models that integrate the human‐natural system to understand the impacts of climate change.