Macroscale water fluxes 3. Effects of land processes on variability of monthly river discharge

P. C. D. Milly; R.T. Wetherald

doi:10.1029/2001wr000761

Home/Archive/USGS

USGS · 70024193

Macroscale water fluxes 3. Effects of land processes on variability of monthly river discharge

P. C. D. Milly·R.T. Wetherald

View original source ↗Export metadata ↓

Abstract

A salient characteristic of river discharge is its temporal variability. The time series of flow at a point on a river can be viewed as the superposition of a smooth seasonal cycle and an irregular, random variation. Viewing the random component in the spectral domain facilitates both its characterization and an interpretation of its major physical controls from a global perspective. The power spectral density functions of monthly flow anomalies of many large rivers worldwide are typified by a “red noise” process: the density is higher at low frequencies (e.g., <1 y −1 ) than at high frequencies, indicating disproportionate (relative to uncorrelated “white noise”) contribution of low frequencies to variability of monthly flow. For many high‐latitude and arid‐region rivers, however, the power is relatively evenly distributed across the frequency spectrum. The power spectrum of monthly flow can be interpreted as the product of the power spectrum of monthly basin total precipitation (which is typically white or slightly red) and several filters that have physical significance. The filters are associated with (1) the conversion of total precipitation (sum of rainfall and snowfall) to effective rainfall (liquid flux to the ground surface from above), (2) the conversion of effective rainfall to soil water excess (runoff), and (3) the conversion of soil water excess to river discharge. Inferences about the roles of each filter can be made through an analysis of observations, complemented by information from a global model of the ocean‐atmosphere‐land system. The first filter causes a snowmelt‐related amplification of high‐frequency variability in those basins that receive substantial snowfall. The second filter causes a relatively constant reduction in variability across all frequencies and can be predicted well by means of a semiempirical water balance relation. The third filter, associated with groundwater and surface water storage in the river basin, causes a strong reduction in high‐frequency variability of many basins. The strength of this reduction can be quantified by an average residence time of water in storage, which is typically on the order of 20–50 days. The residence time is demonstrably influenced by freezing conditions in the basin, fractional cover of the basin by lakes, and runoff ratio (ratio of mean runoff to mean precipitation). Large lake areas enhance storage and can greatly increase total residence times (100 to several hundred days). Freezing conditions appear to cause bypassing of subsurface storage, thus reducing residence times (0–30 days). Small runoff ratios tend to be associated with arid regions, where the water table is deep, and consequently, most of the runoff is produced by processes that bypass the saturated zone, leading to relatively small residence times for such basins (0–40 days).

Explore related subjects

Water Resources Research Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

P. C. D. Milly, R.T. Wetherald. 2002-11-15. Macroscale water fluxes 3. Effects of land processes on variability of monthly river discharge. https://doi.org/10.1029/2001wr000761

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Predicting Minnesota lake ice phenology with deep learning, explainable methods, and a physically based benchmark, 1980-2018

Globally, lakes are losing ice cover, but our understanding of the rates and variability of change are biased toward few lakes with long-term records. Recent works have used remote sensing and predictive modeling to supplement the observational record, but these approaches produce large errors when estimating ice phenology for individual lakes. Accurate, lake-specific estimates of ice formation and breakup through time are important for estimating variability in ice loss across heterogenous lakes and for understanding broader implications. This work explores machine learning (ML) approaches for hindcasting ice phenology using observations from 1980 to 2018 across 625 Minnesota lakes, covering 4359 lake-years of record. We used daily weather and static lake attributes to develop 60 neural networks for hindcasting lake ice time series. We considered LSTMs and attention-based transformers of varying size and initial parameters. We found that the largest LSTM performed most accurately on withheld data, and on a test set of unseen years and lakes, it outperformed a state-of-the-art physically based model in daily accuracy (97% vs. 95%), year-level metrics (RMSE for ice formation = 6.9 vs. 13.2 days; ice breakup = 7.5 vs. 12.2 days; ice duration = 9.6 vs. 13.7 days), and estimating loss of ice cover from 1980 to 2018 (observed = 8.4 days, LSTM = 7.2 days, GLM = 3.4 days). This demonstrates that ML can estimate historical lake ice formation and breakup dates within a week of observed phenology across heterogeneous lakes and recreate broad scale trends in ice cover phenology with limited long-term ice records.

Minnesota↗

USGS2026

A roadmap for identifying and interpreting physical processes and national water model prediction bias associated with baseflow index regimes across the contiguous United States

Understanding how groundwater–surface water interactions shape streamflow variability is critical for diagnosing low flow behavior and prediction bias in continental scale hydrologic models. We present a process informed framework that links observed baseflow (BF) dynamics, watershed attributes, and National Water Model (NWM) performance across the contiguous United States. Using daily observed streamflow from 797 reference quality streamgages, we developed monthly baseflow index (BFI) signatures using a streamgage specific, calibrated digital filter. Hierarchical clustering of these signatures identified seven distinct BFI regimes capturing regional and seasonal variability. We evaluated NWM v3.0 retrospective streamflow performance within each regime using multiple hydrograph and flow duration curve-based metrics. Model skill varied systematically across regimes: mixed flow systems were simulated most accurately, while predominantly BF dominated and quickflow dominated regimes exhibited substantially poorer performance. Across nearly all regimes, the NWM underestimated observed BFI magnitude and frequently failed to reproduce seasonal BF patterns, indicating systematic biases in simulated low flow contributions. To relate these regimes to potential process controls, we trained a Random Forest classifier using static watershed attributes and applied Shapley Additive Explanations to identify features most strongly associated with each regime. Results highlight regionally varying influences, including the dominant role of snow fraction and seasonal runoff timing in snow dominated basins and the importance of evapotranspiration and aridity in quickflow dominated systems. Collectively, these findings demonstrate how hydrologic signatures combined with interpretable machine learning can diagnose regime specific model biases and generate process-based hypotheses about limitations in large scale hydrologic prediction systems.

contiguous United States↗

USGS2026

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research↗