Geology ReportsSearch

SEARCH · Geology Reports

Results for “Michigan Conservation”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Migration patterns and wintering distribution of common loons breeding in the Upper Midwest

Identification of geographic linkages among breeding, migratory and wintering common loon Gavia immer populations is needed to inform regional and national conservation planning efforts and compensation of loons lost during marine oil spill events. Satellite telemetry and archival geolocator tags were used to determine the migration patterns and wintering locations of breeding adult and young of the year juvenile common loons captured and marked on lakes in Minnesota, Wisconsin and Michigan. Adult loons typically traveled from breeding lakes, often via larger staging lakes, to the Great Lakes (primarily Lake Michigan) and then on to wintering areas. Most radiomarked juvenile common loons utilized natal lakes or local lakes through mid-November. Subsequently, the first fall migration of juvenile loons was generally initiated later, and more direct and quicker to wintering areas relative to adults. Among adult (n = 103) and juvenile (n = 23) loons that completed fall migration, most wintered in the Gulf of Mexico (GOM), with smaller proportions wintering off the southern Atlantic Coast or impoundments in the southeastern United States. Spring migration of adults to breeding lakes was less prolonged than fall migration, with adult male loons tending to depart wintering areas earlier than adult females. Juvenile common loons migrated during their first spring from wintering sites in the GOM to summer in the Gulf of St Lawrence/Nova Scotia Coastal region. Juvenile mortality was largely linked to parasitic infection and emaciation; spring appeared to be a survival bottleneck among juvenile loons monitored in our study. Our results identify several areas where common loon conservation efforts could be directed to protect key habitats and minimize stressors during the non-breeding period.

Michigan, Minnesota, Wisconsin

Lake fisheries need lamprey control and research

Since 1921, when the first sea lamprey was recorded from Lake Erie, concern about this parasite in the Great Lakes above Niagara Falls, where previously it had never occurred, grew successively. At first, the concern was shared only in scientific circles, but as the parasite continued its persistent and rapid spread throughout the upper Great Lakes this concern was voiced by state conservation departments, the U.S. Fish and Wildlife Service, and interested fishermen. Catches of lake trout especially, and other species secondarily, began to fall below anything representing normal fluctuations in abundance. The fishing industry on Lake Huron and Lake Michigan became extremely concerned due to the fact that income was diminishing greatly. Producers on Lake Superior were fearful that the same decline in production would soon characterize their fishery.

The Fisherman

Managing forest habitat for conservation-reliant species in a changing climate: The case of the endangered Kirtland’s Warbler

Conservation and recovery of species of concern necessitates evaluating forest habitat conditions under changing climate conditions, especially in the early stages of the delisting process. Managers must weigh implications of near-term habitat management activities within the context of changing environmental conditions and a species’ biological traits that may influence their vulnerability to changing conditions. Here we applied established population-habitat relationships based on decades of monitoring and research-management collaborations for the Kirtland’s Warbler ( Setophaga kirtlandii ) to project potential impacts of changing environmental conditions to breeding habitat distribution, quantity, and quality in the near future. Kirtland’s warblers are habitat-specialists that nest exclusively within dense jack pine ( Pinus banksiana ) forests between ca. 5–20 years of age. Using Random Forests to predict changes in distribution and growth rate of jack pine under future scenarios, results indicate the projected distribution of jack pine will contract considerably (ca. 75%) throughout the Lake States region, U.S.A. in response to projected environmental conditions in 2099 under RCP 4.5 and 8.5 climate scenarios regardless of climate model. Reduced suitability for jack pine regeneration across the Lake States may constrain management options, especially for creating high stem-density plantations nesting habitat. However, conditions remain suitable for jack pine regeneration within their historical and current core breeding range in northern Lower Michigan and several satellite breeding areas. Projected changes in jack pine growth rates varied within the core breeding area, but altered growth rates did not greatly alter the duration that habitat remained suitable for nesting by the Kirtland’s Warblers. These findings contribute to Kirtland’s Warbler conservation by informing habitat spatial planning of plantation management to provide a constant supply of nesting habitat based on the spatial variability of potential loss or gain of lands environmentally suitable for regenerating jack pine in the long-term.

Michigan, Wisconsin

Beach monitoring criteria: reading the fine print

Beach monitoring programs aim to decrease swimming-related illnesses resulting from exposure to harmful microbes in recreational waters, while providing maximum beach access. Managers are advised by the U.S. EPA to estimate microbiological water quality based on a 5-day geometric mean of fecal indicator bacteria (FIB) concentrations or on a jurisdiction-specific single-sample maximum; however, most opt instead to apply a default single-sample maximum to ease application. We examined whether re-evaluation of the U.S. EPA ambient water quality criteria (AWQC) and the epidemiological studies on which they are based could increase public beach access without affecting presumed health risk. Single-sample maxima were calculated using historic monitoring data for 50 beaches along coastal Lake Michigan on various temporal and spatial groupings to assess flexibility in the application of the AWQC. No calculation on either scale was as low as the default maximum (235 CFU/100 mL) that managers typically use, indicating that current applications may be more conservative than the outlined AWQC. It was notable that beaches subject to point source FIB contamination had lower variation, highlighting the bias in the standards for these beaches. Until new water quality standards are promulgated, more site-specific application of the AWQC may benefit beach managers by allowing swimmers greater access to beaches. This issue will be an important consideration in addressing the forthcoming beach monitoring standards.

Environmental Science & Technology

The genetic composition of wild recruits in a recovering lake trout population in Lake Michigan

Strain performance evaluations are vital for developing successful fishery management and restoration strategies. Here, we utilized genotypes from 36 microsatellites to investigate hatchery strain contribution to collections of naturally produced lake trout ( Salvelinus namaycush ) sampled across Lake Michigan. Strain composition varied by area, with recoveries of Seneca Lake strain exceeding expectations based on stocking records in northern Lake Michigan but performing similarly to other strains in southern Lake Michigan. Interstrain hybrids were present at moderate frequencies similar to expectations based on simulations, suggesting that strains are interbreeding randomly. We hypothesize that the superior performance of the Seneca Lake strain in northern Lake Michigan is partially due to adaptive advantages that facilitate increased survival in areas with high mortality from sea lamprey ( Petromyzon marinus ) predation, such as northern Lake Michigan. However, when this selective pressure is lessened, the Seneca Lake strain performs similarly to other strains. Our study demonstrates that strain performance can vary across small spatial scales and illustrates the importance of conducting thorough strain evaluations to inform management and conservation.

Lake Michigan

Tracking phosphorus and sediment sources and transport from fields and channels in Great Lakes Restoration Initiative priority watersheds

A multi-agency, integrated series of studies were initiated in 2017 under the Great Lakes Restoration Initiative (GLRI) by the U.S. Geological Survey, U.S. Forest Service, and the University of Minnesota to quantify the source, downstream travel time, and storage of particulate-bound phosphorus and sediment in agricultural tributaries to the Great Lakes. Of particular interest are contributions at the edge of field, channels, and riparian corridors. Results will be used to help identify upland and stream conservation practices that may reduce phosphorus and sediment inputs to the Great Lakes. The two study watersheds are the 50 km2 Black Creek in the Maumee River basin (Lake Erie) and the 90 km2 Plum Creek in the Lower Fox River basin (Lake Michigan). As part of other GLRI work, Black Creek and Plum Creek have existing, nested, edge-of-field studies in addition to phosphorus and sediment monitoring stations along their mainstems. Sediment-source tracking provides a direct method to quantify suspended sediment, and consequently phosphorus, sources by identifying a minimal set of properties (or fingerprint) that uniquely defines each source of sediment in the basin. This fingerprint can then be used to apportion sources of sediment from agricultural fields as well as other uplands including developed areas, forests, and pastures. These methods can also help distinguish sediment from ditches, ravines and eroding slopes, and streambanks. Multiple tracking methods are being used and adapted for best results in these watersheds, including a suite of trace elements for overall source apportionment in addition to short-term fallout radionuclides beryllium-7 (7Be) and lead-210 (210Pb) for high-flow event-based transport on fields and in stream channels. Tile drain connectivity to the surface is also of interest, especially in the Black Creek watershed. Preliminary results from the overall source apportionment from analyses of streambed sediment and monthly suspended sediment in Plum Creek indicate that the proportion attributed to different land cover varied by season and events. Further data analyses are being conducted for examining event-based pathways on individual fields, while most basin-wide sampling was monthly. Results from both watersheds will help describe the variations in transport of particulate-bound phosphorus across both steep and gentle landscapes representative of the Great Lakes basin. (this is from IP-101450)

Great Lakes

Influence of climate change and postdelisting management on long-term population viability of the conservation-reliant Kirtland's Warbler

Rapid global climate change is resulting in novel abiotic and biotic conditions and interactions. Identifying management strategies that maximize probability of long-term persistence requires an understanding of the vulnerability of species to environmental changes. We sought to quantify the vulnerability of Kirtland's Warbler ( Setophaga kirtlandii ), a rare Neotropical migratory songbird that breeds almost exclusively in the Lower Peninsula of Michigan and winters in the Bahamian Archipelago, to projected environmental changes on the breeding and wintering grounds. We developed a population-level simulation model that incorporates the influence of annual environmental conditions on the breeding and wintering grounds, and parameterized the model using empirical relationships. We simulated independent and additive effects of reduced breeding grounds habitat quantity and quality, and wintering grounds habitat quality, on population viability. Our results indicated the Kirtland's Warbler population is stable under current environmental and management conditions. Reduced breeding grounds habitat quantity resulted in reductions of the stable population size, but did not cause extinction under the scenarios we examined. In contrast, projected large reductions in wintering grounds precipitation caused the population to decline, with risk of extinction magnified when breeding habitat quantity or quality also decreased. Our study indicates that probability of long-term persistence for Kirtland's Warbler will depend on climate change impacts to wintering grounds habitat quality and contributes to the growing literature documenting the importance of considering the full annual cycle for understanding population dynamics of migratory species.

Michigan

Range-wide relative abundance of Pyrgus centaureae wyandot (W.H. Edwards, 1863) (Appalachian Grizzled Skipper)

The US Fish and Wildlife Service has designated the butterfly Pyrgus centaureae wyando t (Appalachian Grizzled Skipper [AGS]) to be at-risk, based on its declining populations and the lack of information on its status. Our objective was to conduct range-wide surveys to locate extant AGS colonies and to quantify the number of AGS observed at each location. From 2021 to 2024, we conducted 284 surveys in 25 unique (i.e., distinct) counties in 8 states (Maryland, Michigan, New York, North Carolina, Ohio, Pennsylvania, Virginia, and West Virginia). We found AGS in only 2 counties: Alleghany County, VA, and Greenbrier County, WV. AGS were observed 180 times in these counties. Our results can inform US Fish and Wildlife decisions about where and how future AGS conservation efforts can be implemented.

Northeastern Naturalists

Modeling brook trout presence and absence from landscape variables using four different analytical methods

As a part of the Great Lakes Regional Aquatic Gap Analysis Project, we evaluated methodologies for modeling associations between fish species and habitat characteristics at a landscape scale. To do this, we created brook trout Salvelinus fontinalis presence and absence models based on four different techniques: multiple linear regression, logistic regression, neural networks, and classification trees. The models were tested in two ways: by application to an independent validation database and cross-validation using the training data, and by visual comparison of statewide distribution maps with historically recorded occurrences from the Michigan Fish Atlas. Although differences in the accuracy of our models were slight, the logistic regression model predicted with the least error, followed by multiple regression, then classification trees, then the neural networks. These models will provide natural resource managers a way to identify habitats requiring protection for the conservation of fish species.

Book chapter

Range-wide relative abundance of the Appalachian grizzled skipper ( Pyrgus centaureae wyandot ) in the Eastern United States

The U.S. Fish and Wildlife Service has designated the Pyrgus centaureae wyandot (Appalachian Grizzled Skipper [AGS]) to be at-risk, based on its declining populations and the lack of information on its status. The objective of this study was to complete range-wide surveys to locate extant AGS colonies and to quantify the number of AGS observed at each location. From 2021–24, 284 surveys were done in 25 unique (that is, distinct) counties in 8 States in the Eastern United States — Maryland, Michigan, New York, North Carolina, Ohio, Pennsylvania, Virginia, and West Virginia. We found AGS in only two counties: Alleghany County, Virginia, and Greenbrier County, West Virginia. AGS were observed 180 times in these two counties. Our results can inform U.S. Fish and Wildlife decisions about where and how future AGS conservation efforts can be implemented.

Maryland, Michigan, New York, North Carolina, Ohio

Farmland in U.S. Conservation Reserve Program has unique floral composition that promotes bee summer foraging

Bee conservation is a topic of global concern, particularly in agroecosystems where their contribution to crop pollination is highly valued. Over a decade ago, bees and other pollinators were made a priority of the Conservation Reserve Program (CRP), a U.S. federal program that pays land owners to establish a conservation cover, typically grassland, on environmentally sensitive farmland. Despite large financial investment in this program, few studies have measured the benefit of CRP to bees, particularly in complex agroecosystems with abundant alternative forage. To determine if CRP land seeded with pollinator-attractive native flowers and/or introduced legumes provides distinct floral composition that attracts more foraging bees than non-CRP habitats, we compared CRP land to paired non-CRP fields and roadsides at 31 sites in Michigan, U.S.A.. We found CRP land had unique floral species community composition, higher floral abundance, greater species richness, more native floral species, and greater inflorescence coverage. Greater inflorescence coverage on CRP land was associated with a greater abundance of both honey bees and wild bees than either non-CRP fields or roadsides, as was native flower abundance for wild bees. Showy native plant species were important forage resources on CRP land: Monarda fistulosa was the most foraged upon species by both honey bees and wild bees, and goldenrod species were important late-summer forage resources for honey bees. These findings demonstrate the benefit of managing CRP land with herbaceous seed mixes to create dense, showy, native plant communities that provide summer-long resources to both bee groups. Insights from this study could be used to enhance the composition of future conservation program investments and management of non-CRP land to benefit pollinators.

Michigan

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming

Conservation and management of fisheries and aquatic communities in Great Lakes connecting channels

The North American Laurentian Great Lakes are linked by a unique series of riverine and lacustrine waters known as the Great Lakes connecting channels that are as integral to the basin's ecology and economies as the lakes themselves. The St. Marys River (SMR) is the northernmost channel and flows from Lake Superior to Lake Huron. Waters from the upper Great Lakes (Lakes Superior, Michigan, and Huron) empty from Lake Huron via the St. Clair–Detroit River system (SCDRS, also known as the Huron–Erie Corridor) into Lake Erie. The SCDRS is composed of the St. Clair River, Lake St. Clair, and the Detroit River. The Niagara River (NR) serves as the outflow from Lake Erie into Lake Ontario. The NR above Niagara Falls is bisected by Grand Island and contains several other islands and man-made embayments whereas the NR below the falls is more linear. The outflow from Lake Ontario, representing the natural outlet of all the Great Lakes, is the St. Lawrence River (SLR) which empties into the Gulf of St. Lawrence in the northwest Atlantic Ocean.

Great Lakes

eDNA metabarcoding outperforms traditional fisheries sampling and reveals fine-scale heterogeneity in a temperate freshwater lake

Understanding biodiversity in aquatic systems is critical to ecological research and conservation efforts, but accurately measuring species richness using traditional methods can be challenging. Environmental DNA (eDNA) metabarcoding, which uses high-throughput sequencing and universal primers to amplify DNA from multiple species present in an environmental sample, has shown great promise for augmenting results from traditional sampling to characterize fish communities in aquatic systems. Few studies, however, have compared exhaustive traditional sampling with eDNA metabarcoding of corresponding water samples at a small spatial scale. We intensively sampled Boardman Lake (1.4 km 2 ) in Michigan, USA, from May to June in 2019 using gill and fyke nets and paired each net set with lake water samples collected in triplicate. We analyzed water samples using eDNA metabarcoding with 12S and 16S fish-specific primers and compared estimates of fish diversity among methods. In total, we set 60 nets and analyzed 180 1 L lake water samples. We captured a total of 12 fish species in our traditional gear and detected 40 taxa in the eDNA water samples, which included all the species observed in nets. The 12S and 16S assays detected a comparable number of taxa, but taxonomic resolution varied between the two genes. In our traditional gear, there was a clear difference in the species selectivity between the two net types, and there were several species commonly detected in the eDNA samples that were not captured in nets. Finally, we detected spatial heterogeneity in fish community composition across relatively small scales in Boardman Lake with eDNA metabarcoding, but not with traditional sampling. Our results demonstrated that eDNA metabarcoding was substantially more efficient than traditional gear for estimating community composition, highlighting the utility of eDNA metabarcoding for assessing species diversity and informing management and conservation.

Michigan

Use of stable isotope signatures to determine mercury sources in the Great Lakes

Sources of mercury (Hg) in Great Lakes sediments were assessed with stable Hg isotope ratios using multicollector inductively coupled plasma mass spectrometry. An isotopic mixing model based on mass-dependent (MDF) and mass-independent fractionation (MIF) (δ 202 Hg and Δ 199 Hg) identified three primary Hg sources for sediments: atmospheric, industrial, and watershed-derived. Results indicate atmospheric sources dominate in Lakes Huron, Superior, and Michigan sediments while watershed-derived and industrial sources dominate in Lakes Erie and Ontario sediments. Anomalous Δ 200 Hg signatures, also apparent in sediments, provided independent validation of the model. Comparison of Δ 200 Hg signatures in predatory fish from three lakes reveals that bioaccumulated Hg is more isotopically similar to atmospherically derived Hg than a lake’s sediment. Previous research suggests Δ 200 Hg is conserved during biogeochemical processing and odd mass-independent fractionation (MIF) is conserved during metabolic processing, so it is suspected even is similarly conserved. Given these assumptions, our data suggest that in some cases, atmospherically derived Hg may be a more important source of MeHg to higher trophic levels than legacy sediments in the Great Lakes.

Lake Erie, Lake Huron, Lake Michigan, Lake Ontario

Classification tree models for predicting distributions of michigan stream fish from landscape variables

Traditionally, fish habitat requirements have been described from local-scale environmental variables. However, recent studies have shown that studying landscape-scale processes improves our understanding of what drives species assemblages and distribution patterns across the landscape. Our goal was to learn more about constraints on the distribution of Michigan stream fish by examining landscape-scale habitat variables. We used classification trees and landscape-scale habitat variables to create and validate presence-absence models and relative abundance models for Michigan stream fishes. We developed 93 presence-absence models that on average were 72% correct in making predictions for an independent data set, and we developed 46 relative abundance models that were 76% correct in making predictions for independent data. The models were used to create statewide predictive distribution and abundance maps that have the potential to be used for a variety of conservation and scientific purposes. ?? Copyright by the American Fisheries Society 2008.

Transactions of the American Fisheries Society

Genetic diversity of lake whitefish in lakes Michigan and Huron: sampling, standardization, and research priorities

We combined data from two laboratories to increase the spatial extent of a genetic data set for lake whitefish Coregonus clupeaformis from lakes Huron and Michigan and saw that genetic diversity was greatest between lakes, but that there was also structuring within lakes. Low diversity among stocks may be a reflection of relatively recent colonization of the Great Lakes, but other factors such as recent population fluctuation and localized stresses such as lamprey predation or heavy exploitation may also have a homogenizing effect. Our data suggested that there is asymmetrical movement of lake whitefish between Lake Huron and Lake Michigan; more genotypes associated with Lake Michigan were observed in Lake Huron. Adding additional collections to the calibrated set will allow further examination of diversity in other Great Lakes, answer questions regarding movement among lakes, and estimate contributions of stocks to commercial yields. As the picture of genetic diversity and population structure of lake whitefish in the Great Lakes region emerges, we need to develop methods to combine data types to help identify important areas for biodiversity and thus conservation. Adding genetic data to existing models will increase the precision of predictions of the impacts of new stresses and changes in existing pressures on an ecologically and commercially important species.

Journal of Great Lakes Research

Hydrologic Unit Map -- 1974, State of Michigan

This map and accompanying table show Hydrologic Units that are basically hydrographic in nature. The Cataloging Units shown will supplant the Cataloging Units previously used by the U.S. Geological Survey in its Catalog of Information on Water Data (1966-72). The previous U.S. Geological Survey Catalog-Indexing System was by map number and letter, such as 49 M. The boundaries shown have been adapted from "The Catalog of Information on Water Data" (1972), "Water Resources Regions and Subregions for the National Assessment of Water and Related Land Resources" by the U.S. Water Resources Council (1970), " River Basins of the United States" by the U.S. Soil Conservation Service (1963, 1970), "River Basin Maps Showing Hydrologic Stations" by the Inter-Agency Committee on Water Resources, Subcommittee on Hydrology (1961), and States planning maps.

Michigan