BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

Characterizing Storm Hydrograph Rise and Fall Dynamics With Stream Stage Data1

Abstract: Storm‐flow transients (i.e., hydrograph rise and fall dynamics) may represent an important aspect of understanding streamflow dynamics. However, little is known about how temporal resolution of transient data and climate variability may color these potential indicators of hydrologic pattern or condition. Warm‐season stream stage and rainfall were monitored continuously (5 min) during the 2002 water year in eight tributaries of the Little Miami River (Ohio), which drain 17‐58 km² catchments. Rise rates generated using 5‐min data were different than those generated with mean daily data [calculated with the Indicators of Hydrologic Alteration (IHA) software], though fall rates were similar for fine and coarse temporal data. This result suggests that data with low temporal resolution may not be adequate to fully represent the dynamics of storm rise rates. Conversely, fall rates based on daily stage data (via IHA) were similar to those based on the 5‐min data, and so daily mean data may be appropriate for characterizing fall rates. We next analyzed the possible correlations between rainfall variability and storm‐flow stage dynamics. We derived rise and recession rates from storm stage hydrographs by assuming exponential rise and decay of a runoff peak. We found that raw rise rates (R_raw) were correlated with both the maximum rainfall rate and the time to the centroid of a rain event. We subsequently removed the trend based on these rainfall characteristics, which yielded new representations of rise rates abbreviated as R_rate and R_tcent, respectively, and that had lower variability than the uncorrected (raw) data. Fall rates were found to be independent of rainfall characteristics. Due to the predominant influence of stream hydrology upon aquatic biota and nutrient fluxes, our work suggests that these stage data analysis protocols can refine or otherwise reduce variability in these indices by accounting for relevant factors such as rainfall forcing. These protocols for derivation of transient indices should be tested for their potential to improve correlations between stream hydrology and temporally aligned biotic data and dissolved nutrient fluxes in streams.     

ANALYSIS AND EXTRACTION OF LOW FLOW RECESSION CHARACTERISTICS1

ABSTRACT: It is important to extract and assess low flow recession characteristics for water resources planning in the upper reaches of streams. However, it is very difficult to express synthetically the low flow recession characteristics for a stream flow. In this paper, first a new method of constructing the master recession curve based on the exponential expression is proposed and applied with the restriction that there are no regulation or diversion structures in the upper reaches above the measurement station. Daily precipitation and stream flow were used for the analysis. Second, analysis for a recession constant was conducted and the relationship between the recession constant and low flow and/or geology was qualitatively examined. In conclusion, the application of the proposed method indicated that it is objective and useful for constructing the master recession curve. It became apparent that the recession constant of a master recession curve may be defined as the total index of low flow characteristics. In addition, it was found that baseflow value increases in the order of Paleozoic, Mesozoic, Tertiary, and Quaternary.     

Analysis of Appropriate Timescales for Water Diversion Permits in Brazil

Agreement on the criteria for granting the right to use water resources between governing bodies represents a significant advance in the process of sharing water use. To aid water resource management agencies in optimizing water use, the impact of using different criteria for permitting water use in the Paracatu river basin, Brazil, was evaluated in this study. The streamflow criteria corresponding to 30 % of the annual Q_7,10 (used by the governing body of Minas Gerais), 70 % of the annual Q₉₅ (used by the governing body of the union), 30 % of the monthly Q_7,10, and 70 % of the monthly Q₉₅ were evaluated. The use of criteria based on the monthly streamflow allows for better management of water use because it allows for greater utilization of this resource in times when there is high water availability and imposes a more realistic restriction during critical periods. Substitution of the annual Q_7,10 for the monthly Q_7,10 significantly increases the streamflow permitted in some months, for example, from December to May. Use of the criterion of 70 % of the annual Q₉₅ involves a high risk of drought in critical months, while the criterion of 70 % of the monthly Q₉₅ minimizes this risk.     

Long‐Term Streamflow Response to Climatic Variability in the Loess Plateau,China1

Abstract: The Loess Plateau region in northwestern China has experienced severe water resource shortages due to the combined impacts of climate and land use changes and water resource exploitation during the past decades. This study was designed to examine the impacts of climatic variability on streamflow characteristics of a 12‐km² watershed near Tianshui City, Gansu Province in northwestern China. Statistic analytical methods including Kendall’s trend test and stepwise regression were used to detect trends in relationship between observed streamflow and climatic variables. Sensitivity analysis based on an evapotranspiration model was used to detect quantitative hydrologic sensitivity to climatic variability. We found that precipitation (P), potential evapotranspiration (PET) and streamflow (Q) were not statistically significantly different (p > 0.05) over the study period between 1982 and 2003. Stepwise regression and sensitivity analysis all indicated that P was more influential than PET in affecting annual streamflow, but the similar relationship existed at the monthly scale. The sensitivity of streamflow response to variations of P and PET increased slightly with the increase in watershed dryness (PET/P) as well as the increase in runoff ratio (Q/P). This study concluded that future changes in climate, precipitation in particular, will significantly impact water resources in the Loess Plateau region an area that is already experiencing a decreasing trend in water yield.     

AUTOMATED METHODS FOR ESTIMATING BASEFLOW AND GROUND WATER RECHARGE FROM STREAMFLOW RECORDS1

ABSTRACT: To quantify and model the natural ground water recharge process, six sites located in the midwest and eastern United States where previous water balance observations had been made were compared to computerized techniques to estimate: (1) base flow and (2) ground water recharge. Results from an existing automated digital filter technique for separating baseflow from daily streamflow records were compared to baseflow estimates made in the six water balance studies. Previous validation of automated baseflow separation techniques consisted only of comparisons with manual techniques. In this study, the automated digital filter technique was found to compare well with measured field estimates yielding a monthly coefficient of determination of 0.86. The recharge algorithm developed in this study is an automated derivation of the Rorabaugh hydrograph recession curve displacement method that utilizes daily streamflow. Comparison of annual recharge from field water balance measurements to those computed with the automated recession curve displacement method had coefficients of determination of 0.76 and predictive efficiencies of 71 percent. Monthly estimates showed more variation and are not advocated for use with this method. These techniques appear to be fast, reproducible methods for estimating baseflow and annual recharge and should be useful in regional modeling efforts and as a quick check on mass balance techniques for shallow water table aquifers.     

GENERALIZED LOW-FLOW FREQUENCY RELATIONSHIPS FOR UNGAGED SITES IN MASSACHUSETTS1

ABSTRACT: Regional hydrologic procedures such as generalized least squares regression and streamflow record augmentation have been advocated for obtaining estimates of both flood-flow and low-flow statistics at ungaged sites. While such procedures are extremely useful in regional flood-flow studies, no evaluation of their merit in regional low-flow estimation has been made using actual streamflow data. This study develops generalized regional regression equations for estimating the d-day, T-year low-flow discharge, Q_{d, t}, at ungaged sites in Massachusetts where d = 3, 7, 14, and 30 days. A two-parameter lognormal distribution is fit to sequences of annual minimum d-day low-flows and the estimated parameters of the lognormal distribution are then related to two drainage basin characteristics: drainage area and relief. The resulting models are general, simple to use, and about as precise as most previous models that only provide estimates of a single statistic such as Q_7,10. Comparisons are provided of the impact of using ordinary least squares regression, generalized least squares regression, and streamflow record augmentation procedures to fit regional low-flow frequency models in Massachusetts.     

Climate Sensitivity of Phosphorus Loadings to an Urban Stream

We investigate the sensitivity of phosphorus loading (mass/time) in an urban stream to variations in climate using nondimensional sensitivity, known as elasticity, methods commonly used by economists and hydrologists. Previous analyses have used bivariate elasticity methods to represent the general relationship between nutrient loading and a variable of interest, but such bivariate relations cannot reflect the complex multivariate nonlinear relationships inherent among nutrients, precipitation, temperature, and streamflow. Using fixed‐effect multivariate regression methods, we obtain two phosphorus models (nonparametric and parametric) for an urban stream with high explanatory power that can both estimate phosphorus loads and the elasticity of phosphorus loading to changes in precipitation, temperature, and streamflow. A case study demonstrates total phosphorus loading depends significantly on season, rainfall, combined sewer overflow events, and flow rate, yet the elasticity of total phosphorus to all these factors remains relatively constant throughout the year. The elasticity estimates reported here can be used to examine how nutrient loads may change under future climate conditions.     

ESTIMATION OF STREAMFLOW,BASE FLOW,AND NITRATE‐NITROGEN LOADS IN IOWA USING MULTIPLE LINEAR REGRESSION MODELS1

Nineteen variables, including precipitation, soils and geology, land use, and basin morphologic characteristics, were evaluated to develop Iowa regression models to predict total streamflow (Q), base flow (Qb), storm flow (Qs) and base flow percentage (%Qb) in gauged and ungauged watersheds in the state. Discharge records from a set of 33 watersheds across the state for the 1980 to 2000 period were separated into Qb and Qs. Multiple linear regression found that 75.5 percent of long term average Q was explained by rainfall, sand content, and row crop percentage variables, whereas 88.5 percent of Qb was explained by these three variables plus permeability and floodplain area variables. Qs was explained by average rainfall and %Qb was a function of row crop percentage, permeability, and basin slope variables. Regional regression models developed for long term average Q and Qb were adapted to annual rainfall and showed good correlation between measured and predicted values. Combining the regression model for Q with an estimate of mean annual nitrate concentration, a map of potential nitrate loads in the state was produced. Results from this study have important implications for understanding geomorphic and land use controls on streamflow and base flow in Iowa watersheds and similar agriculture dominated watersheds in the glaciated Midwest.     

A COMPARISON OF SIX POTENTIAL EVAPOTRANSPIRATION METHODS FOR REGIONAL USE IN THE SOUTHEASTERN UNITED STATES1

ABSTRACT: Potential evapotranspiration (PET) is an important index of hydrologic budgets at different spatial scales and is a critical variable for understanding regional biological processes. It is often an important variable in estimating actual evapotranspiration (AET) in rainfall‐runoff and ecosystem modeling. However, PET is defined in different ways in the literature and quantitative estimation of PET with existing mathematical formulas produces inconsistent results. The objectives of this study are to contrast six commonly used PET methods and quantify the long term annual PET across a physiographic gradient of 36 forested watersheds in the southeastern United States. Three temperature based (Thornthwaite, Hamon, and Hargreaves‐Samani) and three radiation based (Turc, Makkink, and Priestley‐Taylor) PET methods are compared. Long term water balances (precipitation, streamflow, and AET) for 36 forest dominated watersheds from 0.25 to 8213 km² in size were estimated using associated hydrometeorological and land use databases. The study found that PET values calculated from the six methods were highly correlated (Pearson Correlation Coefficient 0.85 to 1.00). Multivariate statistical tests, however, showed that PET values from different methods were significantly different from each other. Greater differences were found among the temperature based PET methods than radiation based PET methods. In general, the Priestley‐Taylor, Turc, and Hamon methods performed better than the other PET methods. Based on the criteria of availability of input data and correlations with AET values, the Priestley‐Taylor, Turc, and Hamon methods are recommended for regional applications in the southeastern United States.     

THE USE OF FLOOD POTENTIAL INDICES FOR FLOOD PEAK ESTIMATION ON UNGAGED WATERSHEDS1

: The construction of a flood peak index map was attempted for use by hydrologists in the simple format of rainfall maps. Since flood peaks are highly dependent on watershed area, the effect of area was removed. By regression analysis flood peaks of 2.33 and 100-year return periods were found to be proportional to watershed area to the 0.8 and 0.7 powers, respectively. Therefore, indices C_{2 33}= Q_{2 33}/_A0.7 were completed at each gage and plotted on a Pennsylvania map. It was attempted to further remove some of the scatter by regression of C with several other watershed parameters like slope, percent forest cover, and watershed shape, but no significant correlation could be found. The index maps, drawn without attenuation of the scatter, can be used by hydrologists to compute flood peaks as Q = CAⁿ (with n = 0.8 and 0.7 for the 2.33 and 100-year flood peaks, respectively). Flood peak safety factors can be based on visual observation of the index variation in the vicinity of the location for which the flood peak estimate is needed.     

EFFECTS OF LOMA PRIETA EARTHQUAKE ON SURFACE WATERS IN WADDELL VALLEY1

ABSTRACT: Waddell Creek is a perennial stream near Santa Cruz, California, which discharges into the Pacific Ocean about 38 km from the epicenter of the Loma Prieta earthquake (October 17, 1989, M 7.1). The earthquake caused no significant changes or damage to structures or surface features in the area, but was followed by an order of magnitude transient increase in discharge of Waddell Creek and nearby springs. The transients approximate a step rise followed by an exponential recession with time constants of from 37 to 59 days. The flow recession pattern of Waddell Creek was obscured by heavy rains after about 50 days; however, the springs maintained an exponential recession with minimal rain interference for periods of several time constants, then abruptly ceased flowing. A comparison of the pre-quake and post-quake chemical signature of one of the springs shows only minor changes in water chemistry indicating a common source for pre- and post-quake waters. Increased permeability of the subsurface ground material appears to be the most probable cause of the hydrologic changes. Post-seismic reduction of permeability has not been observed or inferred.     

Neural Network Input Selection for Hydrological Forecasting Affected by Snowmelt1

Abstract: Snowmelt largely affects runoff in watersheds in Nordic countries. Neural networks (NN) are particularly attractive for streamflow forecasting whereas they rely at least on daily streamflow and precipitation observations. The selection of pertinent model inputs is a major concern in NNs implementation. This study investigates performance of auxiliary NN inputs that allow short‐term streamflow forecasting without resorting to a deterministic snowmelt routine. A case study is presented for the Rivière des Anglais watershed (700 km²) located in Southern Québec, Canada. Streamflow (Q), precipitations (rain R and snow S, or total P), temperature (T) and snow lying (A) observations, combined with climatic and snowmelt proxy data, including snowmelt flow (Q_SM) obtained from a deterministic model, were tested. NN implemented with antecedent Q and R produced the largest gains in performance. Introducing increments of A and T to the NNs further improved the performance. Long‐term averages, seasonal data, and Q_SM failed to improve the networks.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

Automated Algorithms for Heuristic Base‐Flow Separation1

Abstract: The subjective nature of graphical base‐flow separation combined with the many applications of base‐flow time series derived from continuous streamflow data, motivates the development and application of automated algorithms for heuristic base‐flow separation. Base‐flow time series derived from gauged streamflow support diverse applications in engineering hydrology, catchment analysis, hydrogeologic investigations, regional low‐flow analysis, and recharge estimation. Whether based on graphical procedures for recession analysis or analytical expressions derived from fundamental equations of ground‐water flow, the variety of base‐flow separation algorithms belies the array of base‐flow definitions and interpretations that variously refer to dominant process, source, flow path, and characteristic response time. Algorithms that are invariant in their consistent – though heuristic – characterization of base‐flow response are particularly useful for interbasin comparisons of low‐flow characteristics and hydrologic regionalization. More adaptable algorithms provide application‐specific flexibility in allocating flow components like interflow to either quickflow or slowflow. Four widely used algorithms that produce consistent base‐flow time series using only gauged streamflow records are compared and contrasted with a complementary heuristic algorithm that incorporates hydrologic judgment explicitly, through manual parameterization. The utility of these inherently subjective algorithms is illustrated through a simple example of flow phase separation in a two‐component end‐member mixing model of dissolved chlorides in the Cuyahoga River.     

STATISTICAL ANALYSIS OF THE IMPACT OF GROUND WATER PUMPAGE ON LOW-FLOW HYDROLOGY1

ABSTRACT: Ground-water pumpage withdrew 57 cubic feet per second from aquifers beneath the Yahara River Basin in 1970. Forty-six cubic feet per second were exported by the diversion of treated wastewater from the drainage basin. The low-flow hydrology of the upper Yahara River has been impacted by this diversion. Prior to 1959, the wastewater was discharged into the river, augmenting the baseflow during low-flow periods. As much as 85% of streamflow was due to effluent discharge. In 1959 the wastewater was transferred from the river basin. The result was a decrease of about one-third in mean annual streamflow, and a decrease of more than 50% in the 7Q₂ and 7Q₁₀. Regression analysis showed the annual 7-day low-flow and 60-day low-flow have a statistically significant correlation with mean annual flow. Using predictions of future mean annual discharge of the river with increasing interbasin transfers, it is shown that by 1990 there is a significant probability that in some years the 60-day low-flow in the river will be zero.     

THE EFFECTS OF CLIMATE CHANGE AND IRRIGATION ON CRITERION LOW STREAMFLOWS USED FOR DETERMINING TOTAL MAXIMUM DAILY LOADS1

ABSTRACT: This paper addresses the possible impacts of global climate change on low streamflows in the Midwest, both directly, through lower precipitation, and indirectly, by rendering irrigation profitable in areas where it has found little application in the past. In the analysis presented here, streamflow data are altered to represent the effect of climate change and stream-supplied irrigation, and then used to estimate new values for two low-flow criteria, the one- and seven-day-ten-year low flows (₇Q₁₀ and ₁Q₁₀) under 20 climate change and irrigation scenarios. Additionally, the frequencies of violation of these two criteria, and multiple violations in a three-year period, are determined. Results show that the potential impact of the assumed climate change scenarios on low flow standards is substantial. A 25 percent decrease in mean precipitation results in a 63 percent reduction in design flow, even in the absence of irrigation. With irrigation, the reduction can be as much as 100 percent. The frequency of single violations of low flow criteria is found to increase several fold with irrigation. The frequency of multiple violations of low flow criteria in a three-year period is sensitive to climate change, increasing from around 20 percent to nearly 100 percent as the climate change becomes more severe.     

Basic hydrologic studies for assessing impacts of flow diversions on riparian vegetation: Examples from streams of the eastern Sierra Nevada,California, USA

As the number of proposals to divert streamflow for power production has increased in recent years, interest has grown in predicting the impacts of flow reductions on riparian vegetation. Because the extent and density of riparian vegetation depend largely on local geomorphic and hydrologic setting, site-specific geomorphic and hydrologic information is needed. This article describes methods for collecting relevant hydrologic data, and reports the results of such studies on seven stream reaches proposed for hydroelectric development in the eastern Sierra Nevada, California, USA. The methods described are: (a) preparing geomorphic maps from aerial photographs, (b) using well level records to evaluate the influence of streamflow on the riparian water table, (c) taking synoptic flow measurements to identify gaining and losing reaches, and (d) analyzing flow records from an upstream-downstream pair of gages to document seasonal variations in downstream flow losses. In the eastern Sierra Nevada, the geomorphic influences on hydrology and riparian vegetation were pronounced. For example, in a large, U-shaped glacial valley, the width of the riparian strip was highly variable along the study reach and was related to geomorphic controls, whereas the study reaches on alluvial fan deposits had relatively uniform geomorphology and riparian strip width. Flow losses of 20% were typical over reaches on alluvial fans. In a mountain valley, however, one stream gained up to 275% from geomorphically controlled groundwater contributions.     

HYDROLOGIC ANALYSIS OF DISCHARGE SUSTAINABILITY FROM AN ABANDONED UNDERGROUND COAL MINE1

ABSTRACT: Discharge from flooded abandoned subsurface coal mines is considered a potential source for water supplies where other acceptable water sources are not available. The objective of this study was to develop procedures for determining sustainability of mine‐water discharge using rainfall and discharge data for a case study site. The study site is located in southwest Virginia where Late Paleozoic sequences of sandstone, coal, and shale predominate. A rain gauge and a flow rate monitoring system were installed at the site and data were collected for a period of 100 days. The recording period corresponded with one of the driest periods in recent years and, therefore, provided valuable information regarding the flow sustainability during baseflow conditions. From available data on underground mining patterns, geology, and ground water flow regimes, it was determined that a coal mine aquifer exhibits hydraulic characteristics very similar to the extremely heterogeneous systems observed in karst aquifers, and the mine discharge is analogous to springflow. Thus, techniques commonly used in karst‐water systems and springflow analysis were used to develop rainfall/mine‐discharge relationships. Springflow recession analysis was performed on five rainfall recessions and the coefficient for each recession was compared and interpreted in light of known geologic information. It was found that the recession coefficients described the mine discharge adequately and the mine aquifer response to a rainfall pulse was very similar to the response from certain types of karst aquifers. A cross‐correlation analysis was performed to verify the results of the recession analysis and to develop a “black box” statistical model for discharge data. The correlation analysis proved the validity of springflow recession analysis for mine discharge. The recorded data length was not adequate to create a statistical model, however, but a procedure was proposed for a statistical model that could be used with large flow records. For the study site, the mine discharge was found to be sustainable for a prolonged period of time.     

CALIBRATION OF A WATER-BALANCE MODEL FOR SMALL WATERSHEDS IN EASTERN OREGON1

ABSTRACT: The BURP water-balance model was calibrated for 13 small (0.46 to 7.00 mi²), forested watersheds in the Blue Mountains of eastern Oregon where snowmelt is the dominant source of runoff. BURP is the model name and is not an acronym. Six of the 16 parameters in BURP were calibrated. The subsurface recession coefficient and three subsurface water-storage parameters were most sensitive for simulating monthly flow. Calibrated subsurface recession coefficients ranged from 0.988 to 0.998. The subsurface-water storage parameters were calibrated at between 20 to 120 percent of their initial values obtained from a category III soil survey. That reconnaissance-level survey was apparently too broad to accurately reflect subsurface-water storage in small watersheds. Tests of model performance showed BURP is capable of producing accurate simulations of monthly flow for mountainous, snow-dominated watersheds with shallow (< 4 ft) soils when calibrated with 2 to 4 years of streamflow data. A regression of observed versus simulated monthly flows with data from all watersheds combined showed that BURP accounted for 85 percent of the variability in observed flows, which ranged from 0.01 to 20.8 inches, but underpredicted high flow months, with a slope of 1.15 that is significantly different from 1.0 (p = 0.05). Without prior calibration, subsurface-water storage parameters appeared to be the greatest source of potential error.