CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

MONTHLY RUNOFF DISTRIBUTION RELATED TO BASIN GEOGRAPHY1

ABSTRACT: Mean monthly runoff from ungaged drainage basins that have significant snowpacks each year can be estimated quite well by assuming that the time duration between snowfall and snowmelt is the predominant factor in temporal runoff distribution. That time span is related to basin temperatures which are, in turn, functions of basin elevation and latitude. Regional hydrologic analyses of gaged basin data create regression equations for estimating runoff distribution by month. These equations then can be applied to ungaged basins. Basin latitude and mean elevation are two independent variables that can be used in estimating monthly runoff distributions.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

HYDROLOGIC EFFECTS OF THE FLOOD ABATEMENT PROGRAM IN SOUTHWESTERN OKLAHOMA1

ABSTRACT: Winter Creek is a tributary of the Washita River in south-western Oklahoma. The Soil Conservation Service installed floodwater retarding structures which controlled runoff from 56 percent of a 33-square-mile (8550-hectare) gaged drainage area. Application of a hydrologic model to the flood peaks indicated that the structural treatment reduced the flood peaks an average of 61 percent. The Winter Creek channel has narrowed and deepened since the structural treatment was applied. The severe bank erosion occurring before treatment has been arrested and sediment yield from the watershed has been reduced 50 to 60 percent. In some reaches of the channel there has been a dense growth of trees in recent years.     

Simulation of Potential Oyster Density with Variable Freshwater Inflow (1965–2000) to the Caloosahatchee River Estuary, Southwest Florida, USA

Oyster beds are disappearing worldwide through a combination of over-harvesting, diseases, and salinity alterations in the coastal zone. Sensitivity of oysters to variable discharge and salinity is particularly acute in small sub-tropical estuaries subject to regulated freshwater releases. South Florida has sub-tropical estuaries where watershed flood control sometimes results in excessive freshwater inflow to estuaries during the wet season (May–Oct) and reduced discharge and increased salinities in the dry season (Nov–Apr). The potential to reserve freshwater accumulated during the wet season could offer the capacity to regulate freshwater at different temporal scales, thus optimizing salinity conditions for estuarine biota. The goal of this study was to use simulation modeling to explore the effects of freshwater inflows and salinity on adult oyster survival in the Caloosahatchee River Estuary (CRE) in southwest Florida. Water managers derived three different freshwater inflow scenarios for the CRE based on historical and modified watershed attributes for the time period of 1965–2000. Three different salinity time series were generated from the inflow scenarios at each of three sites in the lower CRE and used to conduct nine different oyster simulations. Overall, the predicted densities of adult oysters in the upstream site were 3–4 times greater in seasons that experienced reduced freshwater inflow (e.g., increased salinity) with oyster density in the lower estuary much less influenced by the inflows. Potential storage of freshwater reduced the frequency of extreme flows in the wet season and helped to maintain minimum inflow in the dry season near the estuarine mouth. Analyses of inflows indicated that discharges ranging from 0 to 1,500 cfs could promote favorable salinities of 10–25 in the lower CRE depending on wet versus dry season climatic conditions. This range of inflows is similar to that derived in other studies of the CRE and emphasizes the value of simulation models to help prescribe freshwater releases which benefit estuarine biota.     

LEGAL CHALLENGE OF THE UNGAGED WATERSHED1

The lack of uniform techniques for estimating design discharges in ungaged areas is a source of growing concern in the courts now faced with challenges to floodplain boundaries and culvert design for highway crossings. This paper summarizes a court case in which calculations of the design discharge and the hydraulic backwater effects of a major highway culvert were contested by the plaintiff. Emphasis in the paper is placed on the variation in computed flows and the interpretation of the court in the face of diverse hydrologic methods for the ungaged watershed. The results of a preliminary evaluation of ungaged watershed methods applicable to Pennsylvania are also reported in terms of standard error and bias.     

MODELING NORTH CASCADES ANNUAL RUNOFF DISTRIBUTION USING DRAINAGE BASIN ELEVATION DISTRIBUTION1

ABSTRACT: The annual distribution of flow in a drainage basin within a given region is a function of many factors. These may include annual distribution of rainfall, basin orientation, ground cover, or presence of glaciers. Since the North Cascades region of northern Washington State has little variation in precipitation distribution by month, and the region has significant snowpack, one would predict that in an unregulated basin, basin elevation would be one of the most important factors impacting an annual hydrograph distribution. Such a prediction can be made since the higher a drainage basin is, the larger the portion of runoff that would occur as late spring snowmelt. Given that there is a relationship between elevation distribution and annual hydrograph, the problem becomes one of how to use this relationship to model an ungaged basin's hydrograph. This study concludes that, within the North Cascades region and perhaps within other regions, an effective method of determining annual flow distribution is to model ungaged flows in the same manner as flows from a gaged basin with an elevation distribution similar to that of the subject basin.     

A SIMPLE METHOD FOR ESTIMATING BASEFLOW AT UNGAGED LOCATIONS1

ABSTRACT: Baseflow, or water that enters a stream from slowly varying sources such as ground water, can be critical to humans and ecosystems. We evaluate a simple method for estimating base‐flow parameters at ungaged sites. The method uses one or more baseflow discharge measurements at the ungaged site and longterm streamflow data from a nearby gaged site. A given baseflow parameter, such as the median, is estimated as the product of the corresponding gage site parameter and the geometric mean of the ratios of the measured baseflow discharges and the concurrent discharges at the gage site. If baseflows at gaged and ungaged sites have a bivariate lognormal distribution with high correlation and nearly equal log variances, the estimated baseflow parameters are very accurate. We tested the proposed method using long‐term streamflow data from two watershed pairs in the Driftless Area of southwestern Wisconsin. For one watershed pair, the theoretical assumptions are well met; for the other the log‐variances are substantially different. In the first case, the method performs well for estimating both annual and long‐term baseflow parameters. In the second, the method performs remarkably well for estimating annual mean and annual median baseflow discharge, but less well for estimating the annual lower decile and the long‐term mean, median, and lower decile. In general, the use of four measurements in a year is not substantially better than the use of two.     

ADAPTATION OF FLOOD PEAKS AND DESIGN HYDROGRAPHS FROM GAGED TO NEARBY UNGAGED WATERSHEDS1

ABSTRACT: Equations were developed to transform peak flows and to adapt design hydrographs and unit hydrographs from gaged watersheds to ungaged watersheds with similar hydrologic characteristics. Dimensional analysis was used to develop adjustment equations for peak flow and time base, and these equations were reinforced with results from regional flood frequency research. The authors believe that the use of these transformation equations should yield more reliable flood peak values and hydrogrphs than the common use of empirical flood estimating curves or equations.     

CHANGES IN HYDROLOGIC RESPONSE DUE TO STREAM NETWORK EXTENSION VIA LAND DRAINAGE ACTIVITIES1

ABSTRACT: The objective of this work is to determine the effects of extension of a stream network through land drainage activities during the late 1800s on the hydrologic response of a watershed. The Mackinaw River Basin in Central Illinois was chosen as the focus and the pre‐land and post‐land drainage activity hydrologic responses were obtained through convolution of the hill slope and channel responses and compared. The hill slope response was computed using the kinematic wave model and the channel response was determined using the geomorphologic instantaneous unit hydrograph method. Our hypothesis was that the hydrologic response of the basin would exhibit the characteristic effects of settlement (i.e., increases in peak discharges and decreases in times to peak). This, indeed, is what occurred; however, the increase in peak discharges diminishes as scale increases, leaving only the decrease in times to peak. At larger scales, the dispersive effects of the longer hill slope lengths in the pre‐settlement scenario seem to balance the depressive effects of the longer path lengths in the post‐settlement scenario, thus the pre‐settlement and post‐settlement peak discharges are approximately equivalent. At small scales, the dispersion caused by the hill slope is larger in the pre‐settlement case; thus, the post‐settlement peak discharges are greater than the pre‐settlement.     

AN APPLICATION RUNOFF MODEL STRATEGY FOR UNGAGED WATERSHEDS1

ABSTRACT: Mathematical models for predicting watershed surface flow responses are available, most of which are elaborate nonlinear numerical surface and channel flow models linked with infiltration models. Such models may be used to make predictions for ungaged areas, assuming an acceptable fitting of the model to the topography and roughness of the real system. For some application purposes, these models are impractical because of their complexity and expensive computer solutions. A procedure is developed that uses a complex model of an ungaged area to derive a simpler parametric nonlinear system model for repetitious simulation with input sequences. The predicted flow outputs are obtained with the simpler model at significant savings of money and time. The procedures for constructing a complex kinematic model of a 40 acre (161,880 m²) reference watershed and deriving the simpler system model are outlined. The results of predictions from both models are compared with a selected set of measured events, all having essentially the same initial conditions. Peak discharges ranged from 3 to 118 ft³/sec (0.085 to 3.34 m³/sec), which includes the largest event of record. The inherent limitations of lumped systems models are demonstrated, including the bias caused by their inability to model infiltration losses after rainfall ceases. Computer costs and times for the models were compared. The derived simple model has a cost advantage when repeated use of a model is required. Such an applications hydrologic model has an engineering tradeoff of reduced accuracy, and lumping bias, but is more economical for certain design purposes.     

THE MISUSE OF HYDROLOGIC UNIT MAPS FOR EXTRAPOLATION,REPORTING, AND ECOSYSTEM MANAGEMENT1

ABSTRACT: The use of watersheds to conduct research on land/water relationships has expanded recently to include both extrapolation and reporting of water resource information and ecosystem management. More often than not, hydrologic units (HUs) are used for these purposes, with the implication that hydrologic units are synonymous with watersheds. Whereas true topographic watersheds are areas within which apparent surface water drains to a particular point, generally only 45 percent of all hydrologic units, regardless of their hierarchical level, meet this definition. Because the area contributing to the downstream point in many hydrologic units extends far beyond the unit boundaries, use of the hydrologic unit framework to show regional and national patterns of water quality and other environmental resources can result in incorrect and misleading illustrations. In this paper, the implications of this misuse are demonstrated using four adjacent HUs in central Texas. A more effective way of showing regional patterns in environmental resources is by using data from true watersheds representative of different ecological regions containing particular mosaics of geographical characteristics affecting differences in ecosystems and water quality.     

AN OBJECTWE TEST FOR HYDROLOGIC SCALE1

ABSTRACT: Improving the reliability of parametric hydrologic models (sometimes called cenceptual rainfall-runoff models) in the continuous simulation of runoff from ungaged catchments has been frustrated by difficulties in estimating model parameters from catchment characteristics. An underlying problem is that these models use parameters to represent catchments as a whole, whereas data on catchment characteristics are collected at multiple field locations and are difficult to transform into one measure of collective impact. Subdividing the catchment and calibrating a stochastic parametric model to estimate distributions for the parameters that covered the range of observed streamflow values was found to improve the simulations. This paper presents an optimization of the amount of subdivision to use in simulation with a version of the Stanford Watershed Model using available climatological data. The calibration process assumes that catchment heterogeneity introduces errors that can be reduced by calibrating parameters as spatial distributions rather than single values. Calibrations for three diverse small gaged catchments located in California and in Virginia found the optimal number of subdivisions to range from 4 to 25 and the optimal scale to range from 0.3 to 2.1 mi².     

SEARCH: A GAGED WATERSHED INVENTORY1

ABSTRACT: A method, adapted to an APL interactive terminal, is described which allows the operator to thoroughly search a large set of gaged watersheds in order to find sources of comparable hydrologic data for detailed analysis. Bases of the search - and inventory - include drainage basin size and elevation, and geographical and temporal parameters, and should enhance opportunities for more reliable use of existing data.     

MULTIVARIATE ANALYSIS OF WATER QUALITY IN THE RICHIBUCTO DRAINAGE BASIN (NEW BRUNSWICK,CANADA)1

ABSTRACT: Specific conductivity, pH, dissolved oxygen, carbon, phosphorous, and nitrogen species were measured at 36 stations in the Richibucto River drainage basin, including the estuary, in New Brunswick, Canada, over the six‐year period 1996 through 2001. Each station was sampled between 1 and 26 times (mean = 7.5, standard deviation = 6.0) during the ice free seasons without regard to tide. There was significant variance among stations in most parameters. Principal component analysis (PCA) was used to identify the processes explaining the observed variance in water quality. Because of the high variability in specific conductance, stations were first grouped in a freshwater subset and an estuarine (brackish water) subset. For freshwater stations, most of the variance in water quality was explained by pH and total organic carbon, as well as high nutrient concentrations. These high nutrient concentrations, along with water salinity, which varies with flow and tides, are also important in determining water quality variability in brackish water. It is recommended that water quality parameters that were found to explain most of the variance by PCA be monitored more closely, as they are key elements in understanding the variability in water quality in the Richibucto drainage basin. Cluster analyses showed that high phosphorous and nitrate concentrations were mostly found in areas of peat runoff, tributaries receiving treated municipal effluent, and lentic zones upstream of culverts. Peat runoff was also shown to be acidic, whether it is runoff from a harvested area or a natural bog.     

USE OF HISTORIC FLOOD INFORMATION IN ESTIMATING FLOOD PEAKS ON UNGAGED WATERSHEDS1

Regional procedures to estimate flood magnitudes for ungaged watersheds typically ignore available site-specific historic flood information such as high water marks and the corresponding flow estimates, otherwise referred to as limited site-specific historic (LSSH) flood data. A procedure to construct flood frequency curves on the basis of LSSH flood observations is presented. Simple inverse variance weighting is employed to systematically combine flood estimates obtained from the LSSH data base with those from a regional procedure to obtain improved estimtes of flood peaks on the ungaged watershed. For the region studied, the variance weighted estimates of flow had a lower logarithmic standard error than either the regional or the LSSH flow estimates, when compared to the estimates determined by three standard distributions for gaged watersheds investigated in the development of the methodology. Use of the simple inverse variance weighting procedure is recommended when “reliable” estimates of LSSH floods for the ungaged site are available.     

AVERAGE DISCHARGE,PERENNIAL FLOW INITIATION,AND CHANNEL INITIATION ‐ SMALL SOUTHERN APPALACHIAN BASINS1

ABSTRACT: Regional average evapotranspiration estimates developed by water balance techniques are frequently used to estimate average discharge in ungaged streams. However, the lower stream size range for the validity of these techniques has not been explored. Flow records were collected and evaluated for 16 small streams in the Southern Appalachians to test whether the relationship between average discharge and drainage area in streams draining less than 200 acres was consistent with that of larger basins in the size range (> 10 square miles) typically gaged by the U.S. Geological Survey (USGS). This study was designed to evaluate predictors of average discharge in small ungaged streams for regulatory purposes, since many stream regulations, as well as recommendations for best management practices, are based on measures of stream size, including average discharge. The average discharge/drainage area relationship determined from gages on large streams held true down to the perennial flow initiation point. For the southern Appalachians, basin size corresponding to perennial flow is approximately 19 acres, ranging from 11 to 32 acres. There was a strong linear relationship (R²= 0.85) between average discharge and drainage area for all streams draining between 16 and 200 acres, and the average discharge for these streams was consistent with that predicted by the USGS Unit Area Runoff Map for Georgia. Drainage area was deemed an accurate predictor of average discharge, even in very small streams. Channel morphological features, such as active channel width, cross‐sectional area, and bankfull flow predicted from Manning's equation, were not accurate predictors of average discharge. Monthly baseflow statistics also were poor predictors of average discharge.     

Linking river, floodplain, and vadose zone hydrology to improve restoration of a coastal river affected by saltwater intrusion

Floodplain forests provide unique ecological structure and function, which are often degraded or lost when watershed hydrology is modified. Restoration of damaged ecosystems requires an understanding of surface water, groundwater, and vadose (unsaturated) zone hydrology in the floodplain. Soil moisture and porewater salinity are of particular importance for seed germination and seedling survival in systems affected by saltwater intrusion but are difficult to monitor and often overlooked. This study contributes to the understanding of floodplain hydrology in one of the last bald cypress [Taxodium distichum (L.) Rich.] floodplain swamps in southeast Florida. We investigated soil moisture and porewater salinity dynamics in the floodplain of the Loxahatchee River, where reduced freshwater flow has led to saltwater intrusion and a transition to salt-tolerant, mangrove-dominated communities. Twenty-four dielectric probes measuring soil moisture and porewater salinity every 30 min were installed along two transects-one in an upstream, freshwater location and one in a downstream tidal area. Complemented by surface water, groundwater, and meteorological data, these unique 4-yr datasets quantified the spatial variability and temporal dynamics of vadose zone hydrology. Results showed that soil moisture can be closely predicted based on river stage and topographic elevation (overall Nash-Sutcliffe coefficient of efficiency = 0.83). Porewater salinity rarely exceeded tolerance thresholds (0.3125 S m(-1)) for bald cypress upstream but did so in some downstream areas. This provided an explanation for observed vegetation changes that both surface water and groundwater salinity failed to explain. The results offer a methodological and analytical framework for floodplain monitoring in locations where restoration success depends on vadose zone hydrology and provide relationships for evaluating proposed restoration and management scenarios for the Loxahatchee River.     

TREE-RING RECONSTRUCTION OF UPPER GILA RWER DISCHARGE1

ABSTRACT: Effective planning for use of water resources requires accurate information on hydrologic variability induced by climatic fluctuations. Tree-ring analysis is one method of extending our knowledge of hydrologic variability beyond the relatively short period covered by gaged streamflow records. In this paper, a network of recently developed tree-ring chronologies is used to reconstruct annual river discharge in the upper Gila River drainage in southeastern Arizona and southwestern Arizona since A.D. 1663. The need for data on hydrologic variability for this semi-arid basin is accentuated because water supply is inadequate to meet current demand. A reconstruction based on multiple linear regression (R²=0.66) indicates that 20th century is unusual for clustering of high-discharge years (early 1900s), severity of multiyear drought (1950s), and amplification of low-frequency discharge variations. Periods of low discharge recur at irregular intervals averaging about 20 years. Comparison with other tree-ring reconstructions shows that these low-flow periods are synchronous from the Gila Basin to the southern part of the Upper Colorado River Basin.     

FACTORS AFFECTING FLOODS FROM WATERSHEDS IN HUMID REGIONS OF NORTHEASTERN OHIO1

This paper describes the results of a study of hydrologic factors affecting floods from humid region in northeastern Ohio. Statistical multiple correlation analysis was used to relate floods to hydrologic and basin characteristics. Results of the study emphasize that the characteristics of floods from small and large watersheds are so significantly different that the two problems cannot be combined into one solution. The studies show that the most important hydrologic characteristics in large watersheds were: drainage area size and main channel slope. For small watersheds the most important hydrologic characteristics were: drainage area size, rainfall intensity and soil index. For watershed effect by reservoir storage it was found that: (1) small drainage areas are relatively more affected by storage than large drainage areas; (2) storage of less than 25 acre feet per square mile will not have significant effect on the mean annual flood (for drainage area above 70 square miles).