USGS Science⌕ Search

USGS · sir20065031

A revised logistic regression equation and an automated procedure for mapping the probability of a stream flowing perennially in Massachusetts

Abstract

A revised logistic regression equation and an automated procedure were developed for mapping the probability of a stream flowing perennially in Massachusetts. The equation provides city and town conservation commissions and the Massachusetts Department of Environmental Protection a method for assessing whether streams are intermittent or perennial at a specific site in Massachusetts by estimating the probability of a stream flowing perennially at that site. This information could assist the environmental agencies who administer the Commonwealth of Massachusetts Rivers Protection Act of 1996, which establishes a 200-foot-wide protected riverfront area extending from the mean annual high-water line along each side of a perennial stream, with exceptions for some urban areas. The equation was developed by relating the observed intermittent or perennial status of a stream site to selected basin characteristics of naturally flowing streams (defined as having no regulation by dams, surface-water withdrawals, ground-water withdrawals, diversion, wastewater discharge, and so forth) in Massachusetts. This revised equation differs from the equation developed in a previous U.S. Geological Survey study in that it is solely based on visual observations of the intermittent or perennial status of stream sites across Massachusetts and on the evaluation of several additional basin and land-use characteristics as potential explanatory variables in the logistic regression analysis. The revised equation estimated more accurately the intermittent or perennial status of the observed stream sites than the equation from the previous study. Stream sites used in the analysis were identified as intermittent or perennial based on visual observation during low-flow periods from late July through early September 2001. The database of intermittent and perennial streams included a total of 351 naturally flowing (no regulation) sites, of which 85 were observed to be intermittent and 266 perennial. Stream sites included in the database had drainage areas that ranged from 0.04 to 10.96 square miles. Of the 66 stream sites with drainage areas greater than 2.00 square miles, 2 sites were intermittent and 64 sites were perennial. Thus, stream sites with drainage areas greater than 2.00 square miles were assumed to flow perennially, and the database used to develop the logistic regression equation included only those stream sites with drainage areas less than 2.00 square miles. The database for the equation included 285 stream sites that had drainage areas less than 2.00 square miles, of which 83 sites were intermittent and 202 sites were perennial. Results of the logistic regression analysis indicate that the probability of a stream flowing perennially at a specific site in Massachusetts can be estimated as a function of four explanatory variables: (1) drainage area (natural logarithm), (2) areal percentage of sand and gravel deposits, (3) areal percentage of forest land, and (4) region of the state (eastern region or western region). Although the equation provides an objective means of determining the probability of a stream flowing perennially at a specific site, the reliability of the equation is constrained by the data used in its development. The equation is not recommended for (1) losing stream reaches or (2) streams whose ground-water contributing areas do not coincide with their surface-water drainage areas, such as many streams draining the Southeast Coastal Region-the southern part of the South Coastal Basin, the eastern part of the Buzzards Bay Basin, and the entire area of the Cape Cod and the Islands Basins. If the equation were used on a regulated stream site, the estimated intermittent or perennial status would reflect the natural flow conditions for that site. An automated mapping procedure was developed to determine the intermittent or perennial status of stream sites along reaches throughout a basin. The procedure delineates the drainage area boundaries, determines values for the four explanatory variables, and solves the equation for estimating the probability of a stream flowing perennially at two locations on a headwater (first-order) stream reach-one near its confluence or end point and one near its headwaters or start point. The automated procedure then determines the intermittent or perennial status of the reach on the basis of the calculated probability values and a probability cutpoint (a stream is considered to flow perennially at a cutpoint of 0.56 or greater for this study) for the two locations or continues to loop upstream or downstream between locations less than and greater than the cutpoint of 0.56 to determine the transition point from an intermittent to a perennial stream. If the first-order stream reach is determined to be intermittent, the procedure moves to the next downstream reach and repeats the same process. The automated procedure then moves to the next first-order stream and repeats the process until the entire basin is mapped. A map of the intermittent and perennial stream reaches in the Shawsheen River Basin is provided on a CD-ROM that accompanies this report. The CD-ROM also contains ArcReader 9.0, a freeware product, that allows a user to zoom in and out, set a scale, pan, turn on and off map layers (such as a USGS topographic map), and print a map of the stream site with a scale bar. Maps of the intermittent and perennial stream reaches in Massachusetts will provide city and town conservation commissions and the Massachusetts Department of Environmental Protection with an additional method for assessing the intermittent or perennial status of stream sites.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Gardner C. Bent, Peter A. Steeves. 2006. A revised logistic regression equation and an automated procedure for mapping the probability of a stream flowing perennially in Massachusetts. https://doi.org/10.3133/sir20065031

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

KEEP EXPLORING

Related USGS reports

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.

Scientific Investigations Report↗

Groundwater budget for the Portneuf River Basin, southeastern Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, developed the first groundwater budget for the upper Portneuf River Basin in southeastern Idaho, a tributary of the Snake River. A groundwater budget for the Portneuf River Basin supports water resource and administrative management decisions within the basin and for downstream users of the eastern Snake River Plain aquifer. The analysis quantifies inflows, outflows, and residuals for an aquifer system defined by the Idaho Department of Water Resources under three recent precipitation conditions: an average year (2016), a dry year (2021), and a wet year (2023). The budget was calculated for the basin and for two subbasins delineated by the U.S. Geological Survey streamgage Portneuf River at Topaz, Idaho (13073000): the above Topaz subbasin, which includes the Portneuf-Gem valley, and the below Topaz subbasin, which includes Marsh Valley. Aquifer inflows include mountain-front recharge; areal recharge; recharge from applied water, including recharge from irrigation, canal seepage, and effluent returns; groundwater flow into the basin; and reservoir seepage. Outflows include groundwater pumping for agricultural, municipal, and domestic uses; subsurface outflow through the Portneuf gap; and discharge to streams as base flow. Estimated annual inflows for the full basin ranged from 116,880 acre-feet in 2021 to 154,880 acre-feet in 2023, and outflows ranged from 132,190 acre-feet in 2016 to 158,490 acre-feet in 2023. Basinwide residuals, the change in storage plus uncertainty, were −2,410 acre-feet in 2016, −18,410 acre-feet in 2021, and −3,610 acre-feet in 2023, and each subbasin had a mixture of negative and positive residuals. These basinwide results, interpreted as change in storage plus water-budget uncertainty, indicate minor net losses that increase with dry conditions, consistent with well depth-to-water measurements indicating gradual multi-decadal declines. Patterns in subbasin residuals are driven by irrigation practices and surface-water management; for example, the above Topaz subbasin relies more on groundwater for irrigation, resulting in more negative residuals during dry years. In contrast, the below Topaz subbasin depends more on surface water, producing positive residuals during dry years because of irrigation losses and canal seepage. The analyses for both subbasins indicate decades-long groundwater declines, which agree with the overall negative residuals calculated. The largest sources of uncertainty include climate inputs (precipitation and evapotranspiration) and poorly constrained components such as groundwater inflow and tributary stream infiltration. These findings provide insight into how climatic variables and irrigation practices affect groundwater inflows and outflows within the Portneuf River Basin.

Idaho↗

World minerals outlook—Antimony, copper, graphite, niobium, potash, silicon, silver, and tungsten through 2030

Demand for mineral commodities continues to increase owing to rapid technological advancements and global economic growth. To support this increase in demand, information about current and near-term supplies of these mineral commodities is essential. Such data help evaluate the health and stability of global supply chains and protect against potential supply risks. Current production and production capacity data coupled with information regarding near-term expected changes (such as facility closures, planned expansions, or the completion of new mineral production projects) can provide insight into expected future market developments and expected available supply of mineral commodities. This report within the “World Minerals Outlook” series analyzes data for antimony, copper, natural graphite, niobium, potash, silicon metal, silver, and tungsten—eight mineral commodities listed in the 2025 List of Critical Minerals. World production data from 2019 through 2025, as well as estimated world capacity data from 2023 through 2025 were analyzed. Using these data, near-term (2026–2030) expected world capacities were estimated. For antimony, copper, niobium, silicon metal, silver, and tungsten, modest increases in production capacity of 4 to 10 percent were expected between 2025 and 2030, whereas capacities for natural graphite and potash were expected to increase by about 18 to 20 percent.

Scientific Investigations Report↗