Spatial and Temporal Patterns of Low Streamflow and Precipitation Changes in the Chesapeake Bay Watershed

Spatial and temporal patterns in low streamflows were investigated for 183 streamgages located in the Chesapeake Bay Watershed for the period 1939–2013. Metrics that represent different aspects of the frequency and magnitude of low streamflows were examined for trends: (1) the annual time series of seven‐day average minimum streamflow, (2) the scaled average deficit at or below the 2% mean daily streamflow value relative to a base period between 1939 and 1970, and (3) the annual number of days below the 2% threshold. Trends in these statistics showed spatial cohesion, with increasing low streamflow volume at streamgages located in the northern uplands of the Chesapeake Bay Watershed and decreasing low streamflow volume at streamgages in the southern part of the watershed. For a small subset of streamgages (12%), conflicting trend patterns were observed between the seven‐day average minimum streamflow and the below‐threshold time series and these appear to be related to upstream diversions or the influence of reservoir‐influenced streamflows in their contributing watersheds. Using multivariate classification techniques, mean annual precipitation and fraction of precipitation falling as snow appear to be broad controls of increasing and decreasing low‐flow trends. Further investigation of seasonal precipitation patterns shows summer rainfall patterns, driven by the Atlantic Multidecadal Oscillation, as the main driver of low streamflows in the Chesapeake Bay Watershed.     

USE OF HISTORIC FLOOD INFORMATION IN ESTIMATING FLOOD PEAKS ON UNGAGED WATERSHEDS1

Regional procedures to estimate flood magnitudes for ungaged watersheds typically ignore available site-specific historic flood information such as high water marks and the corresponding flow estimates, otherwise referred to as limited site-specific historic (LSSH) flood data. A procedure to construct flood frequency curves on the basis of LSSH flood observations is presented. Simple inverse variance weighting is employed to systematically combine flood estimates obtained from the LSSH data base with those from a regional procedure to obtain improved estimtes of flood peaks on the ungaged watershed. For the region studied, the variance weighted estimates of flow had a lower logarithmic standard error than either the regional or the LSSH flow estimates, when compared to the estimates determined by three standard distributions for gaged watersheds investigated in the development of the methodology. Use of the simple inverse variance weighting procedure is recommended when “reliable” estimates of LSSH floods for the ungaged site are available.     

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.     

MODELING STREAMFLOW FROM ARTIFICIALLY DRAINED AGRICULTURAL WATERSHEDS IN ILLINOIS1

ABSTRACT: Evaluation of the applicability and validity of hydrologic simulation models for various cropping systems in different hydrogeologic and soil conditions is needed for a range of spatial scales. We calibrated and tested the ADAPT model for simulating streamflow from 552 to 1,985 km² watersheds in central Illinois, where more than 79 percent of the land is used for maize‐soybean production and tile drainage is common. Model calibration was performed with a seven year period (1987 through1993) of measured streamflow from one of the watersheds, and model testing was done using independent weather and measured streamflow data from the two neighboring watersheds for the same seven year period. Simulations of annual streamflow were accurate with a coefficient of determination and Willmott's index of agreement of 0.98 and 0.99, respectively. For simulation of monthly streamflow, Willmott's index of agreement ranged from 0.93 to 0.95. For simulation of daily streamflow, Willmott's index of agreement ranged from 0.84 to 0.85. The daily simulations challenged the temporal and spatial resolution of our measured precipitation data. Discrepancies between simulated and measured data may result from the model's inability to effectively address frozen soils and snowmelt runoff processes and in accurately representing evapotranspiration.     

Evidence for Changing Flood Risk in New England Since the Late 20th Century1

Abstract: Long‐term flow records for watersheds with minimal human influence have shown trends in recent decades toward increasing streamflow at regional and national scales, especially for low flow quantiles like the annual minimum and annual median flows. Trends for high flow quantiles are less clear, despite recent research showing increased precipitation in the conterminous United States over the last century that has been brought about primarily by an increased frequency and intensity of events in the upper 10th percentile of the daily precipitation distribution – particularly in the Northeast. This study investigates trends in 28 long‐term annual flood series for New England watersheds with dominantly natural streamflow. The flood series are an average of 75 years in length and are continuous through 2006. Twenty‐five series show upward trends via the nonparametric Mann‐Kendall test, 40% (10) of which are statistically significant (p < 0.1). Moreover, an average standardized departures series for 23 of the study gages indicates that increasing flood magnitudes in New England occurred as a step change around 1970. The timing of this is broadly synchronous with a phase change in the low frequency variability of the North Atlantic Oscillation, a prominent upper atmospheric circulation pattern that is known to effect climate variability along the United States east coast. Identifiable hydroclimatic shifts should be considered when the affected flow records are used for flood frequency analyses. Special treatment of the flood series can improve the analyses and provide better estimates of flood magnitudes and frequencies under the prevailing hydroclimatic condition.     

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

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

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

MODIS Biophysical States and NEXRAD Precipitation in a Statistical Evaluation of Antecedent Moisture Condition and Streamflow1

Abstract: The potential of remotely sensed time series of biophysical states of landscape to characterize soil moisture condition antecedent to radar estimates of precipitation is assessed in a statistical prediction model of streamflow in a 1,420 km² watershed in south‐central Texas, Moderate Resolution Imaging Spectroradiometer (MODIS) time series biophysical products offer significant opportunities to characterize and quantify hydrologic state variables such as land surface temperature (LST) and vegetation state and status. Together with Next Generation Weather Radar (NEXRAD) precipitation estimates for the period 2002 through 2005, 16 raw and deseasoned time series of LST (day and night), vegetation indices, infrared reflectances, and water stress indices were linearly regressed against observed watershed streamflow on an eight‐day aggregated time period. Time offsets of 0 (synchronous with streamflow event), 8, and 16 days (leading streamflow event) were assessed for each of the 16 parameters to evaluate antecedent effects. The model results indicated a reasonable correlation (r² = 0.67) when precipitation, daytime LST advanced 16 days, and a deseasoned moisture stress index were regressed against log‐transformed streamflow. The estimation model was applied to a validation period from January 2006 through March 2007, a period of 12 months of regional drought and base‐flow conditions followed by three months of above normal rainfall and a flood event. The model resulted in a Nash‐Sutcliffe estimation efficiency (E) of 0.45 for flow series (in log‐space) for the full 15‐month period, ?0.03 for the 2006 drought condition period, and 0.87 for the 2007 wet condition period. The overall model had a relative volume error of ?32%. The contribution of parameter uncertainties to model discrepancy was evaluated.     

Augmenting Watershed Model Calibration with Incorporation of Ancillary Data Sources and Qualitative Soft Data Sources

Watershed simulation models such as the Soil & Water Assessment Tool (SWAT) can be calibrated using “hard data” such as temporal streamflow observations; however, users may find upon examination of model outputs, that the calibrated models may not reflect actual watershed behavior. Thus, it is often advantageous to use “soft data” (i.e., qualitative knowledge such as expected denitrification rates that observed time series do not typically exist) to ensure that the calibrated model is representative of the real world. The primary objective of this study is to evaluate the efficacy of coupling SWAT‐Check (a post‐evaluation framework for SWAT outputs) and IPEAT‐SD (Integrated Parameter Estimation and Uncertainty Analysis Tool‐Soft & hard Data evaluation) to constrain the bounds of soft data during SWAT auto‐calibration. IPEAT‐SD integrates 59 soft data variables to ensure SWAT does not violate physical processes known to occur in watersheds. IPEAT‐SD was evaluated for two case studies where soft data such as denitrification rate, nitrate attributed from subsurface flow to total discharge ratio, and total sediment loading were used to conduct model calibration. Results indicated that SWAT model outputs may not satisfy reasonable soft data responses without providing pre‐defined bounds. IPEAT‐SD provides an efficient and rigorous framework for users to conduct future studies while considering both soft data and traditional hard information measures in watershed modeling.     

APPLICATION OF THE RHESSys MODEL TO A CALIFORNIA SEMIARID SHRUBLAND WATERSHED1

ABSTRACT: Distributed hydrologic models which link seasonal streamflow and soil moisture patterns with spatial patterns of vegetation are important tools for understanding the sensitivity of Mediterranean type ecosystems to future climate and land use change. RHESSys (Regional Hydro‐Ecologic Simulation System) is a coupled spatially distributed hydroecological model that is designed to be able to represent these feedbacks between hydrologic and vegetation carbon and nutrient cycling processes. However, RHESSys has not previously been applied to semiarid shrubland watersheds. In this study, the hydrologic submodel of RHESSys is evaluated by comparing model predictions of monthly and annual streamflow to stream gage data and by comparing RHESSys behavior to that of another hydrologic model of similar complexity, MIKESHE, for a 34 km² watershed near Santa Barbara, California. In model intercomparison, the differences in predictions of temporal patterns in streamflow, sensitivity of model predictions to calibration parameters and landscape representation, and differences in model estimates of soil moisture patterns are explored. Results from this study show that both models adequately predict seasonal patterns of streamflow response relative to observed data, but differ significantly in terms of estimates of soil moisture patterns and sensitivity of those patterns to the scale of landscape tessellation used to derive spatially distributed elements. This sensitivity has implications for implementing RHESSys as a tool to investigate interactions between hydrology and ecosystem processes.     

Using Paleo Reconstructions to Improve Streamflow Forecast Lead Time in the Western United States

In water stressed regions, water managers are exploring new horizons that would help in long‐range streamflow forecasts. Oceanic‐atmospheric oscillations have been shown to influence streamflow variability. In this study, long‐lead time streamflow forecasts are made using a multiclass kernel‐based data‐driven support vector machine (SVM) model. The extended streamflow records based on tree ring reconstructions were used to provide a longer time series data. Reconstructed data were used from 1658 to 1952 and the instrumental record was used from 1953 to 2007. Reconstructions for oceanic‐atmospheric oscillations included the El Niño‐Southern Oscillation, Pacific Decadal Oscillation, Atlantic Multidecadal Oscillation, and North Atlantic Oscillation. Streamflow forecasts using all four oscillations were made with one‐year to five‐year lead times for 21 gages in the western United States. This is the first study that uses both instrumental and reconstructed data of oscillations in SVM model to improve streamflow forecast lead time. SVM model was able to provide “satisfactory” to “very good” forecasts with one‐ to five‐year lead time for the selected gages. The use of all the oscillation indices helped in achieving better predictability compared to using individual oscillations. The SVM modeling results are better when compared with multiple linear regression model forecasts. The findings are statistical in nature and are expected to be useful for long‐term water resources planning and management.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

PARAMETER UNCERTAINTY IN RAINFALL-RUNOFF MODELING FOR POLDER SYSTEM DESIGN1

ABSTRACT: An approach is developed for incorporating the uncertainty of parameters for estimating runoff in the design of polder systems in ungaged watersheds. Monte Carlo Simulation is used to derive a set of realizations of streamflow hydrographs for a given design rainstorm using the U. S. Soil Conservation Service (SCS) unit hydrograph model. The inverse of the SCS curve number, which is a function of the antecedent runoff condition in the SCS model, is the random input in the Monte Carlo Simulation. Monte Carlo realizations of streamfiow hydrographs are used to simulate the performance of a polder flood protection system. From this simulation the probability of occurrence of flood levels for a particular hydraulic design may be used to evaluate its effectiveness. This approach is demonstrated for the Pluit Polder flood protection system for the City of Jakarta, Indonesia. While the results of the application indicate that uncertainty in the antecedent runoff condition is important, the effects of uncertainty in rainfall data, in additional runoff parameters, such as time to peak, in the hydraulic design, and in the rainfall-runoff model selected should also be considered. Although, the SCS model is limited to agricultural conditions, the approach presented herein may be applied to other flood control systems if appropriate storm runoff models are selected.     

Towards Real‐Time Continental Scale Streamflow Simulation in Continuous and Discrete Space

The National Weather Service (NWS) forecasts floods at approximately 3,600 locations across the United States (U.S.). However, the river network, as defined by the 1:100,000 scale National Hydrography Dataset‐Plus (NHDPlus) dataset, consists of 2.7 million river segments. Through the National Flood Interoperability Experiment, a continental scale streamflow simulation and forecast system was implemented and continuously operated through the summer of 2015. This system leveraged the WRF‐Hydro framework, initialized on a 3‐km grid, the Routing Application for the Parallel Computation of Discharge river routing model, operating on the NHDPlus, and real‐time atmospheric forcing to continuously forecast streamflow. Although this system produced forecasts, this paper presents a study of the three‐month nowcast to demonstrate the capacity to seamlessly predict reach scale streamflow at the continental scale. In addition, this paper evaluates the impact of reservoirs, through a case study in Texas. Validation of the uncalibrated model using observed hourly streamflow at 5,701 U.S. Geological Survey gages shows 26% demonstrate PBias ≤ |25%|, 11% demonstrate Nash‐Sutcliffe Efficiency (NSE) ≥ 0.25, and 6% demonstrate both PBias ≤ |25%| and NSE ≥ 0.25. When evaluating the impact of reservoirs, the analysis shows when reservoirs are included, NSE ≥ 0.25 for 56% of the gages downstream while NSE ≥ 0.25 for 11% when they are not. The results presented here provide a benchmark for the evolving hydrology program within the NWS and supports their efforts to develop a reach scale flood forecasting system for the country.     

ACCURACY AND PRECISION OF NRCS MODELS FOR SMALL WATERSHEDS1

ABSTRACT: In the last 30 years, the National Resource Conservation Service's TR‐55 and TR‐20 models have seen a dramatic increase in use for stormwater management purposes. This paper reviews some of the data that were originally used to develop these models and tests how well the models estimate annual series peak runoff rates for the same watersheds using longer historical data record lengths. The paper also explores differences between TR‐55 and TR‐20 peak runoff rate estimates and time of concentration methods. It was found that of the 37 watersheds tested, 25 were either over‐ or under‐predicting the actual historical watershed runoff rates by more than 30 percent. The results of this study indicate that these NRCS models should not be used to model small wooded watersheds less than 20 acres. This would be especially true if the watershed consisted of an area without a clearly defined outlet channel. This study also supports the need for regulators to allow educated hydrologists to alter pre‐packaged model parameters or results more easily than is currently permitted.     

Optimal Reorganization of NASA Earth Science Data for Enhanced Accessibility and Usability for the Hydrology Community

A long‐standing “Digital Divide” in data representation exists between the preferred way of data access by the hydrology community and the common way of data archival by earth science data centers. Typically, in hydrology, earth surface features are expressed as discrete spatial objects (e.g., watersheds), and time‐varying data are contained in associated time series. Data in earth science archives, although stored as discrete values (of satellite swath pixels or geographical grids), represent continuous spatial fields, one file per time step. This Divide has been an obstacle, specifically, between the Consortium of Universities for the Advancement of Hydrologic Science, Inc. and NASA earth science data systems. In essence, the way data are archived is conceptually orthogonal to the desired method of access. Our recent work has shown an optimal method of bridging the Divide, by enabling operational access to long‐time series (e.g., 36 years of hourly data) of selected NASA datasets. These time series, which we have termed “data rods,” are pre‐generated or generated on‐the‐fly. This optimal solution was arrived at after extensive investigations of various approaches, including one based on “data curtains.” The on‐the‐fly generation of data rods uses “data cubes,” NASA Giovanni, and parallel processing. The optimal reorganization of NASA earth science data has significantly enhanced the access to and use of the data for the hydrology user community.     

Neural Network Input Selection for Hydrological Forecasting Affected by Snowmelt1

Abstract: Snowmelt largely affects runoff in watersheds in Nordic countries. Neural networks (NN) are particularly attractive for streamflow forecasting whereas they rely at least on daily streamflow and precipitation observations. The selection of pertinent model inputs is a major concern in NNs implementation. This study investigates performance of auxiliary NN inputs that allow short‐term streamflow forecasting without resorting to a deterministic snowmelt routine. A case study is presented for the Rivière des Anglais watershed (700 km²) located in Southern Québec, Canada. Streamflow (Q), precipitations (rain R and snow S, or total P), temperature (T) and snow lying (A) observations, combined with climatic and snowmelt proxy data, including snowmelt flow (Q_SM) obtained from a deterministic model, were tested. NN implemented with antecedent Q and R produced the largest gains in performance. Introducing increments of A and T to the NNs further improved the performance. Long‐term averages, seasonal data, and Q_SM failed to improve the networks.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

AN EVALUATION OF FLOOD FREQUENCY ESTIMATES BASED ON RAINFALL/RUNOFF MODELING1

ABSTRACT: An evaluation of flood frequency estimates simulated from a rainfall/runoff model is based on (1) computation of the equivalent years of record for regional estimating equations based on 50 small stream sites in Oklahoma and (2) computation of the bias for synthetic flood estimates as compared to observed estimates at 97 small stream sites with at least 20 years of record in eight eastern states. Because of the high intercorrelation of synthetic flood estimates between watersheds, little or no regional (spatial) information may be added to the network as a result of the modeling activity. The equivalent years of record for the regional estimating equations based totally on synthetic flood discharges is shown to be considerably less than the length of rainfall record used to simulate the runoff. Furthermore, the flood estimates from the rainfall/runoff model consistently underestimate the flood discharges based on observed record, particularly for the larger floods. Depending on the way bias is computed, the synthetic estimate of the 100-year flood discharge varies from 11 to 29 percent less than the value based on observed record. In addition, the correlation between observed and synthetic flood frequency estimates at the same site is also investigated. The degree of correlation between these estimates appears to vary with recurrence interval. Unless the correlation between these two estimates is known, it is not possible to compute a weighted estimate with minimum variance.     

Historical Groundwater Trends in Northern New England and Relations with Streamflow and Climatic Variables

Water‐level trends spanning 20, 30, 40, and 50 years were tested using month‐end groundwater levels in 26, 12, 10, and 3 wells in northern New England (Maine, New Hampshire, and Vermont), respectively. Groundwater levels for 77 wells were used in interannual correlations with meteorological and hydrologic variables related to groundwater. Trends in the contemporary groundwater record (20 and 30 years) indicate increases (rises) or no substantial change in groundwater levels in all months for most wells throughout northern New England. The highest percentage of increasing 20‐year trends was in February through March, May through August, and October through November. Forty‐year trend results were mixed, whereas 50‐year trends indicated increasing groundwater levels. Whereas most monthly groundwater levels correlate strongly with the previous month's level, monthly levels also correlate strongly with monthly streamflows in the same month; correlations of levels with monthly precipitation are less frequent and weaker than those with streamflow. Groundwater levels in May through August correlate strongly with annual (water year) streamflow. Correlations of groundwater levels with streamflow data and the relative richness of 50‐ to 100‐year historical streamflow data suggest useful proxies for quantifying historical groundwater levels in light of the relatively short and fragmented groundwater data records presently available.