Use of instantaneous streamflow measurements to improve regression estimates of index flow for the summer month of lowest streamflow in Michigan

David J. Holtschlag

doi:10.3133/sir20105236

Home/Archive/USGS

USGS · sir20105236

Use of instantaneous streamflow measurements to improve regression estimates of index flow for the summer month of lowest streamflow in Michigan

David J. Holtschlag

View original source ↗Export metadata ↓

Abstract

In Michigan, index flow Q50 is a streamflow characteristic defined as the minimum of median flows for July, August, and September. The state of Michigan uses index flow estimates to help regulate large (greater than 100,000 gallons per day) water withdrawals to prevent adverse effects on characteristic fish populations. At sites where long-term streamgages are located, index flows are computed directly from continuous streamflow records as GageQ50. In an earlier study, a multiple-regression equation was developed to estimate index flows IndxQ50 at ungaged sites. The index equation explains about 94 percent of the variability of index flows at 147 (index) streamgages by use of six explanatory variables describing soil type, aquifer transmissivity, land cover, and precipitation characteristics. This report extends the results of the previous study, by use of Monte Carlo simulations, to evaluate alternative flow estimators, DiscQ50, IntgQ50, SiteQ50, and AugmQ50. The Monte Carlo simulations treated each of the available index streamgages, in turn, as a miscellaneous site where streamflow conditions are described by one or more instantaneous measurements of flow. In the simulations, instantaneous flows were approximated by daily mean flows at the corresponding site. All estimators use information that can be obtained from instantaneous flow measurements and contemporaneous daily mean flow data from nearby long-term streamgages. The efficacy of these estimators was evaluated over a set of measurement intensities in which the number of simulated instantaneous flow measurements ranged from 1 to 100 at a site. The discrete measurement estimator DiscQ50 is based on a simple linear regression developed between information on daily mean flows at five or more streamgages near the miscellaneous site and their corresponding GageQ50 index flows. The regression relation then was used to compute a DiscQ50 estimate at the miscellaneous site by use of the simulated instantaneous flow measurement. This process was repeated to develop a set of DiscQ50 estimates for all simulated instantaneous measurements, a weighted DiscQ50 estimate was formed from this set. Results indicated that the expected value of this weighted estimate was more precise than the IndxQ50 estimate for all measurement intensities evaluated. The integrated index-flow estimator, IntgQ50, was formed by computing a weighted average of the index estimate IndxQ50 and the DiscQ50 estimate. Results indicated that the IntgQ50 estimator was more precise than the DiscQ50 estimator at low measurement intensities of one to two measurements. At greater measurement intensities, the precision of the IntgQ50 estimator converges to the DiscQ50 estimator. Neither the DiscQ50 nor the IntgQ50 estimators provided site-specific estimates. In particular, although expected values of DiscQ50 and IntgQ50 estimates converge with increasing measurement intensity, they do not necessarily converge to the site-specific value of Q50. The site estimator of flow, SiteQ50, was developed to facilitate this convergence at higher measurement intensities. This is accomplished by use of the median of simulated instantaneous flow values for each measurement intensity level. A weighted estimate of the median and information associated with the IntgQ50 estimate was used to form the SiteQ50 estimate. Initial simulations indicate that the SiteQ50 estimator generally has greater precision than the IntgQ50 estimator at measurement intensities greater than 3, however, additional analysis is needed to identify streamflow conditions under which instantaneous measurements will produce estimates that generally converge to the index flows. A preliminary augmented index regression equation was developed, which contains the index regression estimate and two additional variables associated with base-flow recession characteristics. When these recession variables were estimated as the medians of recession parameters compute

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 40° to 48° latitude; -90° to -82° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

David J. Holtschlag. 2011. Use of instantaneous streamflow measurements to improve regression estimates of index flow for the summer month of lowest streamflow in Michigan. https://doi.org/10.3133/sir20105236

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Karst hydrologic processes in carbonate fensters of the western Great Smoky Mountains, Tennessee, 2017–18

The geology of the Great Smoky Mountains in Tennessee and North Carolina is dominated by metamorphic strata. Along the western extent of the range, however, several carbonate fensters expose the Ordovician Jonesboro Limestone along the western edge of the Great Smoky Fault. The fensters (Cades Cove, Tuckaleechee Cove, and Wear Cove) create opportunities for allogenic recharge at the contact between the surrounding, relatively insoluble, strata and the underlying soluble formations, resulting in the development of karst systems. Karst environments have a high degree of connectivity between surface streams and groundwater systems; thus, karst aquifers are inherently vulnerable to contamination from surface activities. Karst environments also provide a specific environment for sensitive biota. Protection of the unique resources of karst environments requires an understanding of the behavior of karst aquifers, knowledge of karst recharge areas, and the mapping of groundwater flow paths. A study began in 2017 to characterize the hydrologic behavior of the karst systems in Great Smoky Mountains National Park through hydrologic monitoring, dye tracing, and seepage investigations. Hydrologic monitoring equipment was installed at several karst hydrogeologic features: two in-cave sites, two major springs, Abrams Creek in Cades Cove, and a karst depression wetland. Dye tracing was conducted to determine resurgences for the cave systems, to delineate spring recharge areas, and to estimate karst groundwater velocities. Seepage investigations along Abrams Creek helped determine if streams were sinking into the Jonesboro Limestone, either at the contact with insoluble strata or as streams flowed across the carbonate fenster. Study results indicate considerable differences in groundwater behavior between the conduit-dominated flow systems in Tuckaleechee Cove compared to the lower-gradient, colluvium-covered karst systems in Cades Cove. Study results also indicate that contaminants in the area of Bull sink and Whiteoak Sink could be rapidly transported through karst groundwater towards springs in Tuckaleechee Cove. Karst resources in Cades Cove are most likely to be affected by local events or activities taking place within the cove itself or in upper Abrams Creek upstream of the Great Smoky Fault. This study provides information on the vulnerability and uniqueness of the karst systems in the western Great Smoky Mountains, which are characterized by complex geology, high-gradient subterranean streams, deep caves, and carbonate fensters.

Tennessee↗

USGS2026

Analysis of groundwater geochemistry in the Lansing, Michigan, area

Groundwater resources in the Lansing, Michigan, area are protected by a wellhead protection plan, but previous studies have indicated that groundwater in this area may be susceptible to surface contamination. The U.S. Geological Survey, in cooperation with the Lansing Board of Water & Light, analyzed geochemical and aquifer characteristics to assess the potential for contamination of groundwater resources. Recharge age categories were determined from tritium concentrations for some wells, and by using ternary diagrams and principal component analysis for other wells. These age classifications, combined with hydrogeologic and geochemical data, were analyzed by using statistical methods to identify trends that may indicate areas of increased groundwater susceptibility. Wells with modern recharge tended to have higher sodium, chloride, and chloride:bromide ratios than did wells with mixed recharge, suggesting greater surface connectivity and potential influence from anthropogenic sources such as road salt. In contrast, wells in the mixed age category exhibited more stable geochemistry and lower constituent concentrations. Classification and regression tree models identified aquifer geometry—particularly the thickness of the uppermost bedrock unit and presence of confining layers—as key predictors of recharge age and geochemical indicators. The trends identified from these analyses could help to reveal locations in the Lansing area that are among the most susceptible to contamination.

Michigan↗

USGS2026

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

California, Idaho, Nevada, Oregon, Utah↗