USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30Linked to original sources
USGS2003

Modelling hydrologic responses in a small forested catchment (Panola Mountain, Georgia, USA): A comparison of the original and a new dynamic TOPMODEL

Preliminary modelling results for a new version of the rainfall-runoff model TOPMODEL, dynamic TOPMODEL, are compared with those of the original TOPMODEL formulation for predicting streamflow at the Panola Mountain Research Watershed, Georgia. Dynamic TOPMODEL uses a kinematic wave routing of subsurface flow, which allows for dynamically variable upslope contributing areas, while retaining the concept of hydrological similarity to increase computational efficiency. Model performance in predicting discharge was assessed for the original TOPMODEL and for one landscape unit (LU) and three LU versions of the dynamic TOPMODEL (a bare rock area, hillslope with regolith <1 m, and a riparian zone with regolith ???5 m). All simulations used a 30 min time step for each of three water years. Each 1-LU model underpredicted the peak streamflow, and generally overpredicted recession streamflow during wet periods and underpredicted during dry periods. The difference between predicted recession streamflow generally was less for the dynamic TOPMODEL and smallest for the 3-LU model. Bayesian combination of results for different water years within the GLUE methodology left no behavioural original or 1-LU dynamic models and only 168 (of 96 000 sample parameter sets) for the 3-LU model. The efficiency for the streamflow prediction of the best 3-LU model was 0.83 for an individual year, but the results suggest that further improvements could be made. ?? 2003 John Wiley & Sons, Ltd.

Hydrological Processes↗
USGS2017

Seasonality of stable isotope composition of atmospheric water input at the southern slopes of Mt. Kilimanjaro, Tanzania

To understand the moisture regime at the southern slopes of Mt. Kilimanjaro, we analysed the isotopic variability of oxygen (δ 18 O) and hydrogen (δD) of rainfall, throughfall, and fog from a total of 2,140 samples collected weekly over 2 years at 9 study sites along an elevation transect ranging from 950 to 3,880 m above sea level. Precipitation in the Kilimanjaro tropical rainforests consists of a combination of rainfall, throughfall, and fog. We defined local meteoric water lines for all 3 precipitation types individually and the overall precipitation, δD prec = 7.45 (±0.05) × δ 18 O prec + 13.61 (±0.20), n = 2,140, R 2 = .91, p < .001. We investigated the precipitation-type-specific stable isotope composition and analysed the effects of amount, altitude, and temperature. Aggregated annual mean values revealed isotope composition of rainfall as most depleted and fog water as most enriched in heavy isotopes at the highest elevation research site. We found an altitude effect of δ 18 O rain = −0.11‰ × 100 m −1 , which varied according to precipitation type and season. The relatively weak isotope or altitude gradient may reveal 2 different moisture sources in the research area: (a) local moisture recycling and (b) regional moisture sources. Generally, the seasonality of δ 18 O rain values follows the bimodal rainfall distribution under the influences of south- and north-easterly trade winds. These seasonal patterns of isotopic composition were linked to different regional moisture sources by analysing Hybrid Single Particle Lagrangian Integrated Trajectory backward trajectories. Seasonality of d excess values revealed evidence of enhanced moisture recycling after the onset of the rainy seasons. This comprehensive dataset is essential for further research using stable isotopes as a hydrological tracer of sources of precipitation that contribute to water resources of the Kilimanjaro region.

Mt. Kilimanjaro↗
USGS2017

Stormwater management network effectiveness and implications for urban watershed function: A critical review

Deleterious effects of urban stormwater are widely recognized. In several countries, regulations have been put into place to improve the conditions of receiving water bodies, but planning and engineering of stormwater control is typically carried out at smaller scales. Quantifying cumulative effectiveness of many stormwater control measures on a watershed scale is critical to understanding how small-scale practices translate to urban river health. We review 100 empirical and modelling studies of stormwater management effectiveness at the watershed scale in diverse physiographic settings. Effects of networks with stormwater control measures (SCMs) that promote infiltration and harvest have been more intensively studied than have detention-based SCM networks. Studies of peak flows and flow volumes are common, whereas baseflow, groundwater recharge, and evapotranspiration have received comparatively little attention. Export of nutrients and suspended sediments have been the primary water quality focus in the United States, whereas metals, particularly those associated with sediments, have received greater attention in Europe and Australia. Often, quantifying cumulative effects of stormwater management is complicated by needing to separate its signal from the signal of urbanization itself, innate watershed characteristics that lead to a range of hydrologic and water quality responses, and the varying functions of multiple types of SCMs. Biases in geographic distribution of study areas, and size and impervious surface cover of watersheds studied also limit our understanding of responses. We propose hysteretic trajectories for how watershed function responds to increasing imperviousness and stormwater management. Even where impervious area is treated with SCMs, watershed function may not be restored to its predevelopment condition because of the lack of treatment of all stormwater generated from impervious surfaces; non-additive effects of individual SCMs; and persistence of urban effects beyond impervious surfaces. In most cases, pollutant load decreases largely result from run-off reductions rather than lowered solute or particulate concentrations. Understanding interactions between natural and built landscapes, including stormwater management strategies, is critical for successfully managing detrimental impacts of stormwater at the watershed scale.

Europe↗
USGS2017

Hydrologic metrics for status-and-trends monitoring in urban and urbanizing watersheds

Local governmental agencies are increasingly undertaking potentially costly “status-and-trends” monitoring to evaluate the effectiveness of stormwater control measures and land-use planning strategies, or to satisfy regulatory requirements. Little guidance is presently available for such efforts, and so we have explored the application, interpretation, and temporal limitations of well-established hydrologic metrics of runoff changes from urbanization, making use of an unusually long-duration, high-quality data set from the Pacific Northwest (USA) with direct applicability to urban and urbanizing watersheds. Three metrics previously identified for their utility in identifying hydrologic conditions with biological importance that respond to watershed urbanization—T Qmean (the fraction of time that flows exceed the mean annual discharge), the Richards-Baker Index (characterizing flashiness relative to the mean discharge), and the annual tally of wet-season day-to-day flow reversals (the total number of days that reverse the prior days’ increasing or decreasing trend)—are all successful in stratifying watersheds across a range of urbanization, as measured by total contributing area of urban development. All metrics respond with statistical significance to multi-decadal trends in urbanization, but none detect trends in watershed-scale urbanization over the course of a single decade. This suggests a minimum period over which dependable trends in hydrologic alteration (or improvement) can be detected with confidence. The metrics also prove less well suited to urbanizing watersheds in a semi-arid climate, with only flow reversals showing a response consistent with prior findings from more humid regions. We also explore the use of stage as a surrogate for discharge in calculating these metrics, recognizing potentially significant agency cost savings in data collection with minimal loss of information. This approach is feasible but cannot be implemented under current data-reporting practices, requiring measurement of water-depth values and preservation of the full precision of the original recorded data. With these caveats, however, hydrologic metrics based on stage should prove as or more useful, at least in the context of status-and-trends monitoring, as those based on subsequent calculations of discharge.

Hydrological Processes↗
USGS2017

Estimating daily lake evaporation from biweekly energy‐budget data

Estimates of daily lake evaporation based on energy‐budget data are poor because of large errors associated with quantifying change in lake heat storage over periods of less than about 10 days. Energy‐budget evaporation was determined during approximately biweekly periods at a northern Minnesota, USA, lake for 5 years. Various combinations of shortwave radiation, air temperature, wind speed, lake‐surface temperature, and vapour‐pressure difference were related to energy‐budget evaporation using linear‐regression models in an effort to determine daily evaporation without requiring the heat‐storage term. The model that combined the product of shortwave radiation and air temperature with the product of vapour‐pressure difference and wind speed provided the second best fit based on statistics but provided the best daily data based on comparisons with evaporation determined with the eddy‐covariance method. Best‐model daily values ranged from −0.6 to 7.1 mm/day over a 5‐year period. Daily averages of best‐model evaporation and eddy‐covariance evaporation were nearly identical for all 28 days of comparisons with a standard deviation of the differences between the two methods of 0.68 mm/day. Best‐model daily evaporation also was compared with two other evaporation models, Jensen–Haise and a mass‐transfer model. Best‐model daily values were substantially improved relative to Jensen–Haise and mass‐transfer values when daily values were summed over biweekly energy‐budget periods for comparison with energy‐budget results.

Northern Minnesota↗
USGS2017

Estimating restorable wetland water storage at landscape scales

Globally, hydrologic modifications such as ditching and subsurface drainage have significantly reduced wetland water storage capacity (i.e., volume of surface water a wetland can retain) and consequent wetland functions. While wetland area has been well documented across many landscapes and used to guide restoration efforts, few studies have directly quantified the associated wetland storage capacity. Here, we present a novel raster-based approach to quantify both contemporary and potential (i.e., restorable) storage capacities of individual depressional basins across landscapes. We demonstrate the utility of this method by applying it to the Delmarva Peninsula, a region punctuated by both depressional wetlands and drainage ditches. Across the entire peninsula, we estimated that restoration (i.e., plugging ditches) could increase storage capacity by 80%. Focusing on an individual watershed, we found that over 59% of restorable storage capacity occurs within 20 m of the drainage network, and that 93% occurs within 1 m elevation of the drainage network. Our demonstration highlights widespread ditching in this landscape, spatial patterns of both contemporary and potential storage capacities, and clear opportunities for hydrologic restoration. In Delmarva and more broadly, our novel approach can inform targeted landscape-scale conservation and restoration efforts to optimize hydrologically mediated wetland functions.

Hydrological Processes↗
USGS2018

Modelling surface-water depression storage in a Prairie Pothole Region

In this study, the Precipitation-Runoff Modelling System (PRMS) was used to simulate changes in surface-water depression storage in the 1,126-km 2 Upper Pipestem Creek basin located within the Prairie Pothole Region of North Dakota, USA. The Prairie Pothole Region is characterized by millions of small water bodies (or surface-water depressions) that provide numerous ecosystem services and are considered an important contribution to the hydrologic cycle. The Upper Pipestem PRMS model was extracted from the U.S. Geological Survey's (USGS) National Hydrologic Model (NHM), developed to support consistent hydrologic modelling across the conterminous United States. The Geospatial Fabric database, created for the USGS NHM, contains hydrologic model parameter values derived from datasets that characterize the physical features of the entire conterminous United States for 109,951 hydrologic response units. Each hydrologic response unit in the Geospatial Fabric was parameterized using aggregated surface-water depression area derived from the National Hydrography Dataset Plus, an integrated suite of application-ready geospatial datasets. This paper presents a calibration strategy for the Upper Pipestem PRMS model that uses normalized lake elevation measurements to calibrate the parameters influencing simulated fractional surface-water depression storage. Results indicate that inclusion of measurements that give an indication of the change in surface-water depression storage in the calibration procedure resulted in accurate changes in surface-water depression storage in the water balance. Regionalized parameterization of the USGS NHM will require a proxy for change in surface-storage to accurately parameterize surface-water depression storage within the USGS NHM.

North Dakota↗
USGS2018

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity ( K s ) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K s , and the spatial correlation length scale associated with K s . The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale K s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.

California↗
USGS2018

Comparison of groundwater age models for assessing nitrate loading, transport pathways, and management options in a complex aquifer system

In an aquifer system with complex hydrogeology, mixing of groundwater with different ages could occur associated with various flow pathways. In this study, we applied different groundwater age estimation techniques (lumped parameter model, and numerical model) to characterize groundwater age distributions and the major pathways of nitrate contamination in the Gosan agricultural field, Jeju Island. According to the lumped parameter model, groundwater age in the study area could be explained by the binary mixing of the young groundwater (4-33 years) and the old water component (>60 years). The complex hydrogeologic regimes and local heterogeneity observed in the study area (multi-layered aquifer, well leakage hydraulics) were particularly well reflected in the numerical model. The numerical model predicted that the regional aquifer of Gosan responded to the fertilizer applications more rapidly (mean age: 9.7-22.3 years) than as estimated by other models. Our study results demonstrated that application and comparison of multiple age estimation methods can be useful to understand better the flow regimes and the mixing characteristics of groundwater with different ages (pathways), hence, to reduce the risk of improper groundwater management plan arising from the aquifer heterogeneity.

Jeju Island↗
USGS2003

Characterization of lake water and ground water movement in the littoral zone of Williams Lake, a closed-basin lake in North central Minnesota

Williams Lake, Minnesota is a closed‐basin lake that is a flow‐through system with respect to ground water. Ground‐water input represents half of the annual water input and most of the chemical input to the lake. Chemical budgets indicate that the lake is a sink for calcium, yet surficial sediments contain little calcium carbonate. Sediment pore‐water samplers (peepers) were used to characterize solute fluxes at the lake‐water–ground‐water interface in the littoral zone and resolve the apparent disparity between the chemical budget and sediment data. Pore‐water depth profiles of the stable isotopes δ 18 O and δ 2 H were non‐linear where ground water seeped into the lake, with a sharp transition from lake‐water values to ground‐water values in the top 10 cm of sediment. These data indicate that advective inflow to the lake is the primary mechanism for solute flux from ground water. Linear interstitial velocities determined from δ 2 H profiles (316 to 528 cm/yr) were consistent with velocities determined independently from water budget data and sediment porosity (366 cm/yr). Stable isotope profiles were generally linear where water flowed out of the lake into ground water. However, calcium profiles were not linear in the same area and varied in response to input of calcium carbonate from the littoral zone and subsequent dissolution. The comparison of pore‐water calcium profiles to pore‐water stable isotope profiles indicate calcium is not conservative. Based on the previous understanding that 40–50 % of the calcium in Williams Lake is retained, the pore‐water profiles indicate aquatic plants in the littoral zone are recycling the retained portion of calcium. The difference between the pore‐water depth profiles of calcium and δ 18 O and δ 2 H demonstrate the importance of using stable isotopes to evaluate flow direction and source through the lake‐water–ground‐water interface and evaluate mechanisms controlling the chemical balance of lakes.

Minnesota↗
USGS2002

Interannual variations in snowpack in the Crown of the Continent Ecosystem

Ecosystem changes such as glacier recession and alpine treeline advance have been documented over the previous 150 years in the Rocky Mountains of northern Montana and southern British Columbia and Alberta, a region known as the Crown of the Continent Ecosystem (CCE). Such changes are controlled, at least partially, by variations in snowpack. The CCE consists primarily of public lands, the majority of which is undeveloped or wilderness. Consequently, this region is well suited for an examination of long-term snowpack variation and associated ecosystem change. Data from nine SNOTEL sites provide an indication of the daily accumulation and ablation of snowpack over the period 1977-2001, as well as the relationship between precipitation, temperature and snowpack. 1 April data from 21 snow courses indicated the extent of regional snowpack variation and trends over the period 1950-2001, and 1 May data from three snow courses in Glacier National Park allow this record to be extended back to 1922. SNOTEL data suggest CCE snowpacks are larger and more persistent than in most regions of the western USA, and that water year precipitation is the primary control on 1 April snow water equivalent (SWE). Snow course data indicate that variations in both 1 April and 1 May mean SWE are closely tied to the Pacific decadal oscillation, an El Nino-southern oscillation-like interdecadal pattern of Pacific Ocean climate variability. Despite relatively stable snowpacks and summer temperatures since 1922, the glaciers in Glacier National Park have receded steadily during this period, implying a significant climatic shift between their Little Ice Age glacial maxima (ca 1860) and 1922. Published in 2002 by John Wiley and Sons, Ltd.

Montana↗
USGS2003

An evaluation of sediment rating curves for estimating suspended sediment concentrations for subsequent flux calculations

In the absence of actual suspended sediment concentration (SSC) measurements, hydrologists have used sediment rating (sediment transport) curves to estimate (predict) SSCs for subsequent flux calculations. Various evaluations of the sediment rating-curve method were made using data from long-term, daily sediment-measuring sites within large (>1 000 000 km2), medium (<1 000 000 to > 1000 km2), and small (<1000 km2) river basins in the USA and Europe relative to the estimation of suspended sediment fluxes. The evaluations address such issues as the accuracy of flux estimations for various levels of temporal resolution as well as the impact of sampling frequency on the magnitude of flux estimation errors. The sediment rating-curve method tends to underpredict high, and overpredict low SSCs. As such, the range of errors associated with concomitant flux estimates for relatively short time-frames (e.g. daily, weekly) are likely to be substantially larger than those associated with longer time-frames (e.g. quarterly, annually) because the over- and underpredictions do not have sufficient time to balance each other. Hence, when error limits must be kept under ??20%, temporal resolution probably should be limited to quarterly or greater. The evaluations indicate that over periods of 20 or more years, errors of <1% can be achieved using a single sediment rating curve based on data spanning the entire period. However, somewhat better estimates for the entire period, and markedly better annual estimates within the period, can be obtained if individual annual sediment rating curves are used instead. Relatively accurate (errors <??20%) annual suspended sediment fluxes can be obtained from hydrologically based monthly measurements/samples. For 5-year periods or longer, similar results can be obtained from measurements/samples collected once every 2 months. In either case, hydrologically based sampling, as opposed to calendar-based sampling is likely to limit the magnitude of flux estimation errors. Published in 2003 by John Wiley and Sons, Ltd.

Hydrological Processes↗
USGS2018

Influence of climate on alpine stream chemistry and water sources

The resilience of alpine/subalpine watersheds may be viewed as the resistance of streamflow or stream chemistry to change under varying climatic conditions, which is governed by the relative size (volume) and transit time of surface and subsurface water sources. Here, we use end‐member mixing analysis in Andrews Creek, an alpine stream in Rocky Mountain National Park, Colorado, from water year 1994 to 2015, to explore how the partitioning of water sources and associated hydrologic resilience change in response to climate. Our results indicate that four water sources are significant contributors to Andrews Creek, including snow, rain, soil water, and talus groundwater. Seasonal patterns in source‐water contributions reflected the seasonal hydrologic cycle, which is driven by the accumulation and melting of seasonal snowpack. Flushing of soil water had a large effect on stream chemistry during spring snowmelt, despite making only a small contribution to streamflow volume. Snow had a large influence on stream chemistry as well, contributing large amounts of water with low concentrations of weathering products. Interannual patterns in end‐member contributions reflected responses to drought and wet periods. Moderate and significant correlations exist between annual end‐member contributions and regional‐scale climate indices (the Palmer Drought Severity Index, the Palmer Hydrologic Drought Index, and the Modified Palmer Drought Severity Index). From water year 1994 to 2015, the percent contribution from the talus‐groundwater end member to Andrews Creek increased an average of 0.5% per year (p < 0.0001), whereas the percent contributions from snow plus rain decreased by a similar amount (p = 0.001). Our results show how water and solute sources in alpine environments shift in response to climate variability and highlight the role of talus groundwater and soil water in providing hydrologic resilience to the system.

Colorado↗
USGS2018

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.

Hydrological Processes↗
USGS2018

Thresholds and relations for soil‐hydraulic and soil‐physical properties as a function of burn severity 4 years after the 2011 Las Conchas Fire, New Mexico, USA

Wildfire effects on soil‐physical and ‐hydraulic properties as a function of burn severity are poorly characterized, especially several years after wildfire. A stratified random sampling approach was used in 2015 to sample seven sites representing a spectrum of remotely sensed burn severity in the area impacted by the 2011 Las Conchas Fire in New Mexico, USA. Replicate samples from each site were analysed in the laboratory. Linear and linear indicator regression were used to assess thresholds in soil‐physical and ‐hydraulic properties and functional relations with remotely sensed burn severity. Significant thresholds were present for initial soil‐water content (θ i ) at 0–6 cm depth between the change in the Normalized Burn Ratio ( dNBR ) equal to 618–802, for bulk density (ρ b ) at 3–6 cm between dNBR equal to 416–533, for gravel fraction at 0–1 cm between dNBR equal to 416–533, for fines (the silt + clay fraction) at 0–1 cm for dNBR equal to 416–533, and for fines at 3–6 cm for dNBR equal to 293–416. Significant linear relations with dNB R were present between ρ b at 0–1 cm, loss on ignition (LOI) at 0–1 cm, gravel fraction at 0–1 cm, and the large organic fraction at 1–3 cm. No thresholds or effects on soil‐hydraulic properties of field‐saturated hydraulic conductivity or sorptivity were observed. These results suggest that ρ b and LOI at 0–1 cm have residual direct impacts from the wildfire heat impulse. The θ i threshold is most likely from delayed groundcover/vegetation recovery that increases evaporation at the highest burn severity sites. Gravel and silt + clay thresholds at 0–1 cm at the transition to high burn severity suggest surface gravel lag development from hydraulic erosion. Thresholds in ρ b from 3 to 6 cm and in silt + clay fraction from 3 to 6 cm appear to be the result of soil variability between sites rather than wildfire impacts. This work suggests that gravel‐rich soils may have increased resilience to sustained surface runoff generation and erosion following wildfire, with implications for assessments of postwildfire hydrologic and erosion recovery potential.

New Mexico↗
USGS2018

Do we know how much fluvial sediment reaches the sea? Decreased river monitoring of U.S. coastal rivers

Given the present and future changing climate and human changes to land use and river control, river sediment fluxes to coastal systems are changing and will continue to change in the future. To delineate these changes and their effects, it is increasingly important to document the fluxes of river-borne sediment discharged to the sea. Unfortunately, broad-scale river sediment monitoring programs established more than 50 years ago in the U.S. have diminished substantially and now focus principally on the largest rivers and estuaries. Unless addressed, these data gaps will provide significant challenges in addressing fundamental scientific and management questions about the effects of climate change and sea-level rise in our estuaries and on our coasts.

Hydrological Processes↗
USGS2018

Systematic variation in evapotranspiration trends and drivers across the Northeastern United States

The direction and magnitude of responses of evapotranspiration (ET) to climate change are important to understand, as ET represents a major water and energy flux from terrestrial ecosystems, with consequences that feed back to the climate system. We inferred multidecadal trends in water balance in 11 river basins (1940–2012) and eight smaller watersheds (with records ranging from 18 to 61 years in length) in the Northeastern United States. Trends in river basin actual ET (AET) varied across the region, with an apparent latitudinal pattern: AET increased in the cooler northern part of the region (Maine) but decreased in some warmer regions to the southwest (Pennsylvania–Ohio). Of the four small watersheds with records longer than 45 years, two fit this geographic pattern in AET trends. The differential effects of the warming climate on AET across the region may indicate different mechanisms of change in more‐ vs. less‐energy‐limited watersheds, even though annual precipitation greatly exceeds potential ET across the entire region. Correlations between AET and time series of temperature and precipitation also indicate differences in limiting factors for AET across the Northeastern U.S. climate gradient. At many sites across the climate gradient, water‐year AET correlated with summer precipitation, implying that water limitation is at least transiently important in some years, whereas correlations with temperature indices were more prominent in northern than southern sites within the region.

Hydrological Processes↗