Simulation of nutrient and sediment concentrations and loads in the Delaware inland bays watershed: Extension of the hydrologic and water-quality model to ungaged segments

Rapid population increases, agriculture, and industrial practices have been identified as important sources of excessive nutrients and sediments in the Delaware Inland Bays watershed. The amount and effect of excessive nutrients and sediments in the Inland Bays watershed have been well documented by the Delaware Geological Survey, the Delaware Department of Natural Resources and Environmental Control, the U.S. Environmental Protection Agency’s National Estuary Program, the Delaware Center for Inland Bays, the University of Delaware, and other agencies. This documentation and data previously were used to develop a hydrologic and water-quality model of the Delaware Inland Bays watershed to simulate nutrients and sediment concentrations and loads, and to calibrate the model by comparing concentrations and streamflow data at six stations in the watershed over a limited period of time (October 1998 through April 2000). Although the model predictions of nutrient and sediment concentrations for the calibrated segments were fairly accurate, the predictions for the 28 ungaged segments located near tidal areas, where stream data were not available, were above the range of values measured in the area.

The cooperative study established in 2000 by the Delaware Department of Natural Resources and Environmental Control, the Delaware Geological Survey, and the U.S. Geological Survey was extended to evaluate the model predictions in ungaged segments and to ensure that the model, developed as a planning and management tool, could accurately predict nutrient and sediment concentrations within the measured range of values in the area. The evaluation of the predictions was limited to the period of calibration (1999) of the 2003 model.

To develop estimates on ungaged watersheds, parameter values from calibrated segments are transferred to the ungaged segments; however, accurate predictions are unlikely where parameter transference is subject to error. The unexpected nutrient and sediment concentrations simulated with the 2003 model were likely the result of inappropriate criteria for the transference of parameter values. From a model-simulation perspective, it is a common practice to transfer parameter values based on the similarity of soils or the similarity of land-use proportions between segments. For the Inland Bays model, the similarity of soils between segments was used as the basis to transfer parameter values. An alternative approach, which is documented in this report, is based on the similarity of the spatial distribution of the land use between segments and the similarity of land-use proportions, as these can be important factors for the transference of parameter values in lumped models. Previous work determined that the difference in the variation of runoff due to various spatial distributions of land use within a watershed can cause substantialloss of accuracy in the model predictions.

The incorporation of the spatial distribution of land use to transfer parameter values from calibrated to uncalibrated segments provided more consistent and rational predictions of flow, especially during the summer, and consequently, predictions of lower nutrient concentrations during the same period. For the segments where the similarity of spatial distribution of land use was not clearly established with a calibrated segment, the similarity of the location of the most impervious areas was also used as a criterion for the transference of parameter values.

The model predictions from the 28 ungaged segments were verified through comparison with measured in-stream concentrations from local and nearby streams provided by the Delaware Department of Natural Resources and Environmental Control. Model results indicated that the predicted edge-of-stream total suspended solids loads in the Inland Bays watershed were low in comparison to loads reported for the Eastern Shore of Maryland from the Chesapeake Bay watershed model. The flatness of the terrain and the low annual surface runoff are important factors in determining the amount of detached sediment from the land that is delivered to streams. The highest predicted total suspended solids loads were found in the southern part of the watershed, where the values are associated with high total streamflow and a high surface-runoff component, and related to soil and aquifer permeability and land use. Nutrient loads from model segments in the southern part of the Inland Bays watershed were also higher than those measured in the northern part of the basin, due to relatively high runoff and the substantial amount of available organic fertilizer (animal waste) that results in over-application of organic fertilizer to crops.

Time series of simulated hourly concentrations indicated a seasonal pattern in the simulated base flow for total nitrogen, with the lowest values occurring during the summer and the highest values during the winter months. Total phosphorus and total-suspended-solids concentrations were less seasonal and were more storm-dependent; in general, base-flow concentrations of total phosphorus and total suspended solids were low. During storm events, the total nitrogen concentrations tended to be diluted and total phosphorus concentrations tended to rise sharply. Nitrogen was transported mainly in the aqueous phase and largely through ground water, whereas phosphorus was strongly associated with sediment, which washes off during rainfall events.