DOI: 10.5281/zenodo.17138387
title: Supplement to "Refining satellite-derived river discharge estimates under non-stationary uncertainty in station observations: Bias correction versus space-based rating curve model recalibration"

description: This supplement contains the data and code associated with the manuscript “Refining satellite-derived river discharge estimates under non-stationary uncertainty in station observations: Bias correction versus space-based rating curve model recalibration,” which has been submitted to a journal for possible publication. Please contact Richard Cornelio (rccornelio@up.edu.ph and/or rcalayegcornelio@gmail.com) for any questions or clarifications. 

Code, Data, and Attribute Descriptions

Uploaded are two directories (indicated in bold font, with their contents detailed below) containing input data and various outputs of the analyses performed for our case study on the Pulangi River at Lumayong (Philippines). Please consider this description equivalent to an omnibus README file for the deposited files. Note that we obtained the hydrometric data considered herein from the archives of the Water Projects Division (WPD) of the Philippine Department of Public Works and Highways (DPWH).  





Three R-related files are included in this deposit: 250917_Pln_Q_DFObc.R, the R script file used in the formal analysis; 250917_Pln_Q_DFObc_gg.R, the R script file used to produce the figures in the manuscript and supplementary material; and 250917_Pln_Q_DFObc.RData, which contains stored objects generated from the script based on the last execution. The script consists of admittedly too many lines of code (e.g., for data wrangling, formal analyses, and producing figures used in the manuscript). Please be guided by the outline and the comments and kindly reach out lest issues with the code arise.




Pln_DFO_in_Zenodo contains the following CSV files and a subdirectory; these files were used as inputs to the R script for formal analyses.






Pln-H_arch.csv contains the mean daily stage values [“H_bar_arch”] (m) for every day in the 2004—2020 sub-period [“Date”] (YYYY-MM-DD). The stages in this file are already corrected for gross errors.




Pln-DM-filtered.csv contains the following information on the gaugings performed on the Pulangi at Lumayong (1983—2020): (i) date [“Date”] (YYYY-MM-DD); (ii) stage [“H”] (m); (iii) discharge in L/s [“Q_lps”] and m3/s [“Q_cms”]; (iv) wetted area [“A_sqm”] (m2); (v) mean flow velocity [“Vel_mps”] (m/s); (vi) channel width [“W_m”] (m); (vii) mean flow depth [“D_ave_m”] (m); (viii) location of the measurement cross-section with respect to the staff gauge, with negative values meaning downstream of the gauge and positive values meaning upstream of the gauge [“XS_loc_wrt_gage”]; (ix) maximum flow depth [“Max_Depth_m”] (m); and (x) minimum streambed elevation [“MINSBE”] (m). Note that this file includes only gaugings that passed our filtering step for measurement location consistency.




Pln-DM_ofcl_rating.csv details the validity periods of the official rating curves used by the DPWH. It contains the following information: (i) subperiod [“Sub-period”]; (ii) rating curve [“Rating”]; (iii) start date of the validity period [“Start_date”] (YYYY-MM-DD); and (iv) end date of the validity period (“End_date”] (YYYY-MM-DD). 




Three CSV files corresponding to the rating tables are included: Pln-DM - RatingTable_A.csv, Pln-DM - RatingTable_B.csv, and Pln-DM - RatingTable_C.csv. They have the following columns: (i) stage [“H”] (m); (ii) discharge in L/s [“Q_lps”]; and (iii) discharge in m3/s [“Q_cms”].




Pln_grid_merra2_slist.csv contains the following information on the MERRA-2 grid points bounded by the Upper Pulangi Basin extent: (i) latitude [“LAT”], (ii) longitude [“LON”] in WGS 84; and (iii) the Euclidean distance of the grid point to the basin centroid [“HubDist”].




DFO-RW_v4.5.csv contains the remotely sensed discharge data from the Dartmouth Flood Observatory’s River Watch (DFO-RW) v4.5 for Philippine satellite gauging reaches (SGRs): (i) site code and name of the SGR [“SGR”]; (ii) date [“Date”] (YYYY-MM-DD); (iii) discharge estimate in m3/s [“Q_cms”]; (iv) parameters of the original DFO-RW rating equation (i.e. the slope [“m”] and the intercept [“b”]); and (v) the dimensionless remote sensing signal for passive microwave radiometry [“HR”].




Pln_Q_rated is a subdirectory containing just one CSV file, Pulangi_dpwh_rtd.csv. It has the following columns: year [“YEAR”] (YYYY); day [“DAY”] (DD); and columns for the months of January to December, under which are the officially rated or “processed” discharge estimates in L/s.





Pln_DFO_out_Zenodo contains multiple RDS files of objects generated by the R code. These objects are computationally intensive and time-consuming to generate from scratch. To facilitate reproducibility and ease of use, we provide them here so that users can simply read them into R without rerunning the full set of scripts.