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.     

Long-Term Relationship between Phosphorus Inputs and Wetland Phosphorus Concentrations in a Northern Everglades Marsh

Assessments of long-term relationships between changes innutrient inputs and wetland nutrient concentrations can becomplicated by fluctuations in other environmental factors aswell as by problems typical of long-term monitoring data.Consequently, statistical analysisof these types of data sets requirescareful consideration of environmental covariates, potentialbiases in the monitoring design, and irregularities caused bychanges in field sampling protocols. We evaluated therelationship between anthropogenic phosphorus (P) inputs andwater-column total P (TP) concentrations in a northernEverglades marsh by statistically analyzing available datacollected from several sampling programs over the past 20 years(1978–1997). Canal inputs of agricultural runoff contributemost of the P to the marsh and have produced a zone ofenrichment within the marsh during the past few decades.Regression analyses showed that both canal and marsh TPconcentrations increased during the 1980s and then decreased inthe 1990s. However, the statistical relationship between canal Pinputs and marsh TP, while significant, generally was weakexcept for marsh locations adjacent to the canal. Strongerrelationships existed between marsh TP and hydrologic parameterssuch as marsh water depth, which is controlled by changes inweather patterns and marsh management. In particular, dryconditions during the 1980s may have contributed to observedincreases in marsh P concentrations and the movement of a P`front' further into the marsh. Higher rainfall and water depthsand agricultural best management programs initiated during the1990s have been associated with reduced P concentrations incanal waters entering the marsh. While it is anticipated thatthis reduction eventually will result in lower marsh TPconcentrations, this effect is not yet evident, possibly due tointernal loading of P from enriched marsh soils. Our findingsillustrate some of the environmental factors that can complicateattempts to develop empirical relationships between P inputs andwetland P concentrations and to use such relationships to forecast changesin marsh concentrations based on past monitoring data alone.     

Implications of Conceptual Channel Representation on SWAT Streamflow and Sediment Modeling

Hydrologic modeling outputs are influenced by how a watershed system is represented. Channel routing is a typical example of the mathematical conceptualization of watershed landscape and processes in hydrologic modeling. We investigated the sensitivity of accuracy, equifinality, and uncertainty of Soil and Water Assessment Tool (SWAT) modeling to channel dimensions to demonstrate how a conceptual representation of a watershed system affects streamflow and sediment modeling. Results showed the amount of uncertainty and equifinality strongly responded to channel dimensions. On the other hand, the model performance did not significantly vary with the changes in the channel representation due to the degree of freedom allowed by the conceptual nature of hydrologic modeling in the parameter calibration. Such findings demonstrated good modeling performance statistics do not necessarily mean small output uncertainty, and partial improvements in the watershed representation may neither increase modeling accuracy nor reduce uncertainty. We also showed the equifinality and uncertainty of hydrologic modeling are case‐dependent rather than specific to models or regions, suggesting great caution should be used when attempting to transfer uncertainty analysis results to other modeling studies, especially for ungauged watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Defining the Diurnal Pattern of Snowmelt Using a Beta Distribution Function

Snow is an important component of the hydrologic cycle for many regions worldwide. In addition to vital water resources, snowmelt can be important for forest ecosystem dynamics and flood risk. However, standard design events in the United States lack a design snowmelt event, including only precipitation events, though snowmelt has been shown to be larger than rainfall. In this article, we present a method using hourly snow water equivalent data to develop and test a function for representing the diurnal pattern of snowmelt. A two‐parameter beta distribution function is modified for the purposes of this study and found to fit the pattern of snowmelt well with a root mean squared error of 0.008. Soil moisture sensors were additionally utilized to assess the timing of the snowmelt water outflow from the base of the snowpack that supports the shape of the function, but suggests that the timing of losses recorded on snow pillows lag as much as 3 h. Further testing of the function showed the shape of the function to be accurate. The methods developed and tested in this paper can be applied for design purposes comparing snowmelt and rainfall events or to improve hydrological models investigating processes such as streamflow or groundwater recharge.     

Does Including Soil Moisture Observations Improve Operational Streamflow Forecasts in Snow‐Dominated Watersheds?

Changing climate and growing water demand are increasing the need for robust streamflow forecasts. Historically, operational streamflow forecasts made by the Natural Resources Conservation Service have relied on precipitation and snow water equivalent observations from Snow Telemetry (SNOTEL) sites. We investigate whether also including SNOTEL soil moisture observations improve April‐July streamflow volume forecast accuracy at 0, 1, 2, and 3‐month lead times at 12 watersheds in Utah and California. We found statistically significant improvement in 0 and 3‐month lead time accuracy in 8 of 12 watersheds and 10 of 12 watersheds for 1 and 2‐month lead times. Surprisingly, these improvements were insensitive to soil moisture metrics derived from soil physical properties. Forecasts were made with volumetric water content (VWC) averaged from October 1 to the forecast date. By including VWC at the 0‐month lead time the forecasts explained 7.3% more variability and increased the streamflow volume accuracy by 8.4% on average compared to standard forecasts that already explained an average 77% of the variability. At 1 to 3‐month lead times, the inclusion of soil moisture explained 12.3‐26.3% more variability than the standard forecast on average. Our findings indicate including soil moisture observations increased statistical streamflow forecast accuracy and thus, could potentially improve water supply reliability in regions affected by changing snowpacks.     

Seawater Intrusion and Sustainable Yield of Basal Aquifers1

Liu, Clark C.K. and John J. Dai, 2012. Seawater Intrusion and Sustainable Yield of Basal Aquifers. Journal of the American Water Resources Association (JAWRA) 48(5): 861‐870. DOI: 10.1111/j.1752‐1688.2012.00659.x Abstract: Basal aquifers, in which freshwater floats on top of saltwater, are the major freshwater supply for the Hawaiian Islands, as well as many other coastal regions around the world. Under unexploited or natural conditions, freshwater and the underlying seawater are separated by a relatively sharp interface located below mean sea level at a depth of about 40 times the hydraulic head. With forced draft, the hydraulic head of a basal aquifer would decline and the sharp interface would move up. It is a serious problem of seawater intrusion as huge amounts of freshwater storage is replaced by saltwater. Also, with forced draft, the sharp interface is replaced by a transition zone in which the salinity increases downward from freshwater to saltwater. As pumping continues, the transition zone expands. The desirable source‐water salinity in Hawaii is about 2% of the seawater salinity. Therefore, the transition zone expansion is another serious problem of seawater intrusion. In this study, a robust analytical groundwater flow and salinity transport model (RAM2) was developed. RAM2 has a simple mathematical structure and its model parameters can be determined satisfactorily with the available field monitoring data. The usefulness of RAM2 as a viable management tool for coastal ground water management is demonstrated by applying it to determine the sustainable yield of the Pearl Harbor aquifer, a principal water supply source in Hawaii.     

Impacts of Climate Change on Hydrology and Water Resources in the Boise and Spokane River Basins1

Jin, Xin and Venkataramana Sridhar, 2012. Impacts of Climate Change on Hydrology and Water Resources in the Boise and Spokane River Basins. Journal of the American Water Resources Association (JAWRA) 48(2): 197‐220. DOI: 10.1111/j.1752‐1688.2011.00605.x Abstract: In the Pacific Northwest, warming climate has resulted in a lengthened growing season, declining snowpack, and earlier timing of spring runoff. This study characterizes the impact of climate change in two basins in Idaho, the Spokane River and the Boise River basins. We simulated the basin‐scale hydrology by coupling the downscaled precipitation and temperature outputs from a suite of global climate models and the Soil and Water Assessment Tool (SWAT), between 2010 and 2060 and assess the impacts of climate change on water resources in the region. For the Boise River basin, changes in precipitation ranged from ?3.8 to 36%. Changes in temperature were expected to be between 0.02 and 3.9°C. In the Spokane River region, changes in precipitation were expected to be between ?6.7 and 17.9%. Changes in temperature appeared between 0.1 and 3.5°C over a period of the next five decades between 2010 and 2060. Without bias‐correcting the simulated streamflow, in the Boise River basin, change in peak flows (March through June) was projected to range from ?58 to +106 m³/s and, for the Spokane River basin, the range was expected to be from ?198 to +88 m³/s. Both the basins exhibited substantial variability in precipitation, evapotranspiration, and recharge estimates, and this knowledge of possible hydrologic impacts at the watershed scale can help the stakeholders with possible options in their decision‐making process.     

Effects of Groundwater-Flow Paths On Nitrate Concentrations Across Two Riparian Forest Corridors1

Speiran, Gary K., 2010. Effects of Groundwater-Flow Paths on Nitrate Concentrations Across Two Riparian Forest Corridors. Journal of the American Water Resources Association (JAWRA) 46(2):246-260. DOI: 10.1111/j.1752-1688.2010.00427.x Abstract: Groundwater levels, apparent age, and chemistry from field sites and groundwater-flow modeling of hypothetical aquifers collectively indicate that groundwater-flow paths contribute to differences in nitrate concentrations across riparian corridors. At sites in Virginia (one coastal and one Piedmont), lowland forested wetlands separate upland fields from nearby surface waters (an estuary and a stream). At the coastal site, nitrate concentrations near the water table decreased from more than 10 mg/l beneath fields to 2 mg/l beneath a riparian forest buffer because recharge through the buffer forced water with concentrations greater than 5 mg/l to flow deeper beneath the buffer. Diurnal changes in groundwater levels up to 0.25 meters at the coastal site reflect flow from the water table into unsaturated soil where roots remove water and nitrate dissolved in it. Decreases in aquifer thickness caused by declines in the water table and decreases in horizontal hydraulic gradients from the uplands to the wetlands indicate that more than 95% of the groundwater discharged to the wetlands. Such discharge through organic soil can reduce nitrate concentrations by denitrification. Model simulations are consistent with field results, showing downward flow approaching toe slopes and surface waters to which groundwater discharges. These effects show the importance of buffer placement over use of fixed-width, streamside buffers to control nitrate concentrations.     

Reclaimed Mineland Curve Number Response to Temporal Distribution of Rainfall1

Warner, Richard C., Carmen T. Agouridis, Page T. Vingralek, and Alex W. Fogle, 2010. Reclaimed Mineland Curve Number Response to Temporal Distribution of Rainfall. Journal of the American Water Resources Association (JAWRA) 46(4): 724-732. DOI: 10.1111/j.1752-1688.2010.00444.x Abstract: The curve number (CN) method is a common technique to estimate runoff volume, and it is widely used in coal mining operations such as those in the Appalachian region of Kentucky. However, very little CN data are available for watersheds disturbed by surface mining and then reclaimed using traditional techniques. Furthermore, as the CN method does not readily account for variations in infiltration rates due to varying rainfall distributions, the selection of a single CN value to encompass all temporal rainfall distributions could lead engineers to substantially under- or over-size water detention structures used in mining operations or other land uses such as development. Using rainfall and runoff data from a surface coal mine located in the Cumberland Plateau of eastern Kentucky, CNs were computed for conventionally reclaimed lands. The effects of temporal rainfall distributions on CNs was also examined by classifying storms as intense, steady, multi-interval intense, or multi-interval steady. Results indicate that CNs for such reclaimed lands ranged from 62 to 94 with a mean value of 85. Temporal rainfall distributions were also shown to significantly affect CN values with intense storms having significantly higher CNs than multi-interval storms. These results indicate that a period of recovery is present between rainfall bursts of a multi-interval storm that allows depressional storage and infiltration rates to rebound.     

Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge1

Boggs, Kevin G., Robert W. Van Kirk, Gary S. Johnson, Jerry P. Fairley, and P. Steve Porter, 2010. Analytical Solutions to the Linearized Boussinesq Equation for Assessing the Effects of Recharge on Aquifer Discharge. Journal of the American Water Resources Association (JAWRA) 46(6):1116–1132. DOI: 10.1111/j.1752-1688.2010.00479.x Abstract: There is a need to develop a general understanding of how variations in aquifer recharge are reflected in discharge. Analytical solutions to the linearized Boussinesq equation governing flow in an unconfined aquifer provide a unified mathematical framework to quantify relationships among lag time, attenuation and distance between aquifer recharge and discharge and the effect of an up-gradient no-flow boundary. We applied this framework to three types of recharge: (1) instantaneous, (2) periodic, and (3) constant rate for a finite duration. When the temporal scale of recharge exceeds the diffusive aquifer time scale, recharge will be reflected in discharge quickly and with little attenuation. When aquifer time scale is large, most recharge events are shorter in scale than that of the aquifer, resulting in large attenuation. Attenuation is more sensitive to boundary effects than lag time, and boundary effects increase as recharge time scale increases. Boundary effects can often be ignored when the recharge source is farther than 1/3 of the domain length away from the no-flow boundary. We illustrate analytical results with application to the economically critical Eastern Snake River Plain Aquifer in Idaho. In this aquifer, detectable annual and decadal cycles in discharge can result from recharge no farther than 20 and 60 km away from the discharge point, respectively. The effects of more distant, long-term recharge can be detected only after a time lag of several decades.