Regional regression equations were calculated in Puerto Rico with generalized least squares techniques to estimate flood frequency statistics at ungaged locations using drainage area as the only explanatory variable. The island was divided into 2 regions to minimize residuals. The region division that resulted in lower and more balanced residuals runs primarily north-south near the center of the island, mostly along an 8-digit hydrologic unit code (HUC8) boundary. The division line runs through a HUC8 polygon on the southern end of the island, but care was taken to include entire watersheds and consideration was given where hydrologic and physiographic properties differed. This data release includes geographic information system files that define the polygons for both regions.

Generalized least squares (GLS) techniques were used with at-site statistics, from PeakFQ, and drainage area for streamgages in Puerto Rico to calculate regional regression equations for ungaged locations. The GLS approach assigns weights based on uncertainty of the observations (ie. record length and variance) and accounts for correlated streamflows and time-sampling errors (Farmer and others, 2019). Two separate regions were defined for the use of these equations to minimize residuals, and performance metrics from WREG were analyzed to determine drainage area as the only explanatory variable provided the best results. Values for the correlation smoothing function in region 1 were alpha equal to 0.006 and theta of 0.96 and for region 2 alpha equal to 0.004 and theta of 0.94. This data release includes (1) a PeakFQ folder that contains input files, both .EXP and .PRT (output from PeakFQ), for use in WREG; (2) a SiteInfo folder for use in WREG that contains information about each site used, separated by region; and (3) a Results folder showing the WREG output of each annual exceedance probability flow, separated by regression region where applicable. References Cited: Farmer, W.H., Kiang, J.E., Feaster, T.D., and Eng, K., 2019, Regionalization of surface-water statistics using multiple linear regression: U.S. Geological Survey Techniques and Methods, book 4, chap. A12, 40 p., Accessed August 4, 2020 at https://doi.org/10.3133/tm4A12.

Generalized least squares (GLS) techniques were used with at-site statistics, from PeakFQ, and drainage area for streamgages in Puerto Rico to calculate regional regression equations for ungaged locations. The GLS approach assigns weights based on uncertainty of the observations (ie. record length and variance) and accounts for correlated streamflows and time-sampling errors (Farmer and others, 2019). Two separate regions were defined for the use of these equations to minimize residuals, and performance metrics from WREG were analyzed to determine drainage area as the only explanatory variable provided the best results. Values for the correlation smoothing function in region 1 were alpha equal to 0.006 and theta of 0.96 and for region 2 alpha equal to 0.004 and theta of 0.94. This data release includes (1) a PeakFQ folder that contains input files, both .EXP and .PRT (output from PeakFQ), for use in WREG; (2) a SiteInfo folder for use in WREG that contains information about each site used, separated by region; and (3) a Results folder showing the WREG output of each annual exceedance probability flow, separated by regression region where applicable. References Cited: Farmer, W.H., Kiang, J.E., Feaster, T.D., and Eng, K., 2019, Regionalization of surface-water statistics using multiple linear regression: U.S. Geological Survey Techniques and Methods, book 4, chap. A12, 40 p., Accessed August 4, 2020 at https://doi.org/10.3133/tm4A12.

Flood region boundary for region 4.

Flood region boundary for region 4.

Abstract: This file contains Hydrologic Unit (HU) polygon boundaries for the United States, Puerto Rico, and the U.S. Virgin Islands. The data is a seamless National representation of HU boundaries from 2 to 14 digits compiled from U.S. Geological Survey (USGS) National Hydrography Dataset (NHD) and U.S. Department of Agriculture (USDA) National Resources Conservation Service (NRCS) Watershed Boundary Dataset (WBD) sources. Purpose: This data is intended primarily for geographic display and analysis of regional and national data, and can also be used for illustration purposes at intermediate or small scales (1:250,000 to 1:2,000,000). See https://apps.nationalmap.gov/help/ for assistance with The National Map viewer, download client, services, or metadata.

Flood region boundary for region 2.

Estimates of the magnitude of peak-flows were updated for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent chance exceedance levels for 91 rural, unregulated streamgaging stations on the main island of Puerto Rico. These stations required 10 or more years of annual peak-flow record, using data to 2017, for inclusion in the study. The magnitude and frequency of floods at selected streamgages in Puerto Rico were estimated using the U.S. Geological Survey PeakFQ program and updated methods outlined in Bulletin 17C (England and others, 2018). Regional regression equations were calculated to estimate flood frequency statistics at ungaged locations using selected basin characteristics as explanatory variables. These variables were determined from digital spatial datasets and geographic information systems using the most recent data available, as referenced in the U.S. Geological Survey web application, StreamStats, and published in Kolb and Ryan (2021). A generalized least squares procedure in the U.S. Geological Survey program, WREG, was used to account for cross-correlation of sites and develop the final regional regression equations using drainage area as the only explanatory variable. Two separate regions were defined for regression equation use in this study to minimize residuals. NOTE: All of the data in the previous version can be found in version 1.1. References Cited: England J.F., Jr., Cohn, T.A., Faber, B.A., Stedinger,J.R., Thomas,W.O.,Jr., Veilleux,A.G., Kiang,J.E., and Mason,R.R., Jr., 2018, Guidelines for determining flood flow frequency —Bulletin 17C: U.S. Geological Survey Techniques and Methods, book 4, chap. B5, 148p., https://doi.org/10.3133/tm4B5. Kolb, K.R., and Ryan, P.J., 2021, Basin Characteristic Rasters for Puerto Rico StreamStats, 2021: U.S. Geological Survey data release, https://doi.org/10.5066/P9HK9SSQ.

Estimates of the magnitude of peak-flows were updated for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent chance exceedance levels for 91 rural, unregulated streamgaging stations on the main island of Puerto Rico. These stations required 10 or more years of annual peak-flow record, using data to 2017, for inclusion in the study. The magnitude and frequency of floods at selected streamgages in Puerto Rico were estimated using the U.S. Geological Survey PeakFQ program and updated methods outlined in Bulletin 17C (England and others, 2018). Regional regression equations were calculated to estimate flood frequency statistics at ungaged locations using selected basin characteristics as explanatory variables. These variables were determined from digital spatial datasets and geographic information systems using the most recent data available, as referenced in the U.S. Geological Survey web application, StreamStats, and published in Kolb and Ryan (2021). A generalized least squares procedure in the U.S. Geological Survey program, WREG, was used to account for cross-correlation of sites and develop the final regional regression equations using drainage area as the only explanatory variable. Two separate regions were defined for regression equation use in this study to minimize residuals. NOTE: All of the data in the previous version can be found in version 1.1. References Cited: England J.F., Jr., Cohn, T.A., Faber, B.A., Stedinger,J.R., Thomas,W.O.,Jr., Veilleux,A.G., Kiang,J.E., and Mason,R.R., Jr., 2018, Guidelines for determining flood flow frequency —Bulletin 17C: U.S. Geological Survey Techniques and Methods, book 4, chap. B5, 148p., https://doi.org/10.3133/tm4B5. Kolb, K.R., and Ryan, P.J., 2021, Basin Characteristic Rasters for Puerto Rico StreamStats, 2021: U.S. Geological Survey data release, https://doi.org/10.5066/P9HK9SSQ.

Flow-duration statistics at the 99th, 98th, 95th, 90th, 80th, 70th, 60th, and 50th percent exceedance probabilities and annual n-day low-flow statistics for the 1-, 7-, 14-, and 30-day mean low flows with 2-year (0.5 nonexceedance probability), 5-year (0.2 nonexceedance probability), and 10-year (0.1 nonexceedance probability) recurrence intervals were computed for 28 selected streamflow gaging stations in Puerto Rico. The 28 selected streamflow gaging stations were required to have 10 or more years of daily mean streamflow data through water year 2018. The flow-duration statistics and n-day low-flow frequencies were computed using the U.S. Geological Survey program, SWToolbox. Regional regression equations were developed to estimate flow-duration statistics and n-day low-flow frequencies at ungaged stream locations using selected basin characteristics as explanatory variables. These variables were determined from digital spatial datasets and geographic information systems using the most recent data available, as referenced in the U.S. Geological Survey web application, StreamStats, and published in Kolb and Ryan (2021). An ordinary least-squares procedure in R Studio was used to develop the final regional flow-duration regression equations using drainage area, mean total annual reference evapotranspiration, and minimum basin elevation as the explanatory variables. A generalized least squares procedure in the U.S. Geological Survey program, WREG, was used to account for cross-correlation of sites and develop the final regional n-day low-flow frequency regression equations using drainage area, mean total annual reference evapotranspiration, and minimum basin elevation as the explanatory variables. This data release includes two child pages: Puerto Rico Flow-Duration Regression Files and Puerto Rico N-day Low-Flow Regression Files, a BasinCharacteristics.csv file that contains 47 basin and climatic characteristics considered in the analyses, a BasinCharacteristics_corrrelation_charts folder that contains .pdf files showing correlation matrices, a R_regsubsets_output folder that contains .txt and .pdf files showing results of the "regsubsets" analyses, a Trend_statistics_nday_lowflow_timeseries.csv file that contains the SWToolbox Mann-Kendall tau statistics, and a NWIS_rdb_files folder that contains the .rdb files used in the analyses. References Cited: Kolb, K.R., and Ryan, P.J., 2021, Basin Characteristic Rasters for Puerto Rico StreamStats, 2021: U.S. Geological Survey data release, https://doi.org/10.5066/P9HK9SSQ.