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.     

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.     

Modeling Variability and Trends in Pesticide Concentrations in Streams1

Abstract: A parametric regression model was developed for assessing the variability and long‐term trends in pesticide concentrations in streams. The dependent variable is the logarithm of pesticide concentration and the explanatory variables are a seasonal wave, which represents the seasonal variability of concentration in response to seasonal application rates; a streamflow anomaly, which is the deviation of concurrent daily streamflow from average conditions for the previous 30 days; and a trend, which represents long‐term (inter‐annual) changes in concentration. Application of the model to selected herbicides and insecticides in four diverse streams indicated the model is robust with respect to pesticide type, stream location, and the degree of censoring (proportion of nondetections). An automatic model fitting and selection procedure for the seasonal wave and trend components was found to perform well for the datasets analyzed. Artificial censoring scenarios were used in a Monte Carlo simulation analysis to show that the fitted trends were unbiased and the approximate p‐values were accurate for as few as 10 uncensored concentrations during a three‐year period, assuming a sampling frequency of 15 samples per year. Trend estimates for the full model were compared with a model without the streamflow anomaly and a model in which the seasonality was modeled using standard trigonometric functions, rather than seasonal application rates. Exclusion of the streamflow anomaly resulted in substantial increases in the mean‐squared error and decreases in power for detecting trends. Incorrectly modeling the seasonal structure of the concentration data resulted in substantial estimation bias and moderate increases in mean‐squared error and decreases in power.     

Assessing Retrospective National Water Model Streamflow with Respect to Droughts and Low Flows in the Colorado River Basin

Streamflow monitoring in the Colorado River Basin (CRB) is essential to ensure diverse needs are met, especially during periods of drought or low flow. Existing stream gage networks, however, provide a limited record of past and current streamflow. Modeled streamflow products with more complete spatial and temporal coverage (including the National Water Model [NWM]), have primarily focused on flooding, rather than sustained drought or low flow conditions. Objectives of this study are to (1) evaluate historical performance of the NWM streamflow estimates (particularly with respect to droughts and seasonal low flows) and (2) identify characteristics relevant to model inputs and suitability for future applications. Comparisons of retrospective flows from the NWM to observed flows from the United States Geological Survey stream gage network over 22 years in the CRB reveal a tendency for underestimating low flow frequency, locations with low flows, and the number of years with low flows. We found model performance to be more accurate for the Upper CRB and at sites with higher precipitation, snow percent, baseflow index, and elevations. Underestimation of low flows and variable model performance has important implications for future applications: inaccurate evaluations of historical low flows and droughts, and less reliable performance outside of specific watershed/stream conditions. This highlights characteristics on which to focus future model development efforts.     

Impact of Climate and Land Use Change on Streamflow and Sediment Yield in a Snow‐Dominated Semiarid Mountainous Watershed

This study investigates the impact of climate and land use change on the magnitude and timing of streamflow and sediment yield in a snow‐dominated mountainous watershed in Salt Lake County, Utah using a scenario approach and the Hydrological Simulation Program — FORTRAN model for the 2040s (year 2035–2044) and 2090s (year 2085–2094). The climate scenarios were statistically and dynamically downscaled from global climate models. Land use and land cover (LULC) changes were estimated in two ways — from a regional planning scenario and from a deterministic model. Results indicate the mean daily streamflow in the Jordan River watershed will increase by an amount ranging from 11.2% to 14.5% in the 2040s and from 6.8% to 15.3% in the 2090s. The respective increases in sediment load in the 2040s and 2090s is projected to be 6.7% and 39.7% in the canyons and about 7.4% to 14.2% in the Jordan valley. The historical 50th percentile timing of streamflow and sediment load is projected to be shifted earlier by three to four weeks by mid‐century and four to eight weeks by late‐century. The projected streamflow and sediment load results establish a nonlinear relationship with each other and are highly sensitive to projected climate change. The predicted changes in streamflow and sediment yield will have implications for water supply, flood control and stormwater management.     

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.     

Multi-scale analysis of oxygen demand trends in an urbanizing Oregon watershed, USA

Human alteration of the landscape has an extensive influence on the biogeochemical processes that drive oxygen cycling in streams. We estimated trends from the mid-1990s to 2003, using the seasonal Mann-Kendall's test, for percent saturation dissolved oxygen (DO), chemical oxygen demand (COD), total Kjeldahl nitrogen (TKN), and ammonia-nitrogen (NH(3)-N) for 12 sites in the Rock Creek watershed, northwest Oregon, USA. In order to understand the influence of landscape change, scale, and stormwater runoff management on dissolved oxygen trends, we calculated land cover change through aerial photo interpretation at full-basin, local (near sample point) basin, and 100m stream buffer scales, for the years 1994 and 2000. Significant (p < or = 0.05) trends occurred in DO (increasing at five sites), COD (decreasing at seven sites), TKN (decreasing at five sites, increasing at one site), and NH(3)-N (decreasing at one site, increasing at one site). Significant land cover change occurred in agricultural land cover (-8% for the entire basin area) and residential land cover (+10% for the entire basin area) (p < or = 0.05). Correlation results indicated that: (1) forest cover negatively influenced COD at the full basin scale and positively influences NH(3)-N at local scales, (2) residential land cover influenced oxygen demand variables at local scales, (3) agricultural land cover did not influence oxygen demand, (4) local topography negatively influenced TKN and NH(3)-N, and (5) stormwater runoff management infrastructure correlated positively with COD at the local scale. This study indicates that landscape factors influencing DO conditions for the study streams act at multiple scales, suggesting that better knowledge of scale-process interactions can guide watershed managers' decision making in order to maintain improving water quality conditions.     

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.     

Estimating upper percentiles of surface water monitoring data with sparse samples

The estimation of upper percentiles of chemical concentrations in surface water systems within sites and regions may be necessary for the assessment of potential risk to ecosystems and human health. Limited sample sizes at monitoring sites often limit the use of direct methods to estimate upper percentiles. In such cases, upper percentiles within regions within a time frame may be estimated by pooling data across sites and years, and then deriving percentile estimates from the pooled dataset. The method uses the observations resulting from either a known probability-sampling design or a sampling design treated like one because its observations come close to matching that of a probability-sample. These observations are then weighted to ensure that estimates are representative of a target population across all the sites within the region and the range of years in the time frame. This method of estimating upper percentiles of annual site concentration profiles is demonstrated using atrazine and validated using the monitoring data from both sparsely sampled and high-frequency water monitoring programs, where point and interval estimates of the 90th, 95th, and 99th pooled population percentiles are provided. This method shows that the pooled data from multiple sparse datasets can be used to provide estimates of near-peak concentrations with greater certainty, which are consistent with those generated by high-frequency sampling monitoring programs.     

Climatic Fluctuations and Forecasting of Streamflow in the Lower Colorado River Basin1

Abstract:  Water‐resource managers need to forecast streamflow in the Lower Colorado River Basin to plan for water‐resource projects and to operate reservoirs for water supply. Statistical forecasts of streamflow based on historical records of streamflow can be useful, but statistical assumptions, such as stationarity of flows, need to be evaluated. This study evaluated the relation between climatic fluctuations and stationarity and developed regression equations to forecast streamflow by using climatic fluctuations as explanatory variables. Climatic fluctuations were represented by the Atlantic Multidecadal Oscillation (AMO), Pacific Decadal Oscillation (PDO), and Southern Oscillation Index (SOI). Historical streamflow within the 25‐ to 30‐year positive or negative phases of AMO or PDO was generally stationary. Monotonic trends in annual mean flows were tested at the 21 sites evaluated in this study; 76% of the sites had no significant trends within phases of AMO and 86% of the sites had no significant trends within phases of PDO. As climatic phases shifted in signs, however, many sites had nonstationary flows; 67% of the sites had significant changes in annual mean flow as AMO shifted in signs. The regression equations developed in this study to forecast streamflow incorporate these shifts in climate and streamflow, thus that source of nonstationarity is accounted for. The R² value of regression equations that forecast individual years of annual flow for the central part of the study area ranged from 0.28 to 0.49 and averaged 0.39. AMO was the most significant variable, and a combination of indices from both the Atlantic and Pacific Oceans explained much more variation in flows than only the Pacific Ocean indices. The average R² value for equations with PDO and SOI was 0.15.     

Streamflow Responses to Climate Change: Analysis of Hydrologic Indicators in a New York City Water Supply Watershed

Recent works have indicated that climate change in the northeastern United States is already being observed in the form of shorter winters, higher annual average air temperature, and more frequent extreme heat and precipitation events. These changes could have profound effects on aquatic ecosystems, and the implications of such changes are less understood. The objective of this study was to examine how future changes in precipitation and temperature translate into changes in streamflow using a physically based semidistributed model, and subsequently how changes in streamflow could potentially impact stream ecology. Streamflow parameters were examined in a New York City water supply watershed for changes from model‐simulated baseline conditions to future climate scenarios (2081‐2100) for ecologically relevant factors of streamflow using the Indicators of Hydrologic Alterations tool. Results indicate that earlier snowmelt and reduced snowpack advance the timing and increase the magnitude of discharge in the winter and early spring (November‐March) and greatly decrease monthly streamflow later in the spring in April. Both the rise and fall rates of the hydrograph will increase resulting in increased flashiness and flow reversals primarily due to increased pulses during winter seasons. These shifts in timing of peak flows, changes in seasonal flow regimes, and changes in the magnitudes of low flow can all influence aquatic organisms and have the potential to impact stream ecology.     

Urban Streamflow Response to Imported Water and Water Conservation Policies in Los Angeles,California

Los Angeles has a long history of importing water; however, drought, climate change, and environmental mitigation have forced the City to focus on developing more local water sources (target of 50% local supply by 2035). This study aims to improve understanding of water cycling in Los Angeles, including the impacts of imported water and water conservation policies. We evaluate the influence of local water restrictions on discharge records for 12 years in the Ballona Creek (urban) and Topanga Creek (natural) watersheds. Results show imported water has significantly altered the timing and volume of streamflow in the urban Ballona watershed, resulting in runoff ratios above one (more streamflow than precipitation). Further analysis comparing pre‐ vs. during‐mandatory water conservation periods shows there is a significant decrease in dry season streamflow during‐conservation in Ballona, indicating that prior to conservation efforts, heavy irrigation and other outdoor water use practices were contributing to streamflow. The difference between summer streamflow pre‐ vs. during‐conservation is enough to serve 160,000 customers in Los Angeles. If Los Angeles returns to more watering days, educating the public on proper irrigation rates is critical for ensuring efficient irrigation and conserving water; however, if water restrictions remain in place, the City must take the new flow volumes into account for complying with water quality standards in the region.     

长山河水质评价及变化趋势分析

以2015 ～2019年长山河4个水质监测站的监测数据为基础,采用综合水质标识指数法对其历年水质进行评价,利用季节性kendall检验法对各监测断面的高锰酸盐指数,氨氮,总磷等3个监测指标变化趋势进行分析.结果 表明:长山河各断面的综合水质均达到或优于水功能区目标要求,水质污染总体呈下降趋势,其中长山河大桥和长山闸一号...     

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.     

Estimation of Daily Streamflow of Southeastern Coastal Plain Watersheds by Combining Estimated Magnitude and Sequence

Commonly used methods to predict streamflow at ungauged watersheds implicitly predict streamflow magnitude and temporal sequence concurrently. An alternative approach that has not been fully explored is the conceptualization of streamflow as a composite of two separable components of magnitude and sequence, where each component is estimated separately and then combined. Magnitude is modeled using the flow duration curve (FDC), whereas sequence is modeled by transferring streamflow sequence of gauged watershed(s). This study tests the applicability of the approach on watersheds ranging in size from about 25‐7,226 km² in Southeastern Coastal Plain (U.S.) with substantial surface storage of wetlands. A 19‐point regionalized FDC is developed to estimate streamflow magnitude using the three most selected variables (drainage area, hydrologic soil index, and maximum 24‐h precipitation with a recurrence interval of 100 years) by a greedy‐heuristic search process. The results of validation on four watersheds (Trent River, North Carolina: 02092500; Satilla River, Georgia: 02226500; Black River, South Carolina: 02136000; and Coosawhatchie River, South Carolina: 02176500) yielded Nash‐Sutcliffe efficiency values of 0.86‐0.98 for the predicted magnitude and 0.09‐0.84 for the predicted daily streamflow over a simulation period of 1960‐2010. The prediction accuracy of the method on two headwater watersheds at Santee Experimental Forest in coastal South Carolina was weak, but comparable to simulations by MIKE‐SHE.     

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.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.     

Differences in Dissolved Organic Carbon and Nitrogen Responses to Storm‐Event and Ground‐Water Conditions in a Forested,Glaciated Watershed in Western New York1

Abstract: Differences in the storm‐event responses of dissolved organic carbon (DOC) and nitrogen (DON) in streamflow and ground water were evaluated for a glaciated forested watershed in western New York. Eight to ten storm events with varying rainfall amounts, intensities, and antecedent moisture conditions were studied for three catchments (1.6, 3.4, and 696 ha) over a three‐year period (2003‐2005). Concentrations of DOC in streamflow exiting the catchments were significantly higher for storm events following a dry period, whereas no similar response was observed for DON. Highest DON concentrations in streamflow were typically associated with storm events following wet antecedent moisture conditions. In addition to antecedent moisture conditions, DOC concentrations were also positively correlated with precipitation amounts, while DON did not reveal a consistent pattern. Streamwater and ground‐water concentrations of DOC during storm events were also strongly correlated with riparian ground‐water depths but a similar relationship was not observed for DON. Ground‐water DON concentrations were also more variable than DOC. We hypothesized that the differences in DOC and DON responses stemmed from the differences in catchment sources of these solutes. This study suggests that while DOC and DON are intrinsically linked as components of dissolved organic matter, their dynamics and exports from watersheds may be regulated by a different set of mechanisms and factors. Identifying these differences is critical for developing more reliable and robust models for transport of dissolved organic matter.     

Characterizing Uncertainty in Daily Streamflow Estimates at Ungauged Locations for the Massachusetts Sustainable Yield Estimator

Hydrologic characterization at ungauged locations is one of the quintessential challenges of hydrology. Beyond simulation of historical streamflows, it is similarly important to characterize the level of uncertainty in hydrologic estimates. In tandem with updates to Massachusetts Sustainable Yield Estimator, this work explores the application of global uncertainty estimates to daily streamflow simulations. Expanding on a method developed for deterministic modeling, this approach produces confidence intervals on daily streamflow developed through nonlinear spatial interpolation of daily streamflow using flow duration curves; the 95% confidence is examined. Archived cross‐validations of daily streamflows from 66 watersheds in and around Massachusetts are used to evaluate an approach to uncertainty characterization. Neighboring sites are treated as ungauged, producing relative errors that can be resampled and applied to target sites. The method, with some modification, is found to provide appropriately narrow confidence intervals that contain 95% of the observed streamflows in cross‐validation. Further characterizing uncertainty, multiday means of daily streamflow are evaluated. Working through cross‐validation in Massachusetts, two‐ to three‐month averages of daily streamflow show the best performance. These two approaches to uncertainty characterization inform how streamflow simulation produced for prediction in ungauged basins can be used for water resources management.     

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.