Seasonal and Flood‐Induced Variations in Groundwater–Surface Water Exchange in a Northern Coldwater Fishery

Groundwater upwelling is important to coldwater fisheries survival. This study used stable isotopes to identify upwelling zones within a watershed, then combined isotope analyses with reach‐scale monitoring to measure surface water–groundwater exchange over time. Research focused on Amity Creek, Minnesota, a basin that exemplifies conditions limiting coldwater species survival along Lake Superior's North Shore where shallow bedrock limits groundwater capacity, lowering baseflows and increasing temperatures. Groundwater‐fed reaches were identified through synoptic isotope sampling, with results highlighting the importance of isolated shallow surficial aquifers (glacially derived sands and gravels) for providing cold baseflow waters. In an alluvial reach, monitoring well results show groundwater was stored in two reservoirs: one that reacts quickly to changes in stream levels, and one that remained isotopically isolated under most flow conditions, but which helps sustain summer baseflows for weeks to months. A 500‐year flood demonstrated the capacity of high‐flow events to alter surface water–groundwater connectivity. The previously isolated reservoir was exchanged or mixed during the flood pulse, while incision lowered the water table for years. The results here provide insight for streams that lack substantial groundwater inputs yet maintain coldwater species at risk in a warming climate and an approach for managers seeking to protect cold baseflow sources.     

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.     

Streamflow Response to Climate as Influenced by Geology and Elevation1

Mayer, Timothy D. and Seth W. Naman, 2011. Streamflow Response to Climate as Influenced by Geology and Elevation. Journal of the American Water Resources Association (JAWRA) 47(4):724‐738. DOI: 10.1111/j.1752‐1688.2011.00537.x Abstract: This study examines the regional streamflow response in 25 predominately unregulated basins to warmer winter temperatures and snowpack reductions over the last half century in the Klamath Basin of California and Oregon. Geologic controls of streamflow in the region result in two general stream types: surface‐dominated and groundwater‐dominated basins. Surface‐dominated basins were further differentiated into rain basins and snowmelt basins on the basis of elevation and timing of winter runoff. Streamflow characteristics and response to climate vary with stream type, as discussed in the study. Warmer winter temperatures and snowpack reductions have caused significantly earlier runoff peaks in both snowmelt and groundwater basins in the region. In the groundwater basins, the streamflow response to changes in snowpack is smoothed and delayed and the effects are extended longer in the summer. Our results indicate that absolute decreases in July‐September base flows are significantly greater, by an order of magnitude, in groundwater basins compared to surface‐dominated basins. The declines are important because groundwater basins sustain Upper Klamath Lake inflows and mainstem river flows during the typically dry summers of the area. Upper Klamath Lake April‐September net inflows have decreased an estimated 16% or 84 thousand acre‐feet (103.6 Mm³) since 1961, with the summer months showing proportionately more decline. These changes will exacerbate water supply problems for agriculture and natural resources in the region.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

Controls on Temperature in Salmonid‐Bearing Headwater Streams in Two Common Hydrogeologic Settings,Kenai Peninsula,Alaska

Headwater streams are the most numerous in terms of both number and length in the conterminous United States and play important roles as spawning and rearing grounds for numerous species of anadromous fish. Stream temperature is a controlling variable for many physical, chemical, and biological processes and plays a critical role in the overall health and integrity of a stream. We investigated the controls on stream temperature in salmon‐bearing headwater streams in two common hydrogeologic settings on the Kenai Peninsula, Alaska: (1) drainage‐ways, which are low‐gradient streams that flow through broad valleys; and (2) discharge‐slopes, which are high gradient streams that flow through narrow valleys. We hypothesize local geomorphology strongly influences surface‐water and groundwater interactions, which control streamflow at the network scale and stream temperatures at the reach scale. The results of this study showed significant differences in stream temperatures between the two hydrogeologic settings. Observed stream temperatures were higher in drainage‐way sites than in discharge‐slope sites, and showed strong correlations as a continuous function with the calculated topographic metric flow‐weighted slope. Additionally, modeling results indicated the potential for groundwater discharge to moderate stream temperature is not equal between the two hydrogeologic settings, with groundwater having a greater moderating effect on stream temperature at the drainage‐way sites.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.     

A System Dynamics Model for Conjunctive Management of Water Resources in the Snake River Basin

The Pacific Northwest is expected to witness changes in temperature and precipitation due to climate change. In this study, we enhance the Snake River Planning Model (SRPM) by modeling the feedback loop between incidental recharge and surface water supply resulting from surface water and groundwater extraction for irrigation and provide a case study involving climate change impacts and management scenarios. The new System Dynamics‐Snake River Planning Model (SD‐SRPM) is calibrated to flow at Box Canyon Springs located along a major outlet of the East Snake Plain Aquifer. A calibration of the model to flow at Box Canyon Springs, based on historic diversions (1950‐1995) resulted in an r² value of 0.74 and a validation (1996‐2005) r² value of 0.60. After adding irrigation entities to the model an r² value of 0.91, 0.88, and 0.87 were maintained for modeled vs. observed (1991‐2005) end‐of‐month reservoir content in Jackson Lake, Palisades, and American Falls, the three largest irrigation reservoirs in the system. The scenarios that compared the impacts of climate change were based on ensemble mean precipitation change scenarios and estimated changes to crop evapotranspiration (ET). Increased ET, despite increased precipitation, generally increased surface water shortages and discharge of springs. This study highlights the need to develop and implement models that integrate the human‐natural system to understand the impacts of climate change.     

Nitrogen Subsidies from Hillslope Alder Stands to Streamside Wetlands and Headwater Streams,Kenai Peninsula,Alaska

We examined nitrogen transport and wetland primary production along hydrologic flow paths that link nitrogen‐fixing alder (Alnus spp.) stands to downslope wetlands and streams in the Kenai Lowlands, Alaska. We expected that nitrate concentrations in surface water and groundwater would be higher on flow paths below alder. We further expected that nitrate concentrations would be higher in surface water and groundwater at the base of short flow paths with alder and that streamside wetlands at the base of alder‐near flow paths would be less nitrogen limited than wetlands at the base of long flow paths with alder. Our results showed that groundwater nitrate‐N concentrations were significantly higher at alder‐near sites than at no‐alder sites, but did not differ significantly between alder‐far sites and no‐alder sites or between alder‐far sites and alder‐near sites. A survey of ¹⁵N stable isotope signatures in soils and foliage in alder‐near and no‐alder flow paths indicated the alder‐derived nitrogen evident in soils below alder is quickly integrated downslope. Additionally, there was a significant difference in the relative increase in plant biomass after nitrogen fertilization, with the greatest increase occurring in the no‐alder sites. This study demonstrates that streamside wetlands and streams are connected to the surrounding landscapes through hydrologic flow paths, and flow paths with alder stands are potential “hot spots” for nitrogen subsidies at the hillslope scale.     

Sensitivity of Thermal Habitat of a Trout Stream to Potential Climate Change,Wisconsin, United States1

Abstract: Airborne thermal remote sensing from four flights on a single day from a single‐engine airplane was used to collect thermal infrared data of a 10.47‐km reach of the upper East Branch Pecatonica River in southwest Wisconsin. The study uses a one‐dimensional stream temperature model calibrated with the longitudinal profiles of stream temperature created from the four thermal imaging flights and validated with three days of continuous stream temperature data from instream data loggers on the days surrounding the thermal remote‐sensing campaign. Model simulations were used to quantify the sensitivity of stream thermal habitat to increases in air and groundwater temperature and changes in base flow. The simulations indicate that stream temperatures may reach critical maximum thresholds for brook trout (Salvelinus fontinalis) and brown trout (Salmo trutta) mortality, particularly if both air temperature increases and base flow declines. The approach demonstrates that thermal infrared data can greatly assist stream temperature model validation due to its high spatial resolution, and that this spatially continuous stream temperature data can be used to pinpoint spatial heterogeneity in groundwater inflow to streams. With this spatially distributed data on thermal heterogeneity and base‐flow accretion, stream temperature models considering various climate change scenarios are able to identify thermal refugia that will be critical for fisheries management under a changing climate.     

Estimating Water Budgets and Vertical Leakages for Karst Lakes in North‐Central Florida (United States) Via Hydrological Modeling1

Lin, Zhulu, 2011. Estimating Water Budgets and Vertical Leakages for Karst Lakes in North‐Central Florida (United States) Via Hydrological Modeling. Journal of the American Water Resources Association (JAWRA) 1‐16. DOI: 10.1111/j.1752‐1688.2010.00513.x Abstract: Newnans, Lochloosa, and Orange Lakes are closely hydrologically connected karst lakes located in north‐central Florida, United States. The complex karst hydrology in this region poses a great challenge to the hydrological modeling that is essential to the development of Total Maximum Daily Loads for these lakes. We used a Hydrological Simulation Program – Fortran model coupled with the parallel Parameter ESTimation model calibration and uncertainty analysis software to estimate effectively the hydrological interactions between the lakes and the underlying upper Floridan aquifer and the water budgets for these three lakes. The net results of the lake‐groundwater interactions in Newnans and Orange Lakes are that both lakes recharge the underlying upper Floridan aquifer, with the recharge rate of the latter one magnitude greater than that of the former. However, for Lochloosa Lake, the net lake‐groundwater interaction is that the lake gains water from groundwater in a significant amount, approximately 40% of its total terrestrial water input. The annual average vertical leakages estimated for Newnans, Lochloosa, and Orange Lakes are 6.0 × 10⁶, ?8.9 × 10⁶, and 44.4 × 10⁶ m³, respectively. The average vertical hydraulic conductance (K_v/b) of the units between a lake bottom and the underlying upper Floridan aquifer in this region are also estimated to be from 1.26 × 10^?4 to 1.01 × 10^?3 day^?1.     

Salinity as a Limiting Factor for Biological Condition in Mining‐Influenced Central Appalachian Headwater Streams

Recent studies have found that Appalachian coal mining causes increased surface water salinity, and that benthic macroinvertebrate communities in salinized mining‐influenced streams differ from communities in streams draining unmined areas. Understanding the role of salinity in shaping these communities is challenging because such streams are often influenced by a variety of stressors in addition to salinity. We characterized associations of salinity with biotic condition while isolating salinity from other stressors through rigorous site selection. We used a multimetric index of biotic condition to characterize benthic macroinvertebrate communities in headwater streams in the Central Appalachian Ecoregion of Virginia across a gradient of sulfate‐dominated salinity. We found strong negative seasonal correlations between biotic condition and three salinity measures (specific conductance, total dissolved solids, and SO₄^2? concentration). We found no evidence to suggest stressors other than salinity as significant influences on biotic condition in these streams. Our results confirm negative associations of salinity with benthic macroinvertebrate community condition, as observed in other studies. Thus, our findings demonstrate that elevated salinity is an important limiting factor for biological condition in Central Appalachian headwater streams.     

Groundwater Level Trends and Drivers in Two Northern New England Glacial Aquifers

We evaluated long‐term trends and predictors of groundwater levels by month from two well‐studied northern New England forested headwater glacial aquifers: Sleepers River, Vermont, 44 wells, 1992‐2013; and Hubbard Brook, New Hampshire, 15 wells, 1979‐2004. Based on Kendall Tau tests with Sen slope determination, a surprising number of well‐month combinations had negative trends (decreasing water levels) over the respective periods. Sleepers River had slightly more positive than negative trends overall, but among the significant trends (p < 0.1), negative trends dominated 67 to 40. At Hubbard Brook, negative trends outnumbered positive trends by a nearly 2:1 margin and all seven of the significant trends were negative. The negative trends occurred despite generally increasing trends in monthly and annual precipitation. This counterintuitive pattern may be a result of increased precipitation intensity causing higher runoff at the expense of recharge, such that evapotranspiration demand draws down groundwater storage. We evaluated predictors of month‐end water levels by multiple regression of 18 variables related to climate, streamflow, snowpack, and prior month water level. Monthly flow and prior month water level were the two strongest predictors for most months at both sites. The predictive power and ready availability of streamflow data can be exploited as a proxy to extend limited groundwater level records over longer time periods.     

Impact of Urban Growth and High Residential Irrigation on Streamflow and Groundwater Levels in a Peri‐Urban Semiarid Catchment

The impact of urbanization on groundwater is not simple to understand, as it depends on a variety of factors such as climate, hydrogeology, water management practices, and infrastructure. In semiarid landscapes, the urbanization processes can involve high water consumptions and irrigation increases, which in turn may contribute to groundwater recharge. We assessed the hydrological impacts of urbanization and irrigation rates in an Andean peri‐urban catchment located in Chile, in a semiarid climate. For this purpose, we built and validated a coupled surface–groundwater model that allows the verification of a strong stream–aquifer interaction in areas with shallow groundwater, higher than some sewers and portions of the stream. Moreover, we also identified a significant local recharge associated with pipe leaks and inefficient urban irrigation. From the evaluation of different future scenarios, we found a sustainable water conservation scenario will decrease the current groundwater levels, while the median flow reduces from 408 to 389 L/s, and the low flow (Q_95%) from 43 to 22L/s. Overall, our results show the relevance of integrating the modeling of surface and subsurface water resources at different spatial and temporal scales, when assessing the effect of urban development and the suitability of urban water practices.     

Large‐Scale Fine‐Resolution Hydrological Modeling Using Parameter Regionalization in the Missouri River Basin

This study simulated crop and water yields in the Missouri River Basin (MRB; 1,371,000 km²), one of the largest river basins in the United States, using the Soil and Water Assessment Tool (SWAT) at a fine resolution of 12‐digit Hydrological Unit Codes (HUCs) using the regionalization calibration approach. Very few studies have simulated the entire MRB, and those that have developed were at a coarser resolution of 8‐digit HUCs and were minimally calibrated. The MRB was first divided into three subbasins and was further divided into eleven regions. A “head watershed” was selected in each region and was calibrated for crop and water yields. The parameters from the calibrated head watershed were extrapolated to other subwatersheds in the region to complete comprehensive spatial calibration. The simulated crop yields at the head watersheds were in close agreement with observed crop yields. Spatial validation of the aggregated crop yields resulted in reasonable predictions for all crops except dryland corn in a few regions. Simulated and observed water yields in head watersheds and also in the validation locations were in close agreement in naturalized streams and poor agreement in streams with high groundwater‐surface water interactions and/or reservoirs found upstream of the gauges. Overall, the SWAT model was able to reasonably capture the hydrological and crop growth dynamics occurring in the basin despite some limitations.     

Evaluating the Slope‐Area Method to Accurately Identify Stream Channel Heads in Three Physiographic Regions

Estimation of stream channel heads is an important task since ephemeral channels play a significant role in the transport of sediment and materials to perennial streams. The slope‐area method utilizes digital elevation model (DEM) and related information to develop slope‐area threshold relationships used to estimate the position of channel heads in the watershed. A total of 162 stream channel heads were mapped across the three physiographic regions of Alabama, including the Southwestern Appalachians (51), Piedmont/Ridge and Valley (61), and Coastal Plains (51). Using Geographic Information System and DEM, the local slope and drainage area for each mapped channel head was calculated and region‐specific models were developed and evaluated. Results demonstrated the local slope and drainage area had an inverse and strong correlation in the Piedmont/Ridge and Valley region (r² = 0.71) and the Southwestern Appalachian region (r² = 0.61). Among three physiographic regions, the weakest correlation was observed in the Coastal Plain region (r² = 0.45). By comparing the locations of modeled channel heads to those located in the field, calculated reliability and sensitivity indices indicated model accuracy and reliance were weak to moderate. However, the slope‐area method helped define the upstream boundaries of a more detailed channel network than that derived from the 1:24,000‐scale National Hydrography Dataset, which is commonly used for planning and regulatory purposes.     

Two‐Stage Ditch Floodplains Enhance N‐Removal Capacity and Reduce Turbidity and Dissolved P in Agricultural Streams

Two‐stage ditches represent an emerging management strategy in artificially drained agricultural landscapes that mimics natural floodplains and has the potential to improve water quality. We assessed the potential for the two‐stage ditch to reduce sediment and nutrient export by measuring water column turbidity, nitrate (NO₃^?), ammonium (NH₄⁺), and soluble reactive phosphorus (SRP) concentrations, and denitrification rates. During 2009‐2010, we compared reaches with two‐stage floodplains to upstream reaches with conventional trapezoid design in six agricultural streams. At base flow, these short two‐stage reaches (<600 m) reduced SRP concentrations by 3‐53%, but did not significantly reduce NO₃^? concentrations due to very high NO₃^? loads. The two‐stage also decreased turbidity by 15‐82%, suggesting reduced suspended sediment export during floodplain inundation. Reach‐scale N‐removal increased 3‐24 fold during inundation due to increased bioreactive surface area with high floodplain denitrification rates. Inundation frequency varied with bench height, with lower benches being flooded more frequently, resulting in higher annual N‐removal. We also found both soil organic matter and denitrification rates were higher on older floodplains. Finally, influence of the two‐stage varied among streams and years due to variation in stream discharge, nutrient loads, and denitrification rates, which should be considered during implementation to optimize potential water quality benefits.     

Modeled Summer Background Concentration of Nutrients and Suspended Sediment in the Mid‐Continent (USA) Great Rivers1

Angradi, Ted R., David W. Bolgrien, Matt A. Starry, and Brian H. Hill, 2012. Modeled Summer Background Concentration of Nutrients and Suspended Sediment in the Mid‐Continent (USA) Great Rivers. Journal of the American Water Resources Association (JAWRA) 48(5): 1054‐1070. DOI: 10.1111/j.1752‐1688.2012.00669.x Abstract: We used regression models to predict summer background concentration of total nitrogen (N), total phosphorus (P), and total suspended solids (TSS), in the mid‐continent great rivers: the Upper Mississippi, the Lower Missouri, and the Ohio. From multiple linear regressions of water quality indicators with land use and other stressor variables, we determined the concentration of the indicators when the predictor variables were all set to zero — the y‐intercept. Except for total P on the Upper Mississippi River, we could predict background concentration using regression models. Predicted background concentration of total N was about the same on the Upper Mississippi and Lower Missouri Rivers (430 μg l^?1), which was lower than percentile‐based values, but was similar to concentrations derived from the response of sestonic chlorophyll a to great river total N concentration. Background concentration of total P on the Lower Missouri (65 μg l^?1) was also lower than published and percentile‐based concentrations. Background TSS concentration was higher on the Lower Missouri (40 mg l^?1) than the other rivers. Background TSS concentration on the Upper Mississippi (16 mg l^?1) was below a threshold (30 mg l^?1) designed to protect aquatic vegetation. Our model‐predicted concentrations for the great rivers are an attempt to estimate background concentrations for water quality indicators independent from thresholds based on percentiles or derived from stressor‐response relationships.     

A Hydrologic Drying Bias in Water‐Resource Impact Analyses of Anthropogenic Climate Change

For water‐resource planning, sensitivity of freshwater availability to anthropogenic climate change (ACC) often is analyzed with “offline” hydrologic models that use precipitation and potential evapotranspiration (E_p) as inputs. Because E_p is not a climate‐model output, an intermediary model of E_p must be introduced to connect the climate model to the hydrologic model. Several E_p methods are used. The suitability of each can be assessed by noting a credible E_p method for offline analyses should be able to reproduce climate models’ ACC‐driven changes in actual evapotranspiration in regions and seasons of negligible water stress (E_w). We quantified this ability for seven commonly used E_p methods and for a simple proportionality with available energy (“energy‐only” method). With the exception of the energy‐only method, all methods tend to overestimate substantially the increase in E_p associated with ACC. In an offline hydrologic model, the E_p‐change biases produce excessive increases in actual evapotranspiration (E), whether the system experiences water stress or not, and thence strong negative biases in runoff change, as compared to hydrologic fluxes in the driving climate models. The runoff biases are comparable in magnitude to the ACC‐induced runoff changes themselves. These results suggest future hydrologic drying (wetting) trends likely are being systematically and substantially overestimated (underestimated) in many water‐resource impact analyses.     

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.