RASTER-BASED HYDROLOGIC MODELING OF SPATIALLY-VARIED SURFACE RUNOFF1

ABSTRACT: The proliferation of watershed databases in raster Geographic Information System (GIS) format and the availability of radar-estimated rainfall data foster rapid developments in raster-based surface runoff simulations. The two-dimensional physically-based rainfall-runoff model CASC2D simulates spatially-varied surface runoff while fully utilizing raster GIS and radar-rainfall data. The model uses the Green and Ampt infiltration method, and the diffusive wave formulation for overland and channel flow routing enables overbank flow storage and routing. CASC2D offers unique color capabilities to display the spatio-temporal variability of rainfall, cumulative infiltrated depth, and surface water depth as thunderstorms unfold. The model has been calibrated and independently verified to provide accurate simulations of catchment response to moving rainstorms on watersheds with spatially-varied infiltration. The model can accurately simulate surface runoff from flashfloods caused by intense thunderstorms moving across partial areas of a watershed.     

EFFECTS OF A NORTH CENTRAL WASHINGTON WILDFIRE ON RUNOFF AND SEDIMENT PRODUCTION1

ABSTRACT: Runoff was measured from a 564-ha catchment located on the Entiat Experimental Forest for nine years before a severe wild-fire in 1970 destroyed the mixed conifer vegetation. Runoff records from the Chelan River (2 393 km²) were used as control data for determining changes in water yield during the seven years following the fire. The first post-fire year was a period of transition in which the soil profile retained more water than in previous years and measured runoff was 8.9 cm greater than the predicted value based on pre-fire conditions. Runoff from the burned catchment during subsequent years was much greater than measured values before the fire. Measured minus predicted runoff, based on the pre-fire calibration equation, varied from 10.7 cm during the dry year of 1977 to 47.2 cm during the abnormally wet year of 1972. Flow duration curves indicated that runoff at each percent value after the fire was at least double the comparable pre-fire value. Sediment production increased dramatically after the fire because of increased flow rates, increased overland flow caused by reduced infiltration capacity, and mass soil movement. Sediment yield is beginning to decrease as stream channels become stabilized and vegetation on upper slopes improves infiltration capacity.     

DISCHARGE CHARACTERISTICS OF PERFORATED PIPE FOR USE IN INFILTRATION TRENCHES1

ABSTRACT: Infiltration trenches are an effective stormwater management alternative for the control of urban runoff from small areas. Perforated pipes buried within the gravel of an infiltration trench are used to distribute the inflowing runoff along the length of the trench. Laboratory tests are described that characterize the hydraulics of the orifices in perforated pipes. The results show that the steady-state exfiltration of water from the pipe into a surrounding gravel trench can be described by the orifice equation.     

ESTIMATION OF SURFACE WATER LAG TIME FROM THE KINEMATIC WAVE EQUATIONS1

ABSTRACT Numerous concepts of surface water lag time have been developed and applied in the past. In this report, hydraulic solutions of a lag time derived by Overton [1970] are presented for several idealistic surfaces using the kinematic wave equations. These surfaces are: (1) a uniform plane; (2) hillslope as a cascade of planes; (3) V-shaped watershed; (4) V-shaped watershed with hill-slopes; (5) converging surface; (6) concave surface. The lag times are shown to be related to roughness, length and catchment slope, and the input rate. These relations may be used immediately in predicting lag time as the parameter in a unit response function. A lag relation has been developed for a nonuniform catchment in terms of the lag of a uniform plane and a convergence factor. A numerical procedure is shown whereby the convergence factor can be evaluated for any nonuniform catchment from observed input and output data.     

Mechanistic Simulation of Tree Effects in an Urban Water Balance Model1

Abstract: A semidistributed, physical‐based Urban Forest Effects – Hydrology (UFORE‐Hydro) model was created to simulate and study tree effects on urban hydrology and guide management of urban runoff at the catchment scale. The model simulates hydrological processes of precipitation, interception, evaporation, infiltration, and runoff using data inputs of weather, elevation, and land cover along with nine channel, soil, and vegetation parameters. Weather data are pre‐processed by UFORE using Penman‐Monteith equations to provide potential evaporation terms for open water and vegetation. Canopy interception algorithms modified established routines to better account for variable density urban trees, short vegetation, and seasonal growth phenology. Actual evaporation algorithms allocate potential energy between leaf surface storage and transpiration from soil storage. Infiltration algorithms use a variable rain rate Green‐Ampt formulation and handle both infiltration excess and saturation excess ponding and runoff. Stream discharge is the sum of surface runoff and TOPMODEL‐based subsurface flow equations. Automated calibration routines that use observed discharge has been coupled to the model. Once calibrated, the model can examine how alternative tree management schemes impact urban runoff. UFORE‐Hydro model testing in the urban Dead Run catchment of Baltimore, Maryland, illustrated how trees significantly reduce runoff for low intensity and short duration precipitation events.     

Spatial variability of soil phosphorus in relation to the topographic index and critical source areas: sampling for assessing risk to water quality

A measure of soil P status in agricultural soils is generally required for assisting with prediction of potential P loss from agricultural catchments and assessing risk for water quality. The objectives of this paper are twofold: (i) investigating the soil P status, distribution, and variability, both spatially and with soil depth, of two different first-order catchments; and (ii) determining variation in soil P concentration in relation to catchment topography (quantified as the "topographic index") and critical source areas (CSAs). The soil P measurements showed large spatial variability, not only between fields and land uses, but also within individual fields and in part was thought to be strongly influenced by areas where cattle tended to congregate and areas where manure was most commonly spread. Topographic index alone was not related to the distribution of soil P, and does not seem to provide an adequate indicator for CSAs in the study catchments. However, CSAs may be used in conjunction with soil P data for help in determining a more "effective" catchment soil P status. The difficulties in defining CSAs a priori, particularly for modeling and prediction purposes, however, suggest that other more "integrated" measures of catchment soil P status, such as baseflow P concentrations or streambed sediment P concentrations, might be more useful. Since observed soil P distribution is variable and is also difficult to relate to nationally available soil P data, any assessment of soil P status for determining risk of P loss is uncertain and problematic, given other catchment physicochemical characteristics and the sampling strategy employed.     

SEASONAL PUMPING VARIATION EFFECTS ON WELLHEAD PROTECTION AREA DELINEATION1

ABSTRACT: Delineation of a welihead protection area (WHPA) is the key element in welihead protection programs for drinking water supplies. WHPAs are often delineated under idealized conditions using simple steady-state assumptions, which lead to an incorrect estimation of area and geometry. In this paper, we compare the results from a simple steady-state procedure commonly employed in WHPA delineation with a more complex transient approach that allows consideration of seasonal variation in pumping rates. We also introduce a transient procedure to delineate time-related capture zones using a numerical ground water flow and transport model. Welihead delineation is examined for two municipal wells in Tipton, Oklahoma, using a ten-year time-of-travel criterion. In the steady-state procedure, where we assumed constant pumping rates, we used GPTRAC, a semi-analytical model, and MOC, a numerical model. The capture zone delineated by GPTRAC is comparable in shape with the capture zone delineated by MOC but not in size due to the differences in solution schemes. In the transient procedure, we used MOC and considered the seasonal variation in pumping rates. The capture zones delineated in this procedure were larger than the capture zones delineated by the steady-state procedure using the same model. Further analysis showed that a higher drawdown was predicted in the transient procedure than in the steady-state procedure, which is the reason for larger capture zones.     

APPLICATION OF SCS INFILTRATION MODEL1

ABSTRACT: The SCS infiltration model was applied to the Ralston Creek watershed in eastern Iowa. The criteria to determine the various model parameters were revised to obtain a better agreement between the observed and computed total runoffs. A procedure to calibrate the infiltration model is presented. The infiltration model was used in conjunction with an overland flow model to develop flood hydrographs. The results indicate that SCS infiltration model adequately describe the distribution of losses.     

OBSERVATIONS ON KINEMATIC RESPONSE IN URBAN RUNOFF MODELS1

ABSTRACT: This paper first discusses the results of sensitivity analyses conducted on various parameters of the San Francisco Stormwater Model ta version of WREM) and the Penn State Runoff Model in terms of their impact on outflow hydrographs. The parameters considered within a idealized catchment include: basin shape, imperivous fraction, overland roughness and slope: deterntion depth; infiltration capacity; and hyetograph timing. Second, the results for the hypothetical catchment are extended to the lazzard laboratory surfaces (asphalt, grass, roofing material) as a mean of illustrating the need for changes in model structure, as opposed to continued parameter adjustment Finally the effect of altering the scale of hydraulic representation in the surface runoff and sewer transport calculations are demonstrated for two gaged watersheds in Hamburg, West Germany.     

WATER HARVESTING TREATMENT EVALUATION AT GRANITE REEF1

ABSTRACT: Yearly runoff efficiencies (total runoff/total precipitation), threshold retentions (precipitation needed to initiate runoff), and runoff-efficiencies-after-thresholds were determined for several water-harvesting catchment treatments at the Granite Reef test site. This information was found to be useful for showing (1) overall performance of catchments with time; (2) the distribution of the precipitation among runoff, surface retention, and infiltration; (3) why, how, and when certain treatments weathered and failed; (4) when to repair treatments; and (5) how to design catchments (size, site preparation, material selection, etc.). New impermeable membranes with smooth surfaces yielded nearly 100% of the precipitation. An asphalt-fiberglass treatment continues this high efficiency after 10 years of weathering – polyethylene after 8; efficiency of butyl sheeting was high initially but decreased rapidly after 9 years weathering. A standard roofing treatment retained up to 30% of the precipitation in the gravel covering. A concrete catchment lost as much as 50% of the total precipitation through micropores and surface cracks. Silicone treatments rapidly lost repellancy and efficiency, while paraffin treatments have weathered 5 years with little loss of efficiency.     

Characterizing long-term hydrologic-response and sediment-transport for the R-5 catchment

Recently there have been several calls to establish long-term data collection networks to monitor near-surface hydrologic response and landscape evolution. The focus of this paper is a long-term dataset from the International Hydrologic Decade (1965-1974). The small upland catchment, known as R-5, located near Chickasha, Olahoma, has been the subject of considerable attention within the event-based hydrologic modeling community for more than 30 yr. Here, for the first time, 8 yr of continuous near-surface hydrologic-response and sediment-transport data are analyzed to show trends in the catchment's long-term behavior. The datasets include precipitation, temperature, solar radiation, soil-water content, infiltration, water discharge, and sediment discharge. Potential and actual evapotranspiration rates were estimated and used to calculate an average annual water balance for the catchment. Findings include, for example, that rainfall intensity rarely exceeds the threshold for Horton-type runoff, soil-water content is both spatially and temporally variable, and the water and sediment discharge rates are positively correlated. The R-5 data provide a unique opportunity to test (and refine) process-based models of continuous hydrologic response and sediment transport at the catchment scale for applications in the emerging fields of hydroecology and hydrogeomorphology.     

USING SOIL TEXTURE TO ESTIMATE SATURATED HYDRAULIC CONDUCTIVITY AND THE IMPACT ON RAINFALL-RUNOFF SIMULATIONS1

ABSTRACT: In this paper a new set of soil texture data is used to estimate the spatial distribution of saturated hydraulic conductivity values for a small rangeland catchment. The estimates of conductivity are used to re-excite and re-evaluate a quasi-physically based rainfall-runoff model. The performance of the model is significantly reduced with conductivity estimates gleaned from soil texture data rather than the infiltration data used in our previous efforts.     

MODERNIZATION OF NATIONAL WEATHER SERVICE RIVER FORECASTING TECHNIQUES1

The National Weather Service is nearing the conclusion of a five year period of transition from index type catchment modelling to the use of conceptual hydrologic models. The decision to make this technological change was based on an extensive research project in which various catchment models were tested in a wide variety of basins and their strong and weak points ascertained. This project is described. Some of the problems involved in the changeover, which are discussed, are practical parameter optimization methods, computer requirements for the more complex technology, data requirements, fitting of the catchment model to major river systems, training of personnel and staffing problems.     

A MODEL OF THE BEHAVIOR OF pH DETERMINING PARAMETERS IN IMPOUNDMENTS1

ABSTRACT: A computer model capable of describing the behavior of pH-determining parameters in an impoundment is presented. The steady-state model, based upon principles of water chemistry, applies the complete-mix concept to the routing of the chemical quantities, alkalinity and CO₂ acidity. An iterative procedure is used to account for surface mass transfer of carbon dioxide. A fourth-degree polynomial is solved for the hydrogen ion concentration at each point of routed alkalinity and CO₂ acidity concentrations. Field data required to verify and apply the model are described, including an example application to the Loyalhanna Creek Reservoir.     

EVALUATION OF HYDROLOGIC BENEFITS OF INFILTRATION BASED URBAN STORM WATER MANAGEMENT1

ABSTRACT: As watersheds are urbanized, their surfaces are made less pervious and more channelized, which reduces infiltration and speeds up the removal of excess runoff. Traditional storm water management seeks to remove runoff as quickly as possible, gathering excess runoff in detention basins for peak reduction where necessary. In contrast, more recently developed “low impact” alternatives manage rainfall where it falls, through a combination of enhancing infiltration properties of pervious areas and rerouting impervious runoff across pervious areas to allow an opportunity for infiltration. In this paper, we investigate the potential for reducing the hydrologic impacts of urbanization by using infiltration based, low impact storm water management. We describe a group of preliminary experiments using relatively simple engineering tools to compare three basic scenarios of development: an undeveloped landscape; a fully developed landscape using traditional, high impact storm water management; and a fully developed landscape using infiltration based, low impact design. Based on these experiments, it appears that by manipulating the layout of urbanized landscapes, it is possible to reduce impacts on hydrology relative to traditional, fully connected storm water systems. However, the amount of reduction in impact is sensitive to both rainfall event size and soil texture, with greatest reductions being possible for small, relatively frequent rainfall events and more pervious soil textures. Thus, low impact techniques appear to provide a valuable tool for reducing runoff for the events that see the greatest relative increases from urbanization: those generated by the small, relatively frequent rainfall events that are small enough to produce little or no runoff from pervious surfaces, but produce runoff from impervious areas. However, it is clear that there still needs to be measures in place for flood management for larger, more intense, and relatively rarer storm events, which are capable of producing significant runoff even for undeveloped basins.     

STEADY-STATE ANALYSIS OF INFILTRATION AND OVERLAND FLOW FOR SPATIALLY-VARIED HILLSLOPES1

ABSTRACT: An envelope of steady-state surface runoff response for a hilislope is established in terms of the probability distribution and spatial arrangement of individual point infiltration capacities and the rainfall intensity. Minimum overland flow is shown to occur when point infiltration capacities are ordered with the highest at the slope bottom, while maximum overland flow occurs when the highest point capacities are at the top of the slope. Equations for envelope curves are developed for both continuous distributions and discretely sampled data; examples for each case are given. Use of the analysis as a rainfall-runoff model is also discussed.     

Spatial Variability of Nitrate Concentrations Under Diverse Conditions in Tributaries to a Lake Watershed1

Abstract: Nitrate‐nitrogen (NO₃‐N) concentrations in stream water often respond uniquely to changes in inter‐annual conditions (e.g., biological N uptake and precipitation) in individual catchments. In this paper, we assess (1) how the spatial distribution of NO₃‐N concentrations varies across a dense network of nonnested catchments and (2) how relationships between multiple landscape factors [within whole catchments and hydrologically sensitive areas (HSAs) of the catchments] and stream NO₃‐N are expressed under a variety of annual conditions. Stream NO₃‐N data were collected during two synoptic sampling events across >55 tributaries and two synoptic sampling periods with >11 tributaries during summer low flow periods. Sample tributaries drain mixed land cover watersheds ranging in size from 0.150 to 312 km² and outlet directly to Cayuga Lake, New York. Changes in NO₃‐N concentration ratios between each sampling event suggest a high degree of spatial heterogeneity in catchment response across the Cayuga Lake Watershed, ranging from 0.230 to 61.4. Variations in NO₃‐N concentrations within each of the large synoptic sampling events were also high, ranging from 0.040 to 8.7 mg NO₃‐N/l (March) and 0.090 to 15.5 mg NO₃‐N/l (October). Although Pearson correlation coefficients suggest that this variability is related to multiple landscape factors during all four sampling events, partial correlations suggest percentage of row crops in the catchments as the only similar factor in March and October and catchment area as the only factor during summer low flows. Further, the strength of the relationships is typically lower in the HSAs of catchment. Advancing current understanding of such variations and relationships to landscape factors across multiple catchments – and under a variety of biogeochemical and hydrological conditions – is important, as (1) nitrate continues to be employed as an indicator of regional aquatic ecosystem health and services and (2) a unified framework approach for understanding individual catchment processes is a rapidly evolving focus for catchment‐based science and management.     

Phosphorus load to surface water from bank erosion in a Danish lowland river basin

Phosphorus loss from bank erosion was studied in the catchment of River Odense, a lowland Danish river basin, with the aim of testing the hypothesis of whether stream banks act as major diffuse phosphorus (P) sources at catchment scale. Furthermore, the study aimed at analyzing the impact of different factors influencing bank erosion and P loss such as stream order, anthropogenic disturbances, width of uncultivated buffer strips, and the vegetation of buffer strips. A random stratified procedure in geographical information system (GIS) was used to select two replicate stream reaches covering different stream orders, channelized vs. naturally meandering channels, width of uncultivated buffer strips (≤ 2 m and ≥ 10 m), and buffer strips with different vegetation types. Thirty-six 100-m stream reaches with 180 bank plots and a total of 3000 erosion pins were established in autumn 2006, and readings were conducted during a 3-yr period (2006-2009). The results show that neither stream size nor stream disturbance measured as channelization of channel or the width of uncultivated buffer strip had any significant ( < 0.05) influence on bank erosion and P losses during each of the 3 yr studied. In buffer strips with natural trees bank erosion was significantly ( < 0.05) lower than in buffer strips dominated by grass and herbs. Gross and net P input from bank erosion amounted to 13.8 to 16.5 and 2.4 to 6.3 t P, respectively, in the River Odense catchment during the three study years. The net P input from bank erosion equaled 17 to 29% of the annual total P export and 21 to 62% of the annual export of P from diffuse sources from the River Odense catchment. Most of the exported total P was found to be bioavailable (71.7%) based on a P speciation of monthly suspended sediment samples collected at the outlet of the river basin. The results found in this study have a great importance for managers working with P mitigation and modeling at catchment scale.