Estimation of magnitude and frequency of floods for rural, unregulated streams in and near Virginia and West Virginia

Magnitude and frequency of annual peak streamflows were computed for 813 streamgages on rural, unregulated streams with annual peak streamflow data from 1791 through the 2021 water years in and near Virginia and West Virginia. The study was done in cooperation with the Federal Emergency Management Agency, the West Virginia Department of Transportation, and the Virginia Department of Transportation.

Regression equations were developed for estimating flood frequency and magnitude. Twelve regions with homogeneous flood characteristics were identified. Generalized least squares regression equations relating logarithmic-transformed drainage area and peak streamflow were developed for the 0.5, 0.2, 0.1, 0.04, 0.02, 0.01, 0.005, and 0.002 annual exceedance probabilities (AEPs). Drainage area was the only significant variable for all equations. The range of drainage areas used to develop the equations differed for each region; the smallest drainage area in any region was 0.21 square miles (mi²) and the largest drainage area in any region is 2,966 mi². Pseudo coefficient of determination (pseudo-R²) values for regression equations ranged from 0.481 to 0.995 for all regions and AEPs. Performance metrics and diagnostic plots indicated that equations for 11 of the 12 regions showed generally good performance, with pseudo-R² values ranging from 0.762 to 0.968 for the 0.01 AEP.

The overall average change in at-site 0.01 AEP annual peak streamflows at individual streamgages was 0.5 percent compared to the most recent 2011 Virginia study and 2.3 percent compared to the most recent 2010 West Virginia study. Changes from the previous studies for estimates from regional equations for the 0.01 AEP, solved specifically for a 50 mi² basin, ranged from a 30 percent increase to a 45 percent decrease in areas where the previous regions overlapped with the current regions by 750 mi² or more.

New regional skews were developed using Bayesian weighted least-squares/Bayesian generalized least-squares regression for two skew regions that included the study area. A constant regional skew of 0.50 was computed for streams in Virginia, West Virginia, and Maryland that drain to the Atlantic Ocean. A constant regional skew of 0.048 was computed for streams that drain to the Gulf of America, including streams in Kentucky and Tennessee, most of West Virginia, far southwestern Virginia, and part of western Maryland.

About 12 percent of the 418 streamgages with 30 or more gaged peaks had statistically significant (p-value [significance level] less than or equal to 0.05) trends, with 40 of these exhibiting positive trends and 11 exhibiting negative trends. Streamgages with 30 percent or greater development were excluded from regression analyses.

A regulation index was developed that accounted for storage and drainage area of dams and drainage area at the streamgage; a value of 0.0040 or more for the regulation index indicates regulated peak streamflow. Frequency analyses were done at 86 streamgages on regulated streams.

Regression procedures developed in this study are applicable only to rural, unregulated streams within Virginia and West Virginia with drainage basins that (1) are within the range of drainage areas used to develop the equations for each region, (2) included less than 30 percent of developed area, and (3) had a regulation index less than 0.0040.