Flood region boundary for region 4.

This model archive contains R source code for the Weighted-Multiple Linear Regression Program (WREG), input files, and associated output files needed to recreate regression models that are discussed in the report: Levin, S.B. and Sanocki, C.A., Methods for estimating flood magnitude and frequency for unregulated streams in Wisconsin, U.S. Geological Survey Scientific Investigations Report 2022-5118 (https://doi.org/10.3133/sir20225118). More information and instructions for running the model archive are included in the README.txt file. Information regarding the WREG program can be found at: https://water.usgs.gov/software/WREG/.

Hydrographic and Impairment Statistics (HIS) is a National Park Service (NPS) Water Resources Division (WRD) project established to track certain goals created in response to the Government Performance and Results Act of 1993 (GPRA). One water resources management goal established by the Department of the Interior under GRPA requires NPS to track the percent of its managed surface waters that are meeting Clean Water Act (CWA) water quality standards. This goal requires an accurate inventory that spatially quantifies the surface water hydrography that each bureau manages and a procedure to determine and track which waterbodies are or are not meeting water quality standards as outlined by Section 303(d) of the CWA. This project helps meet this DOI GRPA goal by inventorying and monitoring in a geographic information system for the NPS: (1) CWA 303(d) quality impaired waters and causes; and (2) hydrographic statistics based on the United States Geological Survey (USGS) National Hydrography Dataset (NHD). Hydrographic and 303(d) impairment statistics were evaluated based on a combination of 1:24,000 (NHD) and finer scale data (frequently provided by state GIS layers).

The U.S. Geological Survey (USGS), in cooperation with the Puerto Rico Environmental Quality Board, has compiled a series of geospatial datasets for Puerto Rico to be implemented into the USGS StreamStats application (https://streamstats.usgs.gov/ss/). These geospatial datasets, along with basin characteristics datasets for Puerto Rico published as a separate USGS data release (https://doi.org/10.5066/P9HK9SSQ), were used to delineate watersheds and develop the peak-flow and low-flow regression equations used by StreamStats.

The Watershed Boundary Dataset (WBD) is a comprehensive aggregated collection of hydrologic unit data consistent with the national criteria for delineation and resolution. It defines the areal extent of surface water drainage to a point except in coastal or lake front areas where there could be multiple outlets as stated by the "Federal Standards and Procedures for the National Watershed Boundary Dataset (WBD)" "Standard" (https://pubs.usgs.gov/tm/11/a3/). Watershed boundaries are determined solely upon science-based hydrologic principles, not favoring any administrative boundaries or special projects, nor particular program or agency. This dataset represents the hydrologic unit boundaries to the 12-digit for the entire United States. Some areas may also include additional subdivisions representing the 14- and 16-digit hydrologic unit (HU). At a minimum, the HUs are delineated at 1:24,000-scale in the conterminous United States, 1:25,000-scale in Hawaii, Pacific basin and the Caribbean, and 1:63,360-scale in Alaska, meeting the National Map Accuracy Standards (NMAS). Higher resolution boundaries are being developed where partners and data exist and will be incorporated back into the WBD. WBD data are delivered as a dataset of polygons and corresponding lines that define the boundary of the polygon. WBD polygon attributes include hydrologic unit codes (HUC), size (in the form of acres and square kilometers), name, downstream hydrologic unit code, type of watershed, non-contributing areas, and flow modifications. The HUC describes where the unit is in the country and the level of the unit. WBD line attributes contain the highest level of hydrologic unit for each boundary, line source information and flow modifications. For additional information on NHD, go to https://www.usgs.gov/national-hydrography.

The U.S. Geological Survey (USGS), in cooperation with the New Mexico Department of Transportation, estimated the magnitude and frequency of floods corresponding to the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities (AEPs; otherwise known as the 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year floods, respectively) for 346 selected streamgages in New Mexico and parts of Arizona, Colorado, Oklahoma, Texas, and Utah using data through water year 2020. An updated regional flood skew, -0.145, standard error 0.454, was computed for the study area. Regression equations were developed which can be used to estimate the magnitude and frequency of floods at ungaged locations on unregulated streams in the study area. The methods and results of the study are published in the parent report (Bell and others, 2022, https://doi.org/10.5066/XXXXXXXX). For the 346 selected streamgages, this dataset includes peak-flow (*.pkf) and specification (*.psf), output (*.PRT), and export (*.EXP) files from version 7.3 of USGS PeakFQ software (Veilleux and others, 2014; Flynn and others, 2006). Within PeakFQ software, the Expected Moments Algorithm (EMA) was used to conduct frequency analyses to estimate stream discharges corresponding to the 0.5, 0.2, 0.1, 0.04, 0.02, 0.01, 0.005, and 0.002 AEPs. When appropriate, the updated regional skew was used to weight the at-site skew in the frequency analyses. Results of the frequency analyses were used in generalized least-squares (GLS) regression to generate equations that predict discharges corresponding to selected AEPs at ungaged locations on streams in the study area (Bell and others, 2022).