Abstract

Flooding induced by extreme rainfall events causes tremendous loss of life and property and infrastructure failure. Accurate representation of precipitation, which has high variation in space and time, is critical to hydrologic model simulations and flood analyses. In this study, we examined responses of differently sized United States Geological Survey (USGS) hydrologic units to heavy precipitation using three different data sets. The first consists of rainfall observed at individual meteorological gauges. The second uses the National Centers for Environmental Prediction (NCEP) Environmental Modeling Center (EMC) 4 km gridded radar-estimated precipitation (GRIB) Stage IV data. The third one derives from the method we developed that blends gauge data with the spatial coverage of the Parameter-elevation Relationships on Independent Slopes Model (PRISM) data. We examined how two watersheds in South Carolina respond to the three different representations of heavy rainfall, using the Hydrologic Engineering Center’s Hydrologic Modeling System (HEC-HMS) developed by the U.S. Army Corps of Engineers. We found that the latter two precipitation inputs that consider spatial representation of rainfall yielded similar performance and improved simulated streamflow as compared to simulation using rainfall observed at individual meteorological gauges. The method we developed overcomes the spatial sparsity of rain gauges required for interpolation and extends availability of precipitation surfaces. Our study advances the understanding of advantages and limitations of different precipitation products for flood simulation.

1. Introduction

Flooding is one of the most frequently occurring and costly natural hazards in the world, causing death and injury, displacement of communities, and loss of private and public properties [1]. Hydrological models improve understanding of hydrological processes and flood forecasting in ways that could mitigate impacts [26]. However, the performance of hydrological models depends on the spatial and temporal accuracy of rainfall data [7, 8].

The most widely used source of rainfall data is observations at rain gauges, which provide direct measurement of precipitation intensity and duration at individual points. Rain gauge observations are considered as the most accurate precipitation data source and provide the basis for evaluating and calibrating other precipitation products such as those based on radar data [9, 10]. Data availability at some rain gauges spans several decades, which makes them a valuable data source for the study of historic extreme precipitation and associated flooding [1114]. Despite the advantages of reliability and longtime coverage, rain gauge observations have poor representation of areal precipitation, especially where the gauge network is sparse [15, 16]. The network of rain gauges with high temporal resolution is particularly limited. In the United States, there are over 25,000 daily recording stations and only 8,000 hourly stations [17]. Rain gauge observations inadequately capture the spatial and temporal variability of short-duration storm events, especially in small catchments, thus limiting the accuracy of streamflow simulation [18].

In contrast to the shortcomings of rain gauge observations in representing spatial and temporal structures of precipitation, radar precipitation data provide a better capture of precipitation over spatial and temporal scales that are paramount to hydrologic applications [19]. Radars do not directly measure rainfall amounts but estimate the rainfall rate via an empirical relationship with the radar reflectivity factor. Great efforts have been exerted to improve estimations of the reflectivity factor, calibrate radar data, and increase their usability [2023]. The evolvement of radar data in the United States has experienced four generations [2427]. The most recent generation Stage IV radar data take advantage of manual quality control performed on the Stage III data and have mosaiced Stage III data covering the entire continental United States, providing hourly precipitation at the scale of 4 km polar stereographic grids starting from January 1, 2002 [27]. Nevertheless, the temporal coverage of radar data (i.e., 2002 to present) limits the usage of radar for the hydrological simulation prior to 2002.

The potential and effectiveness of using radar data to improve hydrological simulations has captured the interest of the research and operational communities for some time. Comparisons of hydrological simulations using radar and gauge data as precipitation input have been conducted. Lopez et al. suggested that the Stage IV radar product had better hydrologic skills over other rainfall data sources based on their investigation in the Treja catchment in Italy using a rainfall-runoff model named “geomorphological instantaneous unit hydrograph (GIUH)” [28]. Similar findings were concluded from hydrologic simulation conducted in four river basins in the United States using the distributed Clark model [29]. One particular interest is tests on rainfall gauge data, radar data, and remote-sensed rainfall products (e.g., Tropical Rainfall Measuring Mission (TRMM)) using the Hydrologic Engineering Center’s Hydrologic Modeling System (HEC-HMS) developed by the U.S. Army Corps of Engineers [30], which has been applied to a wide range of geographic areas for solving a variety of hydrological problems. Hamedi and Fuentes examined the effectiveness of Gridded Surface Subsurface Hydrologic Analysis (GSSHA) and HEC-HMS to predict outlet hydrographs in the Park City catchment basin in Utah using next-generation radar (NEXRAD) and gauge data as input and suggested the reliability of using NEXRAD [19]. Investigation over the Brays Bayou watershed in Houston proved simulated outflows from radar data to be accurate, and in some cases more accurate, than those from rain gauges [31]. However, the study in two subbasins of the Cumberland River basin in Middle Tennessee found that radar-driven simulations offered no improvement in predicting the streamflow volume, magnitude, and time to peak over gauge-driven simulations [15].

The comparison of hydrologic simulations using radar and gauge precipitation data involves the spatial and temporal resolution and coverage of the rainfall data, model mechanisms (e.g., simplicity or complexity and parameterization), and watershed characteristics (e.g., size and climatologic and physiographic settings) [15, 32]. South Carolina was impacted by the heavy rainfall and flooding in October 2015 and Hurricane Matthew in October 2016. The former led to $138 million paid losses and the failure of 52 dams [33, 34]. The latter costs $116.9 million federal disaster assistance for the survivors in South Carolina [35]. We assess the suitability of different precipitation data sources that are constrained by spatial and temporal coverage and resolutions of these data in the flood simulation using HEC-HMS for South Carolina. We also propose a new method that integrates merits of precipitation gauge data and the widely used gridded daily Parameter-elevation Relationships on Independent Slopes Model (PRISM) [36] data to overcome the sparsity of rain gauge data and test the effectiveness of the new method.

2. Methods

Flood simulation was conducted in the Gills Creek (hydrologic unit code (HUC): 0305011002; 193 km2) watershed in central South Carolina and Waccamaw (HUC: 03040206; 4277 km2) watershed at the border of South Carolina and North Carolina (Figures 1 and 2) using HEC-HMS [30]. One time period was selected to calibrate the model, and two other time periods were chosen to test the effects of different precipitation sources on model performance. Observed streamflow at United States Geological Survey (USGS) gauges within the investigated watersheds was ranked to obtain the percentiles of streamflow. Time periods with high streamflow that exceeded the 99 percentile were selected for calibration and testing, subject to the availability of streamflow and precipitation data (Tables 1 and 2). HEC-HMS was calibrated using a variety of data sources and methods. Topography was derived from the USGS National Elevation Dataset (NED), which provides gridded elevation values at 10 m spatial resolution. Soil data were obtained through the State Soil Geographic Database (STATSGO) maintained by the National Resource Conservation Service. Land use and land cover information was obtained from the National Land Cover Dataset 2011 created by USGS. The gridded curve number (CN) technique suggested by the Soil Conservation Service (SCS), which enables spatially distributed infiltration calculations, was used to simulate infiltration loss. The CN technique quantifies the infiltration capacity area based on land use, soil and land cover type, and hydrologic soil group [37]. Translation of excess precipitation to runoff was conducted using the Clark unit hydrograph transformation [38]. Baseflow was modeled using an exponential decrease function.

We investigated the effects of three precipitation sources on flood simulation including (1) station observations (hereafter ST), (2) radar data (hereafter RD), and (3) station observations adjusted by the PRISM data which provide daily total precipitation in the continental United States, from 1981 to present (hereafter ST-PRISM) [36]. Station observations in the Automated Surface Observing System (ASOS) obtained from the Iowa Environmental Mesonet [39] are point-based precipitation values, while radar data [27] are areal representation of precipitation. Radar data were not developed fully until 2002, which constrains the data availability of areal representation of precipitation for flood simulation prior to 2002. Although the PRISM data provide areal representation of precipitation for a longer time period than radar data (i.e., 1981 to present), the daily temporal resolution of the PRISM data is not adequate for HEC-HMS which requires hourly precipitation inputs. To overcome this problem, we exploit the high temporal resolution of station observations and adjust these point-based precipitation values using gridded PRISM data with the following procedures: first, we selected candidate stations that are within or surrounding the watershed where flood simulation was conducted. Second, for each day, we identified the station whose daily total precipitation was closest to the areal daily total of PRISM on that day. Third, the hourly precipitation of the identified station was multiplied by the ratio between the areal daily total of PRISM and daily total of the station. In this way, the daily total precipitation of the station matched the areal daily total of PRISM. The approach takes the advantage of areal representation of PRISM data and high temporal resolution of station data (i.e., hourly observations).

Flood simulation on Gills Creek was compared against observed discharge at the USGS station of Gills Creek at Columbia SC (site ID: 02169570) to validate the simulated discharge of Upper Gills Creek subbasin (Figure 1). Precipitation data included station (ST) data from the observation at Columbia Owens Downtown Airport (CUB), radar (RD) data from the National Centers for Environmental Prediction (NCEP) Environmental Modeling Center (EMC) 4 km gridded radar-estimated precipitation (GRIB) Stage IV data [27], and hybrid (ST-PRISM) data from CUB and PRISM.

Flood simulation for the Waccamaw watershed was conducted in two individual subbasins. Streamflows observed at the USGS station of Waccamaw River at Freeland (site ID: 02109500) and Buck Creek near Longs, SC (site ID: 02110400), were used to validate the simulated streamflow of subbasins 1 and 2, respectively (Figure 2). Observations at Buck Creek near Longs, SC (site ID: 02110400), were not available for the second testing period (i.e., Sep 14–25, 1999). The precipitation data source varied for different subbasins and periods. In the second testing period, radar precipitation data were not available, and North Myrtle Beach (CRE) had the only available precipitation gauge. Precipitation data sources for subbasins 1 and 2 were different in the calibration and the first testing periods (Figure 2 and Table 3). The selection of precipitation stations was based on the closeness of stations to subbasins (Figure 2).

HEC-HMS was calibrated separately using point-based representation (i.e., ST) and two areal representations of precipitation (i.e., RD and ST-PRISM), which yielded a set of parameters for each of the three precipitation inputs. In the testing periods, flood simulation was conducted using the three calibrated models with the same precipitation inputs used to calibrate the models. Model performance was assessed using four statistical indices that compare the simulated streamflow against the observed streamflow: Nash–Sutcliffe Efficiency (NSE), R square (R2), p-bias, and root mean square error (RMSE) (see [40, 41] for descriptions). We also compared the total precipitation amount received in a watershed or subbasin and temporal patterns of precipitation during calibration and testing periods from three different data sources to investigate the effects of different precipitation inputs on flood simulation.

3. Results

The total precipitation amount derived from the two areal representations (i.e., RD and ST-PRISM) was close to each other in both watersheds and for different periods of calibration and testing (Figures 3 and 4). The largest difference was 17 mm over 9 days from October 2 to 10, 2015, in subbasin 2 of Waccamaw (Figure 4(d)). In contrast, the difference between the area-based and the point-based representations of precipitation was much larger. The largest difference was 67 mm (ST-PRISM: 153 mm versus ST: 86 mm) in subbasin 1 of Waccamaw from February 3 to 15, 2016 (Figure 4(a)). The peaking time of the point-based representation and the two areal representations of precipitation was similar in most of the situations, except for a few unmatched peaks in the two subbasins of Waccamaw in October 2015 (Figures 4(c) and 4(d)).

Good calibrations were achieved using the three types of precipitation inputs (Figures 3(a), 4(a), and 4(b), and Tables 4 and 5). As for the performance in testing periods, the three precipitation inputs had good performance in subbasin 1 of Waccamaw in October 2015 (Figure 4(c) and Table 5). Overestimation was observed in Gills Creek in December 2014 and in subbasin 2 of Waccamaw in October 2015 for all the three precipitation inputs (Figures 3(b) and 4(d), and Tables 4 and 5). The ST-calibrated model suffered from severe overestimation, while the two areal representations of precipitation had some underestimation in Gills Creek in September 2011 (Figure 3(c) and Table 4). The ST-calibrated model had less underestimation than the model calibrated by the areal representation of precipitation (i.e., ST-PRISM) in subbasin 1 of Waccamaw in September 1999 (Figure 4(e) and Table 3).

4. Discussion and Conclusion

Precipitation amounts affected two aspects of model performance: parameterization during model calibration and the total amount of water that entered a watershed for flood simulation. In Gills Creek, because of the similar amount and temporal patterns of the area-based and the point-based representations of precipitation in the calibration period, calibration using the three precipitation inputs yielded similar parameters and model performance (Figure 3(a)). However, during testing periods, precipitation estimates differed; thus, total precipitation became the most influential factor on model performance. The overestimation by the ST-calibrated model occurred because during testing periods, ST-estimated precipitation was much higher than RD-estimated or ST-PRISM-estimated precipitation (Figures 3(b) and 3(c)). Overall, models calibrated by the two areal representations of precipitation achieved better performance than the ST-calibrated model in Gills Creek.

In Waccamaw, because of much lower station precipitation (ST) compared to the two areal estimates in the calibration period (Figures 4(a) and 4(b)), the potential retention scale factors in subbasins 1 and 2 were set to very low values. This reduced the loss of precipitation and allowed a similar calibration to the RD-calibrated and ST-PRISM-calibrated models for comparisons in testing periods. In subbasin 2, the low value of the potential retention scale factor and the two much larger peaks of ST caused overestimation for the ST-calibrated model (Figure 4(d) and Table 5). In subbasin 1, the potential retention scale factor was particularly set to an extremely low value, although this setting yielded good simulations in the two testing periods for the ST-calibrated model (Figures 4(c) and 4(e)). Others have noted model parameterization differences between calibrations with the radar and gauge precipitation [42]. They underscored the importance of linking calibration parameters to realistic processes. The potential retention scale factor in the ST-calibrated model was set to very low values, which indicated the watershed had little ability to retain precipitation. However, land use and land cover data from the National Land Cover Dataset suggested 57.8% of the watershed is covered by either forest or wetlands which can retain considerable precipitation amounts. Therefore, although the ST-calibrated model achieved good performance in the two testing periods in subbasin 1, parameterization of the ST-calibrated model did not reflect realistic hydrological processes. Overall, models calibrated by the two areal representations of precipitation (RD and ST-PRISM) had similar performance—better than the ST-calibrated model (Figure 4 and Table 5).

Our study underscores the importance of spatial representation of precipitation for flood simulation [7, 18, 43, 44] as observations at a single station (ST) did not reliably represent the precipitation surface. Yet, despite many new spatial precipitation products (e.g., radar and satellite estimates), resulting flood simulations show inconsistent performance. Our findings are similar to those of the studies [29, 45, 46], which suggest precipitation input using radar data produced better hydrological simulation than station observations, but different from those of the study [15] which found no improvement in hydrological simulation using radar data. Several factors could contribute to these disparities, including local errors in radar and rain gauge data, modeling simplicity or complexity, varying watershed sizes, choice of runoff generation mechanisms, or accuracy of watershed parameters and calibration [32].

In terms of the quality of station observations, radar data, and the PRISM data, rain gauge observations are considered as the most accurate precipitation data source. They are often used to evaluate the accuracy of radar data [9, 10]. The PRISM data are the official spatial climate data sets of the U.S. Department of Agriculture and have been produced using methods that consider location, elevation, coastal proximity, topographic facet orientation, vertical atmospheric layer, topographic position, and orographic effectiveness of the terrain and incorporate radar observations. Readers are referred to the study [36] for the uncertainty assessment of PRISM data.

Also unsolved is the optimal method for incorporating gauge data to create a spatial representation of the precipitation surface for flood simulation. Spatial interpolation is a typical way of converting discrete gauge observations into areal representation of precipitation, but it requires adequate density of rain gauges. In addition, interpolation methods and associated parameters must be appropriately chosen [47], since errors and uncertainties associated with spatial interpolation are to be propagated to hydrological models.

When gauge stations are too sparse to conduct spatial interpolation, other means could be used to convert point-based estimates to area-based precipitation estimates. Our blended ST-PRISM method integrates merits of precipitation gauge data and spatial precipitation represented by PRISM. Although it should not be viewed as a replacement for radar data or spatial interpolation applied to rain gauges, it provides an alternative way to capture spatial precipitation, and, in our example, the resulting surface led to model performance as good as that found in the radar-calibrated model. Compared to radar data, ST-PRISM extends the data availably prior to 2002. It is also useful when the density of rain gauges is too low to perform spatial interpolation and eliminates the need to choose and parameterize spatial interpolation methods.

We conducted flood modeling in two watersheds in South Carolina that have experienced severe flooding recently and tested the effects of different precipitation inputs on model performance. Each precipitation source has its own advantages and disadvantages with respect to spatial and temporal coverage and resolution. Our results suggest that there are ways to combine sources to more accurately simulate floods with a simple hydrological model. Specifically, we demonstrate how integrating precipitation gauge data with the widely used PRISM data set extends the temporal availability of an areally representative precipitation data set for flood simulation.

Data Availability

All data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

A part of the manuscript was presented in a session at the American Association of Geographers Annual Conference 2018.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Advanced Support for Innovative Research Excellence (ASPIRE) I, Track 2B, from the Office of the Vice President for Research, University of South Carolina, and the National Oceanic and Atmospheric Administration (NOAA) Climate Program Office (Grant no. NA16OAR4310163) to the Carolinas Integrated Sciences and Assessments (CISA).