PREDICTING THE SUMMER TEMPERATURE OF SMALL STREAMS IN SOUTHWESTERN WISCONSIN1

ABSTRACT: One of the biggest challenges in managing cold water streams in the Midwest is understanding how stream temperature is controlled by the complex interactions among meteorologic processes, channel geometry, and ground water inflow. Inflow of cold ground water, shade provided by riparian vegetation, and channel width are the most important factors controlling summer stream temperatures. A simple screening model was used to quantitatively evaluate the importance of these factors and guide management decisions. The model uses an analytical solution to the heat transport equation to predict steady‐state temperature throughout a stream reach. The model matches field data from four streams in southwestern Wisconsin quite well (typically within 1°C) and helps explain the observed warming and cooling trends along each stream reach. The distribution of ground water inflow throughout a stream reach has an important influence on stream temperature, and springs are especially effective at providing thermal refuge for fish. Although simple, this model provides insight into the importance of ground water and the impact different management strategies, such as planting trees to increase shade, may have on summer stream temperature.     

ESTIMATION OF EPISODIC STREAM ACIDIFICATION BASED ON MONTHLY OR ANNUAL SAMPLING1

ABSTRACT: Programs of monthly or annual stream water sampling will rarely observe the episodic extremes of acidification chemistry that occur during brief, unpredictable runoff events. When viewed in the context of data from several streams, however, baseflow measurements of variables such as acid neutralizing capacity, pH and NO₃· are likely to be highly correlated with the episodic extremes of those variables from the same stream and runoff season. We illustrate these correlations for a water chemistry record, nearly two years in length, obtained from intensive sampling of 13 small Northeastern U.S. streams studied during USEPA's Episodic Response Project. For these streams, simple regression models estimate episodic extremes of acid neutralizing capacity, pH, NO₃·, Ca²⁺, SO₄^2?, and total dissolved Al with good relative accuracy from statistics of monthly or annual index samples. Model performances remain generally stable when episodic extremes in the second year of sampling are predicted from first-year models. Monthly or annual sampling designs, in conjunction with simple empirical models calibrated and maintained through intensive sampling every few years, may estimate episodic extremes of acidification chemistry with economy and reasonable accuracy. Such designs would facilitate sampling a large number of streams, thereby yielding estimates of the prevalence of episodic acidification at regional scales.     

USE OF FIRST-ORDER UNCERTAINTY ANALYSIS TO OPTIMIZE SUCCESSFUL STREAM WATER QUALITY SIMULATION1

ABSTRACT: A first-order uncertainty technique is developed to quantify the relationship between field data collection and a modeling exercise involving both calibration and subsequent verification. A simple statistic (LTOTAL) is used to quantify the total likelihood (probability) of successfully calibrating and verifying the model. Results from the first-order technique are compared with those from a traditional Monte Carlo simulation approach using a simple Streeter-Phelps dissolved oxygen model. The largest single difference is caused by the filtering or removal of unrealistic outcomes within the Monte Carlo framework. The amount of bias inherent in the first-order approach is also a function of the magnitude of input variability and sampling location. The minimum bias of the first-order technique is approximately 20 percent for a case involving relatively large uncertainties. However the bias is well behaved (consistent) so as to allow for correct decision making regarding the relative efficacy of various sampling strategies. The utility of the first-order technique is demonstrated by linking data collection costs with modeling performance. For a simple and inexpensive project, a wise and informed selection resulted in an LTOTAL value of 86 percent, while an uninformed selection could result in an LTOTAL value of only 55 percent.     

ERRORS RELATED TO RANDOM STREAM TEMPERATURE DATA COLLECTION IN UPPER MISSISSIPPI RIVER WATERSHED1

ABSTRACT: Records of hourly water temperatures for two streams in the Upper Mississippi River basin were used to find the error between instantaneous measurements of stream water temperatures and true daily averages. The instantaneous summer water temperature measurements were assumed to be collected during daylight hours, and measurement times were selected randomly. The absolute error at the 95 percent confidence level of randomly collected stream water temperatures was less than 0.9°C for a 1 to 5m deep large river, but as large as 3.6°C for a 0.3 to lm deep small stream. Temperature readings of morning samples were usually below daily average values, and afternoon readings were usually above. Daily mean water temperatures were obtained with less than 0.23°C standard deviation from true daily averages if the daily maximum and minimum water temperatures were averaged. Sample results were obtained for the open water (summer) season only, since diurnal water temperature fluctuations in ice covered streams are usually negligible.     

USING COMPUTED STREAM FLOW IN WATERSHED ANALYSIS1

Modifications in the computed climatic water budget have made it possible to achieve good agreement between computed and measured stream flow on both a monthly and annual basis in basins without appreciable winter snow cover. Comparisons of computed and measured stream flow in 28 basins on the Delmarva peninsula show that for some basins the agreement is excellent (regression line essentially equals unity), for other basins the regression line has a slope of one but it is displaced above or below the y=x line, while for other basins, the slope of the regression line differs appreciably from unity. Study of the basins where agreement between computed and measured values is only fair to poor reveals that the patterns of disagreement can be used to provide information on the water holding capacity in the root zone of the soil, on the quantity of deep aquifer recharge within the basin, or on the effect of human modifications within the basin. The technique should also reveal the quantity of interbasin transfers or other consumptive uses within the basin. The water budget, thus, becomes a useful tool to study hydrologic characteristics or their changes over time within a basin.     

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.     

STREAM TEMPERATURE INCREASES AND LAND USE IN A FORESTED OREGON WATERSHED1

ABSTRACT: The Salmon Creek Watershed drains 325 km² of forested terrain in the Cascade Mountains of western Oregon. Over a 30–year period (from 1955 to 1984) average daily maximum and minimum stream temperatures, calculated from the 10 warmest days of each year, have risen 6°C and 2°C, respectively. In contrast, a small decrease in maximum air temperatures was found over the same period. Regression analysis indicated a highly significant (p < 0.01) relationship between a cumulative index of forest harvesting and maximum stream temperatures. Maximum temperatures also tended to increase for several years following major peak flow events. The interaction between harvest activity (logging and road construction), changing forest and riparian management practices and the occurrence of natural hydrologic events (peak flows and associated mass soil movements) tend to obscure specific cause-and-effect relationships regarding long-term changes in maximum stream temperature.     

MATHEMATICAL GENERALIZATION OF STREAM TEMPERATURE IN CENTRAL NEW ENGLAND1

ABSTRACT: The empirical fit of an annual harmonic function to stream temperature measurements in central New England can be improved by considering a harmonic period of less than 365 days instead of 365 or 366 days. Generalized equations, developed using periodic temperature data from 27 streamflow stations, allow predictions of stream temperature at any site given (1) the mean basin altitude (E), in meters above mean sea level, and (2) station latitude (LAT), in degrees. Stream temperature t, in degrees Celsius, on day number d, in days starting with January 1, is estimated as: in which, M = 31.48 – 0.0025 (E) ? 0.4635 (LAT) with standard error of estimate of 0.62°C, and τ= 1228.88 – 21.01 (LAT) with standard error of estimate of 14.1 days.     

STREAM TEMPERATURE UNDER AN INADEQUATE BUFFER STRIP IN THE SOUTHEAST PIEDMONT1

ABSTRACT: A paired watershed experiment on the southeastern Piedmont to determine the effect of clearcutting loblolly pine on water quantity, quality, and timing has shown that stream water temperatures we increased as much as 20°F even though a partial buffer strip of trees and shrubs were left in place to shade the stream. Water time minimum stream temperatures were lowered as much as 10°F by the same treatment. A stream temperature model now in use did not predict such elevated temperatures. The authors suggest that forest cover reductions in areas of gentle land relief may elevate the temperature of shallow ground water moving to the stream, even with a substantial buffer strip Ln place.     

EFFECTS OF CLEARCUTTING AND SLASHBURNING ON STREAM TEMPERATURE IN SOUTHWESTERN BRITISH COLUMBIA1

A paired watershed study was conducted in western hemlock/western redcedar/Douglas fir forests of southwestern British Columbia to assess the effects of clearcutting and clearcutting plus slash-burning treatments on stream water characteristics. In the case of stream temperatures, both treatments increased summer temperatures as well as summer daily temperature fluctuations. These effects lasted for seven years in the case of the clearcut stream but longer in the case of the clearcut and slashburned stream. Clearcutting increased winter stream temperatures whereas slashburning caused a decrease. These changes lasted less than four years. Clearcutting and slashburning had a greater impact on stream temperatures than did clearcutting alone.     

COMPUTING THE RISKS ASSOCIATED WITH WASTELOAD) ALLOCATION MODELING1

ABSTRACT: The risks associated with a traditional wasteload allocation (WLA) analysis were quantified with data from a recent study of the Upper Trinity River (Texas). Risk is define here as the probability of failing to meet an established in-stream water quality standard. The QUAL-TX dissolved oxygen (DO) water quality model was modified to a Monte Carlo framework. Flow augmentation coding was also modified to allow an exact match to be computed between the predicted and an established DO concentration standard, thereby providing an avenue for linking input parameter uncertainty to the assignment of a wasteload permit (allowable mass loading rate). Monte Carlo simulation techniques were employed to propagate input parameter uncertainties, typically encountered during WLA analysis, to the computed effluent five-day carbonaceous biochemical oxygen demand requirements for a single major wastewater treatment plant (WWTP). The risk of failing to meet an established in-stream DO criterion may be as high as 96 percent. The uncertainty associated with estimation of the future total Kjeldahl nitrogen concentration for a single tributary was found to have the greatest impact on the determination of allowable WWTP loadings.     

MODEL ACCURACY IN SNOWMELT-RUNOFF FORECASTS EXTENDING FROM 1 TO 20 DAYS1

ABSTRACT: This paper examines the performance of snowmelt-runoff models in conditions approximating real-time forecast situations. These tests are one part of an intercomparison of models recently conducted by the World Meteorological Organization (WMO). Daily runoff from the Canadian snowmelt basin Illecille. waet (1155 km², 509–3150 m a.s.l.) was forecast for 1 to 20 days ahead. The performance of models was better than in a previous WMO project, which dealt with runoff simulations from historical data, for the following reasons: (1) conditions for models were more favorable than a real-time forecast situation because measured input data and not meteorological forecast inputs were distributed to the modelers; (2) the selected test basin was relatively easy to handle and familiar from the previous WMO project; and (3) all kinds of updating were allowed so that some models even improved their accuracy towards longer forecast times. Based on this experience, a more realistic follow-up project can be imagined which would include temperature forecasts and quantitative precipitation forecasts instead of measured data.     

STATISTICAL LOW FLOW ESTIMATION USING GIS ANALYSIS IN HUMID MONTANE REGIONS IN PUERTO RICO1

ABSTRACT: Statistical analysis of watershed parameters derived using a Geographical Information system (GIS) was done to develop equations for estimating the 7d–10yr, 30d–10yr, and 7d–2yr low flow for watersheds in humid montane regions of Puerto Rico. Digital elevation models and land use, geology, soils, and stream network coverages were used to evaluate 21 geomorphic, 10 stream channel, 9 relief, 7 geology, 4 climate, and 2 soil parameters for each watershed. To assess which parameters should be used for further investigation, a correlation analysis was used to determine the independence and collinearity among these parameters and their relationship with low flows. Multiple regression analyses using the selected parameters were then performed to develop the statistical models of low flows. The final models were selected in the basis of the Mallow Cp statistic, the adjusted R², the Press statistic, the degree of collinearity, and an analysis of the residuals. In the final models, drainage density, the ratio of length of tributaries to the length of the main channel, the percent of drainage area with northeast aspect, and the average weighted slope of the drainage were the most significant parameters. The final models had adjusted standard errors of 58.7 percent, 59.2 percent, and 48.6 percent for the 7d–10yr, 30d–10yr, and 7d–2yr low flows respectively. For comparison, the best model based on watershed parameters that can be easily measured without a GIS had an adjusted standard error of 82.8 percent.     

STREAM PHOSPHORUS EXPORTS FROM WATERSHEDS WITH CONTRASTING LAND USES IN SOUTHERN ONTARIO1

: The export of dissolved molybdate reactive phosphorus (DMRP) from 22 watersheds in the Duffin Creek drainage basin near Toronto Ontario was measured over a 25-month period. The annual average loss varied from 0.027 to 2.11 kg P/ha. Phosphorus levels in a number of watersheds were strongly influenced by effluent from a sewage treatment plant which contributed about 68 percent of the annual DMRP input to Duffin Creek. An analysis of 12 watersheds which did not contain major point pollution sources revealed that DMRP concentration and losses had a significant positive correlation with crop area and a strong negative association with forest, abandoned farm land, and area of sand + sandy loam soils. The causal relationships underlying these simple correlations are difficult to evaluate because of considerable multicollinearity between land use, soil, and topographic variables. Analysis of a mass balance for the downstream reaches of Duffin Creek indicated that there was considerable retention of phosphorus in the river channel particularly during summer low flows.     

ESTIMATION OF COASTAL WATER CONTAMINATION LONG-TERM TREND USING SEASONAL LINEAR MODELS1

ABSTRACT: Bivalves are used as bioindicators to assess trends of the chemical quality of coastal and marine environments due to their ability to concentrate chemicals. These shellfish are subject to seasonal physiological changes influencing the chemical concentration. Using quarterly data, we model concentration via linear regression with a biologically based seasonal component. This was applied to cadmium concentration measured in the blue mussel (Mytilus edulis) at three sites in the Seine estuary (Normandy, France). In this case we have a high concentration season from January to June and a “low concentration” season from July to December. This season definition was checked a posteriori, using box-and-whisker plots and a statistical test of comparison of pair-wise adjusted least-squares mean differences, and it appears to be very reasonable. We averaged data by season and across sites. Our final model (R²= 0.846 with N= 27 observations) includes highly significant terms: a season effect, which accounts for 45% of the total variability, a linear and a quadratic time term. Outliers were identified by high Studentized residual values and attributed to bias in the temporal sampling schemes. The methodology developed will further be used with other shellfish and/or other trace elements and organic chemicals.