Modeling Streamflow and Water Quality Sensitivity to Climate Change and Urban Development in 20 U.S. Watersheds

Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.     

Adaptation of Climate Model Projections of Streamflow to Account for Upstream Anthropogenic Impairments

A statistical procedure is developed to adjust natural streamflows simulated by dynamical models in downstream reaches, to account for anthropogenic impairments to flow that are not considered in the model. The resulting normalized downstream flows are appropriate for use in assessments of future anthropogenically impaired flows in downstream reaches. The normalization is applied to assess the potential effects of climate change on future water availability on the Rio Grande at a gage just above the major storage reservoir on the river. Model‐simulated streamflow values were normalized using a statistical parameterization based on two constants that relate observed and simulated flows over a 50‐year historical baseline period (1964–2013). The first normalization constant is a ratio of the means, and the second constant is the ratio of interannual standard deviations between annual gaged and simulated flows. This procedure forces the gaged and simulated flows to have the same mean and variance over the baseline period. The normalization constants can be kept fixed for future flows, which effectively assumes that upstream water management does not change in the future, or projected management changes can be parameterized by adjusting the constants. At the gage considered in this study, the effect of the normalization is to reduce simulated historical flow values by an average of 72% over an ensemble of simulations, indicative of the large fraction of natural flow diverted from the river upstream from the gage. A weak tendency for declining flow emerges upon averaging over a large ensemble, with tremendous variability among the simulations. By the end of the 21st Century the higher‐emission scenarios show more pronounced declines in streamflow.     

Decoupling Streamflow Responses to Climate Variability and Land Use/Cover Changes in a Watershed in Northern China

Restored annual streamflow (Q_r) and measured daily streamflow of the Chaohe watershed located in northern China and associated long‐term climate and land use/cover data were used to explore the effects of land use/cover change and climate variability on the streamflow during 1961‐2009. There were no significant changes in annual precipitation (P) and potential evapotranspiration, whereas Q_r decreased significantly by 0.81 mm/yr (p < 0.001) over the study period with a change point in 1999. We used 1961‐1998 as the baseline period (BP) and 1999‐2009 the change period (CP). The mean Q_r during the CP decreased by 39.4 mm compared with that in the BP. From 1979 to 2009, the grassland area declined by 69.6%, and the forest and shrublands increased by 105.4 and 73.1%, respectively. The land use/cover change and climate variability contributed for 58.4 and 41.6% reduction in mean annual Q_r, respectively. Compared with the BP, median and high flows in the CP decreased by 38.8 and up to 75.5%, respectively. The study concludes that large‐scale ecological restoration and watershed management in northern China has greatly decreased water yield and reduced high flows due to the improved land cover by afforestation leading to higher water loss through evapotranspiration. At a large watershed scale, land use/cover change could play as much of an important role as climate variability on water resources.     

Use of Hydrologic Landscape Classification to Diagnose Streamflow Predictability in Oregon

We implement a spatially lumped hydrologic model to predict daily streamflow at 88 catchments within the state of Oregon and analyze its performance using the Oregon Hydrologic Landscape (OHL) classification. OHL is used to identify the physio‐climatic conditions that favor high (or low) streamflow predictability. High prediction catchments (Nash‐Sutcliffe efficiency of (NS) > 0.75) are mainly classified as rain dominated with very wet climate, low aquifer permeability, and low to medium soil permeability. Most of them are located west of the Cascade Mountain Range. Conversely, most low prediction catchments (NS < 0.6) are classified as snow‐dominated with high aquifer permeability and medium to high soil permeability. They are mainly located in the volcano‐influenced High Cascades region. Using a subset of 36 catchments, we further test if class‐specific model parameters can be developed to predict at ungauged catchments. In most catchments, OHL class‐specific parameters provide predictions that are on par with individually calibrated parameters (NS decline < 10%). However, large NS declines are observed in OHL classes where predictability is not high enough. Results suggest higher uncertainty in rain‐to‐snow transition of precipitation phase and external gains/losses of deep groundwater are major factors for low prediction in Oregon. Moreover, regionalized estimation of model parameters is more useful in regions where conditions favor good streamflow predictability.     

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.     

Mid‐Range Streamflow Forecasts Based on Climate Modeling – Statistical Correction and Evaluation1

Abstract: Mid‐range streamflow predictions are extremely important for managing water resources. The ability to provide mid‐range (three to six months) streamflow forecasts enables considerable improvements in water resources system operations. The skill and economic value of such forecasts are of great interest. In this research, output from a general circulation model (GCM) is used to generate hydrologic input for mid‐range streamflow forecasts. Statistical procedures including: (1) transformation, (2) correction, (3) observation of ensemble average, (4) improvement of forecast, and (5) forecast skill test are conducted to minimize the error associated with different spatial resolution between the large‐scale GCM and the finer‐scale hydrologic model and to improve forecast skills. The accuracy of a streamflow forecast generated using a hydrologic model forced with GCM output for the basin was evaluated by forecast skill scores associated with the set of streamflow forecast values in a categorical forecast. Despite the generally low forecast skill score exhibited by the climate forecasting approach, precipitation forecast skill clearly improves when a conditional forecast is performed during the East Asia summer monsoon, June through August.     

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

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

Watershed Evapotranspiration Increased due to Changes in Vegetation Composition and Structure Under a Subtropical Climate1

Abstract: Natural forests in southern China have been severely logged due to high human demand for timber, food, and fuels during the past century, but are recovering in the past decade. The objective of this study was to investigate how vegetation cover changes in composition and structure affected the water budgets of a 9.6‐km² Dakeng watershed located in a humid subtropical mountainous region in southern China. We analyzed 27 years (i.e., 1967‐1993) of streamflow and climate data and associated vegetation cover change in the watershed. Land use/land cover census and Normalized Difference of Vegetation Index (NDVI) data derived from remote sensing were used to construct historic land cover change patterns. We found that over the period of record, annual streamflow (Q) and runoff/precipitation ratio did not change significantly, nor did the climatic variables, including air temperature, Hamon’s potential evapotranspiration (ET), pan evaporation, sunshine hours, and radiation. However, annual ET estimated as the differences between P and Q showed a statistically significant increasing trend. Overall, the NDVI of the watershed had a significant increasing trend in the peak spring growing season. This study concluded that watershed ecosystem ET increased as the vegetation cover shifted from low stock forests to shrub and grasslands that had higher ET rates. A conceptual model was developed for the study watershed to describe the vegetation cover‐streamflow relationships during a 50‐year time frame. This paper highlighted the importance of eco‐physiologically based studies in understanding transitory, nonstationary effects of deforestation or forestation on watershed water balances.     

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.     

Changes in Patterns of Streamflow From Unregulated Watersheds in Idaho,Western Wyoming,and Northern Nevada1

Clark, Gregory M., 2010. Changes in Patterns of Streamflow From Unregulated Watersheds in Idaho, Western Wyoming, and Northern Nevada. Journal of the American Water Resources Association (JAWRA) 46(3):486-497. DOI: 10.1111/j.1752-1688.2009.00416.x Abstract: Recent studies have identified a pattern of earlier spring runoff across much of North America. Earlier spring runoff potentially poses numerous problems, including increased risk of flooding and reduced summer water supply for irrigation, power generation, and migratory fish passage. To identify changing runoff patterns in Idaho streams, streamflow records were analyzed for 26 U.S. Geological Survey gaging stations in Idaho, western Wyoming, and northern Nevada, each with a minimum of 41 years of record. The 26 stations are located on 23 unregulated and relatively pristine streams that drain areas ranging from 28 to >35,000 km². Four runoff parameters were trend tested at each station for both the period of historical record and from 1967 through 2007. Parameters tested were annual mean streamflow, annual minimum daily streamflow, and the dates of the 25th and 50th percentiles of the annual total streamflow. Results of a nonparametric Mann-Kendall trend test revealed a trend toward lower annual mean and annual minimum streamflows at a majority of the stations, as well as a trend toward earlier snowmelt runoff. Significant downward trends over the period of historical record were most prevalent for the annual minimum streamflow (12 stations) and the 50th percentile of streamflow (11 stations). At most stations, trends were more pronounced during the period from 1967 through 2007. A regional Kendall test for water years 1967 through 2007 revealed significant regional trends in the percent change in the annual mean and annual minimum streamflows (0.67% less per year and 0.62% less per year, respectively), the 25th percentile of streamflow (12.3 days earlier), and the 50th percentile of streamflow (11.5 days earlier).     

Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow1

Miller, W. Paul and Thomas C. Piechota, 2011. Trends in Western U.S. Snowpack and Related Upper Colorado River Basin Streamflow. Journal of the American Water Resources Association (JAWRA) 47(6):1197–1210. DOI: 10.1111/j.1752‐1688.2011.00565.x Abstract: Water resource managers in the Western United States (U.S.) are currently faced with the challenge of adapting to unprecedented drought and uncertain impacts of climate change. Recent research has indicated increasing regional temperature and changes to precipitation and streamflow characteristics throughout the Western U.S. As such, there is increased uncertainty in hydroclimatological forecasts, which impact reservoir operations and water availability throughout the Western U.S., particularly in the Colorado River Basin. Previous research by the authors hypothesized a change in the character of precipitation (i.e., the frequency and amount of rainfall and snowfall events) throughout the Colorado River Basin. In the current study, 398 snowpack telemetry stations were investigated for trends in cumulative precipitation, snow water equivalent, and precipitation events. Observations of snow water equivalent characteristics were compared to observations in streamflow characteristics. Results indicate that the timing of the last day of the snow season corresponds well to the volume of runoff observed over the traditional peak flow season (April through July); conversely, the timing of the first day of the snow season does not correspond well to the volume of runoff observed over the peak flow season. This is significant to water resource managers and river forecasters, as snowpack characteristics may be indicative of a productive or unproductive runoff season.     

Long‐Term Streamflow Response to Climatic Variability in the Loess Plateau,China1

Abstract: The Loess Plateau region in northwestern China has experienced severe water resource shortages due to the combined impacts of climate and land use changes and water resource exploitation during the past decades. This study was designed to examine the impacts of climatic variability on streamflow characteristics of a 12‐km² watershed near Tianshui City, Gansu Province in northwestern China. Statistic analytical methods including Kendall’s trend test and stepwise regression were used to detect trends in relationship between observed streamflow and climatic variables. Sensitivity analysis based on an evapotranspiration model was used to detect quantitative hydrologic sensitivity to climatic variability. We found that precipitation (P), potential evapotranspiration (PET) and streamflow (Q) were not statistically significantly different (p > 0.05) over the study period between 1982 and 2003. Stepwise regression and sensitivity analysis all indicated that P was more influential than PET in affecting annual streamflow, but the similar relationship existed at the monthly scale. The sensitivity of streamflow response to variations of P and PET increased slightly with the increase in watershed dryness (PET/P) as well as the increase in runoff ratio (Q/P). This study concluded that future changes in climate, precipitation in particular, will significantly impact water resources in the Loess Plateau region an area that is already experiencing a decreasing trend in water yield.     

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.     

A Comparative Assessment of Runoff Nitrogen from Turf,Forest, Meadow,and Mixed Landuse Watersheds

Landscaping paradigms that encourage high‐input, intensively managed and mono‐culture turf/lawn landscapes have raised concerns about water quality. We conducted a watershed‐scale assessment of landscaping practices that included turf, urban, forest, native meadow, and mixed landuse watersheds with a professional golf course and a parking lot. The turf site was moderately managed and had lower fertilizer inputs than those typically used by homeowners and golf courses. Stream water sampling was performed during base flow and storm events. Highest nitrate and total nitrogen concentrations in runoff were observed for the mixed watershed draining the golf course. In contrast, concentrations in base flow from the turf watershed were lower than expected and were comparable to those measured in the surrounding meadow and forest sites. Total nitrogen concentrations from the turf site increased sharply during the first storms following fertilization, suggesting that despite optimal management there exists a risk for nutrient runoff following fertilization. Overall, this study suggests that turf or lawns, when managed properly, pose minimal water quality risk to surface waters. Rate, timing of application, and the type of fertilizer appear to be the key factors affecting water quality. Better education of homeowners and landscaping professionals with regard to these factors may be a cost‐effective strategy to reduce nonpoint source pollution.     

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.     

Effects of Rainfall and Ground‐Water Pumping on Streamflow in Mākaha,O’ahu,Hawai’i1

Abstract: Land‐use/land‐cover changes in Mākaha valley have included the development of agriculture, residential dwellings, golf courses, potable water supply facilities, and the introduction of alien species. The impact of these changes on surface water and ground water resources in the valley is of concern. In this study, streamflow, rainfall, and ground‐water pumping data for the upper part of the Mākaha valley watershed were evaluated to identify corresponding trends and relationships. The results of this study indicate that streamflow declined during the ground‐water pumping period. Mean and median annual streamflow have declined by 42% (135 mm) and 56% (175 mm), respectively, and the mean number of dry stream days per year has increased from 8 to 125. Rainfall across the study area appears to have also declined though it is not clear whether the reduction in rainfall is responsible for all or part of the observed streamflow decline. Mean annual rainfall at one location in the study area declined by 14% (179 mm) and increased by 2% (48 mm) at a second location. Further study is needed to assess the effect of ground‐water pumping and to characterize the hydrologic cycle with respect to rainfall, infiltration, ground‐water recharge and flow in the study area, and stream base flow and storm flow.     

Estimation of Evapotranspiration Across the Conterminous United States Using a Regression With Climate and Land‐Cover Data1

Sanford, Ward E. and David L. Selnick, 2012. Estimation of Evapotranspiration Across the Conterminous United States Using a Regression with Climate and Land‐Cover Data. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12010 Abstract: Evapotranspiration (ET) is an important quantity for water resource managers to know because it often represents the largest sink for precipitation (P) arriving at the land surface. In order to estimate actual ET across the conterminous United States (U.S.) in this study, a water‐balance method was combined with a climate and land‐cover regression equation. Precipitation and streamflow records were compiled for 838 watersheds for 1971‐2000 across the U.S. to obtain long‐term estimates of actual ET. A regression equation was developed that related the ratio ET/P to climate and land‐cover variables within those watersheds. Precipitation and temperatures were used from the PRISM climate dataset, and land‐cover data were used from the USGS National Land Cover Dataset. Results indicate that ET can be predicted relatively well at a watershed or county scale with readily available climate variables alone, and that land‐cover data can also improve those predictions. Using the climate and land‐cover data at an 800‐m scale and then averaging to the county scale, maps were produced showing estimates of ET and ET/P for the entire conterminous U.S. Using the regression equation, such maps could also be made for more detailed state coverages, or for other areas of the world where climate and land‐cover data are plentiful.     

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.     

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.     

Reconstructions of Columbia River Streamflow from Tree‐Ring Chronologies in the Pacific Northwest,USA

We developed Columbia River streamflow reconstructions using a network of existing, new, and updated tree‐ring records sensitive to the main climatic factors governing discharge. Reconstruction quality is enhanced by incorporating tree‐ring chronologies where high snowpack limits growth, which better represent the contribution of cool‐season precipitation to flow than chronologies from trees positively sensitive to hydroclimate alone. The best performing reconstruction (back to 1609 CE) explains 59% of the historical variability and the longest reconstruction (back to 1502 CE) explains 52% of the variability. Droughts similar to the high‐intensity, long‐duration low flows observed during the 1920s and 1940s are rare, but occurred in the early 1500s and 1630s‐1640s. The lowest Columbia flow events appear to be reflected in chronologies both positively and negatively related to streamflow, implying low snowpack and possibly low warm‐season precipitation. High flows of magnitudes observed in the instrumental record appear to have been relatively common, and high flows from the 1680s to 1740s exceeded the magnitude and duration of observed wet periods in the late‐19th and 20th Century. Comparisons between the Columbia River reconstructions and future projections of streamflow derived from global climate and hydrologic models show the potential for increased hydrologic variability, which could present challenges for managing water in the face of competing demands.