Global Reconstruction of Naturalized River Flows at 2.94 Million Reaches

Undervalued dry riverbeds: A key factor in equating intermittent river CO2 emissions to perennial rivers

Intermittent rivers in semiarid and arid regions, constituting over half of the world's rivers, alternate the carbon cycle interactions among the biosphere, hydrosphere, and atmosphere. Inadequate quantification of flow duration and river water surface area, along with overlooked CO2 emissions from dry riverbeds, result in notable inaccuracies in global carbon cycle assessments. High-resolution remote sensing images combined with intensive field measurements and hydrological modelling were used to estimate and extract the flow duration, river water surface area and dry riverbed area of Huangfuchuan, an intermittent river watershed that acts as a major tributary of the Yellow River in semiarid Northwest China. CO2 emission rates and partial pressures in water and air across the watershed were in-situ measured. In 2018, the flow duration of Huangfuchuan increased from less than 5 days in the first-order tributary to 150 days in the sixth-order mainstream. River water surface area estimated by remote sensing extraction plus the hydrodynamic model simulation varied from 3.9 to 88.6 km2 under 5 %-95 % discharge frequencies. CO2 emissions from the water-air interface and dry riverbed in 2018 were estimated at 582.3 × 103 and 355.2 × 103 ton, respectively. The estimated total annual emission (937.5 × 103 ton) aligns closely with the range of emissions (67.3 × 103-1377.2 × 103 ton) calculated for the water-air interface alone, derived using DEM river length and hydraulic geometry method. This similarity can be attributed to the overestimation of flow duration and flow velocity, as well as the over- or under-estimation of river water surface area and slope. The new method proposed in this study has large potential to be applied in estimating CO2 emissions from data-scarce intermittent rivers located in mountainous regions and provides a standardized solution in the estimation of CO2 emission. Results of this research reveal the spatiotemporal distribution of CO2 emissions along an intermittent river system and highlight the substantial role of dry riverbed in carbon cycle.

Tibetan Plateau Runoff and Evapotranspiration Dataset by an observation-constrained cryosphere-hydrology model

Runoff and evapotranspiration (ET) are pivotal constituents of the water, energy, and carbon cycles. This research presents a 5-km monthly gridded runoff and ET dataset for 1998-2017, encompassing seven headwaters of Tibetan Plateau rivers (Yellow, Yangtze, Mekong, Salween, Brahmaputra, Ganges, and Indus) (hereinafter TPRED). The dataset was generated using the advanced cryosphere-hydrology model WEB-DHM, yielding a Nash coefficient ranging from 0.77 to 0.93 when compared to the observed discharges. The findings indicate that TPRED's monthly runoff notably outperforms existing datasets in capturing hydrological patterns, as evidenced by robust metrics such as the correlation coefficient (CC) (0.944-0.995), Bias (-0.68-0.53), and Root Mean Square Error (5.50-15.59 mm). Additionally, TPRED's monthly ET estimates closely align with expected seasonal fluctuations, as reflected by a CC ranging from 0.94 to 0.98 when contrasted with alternative ET products. Furthermore, TPRED's annual values exhibit commendable concordance with operational products across multiple dimensions. Ultimately, the TPRED will have great application on hydrometeorology, carbon transport, water management, hydrological modeling, and sustainable development of water resources.

Enhancing Wetland Mapping: Integrating Sentinel-1/2, GEDI Data, and Google Earth Engine

Wetlands are amongst Earth's most dynamic and complex ecological resources, serving productive and biodiverse ecosystems. Enhancing the quality of wetland mapping through Earth observation (EO) data is essential for improving effective management and conservation practices. However, the achievement of reliable and accurate wetland mapping faces challenges due to the heterogeneous and fragmented landscape of wetlands, along with spectral similarities among different wetland classes. The present study aims to produce advanced 10 m spatial resolution wetland classification maps for four pilot sites on the Island of Newfoundland in Canada. Employing a comprehensive and multidisciplinary approach, this research leverages the synergistic use of optical, synthetic aperture radar (SAR), and light detection and ranging (LiDAR) data. It focuses on ecological and hydrological interpretation using multi-source and multi-sensor EO data to evaluate their effectiveness in identifying wetland classes. The diverse data sources include Sentinel-1 and -2 satellite imagery, Global Ecosystem Dynamics Investigation (GEDI) LiDAR footprints, the Multi-Error-Removed Improved-Terrain (MERIT) Hydro dataset, and the European ReAnalysis (ERA5) dataset. Elevation data and topographical derivatives, such as slope and aspect, were also included in the analysis. The study evaluates the added value of incorporating these new data sources into wetland mapping. Using the Google Earth Engine (GEE) platform and the Random Forest (RF) model, two main objectives are pursued: (1) integrating the GEDI LiDAR footprint heights with multi-source datasets to generate a 10 m vegetation canopy height (VCH) map and (2) seeking to enhance wetland mapping by utilizing the VCH map as an input predictor. Results highlight the significant role of the VCH variable derived from GEDI samples in enhancing wetland classification accuracy, as it provides a vertical profile of vegetation. Accordingly, VCH reached the highest accuracy with a coefficient of determination (R2) of 0.69, a root-mean-square error (RMSE) of 1.51 m, and a mean absolute error (MAE) of 1.26 m. Leveraging VCH in the classification procedure improved the accuracy, with a maximum overall accuracy of 93.45%, a kappa coefficient of 0.92, and an F1 score of 0.88. This study underscores the importance of multi-source and multi-sensor approaches incorporating diverse EO data to address various factors for effective wetland mapping. The results are expected to benefit future wetland mapping studies.

Remote Sensing-Based Extension of GRDC Discharge Time Series - A Monthly Product with Uncertainty Estimates

The Global Runoff Data Center (GRDC) data set has faced a decline in the number of active gauges since the 1980s, leaving only 14% of gauges active as of 2020. We develop the Remote Sensing-based Extension for the GRDC (RSEG) data set that can ingest legacy gauge discharge and remote sensing observations. We employ a stochastic nonparametric mapping algorithm to extend the monthly discharge time series for inactive GRDC stations, benefiting from satellite imagery- and altimetry-derived river width and water height observations. After a rigorous quality assessment of our estimated discharge, involving statistical validation, tests and visual inspection, results in the extension of discharge records for 3377 out of 6015 GRDC stations. The quality of discharge estimates for the rivers with a large or medium mean discharge is quite satisfactory (average KGE value > 0.5) however for river reaches with a low mean discharge the average KGE value drops to 0.33.The RSEG data set regains monitoring capability for 83% of total river discharge measured by GRDC stations, equivalent to 7895 km3/month.

Modeling the transport of microplastics along river networks

The excessive use of plastics in modern life has led to a significant increase in production and a corresponding rise in plastic waste generation. The slow degradation of plastics results in the introduction and accumulation of microplastics (MP) in the environment, posing environmental and health risks. River networks, acting as conduits between terrestrial and marine environments, play a crucial role in controlling the transport of MP. Predicting the complex processes of MP pathways in these environments is an ongoing challenge. To address this issue, we propose a model that integrates the advection-dispersion equation with anthropogenic MP loads and hydraulic river network characteristics. The validity of the model was assessed using literature data from three river networks worldwide. Model results show a good agreement between predictions and field observations (R2=0.72). Consequently, predicted MP data was used to perform a potential pollution assessment through the pollution load index, revealing in most cases higher MP contamination in headwaters stream and a dilution effect along the river network. The structure of the proposed model allows its further implementation to account for other transport mechanisms, interactions with other emerging contaminants (i.e., pharmaceuticals), and connections with other riverine environments, making it a valuable tool for understanding and mitigating MP pollution.

Comparing likely effectiveness of urban Nature-based Solutions worldwide: The example of riparian tree planting and water quality

Amongst a spectrum of benefits, Nature-based Solutions (NBS) are increasingly being advocated as improving the quality of aquatic environments in urban areas. Of these, a widely adopted measure is tree planting. Yet, because of the local complexities and spatial variability of urban hydrological response, it is difficult to predict to what extent improvements in water quality will arise. To overcome this barrier, a standardised approach to process-based model simulation of urban river quality is described (QUESTOR-YARDSTICK (QUESTOR-YS)). The approach eliminates the influence of point sources of pollution and harmonises the way in which river hydrodynamics and contributory catchment size are represented. Thereby, it focuses on differences in water quality between cities due solely to climate, river discharge and urban diffuse nutrient pollution factors. The relative sensitivity to NBS establishment between urban water bodies in different cities anywhere across the world can also potentially be quantified. The method can be readily extended to include wastewater effluents. The validity of the approach is demonstrated for a small river in Birmingham, UK; and thence demonstrated for the case of 10 km of riparian tree planting in Birmingham, Oslo (Norway) and Aarhus (Denmark). Modelling suggests that riparian tree planting can substantially improve water quality in each example city for three key indicators of water quality in sensitive summer conditions (water temperature, chlorophyll-a and dissolved oxygen). Results show the level of benefit achievable in response to a fixed amount of planting will depend on the existing level of riparian tree occupancy.

Hydrological model-based streamflow reconstruction for Indian sub-continental river basins, 1951-2021

Streamflow is a vital component of the global water cycle. Long-term streamflow observations are required for water resources planning and management, hydroclimatic extremes analysis, and ecological assessment. However, long-term streamflow observations for the Indian-Subcontinental (ISC) river basins are lacking. Using meteorological observations, state-of-the-art hydrological model, and river routing model, we developed hydrological model-simulated monthly streamflow from 1951-2021 for the ISC river basins. We used high-resolution vector-based routing model (mizuRoute) to generate streamflow at 9579 stream reaches in the sub-continental river basins. The model-simulated streamflow showed good performance against the observed flow with coefficient of determination (R2) and Nash-Sutcliffe efficiency (NSE) above 0.70 for more than 60% of the gauge stations. The dataset was used to examine the variability in low, average, and high flow across the streams. Long-term changes in streamflow showed a significant decline in flow in the Ganga basin while an increase in the semi-arid western India and Indus basin. Long-term streamflow can be used for planning water management and climate change adaptation in the Indian sub-continent.

Global methane emissions from rivers and streams

Methane (CH4) is a potent greenhouse gas and its concentrations have tripled in the atmosphere since the industrial revolution. There is evidence that global warming has increased CH4 emissions from freshwater ecosystems1,2, providing positive feedback to the global climate. Yet for rivers and streams, the controls and the magnitude of CH4 emissions remain highly uncertain3,4. Here we report a spatially explicit global estimate of CH4 emissions from running waters, accounting for 27.9 (16.7-39.7) Tg CH4 per year and roughly equal in magnitude to those of other freshwater systems5,6. Riverine CH4 emissions are not strongly temperature dependent, with low average activation energy (EM = 0.14 eV) compared with that of lakes and wetlands (EM = 0.96 eV)1. By contrast, global patterns of emissions are characterized by large fluxes in high- and low-latitude settings as well as in human-dominated environments. These patterns are explained by edaphic and climate features that are linked to anoxia in and near fluvial habitats, including a high supply of organic matter and water saturation in hydrologically connected soils. Our results highlight the importance of land-water connections in regulating CH4 supply to running waters, which is vulnerable not only to direct human modifications but also to several climate change responses on land.

Satellites reveal widespread decline in global lake water storage

Climate change and human activities increasingly threaten lakes that store 87% of Earth's liquid surface fresh water. Yet, recent trends and drivers of lake volume change remain largely unknown globally. Here, we analyze the 1972 largest global lakes using three decades of satellite observations, climate data, and hydrologic models, finding statistically significant storage declines for 53% of these water bodies over the period 1992-2020. The net volume loss in natural lakes is largely attributable to climate warming, increasing evaporative demand, and human water consumption, whereas sedimentation dominates storage losses in reservoirs. We estimate that roughly one-quarter of the world's population resides in a basin of a drying lake, underscoring the necessity of incorporating climate change and sedimentation impacts into sustainable water resources management.

Fate identification and management strategies of non-recyclable plastic waste through the integration of material flow analysis and leakage hotspot modeling

Low priority on waste management has impacted the complex environmental issue of plastic waste pollution, as evident by results of this study where it was found that 24.3% of waste generation in Jakarta and Bandung is emitted into the waterway due to the high intensity of human activity in the urban area. In this study, we investigated the viable integration between material flow analysis and leakage hotspot modeling to improve management strategies for plastic pollution in water systems and open environments. Using a multi-criteria assessment of plastic leakage from current waste management, a material flow analysis was developed on a city-wide scale defining the fate of plastic waste. Geospatial analysis was assigned to develop a calculation for identification and hydrological analysis while identifying the potential amount of plastic leakage to the river system. The results show that 2603 tons of plastic accumulated along the mainstream of the Ciliwung River on an annual basis, and a high-density population like that in Bandung discarded 1547 tons in a one-year period to the Cikapundung River. The methods and results of this study are applicable towards improving the control mechanisms of river rejuvenation from plastic leakage by addressing proper management in concentrated locations.

How Have Global River Widths Changed Over Time?

Changes in a river's width reflect natural and anthropogenic impacts on local and upstream/downstream hydraulic and hydrologic processes. Temporal variation of river width also impacts biogeochemical exchange and reflects geomorphologic evolution. However, while global maps of mean river width and dynamic water surface extent exist, there is currently no standardized global assessment of river widths that documents changes over time. Therefore, we made repeated width measurements from Landsat images for all rivers wider than 90 m collected from 1984 to 2020 (named Global LOng-term river Width, GLOW), which consists of ∼1.2 billion cross-sectional river width measurements, with an average of 3,000 width measurements per 10-km reach. With GLOW, we investigated the temporal variations of global river width, quantified by the interquartile range (IQR) and temporal trend. We found that 85% of global rivers have a width IQR <150 m. We also found that 37% of global river segments show significant temporal trends in width over the past 37 years, and this number is higher (46%) for human-regulated rivers. Further, we leveraged machine learning to identify the most important factors explaining river width variations and revealed that these driving factors are significantly different between free-flowing and non-free-flowing rivers. Specifically, the most important factor driving temporal variations in river width is the climate for free-flowing rivers, and is soil condition for human-impacted rivers. Finally, we anticipate that this study and the public release of GLOW will improve the understanding of river dynamics and catalyze additional interdisciplinary studies.

Deriving River Discharge Using Remotely Sensed Water Surface Characteristics and Satellite Altimetry in the Mississippi River Basin

River discharges are critical for understanding hydrologic and ecological systems, yet in situ data are limited in many regions of the world. While approximating river discharge using satellite-derived water surface characteristics is possible, the key challenges are unknown channel bathymetry and roughness. Here, we present an application for merging mean river-reach characteristics and time-varying altimetry measurements to estimate river discharge for sites within the Mississippi River Basin (USA). This project leverages the Surface Water and Ocean Topography (SWOT) River Database (SWORD) for approximating mean river-reach widths and slopes and altimetry data from JASON-2/3 (2008-Present) and Sentinel-3A/B (2015-Present) obtained from the Hydroweb Theia virtual stations. River discharge is calculated using Manning’s Equation, with optimized parameters for surface roughness, bottom elevation, and channel shape determined using the Kling–Gupta Efficiency (KGE). The results of this study indicate the use of optimized characteristics return 87% of sites with KGE > −0.41, which indicates that the approach provides discharges that outperform using the mean discharge. The use of precipitation to approximate missing flows not observed by satellites results in 66% of sites with KGE > −0.41, while the use of TWSA results in 65% of sites with KGE > −0.41. Future research will focus on extending this application for all available sites in the United States, as well as trying to understand how climate and landscape factors (e.g., precipitation, temperature, soil moisture, landcover) relate to river and watershed characteristics.

Financing national scale energy projects in developing countries - An economy-wide evaluation of Ghana's Bui Dam

Large energy infrastructure can imply special financing arrangements between governments in developing economies and investors or lenders. These arrangements can lead to economy-wide and sector-specific impacts which need to be considered in the project economic evaluation. By considering the case of the Bui Dam in Ghana, we use a macroeconomic approach to determine how the economic performance of critical energy infrastructure manifests during the construction, financing and operation phases. The analysis uses an integrated modelling framework that combines a Computable General Equilibrium (CGE) model of Ghana with a water balance model of the Lower Volta River Basin. The results highlight the importance of including indirect and induced effects, in addition to the direct effects from project operation, as they influence the scale and temporal evolution of the economic impacts. The collateral from the infrastructure loan agreement consisting in cocoa exports to China nearly doubles the project's positive GDP impact and has a significant multiplier effect over urban and rural household income compared to a standard commercial loan. We finish with a discussion of how the proposed investment-oriented modelling framework can contribute to ex-ante strategic assessments of proposed energy infrastructure in developing countries.

Hydrological Impact of the New ECMWF Multi-Layer Snow Scheme

The representation of snow is a crucial aspect of land-surface modelling, as it has a strong influence on energy and water balances. Snow schemes with multiple layers have been shown to better describe the snowpack evolution and bring improvements to soil freezing and some hydrological processes. In this paper, the wider hydrological impact of the multi-layer snow scheme, implemented in the ECLand model, was analyzed globally on hundreds of catchments. ERA5-forced reanalysis simulations of ECLand were coupled to CaMa-Flood, as the hydrodynamic model to produce river discharge. Different sensitivity experiments were conducted to evaluate the impact of the ECLand snow and soil freezing scheme changes on the terrestrial hydrological processes, with particular focus on permafrost. It was found that the default multi-layer snow scheme can generally improve the river discharge simulation, with the exception of permafrost catchments, where snowmelt-driven floods are largely underestimated, due to the lack of surface runoff. It was also found that appropriate changes in the snow vertical discretization, destructive metamorphism, snow-soil thermal conductivity and soil freeze temperature could lead to large river discharge improvements in permafrost by adjusting the evolution of soil temperature, infiltration and the partitioning between surface and subsurface runoff.

Assessing placement bias of the global river gauge network

Knowing where and when rivers flow is paramount to managing freshwater ecosystems. Yet stream gauging stations are distributed sparsely across rivers globally and may not capture the diversity of fluvial network properties and anthropogenic influences. Here we evaluate the placement bias of a global stream gauge dataset on its representation of socioecological, hydrologic, climatic and physiographic diversity of rivers. We find that gauges are located disproportionally in large, perennial rivers draining more human-occupied watersheds. Gauges are sparsely distributed in protected areas and rivers characterized by non-perennial flow regimes, both of which are critical to freshwater conservation and water security concerns. Disparities between the geography of the global gauging network and the broad diversity of streams and rivers weakens our ability to understand critical hydrologic processes and make informed water-management and policy decisions. Our findings underscore the need to address current gauge placement biases by investing in and prioritizing the installation of new gauging stations, embracing alternative water-monitoring strategies, advancing innovation in hydrologic modelling, and increasing accessibility of local and regional gauging data to support human responses to water challenges, both today and in the future.

High-quality reconstruction of China's natural streamflow

Reconstruction of natural streamflow is fundamental to the sustainable management of water resources. In China, previous reconstructions from sparse and poor-quality gauge measurements have led to large biases in simulation of the interannual and seasonal variability of natural flows. Here we use a well-trained and tested land surface model coupled to a routing model with flow direction correction to reconstruct the first high-quality gauge-based natural streamflow dataset for China, covering all its 330 catchments during the period from 1961 to 2018. A stronger positive linear relationship holds between upstream routing cells and drainage areas, after flow direction correction to 330 catchments. We also introduce a parameter-uncertainty analysis framework including sensitivity analysis, optimization, and regionalization, which further minimizes biases between modeled and inferred natural streamflow from natural or near-natural gauges. The resulting behavior of the natural hydrological system is represented properly by the model which achieves high skill metric values of the monthly streamflow, with about 83% of the 330 catchments having Nash-Sutcliffe efficiency coefficient (NSE) > 0.7, and about 56% of the 330 catchments having Kling-Gupta efficiency coefficient (KGE) > 0.7. The proposed construction scheme has important implications for similar simulation studies in other regions, and the developed low bias long-term national datasets by statistical postprocessing should be useful in supporting river management activities in China.

The importance of hydrology in routing terrestrial carbon to the atmosphere via global streams and rivers

Stream/river carbon dioxide (CO₂) emission has significant spatial and seasonal variations critical for understanding its macroecosystem controls and plumbing of the terrestrial carbon budget. We relied on direct fluvial CO₂ partial pressure measurements and seasonally varying gas transfer velocity and river network surface area estimates to resolve reach-level seasonal variations of the flux at the global scale. The percentage of terrestrial primary production (GPP) shunted into rivers that ultimately contributes to CO₂ evasion increases with discharge across regions, due to a stronger response in fluvial CO₂ evasion to discharge than GPP. This highlights the importance of hydrology, in particular water throughput, in terrestrial–fluvial carbon transfers and the need to account for this effect in plumbing the terrestrial carbon budget.

The magnitude of stream and river carbon dioxide (CO₂) emission is affected by seasonal changes in watershed biogeochemistry and hydrology. Global estimates of this flux are, however, uncertain, relying on calculated values for CO₂ and lacking spatial accuracy or seasonal variations critical for understanding macroecosystem controls of the flux. Here, we compiled 5,910 direct measurements of fluvial CO₂ partial pressure and modeled them against watershed properties to resolve reach-scale monthly variations of the flux. The direct measurements were then combined with seasonally resolved gas transfer velocity and river surface area estimates from a recent global hydrography dataset to constrain the flux at the monthly scale. Globally, fluvial CO₂ emission varies between 112 and 209 Tg of carbon per month. The monthly flux varies much more in Arctic and northern temperate rivers than in tropical and southern temperate rivers (coefficient of variation: 46 to 95 vs. 6 to 12%). Annual fluvial CO₂ emission to terrestrial gross primary production (GPP) ratio is highly variable across regions, ranging from negligible (<0.2%) to 18%. Nonlinear regressions suggest a saturating increase in GPP and a nonsaturating, steeper increase in fluvial CO₂ emission with discharge across regions, which leads to higher percentages of GPP being shunted into rivers for evasion in wetter regions. This highlights the importance of hydrology, in particular water throughput, in routing terrestrial carbon to the atmosphere via the global drainage networks. Our results suggest the need to account for the differential hydrological responses of terrestrial–atmospheric vs. fluvial–atmospheric carbon exchanges in plumbing the terrestrial carbon budget.

Recent changes to Arctic river discharge

Arctic rivers drain ~15% of the global land surface and significantly influence local communities and economies, freshwater and marine ecosystems, and global climate. However, trusted and public knowledge of pan-Arctic rivers is inadequate, especially for small rivers and across Eurasia, inhibiting understanding of the Arctic response to climate change. Here, we calculate daily streamflow in 486,493 pan-Arctic river reaches from 1984-2018 by assimilating 9.18 million river discharge estimates made from 155,710 satellite images into hydrologic model simulations. We reveal larger and more heterogenous total water export (3-17% greater) and water export acceleration (factor of 1.2-3.3 larger) than previously reported, with substantial differences across basins, ecoregions, stream orders, human regulation, and permafrost regimes. We also find significant changes in the spring freshet and summer stream intermittency. Ultimately, our results represent an updated, publicly available, and more accurate daily understanding of Arctic rivers uniquely enabled by recent advances in hydrologic modeling and remote sensing.

A near-global, high resolution land surface parameter dataset for the variable infiltration capacity model

Hydrologic models predict the spatial and temporal distribution of water and energy at the land surface. Currently, parameter availability limits global-scale hydrologic modelling to very coarse resolution, hindering researchers from resolving fine-scale variability. With the aim of addressing this problem, we present a set of globally consistent soil and vegetation parameters for the Variable Infiltration Capacity (VIC) model at 1/16° resolution (approximately 6 km at the equator), with spatial coverage from 60°S to 85°N. Soil parameters derived from interpolated soil profiles and vegetation parameters estimated from space-based MODIS measurements have been compiled into input files for both the Classic and Image drivers of the VIC model, version 5. Geographical subsetting codes are provided, as well. Our dataset provides all necessary land surface parameters to run the VIC model at regional to global scale. We evaluate VICGlobal's ability to simulate the water balance in the Upper Colorado River basin and 12 smaller basins in the CONUS, and their ability to simulate the radiation budget at six SURFRAD stations in the CONUS.

Global riverine nitrous oxide emissions: The role of small streams and large rivers

Nitrous oxide, N₂O, is the leading cause of stratospheric ozone depletion and one of the most potent greenhouse gases (GHG). Its concentration in the atmosphere has been rapidly increasing since the green revolution in the 1950s and 1960s. Riverine systems have been suggested to be an important source of N₂O, although their quantitative contribution has been estimated with poor precision, ranging between 32.2 and 2100 GgN₂O - N/yr. Here, we quantify reach scale N₂O emissions by integrating a data-driven machine learning model with a physically-based upscaling model. The application of this hybrid modeling approach reveals that small streams (those with widths less than 10 m) are the primary sources of riverine N₂O emissions to the atmosphere. They contribute nearly 36 GgN₂O - N/yr; almost 50% of the entire N₂O emissions from riverine systems (72.8 Gg₂O - N/yr), although they account for only 13% of the total riverine surface area worldwide. Large rivers (widths wider than 175 m), such as the main stems of the Amazon River (~ 6 GgN₂O - N/yr), the Mississippi River (~ 2 GgN₂O - N/yr), the Congo River (~ 1 GgN₂O - N/yr) and the Yang Tze River (~ 0.7 GgN₂O - N/yr), only contribute 26% of global N₂O emissions, which primarily originate from their water column. This study identifies, for the first time, near-global N₂O emission and NO₃ removal hot spots within watersheds and thus can aid the development of local- to global-scale management and mitigation strategies for riverine systems with respect to N₂O emissions. The presented framework can be extended to quantified biogeochemical, besides N₂O emissions, processes at the global scale.

Satellite Flood Inundation Assessment and Forecast Using SMAP and Landsat

The capability and synergistic use of multisource satellite observations for flood monitoring and forecasts is crucial for improving disaster preparedness and mitigation. Here, surface fractional water cover (FW) retrievals derived from Soil Moisture Active Passive (SMAP) L-band (1.4 GHz) brightness temperatures were used for flood assessment over southeast Africa during the Cyclone Idai event. We then focused on five subcatchments of the Pungwe basin and developed a machine learning based approach with the support of Google Earth Engine for daily (24-h) forecasting of FW and 30-m inundation downscaling and mapping. The Classification and Regression Trees model was selected and trained using retrievals derived from SMAP and Landsat coupled with rainfall forecasts from the NOAA Global Forecast System. Independent validation showed that FW predictions over randomly selected dates are highly correlated (R = 0.87) with the Landsat observations. The forecast results captured the flood temporal dynamics from the Idai event; and the associated 30-m downscaling results showed inundation spatial patterns consistent with independent satellite synthetic aperture radar observations. The data-driven approach provides new capacity for flood monitoring and forecasts leveraging synergistic satellite observations and big data analysis, which is particularly valuable for data sparse regions.

Spatio-Temporal Distribution of Environmental Health Risk of Heavy Metals in Industrial Wastewater of China during 1999-2018

In recent decades, environmental health risk caused by heavy metals in industrial wastewater (EHR-IHM) has become a serious issue globally, especially for China. Given the spatial difference of heavy metal emissions, hydrogeography, population distribution, etc., it is essential to estimate China's EHR-IHM from a high-resolution perspective. Based on the framework of USEtox, this study constructs an environmental health risk assessment method for heavy metals discharged from industrial wastewater by coupling the Pollutant Accumulation Model (PAM). This method also considers the process of heavy metal flows between upstream and downstream areas. Based on this constructed method, we investigate the spatio-temporal distribution of EHR-IHM of As, Cd, Cr(VI), Hg, and Pb in China from 1999 to 2018. Results showed that the EHR-IHM in China increased rapidly during 1999-2007 and decreased gradually during 2007-2018, with the highest Damage Level (DL) of 6.8 × 10⁴ disability-adjusted life years (DALY). As and Cr(VI) were the major heavy metal pollutants, which induced 58.9-70.6% and 23.9-36.2% of the total EHR-IHM, respectively. Intake of aquatic products was the dominant exposure route, accounting for over 84.1% of national EHR-IHM, followed by drinking water intake, accounting for 9.5-15.8%. Regarding spatial distribution, the regions with high EHR-IHM are mainly distributed in the middle-lower reaches of the Yangtze River, southeast coastal cities, Bohai Rim, etc.

Analyzing the Suitability of Remotely Sensed ET for Calibrating a Watershed Model of a Mediterranean Montane Forest

The ability to spatially characterize runoff generation and forest health depends partly on the accuracy and resolution of evapotranspiration (ET) simulated by numerical models. A possible strategy to increase the accuracy and resolution of numerically modeled ET is the use of remotely sensed ET products as an observational basis for parameter estimation (model calibration) of those numerical models. However, the extent to which that calibration strategy leads to a realistic representation of ET, relative to ground conditions, is not well understood. We examined this by comparing the spatiotemporal accuracy of ET from a remote sensing product, MODIS MOD16A2, to that from a watershed model (SWAT) calibrated to flow measured at an outlet streamgage. We examined this in the upper Kings River watershed (3999 km2) of California’s Sierra Nevada, a snow-influenced watershed in a Mediterranean climate. We assessed ET accuracies against observations from three eddy-covariance flux towers at elevations of 1160–2700 m. The accuracy of ET from the stream-calibrated watershed model surpassed that of MODIS in terms of Nash-Sutcliffe efficiency (+0.36 versus −0.43) and error in elevational trend (+7.7% versus +81%). These results indicate that for this particular experiment, an outlet streamgage would provide a more effective observational basis than remotely sensed ET product for watershed-model parameter estimation. Based on analysis of ET-weather relationships, the relatively large errors we found in MODIS ET may be related to weather-based corrections to water limitation not representative of the hydrology of this snow-influenced, Mediterranean-climate area.

A new vector-based global river network dataset accounting for variable drainage density

Spatial variability of river network drainage density (Dd) is a key feature of river systems, yet few existing global hydrography datasets have properly accounted for it. Here, we present a new vector-based global hydrography that reasonably estimates the spatial variability of Dd worldwide. It is built by delineating channels from the latest 90-m Multi-Error-Removed Improved Terrain (MERIT) digital elevation model and flow direction/accumulation. A machine learning approach is developed to estimate Dd based on the global watershed-level climatic, topographic, hydrologic, and geologic conditions, where relationships between hydroclimate factors and Dd are trained using the high-quality National Hydrography Dataset Plus (NHDPlusV2) data. By benchmarking our dataset against HydroSHEDS and several regional hydrography datasets, we show the new river flowlines are in much better agreement with Landsat-derived centerlines, and improved Dd patterns of river networks (totaling ~75 million kilometers in length) are obtained. Basins and estimates of intermittent stream fraction are also delineated to support water resources management. This new dataset (MERIT Hydro-Vector) should enable full global modeling of river system processes at fine spatial resolutions.

A Model Tree Generator (MTG) Framework for Simulating Hydrologic Systems: Application to Reservoir Routing

Data-driven algorithms have been widely used as effective tools to mimic hydrologic systems. Unlike black-box models, decision tree algorithms offer transparent representations of systems and reveal useful information about the underlying process. A popular class of decision tree models is model tree (MT), which is designed for predicting continuous variables. Most MT algorithms employ an exhaustive search mechanism and a pre-defined splitting criterion to generate a piecewise linear model. However, this approach is computationally intensive, and the selection of the splitting criterion can significantly affect the performance of the generated model. These drawbacks can limit the application of MTs to large datasets. To overcome these shortcomings, a new flexible Model Tree Generator (MTG) framework is introduced here. MTG is equipped with several modules to provide a flexible, efficient, and effective tool for generating MTs. The application of the algorithm is demonstrated through simulation of controlled discharge from several reservoirs across the Contiguous United States (CONUS).

Estimating River Sediment Discharge in the Upper Mississippi River Using Landsat Imagery

With the decline of operational river gauges monitoring sediments, a viable means of quantifying sediment transport is needed. In this study, we address this issue by applying relationships between hydraulic geometry of river channels, water discharge, water-leaving surface reflectance (SR), and suspended sediment concentration (SSC) to quantify sediment discharge with the aid of space-based observations. We examined 5490 Landsat scenes to estimate water discharge, SSC, and sediment discharge for the period from 1984 to 2017 at nine gauging sites along the Upper Mississippi River. We used recent advances in remote sensing of fluvial systems, such as automated river width extraction, Bayesian discharge inference with at-many-stations hydraulic geometry (AMHG), and SSC-SR regression models. With 621 Landsat scenes available from all the gauging sites, the results showed that the water discharge and SSC retrieval from Landsat imagery can yield reasonable sediment discharge estimates along the Upper Mississippi River. An overall relative bias of −25.4, mean absolute error (MAE) of 6.24 × 104 tonne/day, relative root mean square error (RRMSE) of 1.21, and Nash–Sutcliffe Efficiency (NSE) of 0.49 were obtained for the sediment discharge estimation. Based on these statistical metrics, we identified three of the nine gauging sites (St. Louis, MO; Chester, IL; and Thebes, IL), which were in the downstream portion of the river, to be the best locations for estimating water and sediment discharge using Landsat imagery.

The Reliability of Global Remote Sensing Evapotranspiration Products over Amazon

As a key component of terrestrial water cycle, evapotranspiration (ET), specifically over the Amazon River basin, is of high scientific significance. However, due to the sparse observation network and relatively short observational period of eddy covariance data, large uncertainties remain in the spatial-temporal characteristics of ET over the Amazon. Recently, a great number of long-term global remotely sensed ET products have been developed to fill the observation gap. However, the reliabilities of these global ET products over the Amazon are unknown. In this study, we assessed the consistency of the magnitude, trend and spatial pattern of Amazon ET among five global remotely sensed ET reconstructions. The magnitudes of these products are similar but the long-term trends from 1982 to 2011 are completely divergent. Validation from the eddy covariance data and water balance method proves a better performance of a product grounded on local measurements, highlighting the importance of local measurements in the ET reconstruction. We also examined four hypotheses dealing with the response of ET to brightening, warming, greening and deforestation, which shows that in general, these ET products respond better to warming and greening than to brightening and deforestation. This large uncertainty highlights the need for future studies focusing on ET issues over the Amazon.

Timing of Landsat Overpasses Effectively Captures Flow Conditions of Large Rivers

Satellites provide a temporally discontinuous record of hydrological conditions along Earth’s rivers (e.g., river width, height, water quality). The degree to which archived satellite data effectively capture the overall population of river flow frequency is unknown. Here, we use the entire archives of Landsat 5, 7, and 8 to determine when a cloud-free image is available over the United States Geological Survey (USGS) river gauges located on Landsat-observable rivers. We compare the flow frequency distribution derived from the daily gauge record to the flow frequency distribution derived from ideally sampling gauged discharge based on the timing of cloud-free Landsat overpasses. Examining the patterns of flow frequency across multiple gauges, we find that there is not a statistically significant difference between the flow frequency distribution associated with observations contained within the Landsat archive and the flow frequency distribution derived from the daily gauge data (α = 0.05), except for hydrological extremes like maximum and minimum flow. At individual gauges, we find that Landsat observations span a wide range of hydrological conditions (97% of total flow variability observed in 90% of the study gauges) but the degree to which the Landsat sample can represent flow frequency distribution varies from location to location and depends on sample size. The results of this study indicate that the Landsat archive is, on average, representative of the temporal frequencies of hydrological conditions present along Earth’s large rivers with broad utility for hydrological, ecologic and biogeochemical evaluations of river systems.

Aridity Trends in Central America: A Spatial Correlation Analysis

Trend analyses are common in several types of climate change studies. In many cases, finding evidence that the trends are different from zero in hydroclimate variables is of particular interest. However, when estimating the confidence interval of a set of hydroclimate stations or gridded data the spatial correlation between can affect the significance assessment using for example traditional non-parametric and parametric methods. For this reason, Monte Carlo simulations are needed in order to generate maps of corrected trend significance. In this article, we determined the significance of trends in aridity, modeled runoff using the Variable Infiltration Capacity Macroscale Hydrological model, Hagreaves potential evapotranspiration (PET) and near-surface temperature in Central America. Linear-regression models were fitted considering that the predictor variable is the time variable (years from 1970 to 1999) and predictand variable corresponds to each of the previously mentioned hydroclimate variables. In order to establish if the temporal trends were significantly different from zero, a Mann Kendall and a Monte Carlo test were used. The spatial correlation was calculated first to correct the variance of each trend. It was assumed in this case that the trends form a spatial stochastic process that can be modeled as such. Results show that the analysis considering the spatial correlation proposed here can be used for identifying those extreme trends. However, a set of variables with strong spatial correlation such as temperature can have robust and widespread significant trends assuming independence, but the vast majority of the stations can still fail the Monte Carlo test. We must be vigilant of the statistically robust changes in key primary parameters such as temperature and precipitation, which are the driving sources of hydrological alterations that may affect social and environmental systems in the future.

Remote Sensing of River Discharge: A Review and a Framing for the Discipline

Remote sensing of river discharge (RSQ) is a burgeoning field rife with innovation. This innovation has resulted in a highly non-cohesive subfield of hydrology advancing at a rapid pace, and as a result misconceptions, mis-citations, and confusion are apparent among authors, readers, editors, and reviewers. While the intellectually diverse subfield of RSQ practitioners can parse this confusion, the broader hydrology community views RSQ as a monolith and such confusion can be damaging. RSQ has not been comprehensively summarized over the past decade, and we believe that a summary of the recent literature has a potential to provide clarity to practitioners and general hydrologists alike. Therefore, we here summarize a broad swath of the literature, and find after our reading that the most appropriate way to summarize this literature is first by application area (into methods appropriate for gauged, semi-gauged, regionally gauged, politically ungauged, and totally ungauged basins) and next by methodology. We do not find categorizing by sensor useful, and everything from un-crewed aerial vehicles (UAVs) to satellites are considered here. Perhaps the most cogent theme to emerge from our reading is the need for context. All RSQ is employed in the service of furthering hydrologic understanding, and we argue that nearly all RSQ is useful in this pursuit provided it is properly contextualized. We argue that if authors place each new work into the correct application context, much confusion can be avoided, and we suggest a framework for such context here. Specifically, we define which RSQ techniques are and are not appropriate for ungauged basins, and further define what it means to be ‘ungauged’ in the context of RSQ. We also include political and economic realities of RSQ, as the objective of the field is sometimes to provide data purposefully cloistered by specific political decisions. This framing can enable RSQ to respond to hydrology at large with confidence and cohesion even in the face of methodological and application diversity evident within the literature. Finally, we embrace the intellectual diversity of RSQ and suggest the field is best served by a continuation of methodological proliferation rather than by a move toward orthodoxy and standardization.

A global streamflow reanalysis for 1980-2018

Global and continental scale hydrological reanalysis datasets receive growing attention due to their increasing number of applications, ranging from water resources management, climate change studies, water related hazards and policy support. Until recently, their use was mostly limited to qualitative assessments, due to their coarse spatial and temporal resolution, large uncertainty and bias in the model output, and limited extent of the dataset in space and time. This research reports on the setup of a gridded hydrological model with quasi-global coverage, able to reproduce a seamless 39-year streamflow simulation in all world's medium to large river basins. The model was calibrated at 1226 river sections with a total drainage area of 51 million km2 within 66 countries, using ECMWF's latest atmospheric reanalysis ERA5. A performance assessment revealed large improvements in reproducing past discharge observations, in comparison to the calibration used in the current operational setup of the hydrological model as part of the Copernicus - Global Flood Awareness System (GloFAS, www.globalfloods.eu), with median scores of Kling-Gupta Efficiency KGE = 0.67 and correlation r = 0.8. The simulation bias was also dramatically reduced and narrowed around zero, with more than 60% of stations showing percent bias within ±20%. Pronounced regional differences in the simulation results remain, pointing out the need for detailed investigation of the hydrological processes in specific regions, including parts of Africa and South Asia. In addition, observed discharges with high data quality is key to achieving skillful model output. The new calibrated model will become part of the operational runs of GloFAS in the next system release foreseen for Spring 2020, together with a near real time extension of the streamflow reanalysis.