Supplement to "Feedback between channel conveyance and flood frequency under stage–discharge rating curve uncertainty: Insights from the Pulangi River (Philippines)"

README Updated 2 March 2026 This supplement contains the data and code associated with the manuscript “Feedback between channel conveyance and flood frequency under stage–discharge rating curve uncertainty: Insights from the Pulangi River (Philippines),” which has been submitted to a journal for possible publication. Any questions or clarifications are welcome. 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 rescued the hydrometric data considered herein from the archives of the Water Projects Division (WPD) of the Philippine Department of Public Works and Highways (DPWH). Two main R files are included in this deposit: Pln_RC_upd.R, the R script file, and Pln_RC_glbenv.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 if any issues with the code arise. Three R script files are uploaded as supplements: Pln_RC_Bopha.R contains the additional analysis performed to quantify the effect of Tropical Cyclone Bopha (2012) on estimates of trends in channel capacity. Pln_RC_H9xSA.R contains the code for analyzing the sensitivity of the flood frequency discharge and return period estimates to variations in the flood-stage definition (i.e. the percentile of the daily stage distribution). PlnSgn_RC.R contains the comparative analysis of rating curve models for two gauged Philippine rivers, the Pulangi River and the Saguing River. The discussion on their differences is provided in Appendix B of the study's supplementary information. Pln_RC_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-XS-csv.csv contains data on the gauging transects: gauging ID [“GaugingID”]; date [“Date”] (YYYY-MM-DD); lateral distance from a fixed initial point [“Lat_distance”] (m); width of the gauging vertical [“Width_vert..m.”] (m); depth relative to the water surface [“Depth..m.”] (m); stage at the gauging vertical [“H_m_vert..m.”] (m); and elevation with respect to a fixed arbitrary datum at the station [“Elev..m.”] (m). 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. It has the following columns: (i) stage [“H”] (m); (ii) discharge in L/s [“Q_lps”]; and (iii) discharge in m3/s [“Q_cms”]. Sgn-DM-full.csv contains the following information on the gaugings performed on the Saguing River (2003-2021): (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 [“V_avg”] (m/s); (vi) channel width [“W_m”] (m); (vii) 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”]; (viii) gauging method ["Method"], which is either "Bridge" or "Wading"; and (ix) a column indicating whether or not there had a datum change for each gauging ["DatumChange"]. 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_RC_out_Zenodo contains only one RDS file, Pln_DM_BDRC_tour.rds, which is the saved R object for the Bayesian rating curve modeling results. Generating said object may take time if re-run from scratch, so reading in this RDS can save on computational cost.README 更新于2026年3月2日 本补充材料包含与已提交至某期刊等待发表的手稿《水位-流量关系曲线（stage–discharge rating curve）不确定性下河道输水与洪水频率的反馈：来自菲律宾普朗吉河（Pulangi River）的启示》相关的数据与代码。如有任何疑问或需要澄清之处，欢迎联系我们。 代码、数据与属性说明 本次上传包含两个目录（以粗体标注，内容详述如下），涵盖针对菲律宾 Lumayong 河段普朗吉河（Pulangi River）的案例研究中，正式分析所用的输入数据与各类输出结果。本说明等同于存档文件的综合README文件。需注意，本文所用的水文测量数据系从菲律宾公共工程与公路部（Department of Public Works and Highways, DPWH）水利项目司（Water Projects Division, WPD）的档案中抢救获取。 本次存档包含两个主要的R脚本文件：Pln_RC_upd.R，即R脚本文件；以及Pln_RC_glbenv.RData，其中存储了该脚本末次运行时生成的R对象。本脚本代码行数较多，涵盖数据整理、正式分析以及手稿配图生成等多个环节。请参照脚本大纲与注释进行操作，若代码运行出现问题，请随时与我们联系。 本次上传还包含三个补充R脚本文件： 1. Pln_RC_Bopha.R：包含用于量化热带气旋宝霞（Tropical Cyclone Bopha, 2012）对河道输水能力趋势估计影响的额外分析代码。 2. Pln_RC_H9xSA.R：包含用于分析洪水频率流量与重现期估计结果，对洪水水位定义（即日水位分布的百分位数）变化的敏感性的代码。 3. PlnSgn_RC.R：包含针对菲律宾两条测站河流——普朗吉河（Pulangi River）与萨金河（Saguing River）的水位流量关系曲线模型的对比分析，相关差异讨论详见本研究补充材料的附录B。 Pln_RC_in_Zenodo 目录包含以下CSV文件与一个子目录，这些文件均作为正式分析R脚本的输入数据： Pln-H_arch.csv：包含2004—2020年时段内每日平均水位值["H_bar_arch"]（单位：米），日期格式为YYYY-MM-DD。本文件中的水位数据已完成明显误差修正。 Pln-DM-filtered.csv：包含普朗吉河Lumayong河段1983—2020年的测流数据，各字段信息如下：(i) 日期["Date"]（格式为YYYY-MM-DD）；(ii) 水位["H"]（单位：米）；(iii) 流量单位分别为升每秒["Q_lps"]与立方米每秒["Q_cms"]；(iv) 过水面积["A_sqm"]（单位：平方米）；(v) 平均流速["Vel_mps"]（单位：米每秒）；(vi) 河道宽度["W_m"]（单位：米）；(vii) 平均水深["D_ave_m"]（单位：米）；(viii) 测量断面相对于水位标尺的位置（负值表示标尺下游，正值表示标尺上游）["XS_loc_wrt_gage"]；(ix) 最大水深["Max_Depth_m"]（单位：米）；(x) 最低河床高程["MINSBE"]（单位：米）。注：本文件仅包含通过测量位置一致性筛选步骤的测流数据。 Pln-XS-csv.csv：包含测流断面数据，各字段信息如下：测流ID["GaugingID"]；日期["Date"]（格式为YYYY-MM-DD）；固定起始点的横向距离["Lat_distance"]（单位：米）；测流垂线宽度["Width_vert..m."]（单位：米）；相对于水面的水深["Depth..m."]（单位：米）；测流垂线处水位["H_m_vert..m."]（单位：米）；相对于测站固定任意基准面的高程["Elev..m."]（单位：米）。 Pln-DM_ofcl_rating.csv：详细列出了菲律宾公共工程与公路部（DPWH）所使用的官方水位流量关系曲线的有效时段，各字段信息如下：(i) 子时段["Sub-period"]；(ii) 水位流量关系曲线["Rating"]；(iii) 有效时段开始日期["Start_date"]（格式为YYYY-MM-DD）；(iv) 有效时段结束日期["End_date"]（格式为YYYY-MM-DD）。 本次存档还包含三份对应水位流量关系表的CSV文件：Pln-DM - RatingTable_A.csv、Pln-DM - RatingTable_B.csv与Pln-DM - RatingTable_C.csv，各字段信息如下：(i) 水位["H"]（单位：米）；(ii) 流量（升每秒）["Q_lps"]；(iii) 流量（立方米每秒）["Q_cms"]。 Sgn-DM-full.csv：包含萨金河（Saguing River）2003—2021年的测流数据，各字段信息如下：(i) 日期["Date"]（格式为YYYY-MM-DD）；(ii) 水位["H"]（单位：米）；(iii) 流量单位分别为升每秒["Q_lps"]与立方米每秒["Q_cms"]；(iv) 过水面积["A_sqm"]（单位：平方米）；(v) 平均流速["V_avg"]（单位：米每秒）；(vi) 河道宽度["W_m"]（单位：米）；(vii) 测量断面相对于水位标尺的位置（负值表示标尺下游，正值表示标尺上游）["XS_loc"]；(viii) 测流方法["Method"]，可选值为"Bridge（桥测）"或"Wading（涉水测量）"；(ix) 用于标注每次测流是否发生基准面变动的字段["DatumChange"]。 Pln_Q_rated 是一个子目录，仅包含一个CSV文件Pulangi_dpwh_rtd.csv，各字段信息如下：年份["YEAR"]（格式为YYYY）；日期["DAY"]（格式为DD）；以及1—12月各月列，列中存储了官方评级或经"处理"的流量估计值（单位：升每秒）。 Pln_RC_out_Zenodo 目录仅包含一个RDS文件Pln_DM_BDRC_tour.rds，即贝叶斯水位流量关系曲线建模结果的保存R对象。若从头重新运行以生成该对象可能耗时较长，直接读取该RDS文件可有效节省计算成本。