The U.S. Geological Survey (USGS) in cooperation with the Fountain Creek Watershed, Flood Control and Greenway District assessed remediation scenarios to attenuate peak flows and reduce sediment loads in the Fountain Creek watershed. To evaluate these strategies, the U.S. Army Corps of Engineers Hydrologic Engineering Center (HEC) hydrologic and hydraulic models were employed.

The U.S. Army Corps of Engineers modeling system HEC-HMS (Hydrologic Modeling System) version 3.5 was used to simulate runoff in the Fountain Creek watershed, Colorado, associated with storms of varying magnitude and duration. Rain-gage precipitation data and radar-based precipitation data from the April 28–30, 1999, and September 14–15, 2011, storm events were used in the calibration process for the HEC-HMS model. The curve number and lag time for each subwatershed and Manning's roughness coefficients for each channel reach were adjusted within an acceptable range so that the simulated and measured streamflow hydrographs for each of the 12 USGS streamgages approximated each other.

The U.S. Army Corps of Engineers modeling system HEC-RAS (River Analysis System) versions 4.1 and 4.2 were used to simulate streamflow and sediment transport, respectively, for the Fountain Creek watershed generated by a particular storm event. Data from 15 USGS streamgages were used for model calibration and 7 of those USGS streamgages were used for model validation. The calibration process consisted of comparing the simulated water-surface elevations and the cross-section-averaged velocities from the model with those surveyed in the field at the cross section at the corresponding 15 and 7 streamgages, respectively. The final Manning’s roughness coefficients were adjusted between –30 and 30 percent at the 15 calibration streamgages from the original left, right, and channel-averaged Manning's roughness coefficients upon completion of calibration.

The U.S. Army Corps of Engineers modeling system HEC-RAS version 4.2 was used to simulate streamflow and sediment transport for the Fountain Creek watershed generated by a design-storm event. The Laursen-Copeland sediment-transport function was used in conjunction with the Exner 5 sorting method and the Ruby fall-velocity method to predict sediment transport. Six USGS streamgages equipped with suspended-sediment samplers were used to develop sediment-flow rating curves for the sediment-transport-model calibration. The critical Shields number in the Laursen-Copeland sediment-transport function and the volume of sediment available at a given cross section were adjusted during the HEC-RAS sediment-model calibration process.

HEC-RAS model simulations used to evaluate the 14 remediation scenarios were based on unsteady-state streamflows associated with a 24-hour, 1-percent annual exceedance probability (100-year) National Oceanic and Atmospheric Administration Type II precipitation event. Scenario 0 represents the baseline or current conditions in the watershed and was used to compare the remaining 13 scenarios. Scenarios 1–8 and 12 rely on side-detention facilities to reduce peak flows and sediment transport. Scenario 9 has a diversion channel, and scenario 10 has a reservoir. Scenarios 11 and 13 incorporate channel armoring and channel widening, respectively. Scenarios 8 and 10, the scenario with the most side-detention facilities, and the scenario with the reservoir, respectively, were the most effective at reducing sediment transport and peak flow at the Pueblo, Colorado, streamgage. Scenarios 8 and 10 altered the peak flow by –58.9 and –56.4 percent, respectively. In turn, scenarios 8 and 10 altered the sediment transport by –17.7 and –62.1 percent, respectively.