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.     

Discharge From the Edwards Aquifer Through the Leona River Floodplain,Uvalde, Texas1

Abstract: Analysis of results from an electrical resistivity survey, a magnetic survey, and an aquifer test performed on the Leona River floodplain in south‐central Texas indicates that ground‐water discharge from the Edwards Aquifer through the Leona River floodplain may be as great as 91.7 Mm³/year. When combined with an estimate of 8.8 Mm³/year for surface flow in the Leona River, as much as 100.5 Mm³/year could be discharged from the Edwards Aquifer through the Leona River floodplain. A value of 11,200 acre‐ft/year (13.82 Mm³/year) has been used as the calibration target in existing ground‐water models for total discharge from Leona Springs and the Leona River. Including ground water or underflow discharge would significantly increase the calibration target in future models. This refinement would improve the conceptualization of ground‐water flow in the western portion of the San Antonio segment of the Edwards Aquifer and would thereby allow for more accurate assessment and management of the ground‐water resources provided by the Edwards Aquifer.     

Decreased Runoff Response to Precipitation,Little Missouri River Basin,Northern Great Plains,USA

High variability in precipitation and streamflow in the semiarid northern Great Plains causes large uncertainty in water availability. This uncertainty is compounded by potential effects of future climate change. We examined historical variability in annual and growing season precipitation, temperature, and streamflow within the Little Missouri River Basin and identified differences in the runoff response to precipitation for the period 1976‐2012 compared to 1939‐1975 (n = 37 years in both cases). Computed mean values for the second half of the record showed little change (<5%) in annual or growing season precipitation, but average annual runoff at the basin outlet decreased by 22%, with 66% of the reduction in flow occurring during the growing season. Our results show a statistically significant (p < 0.10) 27% decrease in the annual runoff response to precipitation (runoff ratio). Surface‐water withdrawals for various uses appear to account for <12% of the reduction in average annual flow volume, and we found no published or reported evidence of substantial flow reduction caused by groundwater pumping in this basin. Results of our analysis suggest that increases in monthly average maximum and minimum temperatures, including >1°C increases in January through March, are the dominant driver of the observed decrease in runoff response to precipitation in the Little Missouri River Basin.     

Integrating Temporal Inequality into Conservation Planning to Improve Practice Design and Efficacy

In contrast to spatial inequality, there are currently no methods for leveraging information on temporal inequality to improve conservation efficacy. The objective of this study was to use Lorenz curves to quantify temporal inequality in surface runoff and tile drainage, identify controls on nutrient loading in these flowpaths, and develop design flows for structural conservation practices. Surface runoff (n = 94 site‐years) and tile drainage (n = 90 site‐years) were monitored on 40 fields in Ohio. Results showed, on average, 80% of nitrate‐nitrogen, soluble reactive phosphorus (P), and total P loads occurred between 7 and 12 days per year in surface runoff and between 32 and 58 days per year in tile drainage. Similar temporal inequality between discharge and load provided evidence that loading was transport‐limited and highlighted the critical role hydrologic connectivity plays in nutrient delivery from tile‐drained fields. Design flow criterion for sizing structural practices based on load reduction goals was developed by combining Lorenz curves and flow duration curves. Comparing temporal inequality between fields and the Maumee River, the largest tributary to the western Lake Erie Basin, revealed challenges associated with achieving watershed load reduction goals with field‐scale conservation. In‐field (i.e., improved nutrient and water management), edge‐of‐field (i.e., structural practices), and instream practices will all be required to meet nutrient reduction goals from tile‐drained watersheds.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

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.     

Field Evaluation of DRAINMOD 5.1 Under a Cold Climate: Simulation of Daily Midspan Water Table Depths and Drain Outflows1

Abstract: The hydrologic performance of DRAINMOD 5.1 was assessed for the southern Quebec region considering freezing/thawing conditions. A tile drained agricultural field in the Pike River watershed was instrumented to measure tile drainage volumes. The model was calibrated using water table depth and subsurface flow data over a two‐year period, while another two‐year dataset served to validate the model. DRAINMOD 5.1 accurately simulated the timing and magnitude of subsurface drainage events. The model also simulated the pattern of water table fluctuations with a good degree of accuracy. The R² between the observed and simulated daily WTD for calibration was >0.78, and that for validation was 0.93. The corresponding coefficients of efficiency (E) were >0.74 and 0.31. The R² and E values for calibration/validation of subsurface flow were 0.73/0.48 and 0.72/0.40, respectively. DRAINMOD simulated monthly subsurface flow quite accurately (E > 0.82 and R² > 0.84). The model precisely simulated daily/monthly drain flow over the entire year, including the winter months. Thus DRAINMOD 5.1 performed well in simulating the hydrology of a cold region.     

Effects of Land Use and Sample Location on Nitrate‐Stream Flow Hysteresis Descriptors during Storm Events

The U.S. Geological Survey's New Jersey and Iowa Water Science Centers deployed ultraviolet‐visible spectrophotometric sensors at water‐quality monitoring sites on the Passaic and Pompton Rivers at Two Bridges, New Jersey, on Toms River at Toms River, New Jersey, and on the North Raccoon River near Jefferson, Iowa to continuously measure in‐stream nitrate plus nitrite as nitrogen (NO₃ + NO₂) concentrations in conjunction with continuous stream flow measurements. Statistical analysis of NO₃ + NO₂ vs. stream discharge during storm events found statistically significant links between land use types and sampling site with the normalized area and rotational direction of NO₃ + NO₂‐stream discharge (N‐Q) hysteresis patterns. Statistically significant relations were also found between the normalized area of a hysteresis pattern and several flow parameters as well as the normalized area adjusted for rotational direction and minimum NO₃ + NO₂ concentrations. The mean normalized hysteresis area for forested land use was smaller than that of urban and agricultural land uses. The hysteresis rotational direction of the agricultural land use was opposite of that of the urban and undeveloped land uses. An r² of 0.81 for the relation between the minimum normalized NO₃ + NO₂ concentration during a storm vs. the normalized NO₃ + NO₂ concentration at peak flow suggested that dilution was the dominant process controlling NO₃ + NO₂ concentrations over the course of most storm events.     

The 2007‐2009 Drought in Athens,Georgia, United States: A Climatological Analysis and an Assessment of Future Water Availability1

Campana, Pete, John Knox, Andrew Grundstein, and John Dowd, 2012. The 2007‐2009 Drought in Athens, Georgia, United States: A Climatological Analysis and an Assessment of Future Water Availability. Journal of the American Water Resources Association (JAWRA) 48(2): 379‐390. DOI: 10.1111/j.1752‐1688.2011.00619.x Abstract: Population growth and development in many regions of the world increase the demand for water and vulnerability to water shortages. Our research provides a case study of how population growth can augment the severity of a drought. During 2007‐2009, a drought event that caused extreme societal impacts occurred in the Athens, Georgia region (defined as Clarke, Barrow, Oconee, and Jackson counties). An examination of drought indices and precipitation records indicates that conditions were severe, but not worse than during the 1925‐1927, 1954‐1956, and 1985‐1987 drought events. A drought of similar length to the 2007‐2009 drought would be expected to occur approximately every 25 years. Streamflow analysis shows that discharge levels in area streams were at a record low during 2007 before water restrictions were implemented, because of greater water usage caused by recent population increases. These population increases, combined with a lack of water conservation, led to severe water shortages in the Athens region during late 2007. Only after per capita usage decreased did water resources last despite continuing drought conditions through 2009. Retaining mitigation strategies and withdrawal levels such as seen during the height of the drought will be an essential strategy to prevent water shortages during future extreme drought events. The key mitigation strategy, independent local action to restrict water use in advance of state‐level restrictions, is now prohibited by Georgia State Law.     

On the Role of Spatial,Temporal, and Climatic Forces on Stream Sediment Loading from Rural and Urban Ecosystems

In this study, we characterize the greatest sediment loading events by their sediment delivery behavior; dominant climate, watershed, and antecedent conditions; and their seasonal distribution for rural and urban land uses. The study area is Paradise Creek Watershed, a mixed land use watershed in northern Idaho dominated by saturation excess processes in the upstream rural area and infiltration excess in the downstream urban area. We analyzed 12 years of continuous streamflow, precipitation, and watershed data at two monitoring stations. We identified 137 sediment loading events in the upstream rural section of the watershed and 191 events in the downstream urban section. During the majority of these events conditions were transport limited and the sediment flush occurred early in the event, generally in the first 20% of elapsed event time. Statistical analysis including two dozen explanatory variables showed peak discharge, event duration, and antecedent baseflow explained most of the variation in event sediment load at both stations and for the watershed as a whole (R² = 0.73‐0.78). In the rural area, saturated soils combined with spring snowmelt in March led to the greatest loading events. The urban area load contribution peaked in January, which could be a re‐suspension of streambed sediments from the previous water year. Throughout the study period, one event contributed, on average, 33% of the annual sediment load but only accounted for 2% of the time in a year.     

Mitigating Impact of Devils Lake Flooding on the Sheyenne River Sulfate Concentration

Devils Lake is a terminal lake located in northeast North Dakota. Because of its glacial origin and accumulated salts from evaporation, the lake has a high concentration of sulfate compared to the surrounding water bodies. From 1993 to 2011, Devils Lake water levels rose by ~10 m, which flooded surrounding communities and increased the chance of an overspill to the Sheyenne River. To control the flooding, the State of North Dakota constructed two outlets to pump the lake water to the river. However, the pumped water has raised concerns about of water quality degradation and potential flooding risk of the Sheyenne River. To investigate these perceived impacts, a Soil and Water Assessment Tool (SWAT) model was developed for the Sheyenne River and it was linked to a coupled SWAT and CE‐QUAL‐W2 model that was developed for Devils Lake in a previous study. While the current outlet schedule has attempted to maintain the total river discharge within the confines of a two‐year flood (36 m³/s), our simulation from 2012 to 2018 revealed that the diversion increased the Sheyenne River sulfate concentration from an average of 125 to >750 mg/L. Furthermore, a conceptual optimization model was developed with a goal of better preserving the water quality of the Sheyenne River while effectively mitigating the flooding of Devils Lake. The optimal solution provides a “win–win” outlet management that maintains the efficiency of the outlets while reducing the Sheyenne River sulfate concentration to ≤600 mg/L.     

Modeling Water Quantity and Sulfate Concentrations in the Devils Lake Watershed Using Coupled SWAT and CE‐QUAL‐W2

Devils Lake is an endorheic lake in the Red River of the North basin in northeastern North Dakota. During the last two decades, the lake water level has risen by nearly 10 m, causing floods that have cost more than 1 billion USD in mitigation measures. Another increase of approximately 1.5 m in the lake water level would cause spillage into the Sheyenne River. To alleviate this potentially catastrophic spillage, two artificial outlets were constructed. However, the artificial drainage of water into the Sheyenne River raises water quality concerns because the Devils Lake water contains significantly higher concentrations of dissolved solids, particularly sulfate. In this study, the Soil and Water Assessment Tool (SWAT) was coupled with the CE‐QUAL‐W2 model to simulate both water balance and sulfate concentrations in the lake. The SWAT model performed well in simulating daily flow in tributaries with E_NS > 0.5 and |PBIAS| < 25%, and reproduced the lake water level with a root mean square error of 0.35 m for the study period from 1995 to 2014. The water temperature and sulfate concentrations simulated by CE‐QUAL‐W2 for the lake are in general agreement with the field observations. The model results show that the operation of the two outlets since August 2005 has lowered the lake level by 0.70 m. Furthermore, the models show pumping water from the two outlets raises sulfate concentrations in the Sheyenne River from ~100 to >500 mg/L. 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.     

Observed Changes in Climate and Streamflow in the Upper Rio Grande Basin

Observed streamflow and climate data are used to test the hypothesis that climate change is already affecting Rio Grande streamflow volume derived from snowmelt runoff in ways consistent with model‐based projections of 21st‐Century streamflow. Annual and monthly changes in streamflow volume and surface climate variables on the Upper Rio Grande, near its headwaters in southern Colorado, are assessed for water years 1958–2015. Results indicate winter and spring season temperatures in the basin have increased significantly, April 1 snow water equivalent (SWE) has decreased by approximately 25%, and streamflow has declined slightly in the April–July snowmelt runoff season. Small increases in precipitation have reduced the impact of declining snowpack on trends in streamflow. Changes in the snowpack–runoff relationship are noticeable in hydrographs of mean monthly streamflow, but are most apparent in the changing ratios of precipitation (rain + snow, and SWE) to streamflow and in the declining fraction of runoff attributable to snowpack or winter precipitation. The observed changes provide observational confirmation for model projections of decreasing runoff attributable to snowpack, and demonstrate the decreasing utility of snowpack for predicting subsequent streamflow on a seasonal basis in the Upper Rio Grande Basin.     

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.     

AI-based runoff simulation based on remote sensing observations: A case study of two river basins in the United States and Canada

Data-driven techniques are used extensively for hydrologic time-series prediction. We created various data-driven models (DDMs) based on machine learning: long short-term memory (LSTM), support vector regression (SVR), extreme learning machines, and an artificial neural network with backpropagation, to define the optimal approach to predicting streamflow time series in the Carson River (California, USA) and Montmorency (Canada) catchments. The moderate resolution imaging spectroradiometer (MODIS) snow-coverage dataset was applied to improve the streamflow estimate. In addition to the DDMs, the conceptual snowmelt runoff model was applied to simulate and forecast daily streamflow. The four main predictor variables, namely snow-coverage (S-C), precipitation (P), maximum temperature (T_max), and minimum temperature (T_min), and their corresponding values for each river basin, were obtained from National Climatic Data Center and National Snow and Ice Data Center to develop the model. The most relevant predictor variable was chosen using the support vector machine-recursive feature elimination feature selection approach. The results show that incorporating the MODIS snow-coverage dataset improves the models' prediction accuracies in the snowmelt-dominated basin. SVR and LSTM exhibited the best performances (root mean square error = 8.63 and 9.80) using monthly and daily snowmelt time series, respectively. In summary, machine learning is a reliable method to forecast runoff as it can be employed in global climate forecasts that require high-volume data processing.     

AN ORGANIZED SIGNAL IN SNOWMELT RUNOFF OVER THE WESTERN UNITED STATES1

ABSTRACT: Daily‐to‐weekly discharge during the snowmelt season is highly correlated among river basins in the upper elevations of the central and southern Sierra Nevada (Carson, Walker, Tuolumne, Merced, San Joaquin, Kings, and Kern Rivers). In many cases, the upper Sierra Nevada watershed operates in a single mode (with varying catchment amplitudes). In some years, with appropriate lags, this mode extends to distant mountains. A reason for this coherence is the broad scale nature of synoptic features in atmospheric circulation, which provide anomalous insolation and temperature forcing that span a large region, sometimes the entire western U.S. These correlations may fall off dramatically, however, in dry years when the snowpack is spatially patchy.     

Field Testing the Riparian Ecosystem Management Model on a Riparian Buffer in the North Carolina Upper Coastal Plain

The riparian ecosystem management model (REMM) was field tested using five years (2005‐2009) of measured hydrologic and water quality data on a riparian buffer located in the Tar‐Pamlico River Basin, North Carolina. The buffer site received NO₃‐N loading from an agricultural field that was fertilized with inorganic fertilizer. Field results showed the buffer reduced groundwater NO₃‐N concentration moving to the stream over a five‐year period. REMM was calibrated hydrologically using daily field‐measured water table depths (WTDs), and with monthly NO₃‐N concentrations in groundwater wells. Results showed simulated WTDs and NO₃‐N concentrations in good agreement with measured values. The mean absolute error and Willmott's index of agreement for WTDs varied from 13‐45 cm and 0.72‐0.92, respectively, while the root mean square error and Willmott's index of agreement for NO₃‐N concentrations ranged from 1.04‐5.92 mg/l and 0.1‐0.86, respectively, over the five‐year period. REMM predicted plant nitrogen (N) uptake and denitrification were within ranges reported in other riparian buffer field studies. The calibrated and validated REMM was used to simulate 33 years of buffer performance at the site. Results showed that on average the buffer reduced NO₃‐N concentrations from 12 mg/l at the field edge to 0.7 mg/l at the stream edge over the simulation period, while the total N and NO₃‐N load reductions from the field edge to the stream were 77 and 82%, respectively.     

Modeling the Effects of Extreme Drought on Pollutant Transport Processes in the Yangtze River Estuary

Water resources in the Yangtze River Estuary (YRE), which is the vital water supply for Shanghai, decreased by approximately 2.45 Gm³ in 2006, the second‐worst recorded drought year. A numerical model was developed to investigate the effects of this extreme drought on pollutant transport processes in the YRE. The model was calibrated against observations and displayed good agreement. Residence time, a critical hydrodynamic indicator, was implemented to indicate pollutant transport processes. Numerical experiments were conducted to examine the possibly drought‐induced influences. The model results demonstrated that the influences on pollutant transport processes varied spatially and temporally, and these influences could partly explain the observed temporal and spatial variations of total nitrate in 2006. The area most susceptible to drought is in the north branch with 2‐11 days' extension of residence time. As the drought occurred in both the high and normal water periods, its influences were more significant during the normal water period with saltwater intrusion into the north branch. The drought also introduced a pollutant transport lag in timescale of approximately five days by diminishing the seaward advection flux with freshwater discharge. In 2006, the magnified tidal influence during the drought contributed more than usual to structuring pollutant transport, as the pollutant transport processes were intensely associated with tidal flow and tidal cycle.     

Estimating the Frequency and Extent of Overbank Flow of a Baseflow‐Dominated Stream Flowing through a Wetland

River floodplains provide critical habitat for a wide range of animal and plant species and reduce phosphorus and nitrogen loads in streams. It has been observed that baseflow‐dominated streams flowing through wetlands are commonly at or near bankfull and overflow their banks much more frequently than other streams. However, there is very little published quantitative support for this observation. The study focuses on a 1‐km reach of Black Earth Creek, a stream in the Midwestern United States (U.S.). We used one‐dimensional hydraulic modeling to estimate bankfull discharge at evenly spaced stream cross sections, and two‐dimensional modeling to quantitate the extent of wetland inundation as a function of discharge. We then used historical streamflow data from two U.S. Geological Survey gaging stations to quantitate the frequency of wetland inundation. For the with‐sediment case, the frequency of overbank conditions at the 38 cross sections in the wetland ranged from 3 to 85 days per year and averaged 43 days per year. Ten percent of the wetland was inundated for an average of 35 days per year. For the without‐sediment case, the frequency of overbank conditions ranged from 2.6 to 48 days per year and averaged 14 days per year. Also, 10% of the wetland was inundated for an average of 25 days per year. These unusually high rates of floodplain inundation are likely due in part to the very low stream gradient and shallow depths of overbank flow.     

The South‐to‐North Water Transfer Project of China: Environmental Implications and Monitoring Strategy1

Abstract: In 2002, China launched the South‐to‐North Water Transfer Project after completing a 50‐year feasibility study. By 2050, the three‐route (i.e., East, Middle, and West) project will be capable of transferring 44.8 billion m³/year of water from the water rich Yangtze River to the arid north to alleviate water shortage and help secure a balanced social and economic development across the nation. However, diversion of such a large quantity of water could profoundly change the riverine environment of the upper Yellow River and the lower reach of the Han River, a tributary of the Yangtze River and the water supplying area of the project’s Middle Route, because of changes in the annual discharge. Secondary salinization seems inevitable in the water receiving areas of the North China Plain, and decrease in the discharge of the Yangtze River will result in seawater intrusion into the Yangtze Delta. This paper describes the project and discusses its environmental implications. Additionally, a long‐term monitoring strategy under the umbrella of the Chinese Ecological Research Network is proposed for environmental monitoring.