DOI: 10.5281/zenodo.10013744
title: RRR input and output files corresponding to "Global patterns in river water storage dependent on residence time"

description: Corresponding peer-reviewed publication

This dataset corresponds to all the RRR input and output files that were used in the study reported in:



Collins, E. L., C. H. David, R. Riggs, G. H. Allen, T. M. Pavelsky, P. Lin, M. Pan, D. Yamazaki, R. K. Meentemeyer, and G. M. Sanchez (2024), Global patterns in river water storage dependent on residence time.


When making use of any of the files in this dataset, please cite both the aforementioned article and the dataset herein.

Input datasets

The dataset herein benefited from and was built upon other datasets:



MERIT-Basins (version 1.0) derived from MERIT-Hydro (version 0.7) available under a CC BY-NC-SA 4.0 license. https://www.reachhydro.org/home/params/merit-basins

GLDAS VIC Land Surface Model L4 monthly 1.0 x 1.0 degree (V2.0) available under a CC0 license. https://science.nasa.gov/spd-41/. DOI: 10.5067/ZRIHVF29X43C

GLDAS Noah Land Surface Model L4 monthly 1.0 x 1.0 degree (V2.0) available under a CC0 license. https://science.nasa.gov/spd-41/. DOI: 10.5067/QN80TO7ZHFJZ

GLDAS Catchment Land Surface Model L4 monthly 1.0 x 1.0 degree (V2.0) available under a CC0 license. https://science.nasa.gov/spd-41/. DOI: 10.5067/SGSL3LNKGJWW


Known bugs and limitations in this dataset or the associated manuscript

In the original files for MERIT-Basins (version 1.0) derived from MERIT-Hydro (version 0.7), as obtained from the pfaf_level_02 zip files, two river reaches do not have a corresponding catchment: COMID=31000001 and COMID=61000003. These two river reaches were removed from our analysis, hence bringing the total number of river reaches used from 2,938,143 to 2,938,141, the same value as the total number of catchments.

In the original files for MERIT-Basins (version 1.0) derived from MERIT-Hydro (version 0.7), as obtained from the pfaf_level_02 zip files, the attribute table for the river reach with COMID=67072432 contains only four upstream reaches (COMID=64072433, COMID=64072735, COMID=64073772, and COMID=64073786), while the reach actually has five upstream reaches on the map (those above and COMID=64073793). This was modified herein because all connectivity in this dataset was created by associating upstream and downstream nodes rather than using the attribute table.







Our monthly dataset of river gauges with 95% daily data availability for 1980–2009 and with average discharge greater than or equal to 100 m3/s contains 1,148 stations. After mapping gauges to the MERIT Hydro river network using a joint criteria of 1) within a distance of 0.05 degrees and 2) uncorrected 30-yr average simulated discharge of the river reach within one order of magnitude relative to the observed average discharge of the gauge; the number of gauges reduces to 1129. In cases of multiple gauges per reach, we selected the gauge that had the closest observed average discharge to the uncorrected 30-yr averages simulated, resulting in 1,001 gauges. The dataset was further split in calibration gauges (702) and validation gauges (299). Three specific gauges are located in places with 0 accumulation of runoff in their specific sub-basin, they are located at river reaches with ID 41003992, 41001545, and 24005996. These three gauges were removed from the overall analysis at all gauges, resulting in 998 total gauges. All three gauges are part of the 299 validation gauges, not the 702 calibration gauges. They were retained in the analysis of 299 validation gauges. 

 




Files included in this version

All monthly discharge observation files for 1148 gauges, in shapefile and CSV formats:



sites_1980-01_2009-12_100cms_095p.zip

timeseries_obs_1980-01_2009-12_100cms_095p_monthly.zip


All CSV files corresponding to the hydrography of all pfaf_level_02 regions of MERIT-Hydro v0.7 Basins v1.0 (i.e. the 61 values of ii):



coords_pfaf_ii.zip

kfac_pfaf_ii_1km_hour.zip

k_pfaf_ii_low.zip

k_pfaf_ii_nrm.zip

k_pfaf_ii_hig.zip

rapid_catchment_pfaf_ii.zip

rapid_connect_pfaf_ii.zip

rapid_coupling_pfaf_ii_GLDAS.zip

riv_bas_id_pfaf_ii_topo.zip

sort_pfaf_ii_topo.zip

xfac_pfaf_ii_0.1.zip

x_pfaf_ii_low.zip

x_pfaf_ii_nrm.zip

x_pfaf_ii_hig.zip


All netCDF4 files with monthly surface and subsurface runoff that were combined from GLDAS 2.0 and corresponding to the following land surface models for the years 1980 to 2009 (both included):



GLDAS_CLSM_M_1980-01_2009-12_utc.zip

GLDAS_NOAH_M_1980-01_2009-12_utc.zip

GLDAS_VIC_M_1980-01_2009-12_utc.zip

GLDAS_ENS_M_1980-01_2009-12_utc.zip (the 3-model ensemble average based on the three above models).


All netCDF4 files with monthly lateral inflows (based on the GLDAS data above), corresponding to each one of the 61 values of ii, the pfaf_level_02 code:



m3_riv_pfaf_ii_GLDAS_CLSM_M_1980-01_2009-12_utc.zip

m3_riv_pfaf_ii_GLDAS_NOAH_M_1980-01_2009-12_utc.zip

m3_riv_pfaf_ii_GLDAS_VIC_M_1980-01_2009-12_utc.zip

m3_riv_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc.zip


All netCDF4 files with monthly discharge (based on the m3_riv data above), corresponding to each one of the 61 values of ii, the pfaf_level_02 code: 



Qout_pfaf_ii_GLDAS_CLSM_M_1980-01_2009-12_utc.zip

Qout_pfaf_ii_GLDAS_NOAH_M_1980-01_2009-12_utc.zip

Qout_pfaf_ii_GLDAS_VIC_M_1980-01_2009-12_utc.zip

Qout_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc.zip


All netCDF4 files with monthly storage (based on the Qout_pfaf_ii_GLDAS_ENS data above), corresponding to each one of the 61 values of ii, the pfaf_level_02 code, and to each characteristic value of residence time (low, nrm, hig):



V_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc_low.zip

V_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc_nrm.zip

V_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc_hig.zip


All shapefiles with MERIT Hydro (v0.7) Basins (v1.0) river network on which the 30-year mean of lumped discharge was appended (based on the Qout_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc.zip data above):



riv_pfaf_ii_MERIT_Hydro_v07_Basins_v01_GLDAS_ENS.zip


All monthly discharge observation files for 1148 gauges, in shapefile format, with observed monthly mean appended:



sites_1980-01_2009-12_100cms_095p_meanQ.zip


All monthly discharge observation files for 1001 gauges snapped on the MERIT Hydro (v0.7) Basins (v1.0) river network, in shapefile and CSV formats, corresponding to each one of the 61 values of ii, the pfaf_level_02 code: 



sites_1980-01_2009-12_100cms_095p_meanQ_COR_pfaf_ii.zip

sites_1980-01_2009-12_100cms_095p_meanQ_COR.zip

obs_tot_id_1980-01_2009-12_100cms_095p_meanQ_COR_pfaf_ii.zip

Qobs_1980-01_2009-12_100cms_095p_meanQ_COR_pfaf_ii.zip


All monthly discharge observation files for 1001 gauges snapped on the MERIT Hydro (v0.7) Basins (v1.0) river network, separated into calibration (CAL) and validation (VAL) gauges, in shapefile and CSV formats, corresponding to each one of the 61 values of ii, the pfaf_level_02 code: 



sites_1980-01_2009-12_100cms_095p_meanQ_CAL.zip

obs_bas_id_1980-01_2009-12_100cms_095p_meanQ_CAL_pfaf_ii.zip

sites_1980-01_2009-12_100cms_095p_meanQ_VAL.zip


All netCDF4 files with monthly lateral inflows (based on the m3_riv_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc.zip data above) corrected using Long-Term Inverse Routing for calibration (CAL) gauges only and associated monthly discharge data, corresponding to each one of the 61 values of ii, the pfaf_level_02 code:



m3_riv_pfaf_ii_GLDAS_CAL_M_1980-01_2009-12_utc.zip

Qout_pfaf_ii_GLDAS_CAL_M_1980-01_2009-12_utc.zip


All netCDF4 files with monthly lateral inflows (based on the m3_riv_pfaf_ii_GLDAS_ENS_M_1980-01_2009-12_utc.zip data above) corrected using Long-Term Inverse Routing for all correction (COR) gauges only and associated monthly discharge data, corresponding to each one of the 61 values of ii, the pfaf_level_02 code:



m3_riv_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc.zip

Qout_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc.zip


All netCDF4 files with monthly storage (based on the Qout_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc.zip data above), corresponding to each one of the 61 values of ii, the pfaf_level_02 code, and to each characteristic value of residence time (low, nrm, hig):



V_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc_low.zip

V_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc_nrm.zip

V_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc_hig.zip


All shapefiles with MERIT Hydro (v0.7) Basins (v1.0) river network on which the 30-year mean of lumped discharge was appended (based on the Qout_pfaf_ii_GLDAS_COR_M_1980-01_2009-12_utc.zip data above):



riv_pfaf_ii_MERIT_Hydro_v07_Basins_v01_GLDAS_COR.zip


All shapefiles with MERIT Hydro (v0.7) Basins (v1.0) catchments that have been dissolved and of which the perimeter was extracted, corresponding to each one of the 61 values of ii, the pfaf_level_02 code; as well as global combined files:



cat_pfaf_ii_MERIT_Hydro_v07_Basins_v01_disso.zip

cat_ MERIT_Hydro_v07_Basins_v01_disso.zip

cat_pfaf_ii_MERIT_Hydro_v07_Basins_v01_perim.zip

cat_ MERIT_Hydro_v07_Basins_v01_perim.zip


All shapefiles with MERIT Hydro (v0.7) Basins (v1.0) river network retaining only those reaches that flow to the global coast, corresponding to each one of the 61 values of ii, the pfaf_level_02 code:



riv_pfaf_ii_MERIT_Hydro_v07_Basins_v01_coast.zip