Geology ReportsSearch

USGS · b1077

Geology of the Lake Mary quadrangle, Iron County, Michigan

Abstract

The Lake Mary quadrangle is in eastern Iron County, in the west part of the Upper Peninsula of Michigan. The quadrangle is underlain by Lower and Middle Precambrian rocks, formerly designated Archean and Algonkian rocks, and is extensively covered by Pleistocene glacial deposits. A few Upper Precambrian (Keweenawan) diabase dikes and two remnants of sandstone and dolomite of early Paleozoic age are also found in the area. The major structural feature is the Holmes Lake anticline, the axis of which strikes northwest through the northeast part of the quadrangle. Most of the quadrangle, therefore, is underlain by rock of the west limb of the anticline. To the northwest along the fold axis, the Holmes Lake anticline is separated from the Amasa oval by a saddle of transverse folds in the vicinity of Michigamme Mountain in the Kiernan quadrangle. The Lower Precambrian rocks are represented by the Dickinson group and by porphyritic red granite whose relation to the Dickinson group is uncertain, but which may be older. The rocks of the Dickinson group are chiefly green to black metavolcanic schist and red felsite, some of the latter metarhyolite. The dark schist is commonly magnetic. The Dickinson group underlies the core area of the Holmes Lake anticline, which is flanked by steeply dipping Middle Precambrian formations of the Animikie series. A major unconformity separates the Lower Precambrian rocks from the overlying Middle Precambrian rocks. In ascending order the formations of the Middle Precambrian are the Randville dolomite, the Hemlock formation, which includes the Mansfield iron-bearing slate member, and the Michigamme slate. An unconformity occurs between the Hemlock formation and Michigamme slate. The post-Hemlock unconformity is thought to be represented in the Lake Mary quadrangle by the absence of iron-formation of the Amasa formation, which is known to lie between the Hemlock and the Michigamme to the northwest of the Lake Mary quadrangle in the Crystal Falls quadrangle. Post-Hemlock erosion may account also for the absence of iron-formation of the Fence River formation on the east limb of the Holmes Lake anticline within the Lake Mary quadrangle. The Randville dolomite is not exposed and is known only from diamond drilling in the northeast part of the area where it occurs in the east and west limbs of the Holmes Lake anticline. The formation has a maximum thickness of about 2,100 feet; this includes a lower arkosic phase, some of which is quartz pebble conglomerate, a medial dolomitic phase, and an upper slate phase. The triad is gradational. Included within the formation are a few beds of chloritic schist thought to be of volcanic origin. An unconformity between the Randville and the succeeding Hemlock is not indicated in the quadrangle, but is probably present. The Hemlock formation is best exposed in the northwest and south-central parts of the area. The apparent thickness of the formation is 10,000- 17,000 feet. It is composed mainly of mafic metavolcanic rocks and intercalated slate and iron-formation. In the north part of the quadrangle the volcanic rocks are greenstone, which includes altered basaltic flow rocks, volcanic breccia, tuff, and slate. Pillow structures are common in the metabasalt. It is not certain if any Hemlock rocks are present in the east limb of the Holmes Lake anticline. In the south part of the quadrangle, the rocks of the Hemlock are chiefly chlorite and hornblende schist and hornfels. Pyroxene hornfels is sparingly present. At least two sedimentary slate belts are included in the Hemlock formation. One of these, the Mansfield iron-bearing slate member, includes in its upper part an altered chert-siderite iron-formation 30 to over 150 feet thick from which iron ore has been mined at the Mansfield location. The position of the iron-bearing rocks has been determined magnetically, and past explorations for iron ore are discussed. Though probably; unconformable, the contact between the Hemlock and the Michigamme formations appears conformable. The Michigamme slate consists of at least 4,000 feet of interbedded mica schist and granulite, the altered equivalents of the slate and graywacke characteristic of the Michigamme in adjacent areas. The Michigamme rocks are best exposed in the south part of the quadrangle in the vicinity of Peavy Pond. Two periods of regional metamorphism have resulted in the alteration of almost all of the rocks of the quadrangle. The Lower Precambrian rocks underwent at least one period of metamorphism, uplift, and erosion before the deposition of the Randville dolomite. After the deposition of the Michigamme slate, a post-Middle Precambrian period of regional metamorphism occurred with attending deformation and igneous intrusion. The grade of metamorphism rises toward the south in the area. The rocks in the northern two-thirds of the quadrangle are representative of greenschist facies of regional metamorphism, whereas the rocks in the southern onethird of the quadrangle are representative of the albite-epidote-amphibolite, the amphibolite, and the pyroxene hornfels facies, the metamorphic node centering about the intrusive Peavy Pond complex in the Peavy Pond area. The Precambrian sedimentary and volcanic rocks are cut by intrusive igneous rocks of different types and several different ages. Gabbroic sills and dikes invaded the Hemlock rocks at some time after the Hemlock was deposited and before the post-Middle Precambrian orogeny and metamorphism. Some contact metamorphism attended the intrusion of the major sills. One of the sills, the West Kiernan sill, is well differentiated. A syntectonic igneous body, composed of gabbro and minor ultramafic parts and fringed with intermediate and felsic differentiates and hybrids, the Peavy; Pond complex, was intruded into the Hemlock and Michigamme formations during the post-Middle Precambrian orogeny. The complex is situated in the Peavy Pond area at the crest of the regional metamorphic node. Contact-altered sedimentary and volcanic rocks margin the complex. The effects of regional metamorphism have been superposed on the contact metamorphic rocks peripheral to the complex and on the igneous rocks of the complex as well. The mafic augite-bearing rocks of the complex emplaced early in the orogeny were deformed by granulation at the peak of the deformation and subsequently metamorphosed to hornblende rocks. Some of the intermediate and felsic rocks of the complex were foliated by the deformation, while the more fluid, felsic parts of the complex were intruded under orogenic stress and crystallized after the peak of deformation. The deformation culminated in major faulting during which the formations were dislocated, and some of the granite of the complex was extremely brecciated. A few diabase dikes, probably of Keweenawan age, have intruded the deformed and altered Animikie rocks. The only known metallic resource is iron ore. The Mansfield mine produced 1¥2 million tons of high-grade iron ore between the years 1890 and 1913. Sporadic exploration since 1913 has failed to reveal other ore deposits of economic importance.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 46° to 46.125° latitude; -88.25° to -88.125° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Richard W. Bayley. 1959. Geology of the Lake Mary quadrangle, Iron County, Michigan. https://doi.org/10.3133/b1077

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

KEEP EXPLORING

Related USGS reports

Groundwater flow model for the Little Plover River basin in Wisconsin’s Central Sands

The Little Plover River is a groundwater-fed stream in the sand plains region of central Wisconsin. In this region, sandy sediment deposited during or soon after the last glaciation forms an important unconfined sand and gravel aquifer. This aquifer supplies water for numerous high-capacity irrigation, municipal, and industrial wells that support a thriving agricultural industry. In recent years the addition of many new wells, combined with observed diminished flows in the Little Plover and other nearby rivers, has raised concerns about the impacts of the wells on groundwater levels and on water levels and flows in nearby lakes, streams, and wetlands. Diverse stakeholder groups, including well operators, Growers, environmentalists, local land owners, and regulatory and government officials have sought a better understanding of the local groundwater-surface water system and have a shared desire to balance the water needs of the he liagricultural, industrial, and urban users with the maintenance and protection of groundwater-dependent natural resources. To help address these issues, the Wisconsin Department of Natural Resources requested that the Wisconsin Geological and Natural History Survey and U.S. Geological Survey cooperatively develop a groundwater flow model that could be used to demonstrate the relationships among groundwater, surface water, and well withdrawals and also be a tool for testing and evaluating alternative water management strategies for the central sands region. Because of an abundance of previous studies, data availability, local interest, and existing regulatory constraints the model focuses on the Little Plover River watershed, but the modeling methodology developed during this study can apply to much of the larger central sands of Wisconsin. The Little Plover River groundwater flow model simulates three-dimensional groundwater movement in and around the Little Plover River basin under steady-state and transient conditions. This model explicitly includes all high-capacity wells in the model domain and simulates seasonal variations in recharge and well pumping. The model represents the Little Plover River, and other significant streams and drainage ditches in the model domain, as fully connected to the groundwater system, computes stream base flow resulting from groundwater discharge, and routes the flow along the stream channel. A separate soil-water-balance (SWB) model was used to develop groundwater recharge arrays as input for the groundwater flow model. The SWB model uses topography, soils, land use, and climatic data to estimate recharge as deep drainage from the soil zone. The SWB model explicitly includes recharge originating as irrigation water, and computes irrigation using techniques similar to those used by local irrigation operators. The groundwater flow model uses the U.S. Geological Survey’s MODFLOW modeling code which is freely available, widely accepted, and commonly used by the groundwater community. The groundwater flow model and the SWB model use identical high-resolution numerical grids having model cells 100 feet on a side, with physical properties assigned to each grid cell. This grid allows accurate geographic placement of wells, streams, and other model features. The 3-dimensional grid has three layers; layers 1 and 2 represent the sand and gravel aquifer and layer 3 represents the underlying sandstone. The distribution of material properties in the model (hydraulic conductivity, aquifer thickness, etc.) comes from previous published geologic studies of the region, updated by calibration to recent streamflow and groundwater level data. The SWB model operates on a daily time step. The groundwater flow model was calibrated to monthly stress periods with time steps ranging from 1 to 16 days. More detailed time discretization is possible. The groundwater model was calibrated to water-level and streamflow data collected during 2013 and 2014 by adjusting model parameters (primarily hydraulic conductivity, storage, and recharge) until the model produced a conditionally optimal fit between field observations and model output, subject to consistency with previously published geologic studies. Calibration was performed under both steady and transient conditions, and used a sophisticated parameter-estimation procedure (PEST) for the calibration process and to identify important model parameters. For the Little Plover River, the two most important parameters are the global recharge multiplier and the hydraulic conductivity of the stream bed. The calibrated model produces water-level and mass-balance results that are consistent with field observations and previous studies of the area. The completed model is a powerful tool for testing and demonstrating alternative water-management scenarios. Example model applications described in this report include simulating how the cumulative impacts of pumping and land-use change have affected average baseflow in the Little Plover River. Depletion-potential mapping represents a method for predicting which wells and well locations have the greatest impact on nearby surface-water resources. The completed model is publicly available, along with a companion user’s guide to assist with its operation, at http://wgnhs.org/littleplover- river-groundwater-model.

Wisconsin

The 2016 groundwater flow model for Dane County, Wisconsin

A new groundwater flow model for Dane County, Wisconsin, replaces an earlier model developed in the 1990s by the Wisconsin Geological and Natural History Survey (WGNHS) and the U.S. Geological Survey (USGS). This modeling study was conducted cooperatively by the WGNHS and the USGS with funding from the Capital Area Regional Planning Commission (CARPC). Although the overall conceptual model of the groundwater system remains largely unchanged, the incorporation of newly acquired high-quality datasets, recent research findings, and improved modeling and calibration techniques have led to the development of a more detailed and sophisticated model representation of the groundwater system. The new model is three-dimensional and transient, and conceptualizes the county’s hydrogeology as a 12-layer system including all major unlithified and bedrock hydrostratigraphic units and two high-conductivity horizontal fracture zones. Beginning from the surface down, the model represents the unlithified deposits as two distinct model layers (1 and 2). A single layer (3) simulates the Ordovician sandstone and dolomite of the Sinnipee, Ancell, and Prairie du Chien Groups. Sandstone of the Jordan Formation (layer 4) and silty dolostone of the St. Lawrence Formation (layer 5) each comprise separate model layers. The underlying glauconitic sandstone of the Tunnel City Group makes up three distinct layers: an upper aquifer (layer 6), a fracture feature (layer 7), and a lower aquifer (layer 8). The fracture layer represents a network of horizontal bedding-plane fractures that serve as a preferential pathway for groundwater flow. The model simulates the sandstone of the Wonewoc Formation as an upper aquifer (layer 9) with a bedding-plane fracture feature (layer 10) at its base. The Eau Claire aquitard (layer 11) includes shale beds within the upper portion of the Eau Claire Formation. This layer, along with overlying bedrock units, is mostly absent in the preglacially eroded valleys along the Yahara River valley and in northeastern Dane County. Layer 12 represents the Mount Simon sandstone as the lowermost model layer. It directly overlies the Precambrian crystalline basement rock, whose top surface forms the lower boundary of the model. The model uses the USGS MODFLOW-NWT finite-difference code, a standalone version of MODFLOW-2005 that incorporates the Newton (NWT) solver. MODFLOW-NWT improves the handling of unconfined conditions by smoothing the transition from wet to dry cells. The model explicitly simulates groundwater–surface-water interaction with streamflow routing and lake-level fluctuation. Model input included published and unpublished hydrogeologic data from recent estimates of aquifer hydraulic conductivities. A spatial groundwater recharge distribution was obtained from a recent GIS-based, soil-water-balance model for Dane County. Groundwater withdrawals from pumping were simulated for 572 wells across the entire model domain, which includes Dane County and portions of seven neighboring counties—Columbia, Dodge, Green, Iowa, Jefferson, Lafayette, and Rock. These wells withdrew an average of 60 million gallons per day (mgd) over the 5-year period from 2006 through 2010. Within Dane County, 385 wells were simulated with an average withdrawal rate of 52 mgd. Model calibration used the parameter estimation code PEST, and calibration targets included heads, stream and spring flows, lake levels, and borehole flows. Steady-state calibration focused on the period 2006 through 2010; the transient calibration focused on the 7-week drought period from late May through July 2012. This model represents a significant step forward from previous work because of its finer grid resolution, improved hydrostratigraphic discretization, transient capabilities, and more sophisticated representation of surface-water features and multi-aquifer wells. Potential applications of the model include evaluation of potential sites for and impacts of new high-capacity wells, development of wellhead protection plans, evaluating the effects of changing land use and climate on groundwater, and quantifying the relationships between groundwater and surface water.

Wisconsin

Geologic framework of the Alaska Peninsula, southwest Alaska, and the Alaska Peninsula terrane

The Alaska Peninsula is composed of the late Paleozoic to Quaternary sedimentary, igneous, and minor metamorphic rocks that record the history of a number of magmatic arcs. These magmatic arcs include an unnamed Late Triassic(?) and Early Jurassic island arc, the early Cenozoic Meshik arc, and the late Cenozoic Aleutian arc. Also found on the Alaska Peninsula is one of the most complete nonmetamorphosed, fossiliferous, marine Jurassic sedimentary sections known. As much as 8,500 m of section of Mesozoic sedimentary rocks record the growth and erosion of the Early Jurassic island arc. A thinner, but still thick (as much as 5,400 m), sequence of Tertiary sedimentary rocks that are predominantly continental overlies the Mesozoic section. A brief regression in early Tertiary time on the Alaska Peninsula and granodiorite plutonism in the Shumagin, Semidi, and Sanak Islands was followed by deposition of fluvial and minor marine clastic strata. This was followed by deposition of transgressive marine clastic strata and initiation of the Meshik arc, shown by an areally extensive outpouring of volcanic and volcaniclastic rocks and debris between late Eocene and earliest Miocene time. Late Miocene time was marked by another brief transgression and northwest- to southeast-directed compression, followed by renewed volcanism and plutonism which initiated the modern Aleutian magmatic arc. Extensive glacial and glaciomarine deposits of late Pleistocene age create an extensive lowland physiographic province on the northwest side of the Alaska Peninsula and join isolated mountain masses to the Alaska Peninsula on the southwest. Multiple active volcanoes and volcanic peaks dominate the skyline of the Alaska Peninsula and represent the continuation of magmatic activity that has formed the Aleutian arc since late Miocene time. The Alaska Peninsula has had a long and involved history since Paleozoic time. We propose that the Paleozoic and Mesozoic rocks that constitute much of the Alaska Peninsula be called the Alaska Peninsula terrane. Using the concept of subterranes, we divide the terrane into two distinct but tectonically related subterranes: the Chignik and Iliamna subterranes, which share a limited common geologic history. The Iliamna subterrane has served at most times as a source area for the Chignik subterrane; however, some rock units are in common across the subterranes. The Iliamna and Chignik subterranes are in part separated by the Bruin Bay fault system. The Iliamna subterrane is composed of moderately deformed early Mesozoic marine sedimentary and volcanic rocks and schist, gneiss, and marble of Paleozoic(?) and Mesozoic age, and plutonic rocks of the Alaska-Aleutian Range batholith. Characteristic of the Chignik subterrane are little-deformed, shallow-marine to continental clastic sedimentary rocks ranging in age from Permian to latest Cretaceous. However, deep-marine, volcaniclastic, and calcareous rocks form important components of the older rocks in the subterrane. The two subterranes of the Alaska Peninsula terrane are characterized by radically different structural and metamorphic styles. The nonplutonic rocks of the Iliamna subterrane are characterized by metamorphism up to amphibolite-facies grade and intense folding. In the Chignik subterrane, the structural style is dominated by large, open, en echelon anticlinal structures, normal faulting, and thrust and high-angle reverse faults that have minor displacement in a northwest to southeast direction. In the Outer Shumagin and Sanak Islands, rocks assigned to the Chugach terrane are characterized structurally by tight, generally northeast-trending folds. Dips in these rocks tend to be steep, rarely less than 35°, and overturned beds are locally common. The boundaries separating the Alaska Peninsula terrane from other terranes are commonly indistinct or poorly defined. A few boundaries have been defined at major faults, although the extensions of these faults are speculative through some areas. The west side of the Alaska Peninsula terrane is overlapped by Tertiary sedimentary and volcanic rocks and Quaternary deposits.

Alaska