Water balance characterization of the early 21st century drought in the western United States

Monthly temperature and precipitation data for 923 United States Geological Survey 8-digit hydrologic units are used as inputs to a monthly water balance model to compute monthly actual evapotranspiration, soil moisture storage, and runoff across the western United States (U.S.) for the period 1900 through 2020. Time series of these water balance variables are examined to characterize and explain the dry conditions across the western U.S. since the year 2000. Results indicate that although precipitation deficits account for most of the changes in actual evapotranspiration and runoff, increases in temperature primarily explain decreases in soil moisture storage. Specifically, temperature has been particularly impactful on the magnitude of negative departures of soil moisture storage during the spring (April through June) and summer (July through September) seasons. These effects on soil moisture may be particularly detrimental to agriculture in regions already stressed by drought such as the western U.S.     

SPATIAL VARIABILITY IN WATER‐BALANCE MODEL PERFORMANCE IN THE CONTERMINOUS UNITED STATES1

ABSTRACT: A monthly water‐balance (WB) model was tested in 44 river basins from diverse physiographic and climatic regions across the conterminous United States (U.S.). The WB model includes the concepts of climatic water supply and climatic water demand, seasonality in climatic water supply and demand, and soil‐moisture storage. Exhaustive search techniques were employed to determine the optimal set of precipitation and temperature stations, and the optimal set of WB model parameters to use for each basin. It was found that the WB model worked best for basins with: (1) a mean elevation less than 450 meters or greater than 2000 meters, and/or (2) monthly runoff that is greater than 5 millimeters (mm) more than 80 percent of the time. In a separate analysis, a multiple linear regression (MLR) was computed using the adjusted R‐square values obtained by comparing measured and estimated monthly runoff of the original 44 river basins as the dependent variable, and combinations of various independent variables [streamflow gauge latitude, longitude, and elevation; basin area, the long‐term mean and standard deviation of annual precipitation; temperature and runoff; and low‐flow statistics (i.e., the percentage of months with monthly runoff that is less than 5 mm)]. Results from the MLR study showed that the reliability of a WB model for application in a specific region can be estimated from mean basin elevation and the percentage of months with gauged runoff less than 5 mm. The MLR equations were subsequently used to estimate adjusted R‐square values for 1,646 gauging stations across the conterminous U.S. Results of this study indicate that WB models can be used reliably to estimate monthly runoff in the eastern U.S., mountainous areas of the western U.S., and the Pacific Northwest. Applications of monthly WB models in the central U.S. can lead to uncertain estimates of runoff.     

Operational Evapotranspiration Mapping Using Remote Sensing and Weather Datasets: A New Parameterization for the SSEB Approach

The increasing availability of multi‐scale remotely sensed data and global weather datasets is allowing the estimation of evapotranspiration (ET) at multiple scales. We present a simple but robust method that uses remotely sensed thermal data and model‐assimilated weather fields to produce ET for the contiguous United States (CONUS) at monthly and seasonal time scales. The method is based on the Simplified Surface Energy Balance (SSEB) model, which is now parameterized for operational applications, renamed as SSEBop. The innovative aspect of the SSEBop is that it uses predefined boundary conditions that are unique to each pixel for the “hot” and “cold” reference conditions. The SSEBop model was used for computing ET for 12 years (2000‐2011) using the MODIS and Global Data Assimilation System (GDAS) data streams. SSEBop ET results compared reasonably well with monthly eddy covariance ET data explaining 64% of the observed variability across diverse ecosystems in the CONUS during 2005. Twelve annual ET anomalies (2000‐2011) depicted the spatial extent and severity of the commonly known drought years in the CONUS. More research is required to improve the representation of the predefined boundary conditions in complex terrain at small spatial scales. SSEBop model was found to be a promising approach to conduct water use studies in the CONUS, with a similar opportunity in other parts of the world. The approach can also be applied with other thermal sensors such as Landsat.     

Water Footprint of 65 Mid‐ to Large‐Sized U.S. Cities and Their Metropolitan Areas

Conventional indicators of water use for urban areas account primarily for direct water use. In contrast, our objective here is to employ the water footprint (WF) concept and methodology to include the virtual or indirect water use to assess the production‐side and consumption‐side WF of 65 United States (U.S.) cities. The 65 cities include the largest metropolitan areas and some of the major mid‐sized cities in the U.S. We use metropolitan areas to define our city boundaries as this is the native spatial resolution of the main datasets used. To estimate the urban WFs, we integrated large and disparate datasets, including commodity flow (agricultural, livestock, and industrial commodities), water use, and socioeconomic data. By analyzing the estimated WF values, we found indirect water use accounts on average for 66% of the WF of consumption. We found some cities are net virtual water exporters (11 of 65) because they rely heavily on direct water uses or are heavy producers of industrial commodities. Also, WF patterns vary widely across the U.S. but regional patterns seem to emerge. For example, the dense cities of the U.S. northeast megaregion have a significantly low per capita WF relative to the other cities, while cities in the Gulf Coast megaregion have a significantly higher industrial WF of production and consumption. Furthermore, there is inequality in the WF of consumption where a few cities account for a disproportionate share of the total U.S. urban water uses.     

The Water‐Year Water Balance of the Colorado River Basin

Model‐estimated monthly water balance components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff (R)) for 146 United States (U.S.) Geological Survey 8‐digit hydrologic units located in the Colorado River Basin (CRB) are used to examine the temporal and spatial variability of the CRB water balance for water years 1901 through 2014 (a water year is the period from October 1 of one year through September 30 of the following year). Results indicate that the CRB can be divided into six subregions with similar temporal variability in monthly R. The water balance analyses indicated that approximately 75% of total water‐year R is generated by just one CRB subregion and that most of the R in the basin is derived from surplus (S) water generated during the months of October through April. Furthermore, the analyses show that temporal variability in S is largely controlled by the occurrence of negative atmospheric pressure anomalies over the northwestern conterminous U.S. (CONUS) and positive atmospheric pressure anomalies over the southeastern CONUS. This combination of atmospheric pressure anomalies results in an anomalous flow of moist air from the North Pacific Ocean into the CRB, particularly the Upper CRB. Additionally, the occurrence of extreme dry and wet periods in the CRB appears to be related to variability of the Atlantic Multidecadal Oscillation and the Pacific Decadal Oscillation.     

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.     

REGIONAL HYDROLOGIC RESPONSE OF LOBLOLLY PINE TO AIR TEMPERATURE AND PRECIPITATION CHANGES1

ABSTRACT: Large deviations in average annual air temperatures and total annual precipitation were observed across the southern United States during the last 50 years, and these fluctuations could become even larger during the next century. We used PnET-IIS, a monthly time-step forest process model that uses soil, vegetation, and climate inputs to assess the influence of changing climate on southern U.S. pine forest water use. After model predictions of historic drainage were validated, the potential influences of climate change on loblolly pine forest water use was assessed across the region using historic (1951 to 1984) monthly precipitation and air temperature which were modified by two general circulation models (GCMs). The GCMs predicted a 3.2°C to 7.2°C increase in average monthly air temperature, a -24 percent to + 31 percent change in monthly precipitation and a -1 percent to + 3 percent change in annual precipitation. As a comparison to the GCMs, a minimum climate change scenario using a constant 2°C increase in monthly air temperature and a 20 percent increase in monthly precipitation was run in conjunction with historic climate data. Predicted changes in forest water drainage were highly dependent on the GCM used. PnET-IIS predicted that along the northern range of loblolly pine, water yield would decrease with increasing leaf area, total evapotranspiration and soil water stress. However, across most of the southern U.S., PnET-IIS predicted decreased leaf area, total evapotranspiration, and soil water stress with an associated increase in water yield. Depending on the GCM and geographic location, predicted leaf area decreased to a point which would no longer sustain loblolly pine forests, and thus indicated a decrease in the southern most range of the species within the region. These results should be evaluated in relation to other changing environmental factors (i.e., CO₂ and O₃) which are not present in the current model.     

Spatiotemporal Evaluation of Simulated Evapotranspiration and Streamflow over Texas Using the WRF‐Hydro‐RAPID Modeling Framework

This study assesses a large‐scale hydrologic modeling framework (WRF‐Hydro‐RAPID) in terms of its high‐resolution simulation of evapotranspiration (ET) and streamflow over Texas (drainage area: 464,135 km²). The reference observations used include eight‐day ET data from MODIS and FLUXNET, and daily river discharge data from 271 U.S. Geological Survey gauges located across a climate gradient. A recursive digital filter is applied to decompose the river discharge into surface runoff and base flow for comparison with the model counterparts. While the routing component of the model is pre‐calibrated, the land component is uncalibrated. Results show the model performance for ET and runoff is aridity‐dependent. ET is better predicted in a wet year than in a dry year. Streamflow is better predicted in wet regions with the highest efficiency ~0.7. In comparison, streamflow is most poorly predicted in dry regions with a large positive bias. Modeled ET bias is more strongly correlated with the base flow bias than surface runoff bias. These results complement previous evaluations by incorporating more spatial details. They also help identify potential processes for future model improvements. Indeed, improving the dry region streamflow simulation would require synergistic enhancements of ET, soil moisture and groundwater parameterizations in the current model configuration. Our assessments are important preliminary steps towards accurate large‐scale hydrologic forecasts.     

Water Supply and Stormwater Management Benefits of Residential Rainwater Harvesting in U.S. Cities

This article presents an analysis of the projected performance of urban residential rainwater harvesting systems in the United States (U.S.). The objectives are to quantify for 23 cities in seven climatic regions (1) water supply provided from rainwater harvested at a residential parcel and (2) stormwater runoff reduction from a residential drainage catchment. Water‐saving efficiency is determined using a water‐balance approach applied at a daily time step for a range of rainwater cistern sizes. The results show that performance is a function of cistern size and climatic pattern. A single rain barrel (190 l [50 gal]) installed at a residential parcel is able to provide approximately 50% water‐saving efficiency for the nonpotable indoor water demand scenario in cities of the East Coast, Southeast, Midwest, and Pacific Northwest, but <30% water‐saving efficiency in cities of the Mountain West, Southwest, and most of California. Stormwater management benefits are quantified using the U.S. Environmental Protection Agency Storm Water Management Model. The results indicate that rainwater harvesting can reduce stormwater runoff volume up to 20% in semiarid regions, and less in regions receiving greater rainfall amounts for a long‐term simulation. Overall, the results suggest that U.S. cities and individual residents can benefit from implementing rainwater harvesting as a stormwater control measure and as an alternative source of water.     

Hydrologic Model Framework for Water Resources Planning in the Santa Cruz River,Southern Arizona1

Abstract: The authors develop a model framework that includes a set of hydrologic modules as a water resources management and planning tool for the upper Santa Cruz River near the Mexican border, Southern Arizona. The modules consist of: (1) stochastic generation of hourly precipitation scenarios that maintain the characteristics and variability of a 45‐year hourly precipitation record from a nearby rain gauge; (2) conceptual transformation of generated precipitation into daily streamflow using varied infiltration rates and estimates of the basin antecedent moisture conditions; and (3) surface‐water to ground‐water interaction for four downstream microbasins that accounts for alluvial ground‐water recharge, and ET and pumping losses. To maintain the large inter‐annual variability of streamflow as prevails in Southern Arizona, the model framework is constructed to produce three types of seasonal winter and summer categories of streamflow (i.e., wet, medium, or dry). Long‐term (i.e., 100 years) realizations (ensembles) are generated by the above described model framework that reflects two different regimes of inter annual variability. The first regime is that of the historic streamflow gauge record. The second regime is that of the tree ring reconstructed precipitation, which was derived for the study location. Generated flow ensembles for these two regimes are used to evaluate the risk that the regional four ground‐water microbasins decline below a preset storage threshold under different operational water utilization scenarios.     

Modeling Streamflow and Water Quality Sensitivity to Climate Change and Urban Development in 20 U.S. Watersheds

Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.     

Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.     

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

Precipitation Amount and Intensity Trends Across Southwest Saudi Arabia

This study examines precipitation accumulation and intensity trends across a region in southwest Saudi Arabia characterized by distinct seasonal weather patterns and mountainous terrain. The region is an example of an arid/semiarid area faced with maintaining sustainable water resources with a growing population. Annual and seasonal trends in precipitation amount were examined from 29 rain gages divided among four geographically unique regions from 1945/1946 to 2009. Two of the regions displayed significantly declining annual trends over the time series using a Mann‐Kendall test modified for autocorrelation (α < 0.05). Seasonal analysis revealed insignificant declining trends in at least two of the regions during each season. The largest and most consistent declining trends occurred during wintertime where all regions experienced negative trends. Several intensity metrics were examined in the study area from four additional stations containing daily data from 1985 to 2011. Intensity metrics included total precipitation, wet day count, simple intensity index, maximum daily annual rainfall, and upper/lower precipitation distribution changes. In general, no coherent trends were found among the daily stations suggesting precipitation is intensifying across the study area. The work represents the first of its size in the study area, and one of few in the region due to the lack of available long‐term data needed to properly examine precipitation changes.     

Estimation of Evapotranspiration Across the Conterminous United States Using a Regression With Climate and Land‐Cover Data1

Sanford, Ward E. and David L. Selnick, 2012. Estimation of Evapotranspiration Across the Conterminous United States Using a Regression with Climate and Land‐Cover Data. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12010 Abstract: Evapotranspiration (ET) is an important quantity for water resource managers to know because it often represents the largest sink for precipitation (P) arriving at the land surface. In order to estimate actual ET across the conterminous United States (U.S.) in this study, a water‐balance method was combined with a climate and land‐cover regression equation. Precipitation and streamflow records were compiled for 838 watersheds for 1971‐2000 across the U.S. to obtain long‐term estimates of actual ET. A regression equation was developed that related the ratio ET/P to climate and land‐cover variables within those watersheds. Precipitation and temperatures were used from the PRISM climate dataset, and land‐cover data were used from the USGS National Land Cover Dataset. Results indicate that ET can be predicted relatively well at a watershed or county scale with readily available climate variables alone, and that land‐cover data can also improve those predictions. Using the climate and land‐cover data at an 800‐m scale and then averaging to the county scale, maps were produced showing estimates of ET and ET/P for the entire conterminous U.S. Using the regression equation, such maps could also be made for more detailed state coverages, or for other areas of the world where climate and land‐cover data are plentiful.     

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.     

Lake Water Levels and Associated Hydrologic Characteristics in the Conterminous U.S.

Establishing baseline hydrologic characteristics for lakes in the United States (U.S.) is critical to evaluate changes to lake hydrology. We used the U.S. Environmental Protection Agency National Lakes Assessment 2007 and 2012 surveys to assess hydrologic characteristics of a population of ~45,000 lakes in the conterminous U.S. based on probability samples of ~1,000 lakes/yr distributed across nine ecoregions. Lake hydrologic study variables include water‐level drawdown (i.e., vertical decline and horizontal littoral exposure) and two water stable isotope‐derived parameters: evaporation‐to‐inflow (E:I) and water residence time. We present (1) national and regional distributions of the study variables for both natural and man‐made lakes and (2) differences in these characteristics between 2007 and 2012. In 2007, 59% of the population of U.S. lakes had Greater than normal or Excessive drawdown relative to water levels in ecoregional reference lakes with minimal human disturbances; whereas in 2012, only 20% of lakes were significantly drawn down beyond normal ranges. Water isotope‐derived variables did not differ significantly between survey years in contrast to drawdown. Median E:I was 20% indicating that flow‐through processes dominated lake water regimes. For 75% of U.S. lakes, water residence time was less than one year and was longer in natural vs. man‐made lakes. Our study provides baseline ranges to assess local and regional lake hydrologic status and inform management decisions in changing environmental conditions.     

A MONTHLY TIME SERIES MODEL OF MUNICIPAL WATER DEMAND1

A multivariate time series model is formulated to study monthly variations in municipal water demand. The left hand side variable in the multivariate regression model is municipal water demand (gallons per connection per day) and the right hand side contains (explanatory) variables which include price (constant dollars), average temperature, total precipitation, and percentage of daylight hours. The application of the regression model to Salt Lake City Water Department data produced a high multiple correlation coefficient and F-statistic. The regression coefficients for the right hand side variables all have the appropriate sign. In an ex post forecast, the model accurately predicts monthly variations in municipal water demand. The proposed monthly multivariate model is not only found useful for forecasting water demand, but also useful for predicting and studying the impact of nonstructural management decisions such as the effect of price changes, peak load pricing methods, and other water conservation programs.     

Hydroclimatology of the Mississippi River Basin

Model estimated monthly water balance (WB) components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff [R]) for 848 United States (U.S.) Geological Survey 8‐digit hydrologic units located in the Mississippi River Basin (MRB) are used to examine the temporal and spatial variability of the MRB WB for water years 1901 through 2014. Results indicate the MRB can be divided into nine subregions with similar temporal variability in R. The WB analyses indicated ~79% of total water‐year MRB runoff is generated by four of the nine subregions and most of the R in the basin is derived from surplus (S) water during the months of December through May. Furthermore, the analyses showed temporal variability in S is largely controlled by the occurrence of negative atmospheric pressure anomalies over the western U.S. and positive atmospheric pressure anomalies over the eastern U.S. coast. This combination of atmospheric pressure anomalies results in an anomalous flow of moist air from the Gulf of Mexico into the MRB. In the context of paleo‐climate reconstructions of the Palmer Drought Severity Index, since about 1900 the MRB has experienced wetter conditions than were experienced during the previous 500 years.     

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.