Simulation of regional ground-water flow in the Cambrian-Ordovician aquifer system in the northern Midwest, United States: in Regional aquifer-system analysis

R.J. Mandle; A.L. Kontis

doi:10.3133/pp1405c

Home/Archive/USGS

USGS · pp1405C

Simulation of regional ground-water flow in the Cambrian-Ordovician aquifer system in the northern Midwest, United States: in Regional aquifer-system analysis

R.J. Mandle·A.L. Kontis

View original source ↗Read source document ↗Export metadata ↓

Abstract

A six-State area in the northern Midwest of the United States has become increasingly dependent on ground water from the Cambrian-Ordovician aquifer system, which consists of a sequence of sandstones, carbonate rocks, and shales. Ground-water withdrawals from the aquifer system began in the late 1800's and increased to approximately 684 million gallons per day (1,058 cubic feet per second) by 1980. This withdrawal has caused more than 900 feet of decline in the potentiometric surface of the aquifer system in parts of the Chicago, Ill., area. Pumping in Minneapolis-St. Paul, Minn., Milwaukee, Wis., and central Iowa has produced a few hundreds of feet of water-level decline. A quasi-three-dimensional ground-water flow model was developed to improve understanding of the regional ground-water flow system by simulating the aquifer system under the conditions that existed before and during ground-water development. The Cambrian-Ordovician aquifer system and the overlying rocks were incorporated in the conceptual model as five aquifer layers with four intervening confining layers. The aquifer layers, from top to bottom, are as follows: Quaternary deposits and Cretaceous rocks (aquifer layer 5); basal Devonian carbonate rocks and underlying Silurian carbonate rocks (aquifer layer 4); Middle Ordovician St. Peter Sandstone, Lower Ordovician Prairie du Chien Group, and Upper Cambrian Jordan Sandstone (aquifer layer 3); Upper Cambrian Ironton and Galesville Sandstones (aquifer layer 2); and Upper Cambrian Mount Simon Sandstone and Precambrian Hinckley Sandstone (aquifer layer 1). The effects on the flow system of ground water having variable density and the individual aquifer layer contribution of flow to wells open to several aquifer layers were simulated by incorporating appropriate terms in the ground-water flow equation and corresponding modifications in the ground-water flow model. Results of steady-state simulations are shown as maps of freshwater head and as flow-rate vectors. A comparison was made of available predevelopment head data and simulated predevelopment heads for the Mount Simon, St. Peter-Prairie du Chien-Jordan, and Silurian-Devonian aquifers. The root mean square of the differences between observed and simulated heads for these aquifers was 53, 36, and 44 feet, respectively. Steady-state model simulation indicates that regional recharge areas are located in northwestern Iowa, southeastern Minnesota, much of Wisconsin, northern Illinois, and central Missouri. Regional discharge areas are located near the Mississippi River and its tributaries, the Missouri River, Lake Michigan, and the Illinois basin. Results of a transient simulation of the period 1861 to 1980 are shown as maps of freshwater head and freshwater-head decline. The root mean square of the differences between observed and simulated heads for the St. Peter-Prairie du Chien-Jordan aquifer was 63 feet. In addition, hydrographs of simulated hydraulic head were compared with hydrographs of observed head in wells open to various combinations of aquifers. The numerous head measurements within a model node display a wide range in magnitude; however, their historical trend generally follows the trend of the simulated heads. The simulated recharge from the glacial drift to the immediately underlying bedrock aquifers averages 0.03, 0.06, 0.24, and 0.02 inch per year, respectively, to the Mount Simon, Ironton-Galesville, St. Peter- Prairie du Chien-Jordan, and Silurian-Devonian aquifers for predevelopment conditions and 0.03, 0.07, 0.45, and 0.07 inch per year, respectively, to the same aquifers for 1976-80. These recharge rates are less than 1.5 percent of average annual precipitation of about 30 inches per year. Most of the recharge from precipitation discharges to streams as base flow through local and intermediate ground-water flow systems. Only a small fraction of the precipitation recharges the deeper, regional flow system. The simulated predevelopment recharge of 571 cubic feet per second is balanced by an equivalent discharge. For the period 1976-80, simulated recharge increased to 1,398 cubic feet per second. Total discharge, including pumpage from the four bedrock aquifers, increased to 1,619 cubic feet per second. The difference in recharge and discharge during this period is from aquifer storage. Results of variable-density simulations indicate that the rate of ground-water movement is small in areas where ground water is highly mineralized. The rates and directions are controlled by the intrinsic permeability of the rock, freshwater head gradients, and gravitational force.

Explore related subjects

Professional Paper Report

Source-reported bounding extent: 37.54457732085582° to 46.34692761055676° latitude; -96.7236328125° to -85.78125° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

R.J. Mandle, A.L. Kontis. 1992. Simulation of regional ground-water flow in the Cambrian-Ordovician aquifer system in the northern Midwest, United States: in Regional aquifer-system analysis. https://doi.org/10.3133/pp1405c

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

KEEP EXPLORING

Related USGS reports

USGS2026

Allostratigraphy and formalizing the Cold Harbor Alloformation, Virginia Coastal Plain, Eastern Virginia

Introduction Distinct morphologic features in the Virginia part of the Coastal Plain Province (herein referred to as the Virginia Coastal Plain) were noted in one of the earliest geologic surveys in 1839 (Rogers, 1884). South of the James River, Rogers (1884, p. 252) described a western bench that ranged from “a low rocky coast…of the primeval sea” to an eastern facing scarp east of the City of Suffolk, Virginia. His eastern bench extended from the scarp at Suffolk, which he suggested was a former shoreline, to the present Atlantic Ocean. Notably, the eastern plain was described as less dissected than the western plain. In Virginia, creation and publication of topographic maps by the U. S. Geological Survey began in the mid-1880s. By 1905, many areas in eastern Virginia were mapped as 30-minute and 15-minute series with 10-foot (ft) contour interval accuracy (Virginia Division of Mineral Resources [VDMR] , 2008). After the continued release of these early maps that displayed Virginia Coastal Plain morphology in increasingly greater detail, several geological reports were published where scarps and plains were traced, named, and origins proposed (Shaler, 1890; Shattuck, 1901; Clark and Miller, 1906, 1912; Wentworth, 1930). Shattuck (1901) recognized three continuous scarps in Maryland and most importantly suggested that there was a correlation between surficial deposits (geologic formations) and morphologic terraces; this defined the terrace-formation concept that became the basis of geologic mapping of the Upper Pliocene and Pleistocene sediments in the Virginia Coastal Plain. Clark and Miller (1912) described three Pleistocene formations and supported the terrace-formation concept because the sediments within each formation were lithologically similar to one another, and differentiating them by comparing key vertical sections was not helpful. In Virginia, Wentworth (1930) defined and mapped five Pleistocene terraces on the basis of correlating morphologic terraces with their underlying formation sediments. Similarly, Cooke (1931) identified and traced seven marine shorelines and proposed corresponding terrace names from Maryland to Florida. Cooke (1931) also pointed out that the various elevations of scarps were related to interglacial sea-level highstands whereas prior workers called on vertical land movement of emergence and submergence to create shoreline scarps. Cooke’s (1936) geologic map of eastern South Carolina showing the Pleistocene sediments is also based on the terrace-formation concept. Detailed geologic mapping of the Virginia Coastal Plain by Virginia Division of Mineral Resources (VDMR) and U.S. Geological Survey (USGS) followed lithostratigraphic convention for Early Pliocene and older marine sediments. However, geologic mapping of Upper Pliocene and Pleistocene sediments used the concept of morphologic setting and elevations of scarp toe and flat-plain, or the terrace-formation concept. The 1983 version of the North American Stratigraphic Code (hereafter the Code) (North American Commission on Stratigraphic Nomenclature [NACSN] , 1983) introduced allostratigraphy that gave credence to the terrace-formation concept. According to the current Code (NACSN, 2021; Article 58), “an allostratigraphic unit is a mappable body of rock that is defined and identified on the basis of its bounding discontinuities,” and, “formal allostratigraphic units may be defined or distinguish between different (1) superposed discontinuity-bounded deposits of similar lithology, (2) contiguous discontinuity-bounded deposits of similar lithology, or (3) geographically separated discontinuity-bounded units of similar lithology…The lithology of an allostratigraphic unit plays no part in its definition.” However, definition and placement of bounding discontinuities relies on lateral and vertical knowledge of lithology and is so recognized in Article 58e, “A type locality and type area must be designated; a composite stratotype or a type section and several reference sections are desirable.” We understand that each of the terrace formations of Late Pliocene and Pleistocene age were deposited during a marine transgressive event tied to an interglacial sea-level highstand (scarp above the Poquoson Member of the Tabb Formation may be the result of a still-stand as sea level was falling after deposition of the Lynnhaven Member, refer to fig. 1). The base of each map unit is an unconformity (or is at least in part for the base of the Poquoson) and nomenclature for each unit might be changed from “formation” to “alloformation,” as suggested by Doar and Berquist (2009). Figure 1 shows the relative placement of these mapped and named alloformations as used in the Virginia Coastal Plain and Fall Zone (a zone of transition between the Piedmont and Coastal Plain Provinces) (Berquist, 2021). Except perhaps for the Poquoson, many of these proposed alloformations are unlikely to interfinger or grade laterally seaward into extensive deep marine sediments. After their deposition, erosion during sea-level fall and then rise for the next depositional event would seemingly obliterate any marine connection without a basal unconformity. All these alloformations contain sediments deposited in fluvial, estuarine, and shallow-marine environments and are bound everywhere by unconformities. This paper documents the Cold Harbor Alloformation as a formal stratigraphic unit and the rationale for mapping the sediments using an allostratigraphic classification. The Cold Harbor Alloformation is composed of sediments from a variety of environments such as gravel beds, shoreline beach sands, marsh and estuarine sand and mud, and shallow marine sands. These sediments are found at elevations up to 240 ft and lie unconformably above older sediments and rocks. Our work relied on observations of outcrops and many auger borings to identify basal unconformities.

Professional Paper↗

USGS2026

Formalization of the Elsing Green Alloformation, Virginia Coastal Plain, Eastern Virginia

The Elsing Green Alloformation is herein formalized as a new Pleistocene alloformation in the Virginia Coastal Plain. This allostratigraphic unit is characterized by an unconformity-bounded terrace deposit occupying a morphologic flat with a mapped elevation range of about 26 to 38 feet (7.9 to 11.6 meters); its landward extent is commonly marked at the base of an erosional scarp where the toe is typically at 38 feet (11.6 meters) and is distinguished by its topographic expression and stratigraphic position. Sediments of this unit were previously mapped as part of the Shirley Formation or Sedgefield Member of the Tabb Formation but differ significantly in their bounding scarps and terrace morphology. The type area is within the Tunstall 7.5-minute quadrangle near Elsing Green Plantation, King William County, Virginia. The reference section is provided by boring 2012270-2. Correlation is proposed with other mid-Pleistocene terrace deposits in the mid-Atlantic Coastal Plain, including the younger Lynch Heights Formation of Delaware and the Butlers Bluff Member of the Nassawadox Formation of the Eastern Shore of Virginia. The age of the Elsing Green Alloformation is interpreted to be Marine Isotope Stage 9e (about 330,000 years).

Professional Paper↗

USGS2026

Resurvey of the Marble Canyon and Bridge Canyon dam sites in Grand Canyon National Park—Changes in sediment storage and evidence supporting the occurrence of bedrock incision through the mid-20th century

The Bureau of Reclamation developed an extensive plan for a network of dams, water tunnels, and hydropower plants in and around Grand Canyon, Arizona, in the 1940s through 1960s. The two largest of these planned dams were the Marble Canyon and Bridge Canyon Dams on the Colorado River. Though these dams were ultimately never built, Reclamation conducted extensive topographic, bathymetric, and subsurface exploration work at the sites proposed for these dams in the 1940s and 1950s. Resurveys of these dam sites were conducted between 1998 and 2021 to determine the changes in sediment storage at these dam sites caused by the upstream construction and operation of Glen Canyon Dam and by the recession of Lake Mead, the reservoir impounded by Hoover Dam. The resurveys of the Marble Canyon dam sites indicate that the post-1950s changes in sediment storage at these dam sites are broadly consistent with flux-based estimates of voluminous sand erosion from Marble Canyon since the 1963 closure of Glen Canyon Dam. These resurveys also suggest that the pre-dam longitudinal variation in sediment thickness over bedrock played a key role in determining the locations of the sand erosion induced by Glen Canyon Dam; more sand eroded from locations where more sand was present in the 1950s. The resurvey of the Bridge Canyon dam sites indicates that the Colorado River’s incision of the Lake Mead delta is regulated both by bed-sediment grain size and downstream hydraulic controls. Finally, analyses of bed-sediment thickness and sedimentological data at the dam sites, and observations of bed scour and gravel transport, suggest that sufficient bedrock was exposed to allow bedrock incision during commonly recurring pre-dam snowmelt floods that entrained small boulders into transport.

Arizona↗