SEARCH · USGS Science

Results for “Northeastern Geology”

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

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55Linked to original sources

USGS2012

Water-level altitudes 2012 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973-2011 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

Most of the subsidence in the Houston–Galveston region, Texas, has occurred as a direct result of groundwater withdrawals for municipal supply, commercial and industrial use, and irrigation that depressured and dewatered the Chicot and Evangeline aquifers and caused compaction of the clay layers of the aquifer sediments. This report—prepared by the U.S. Geological Survey in cooperation with the Harris– Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District—is one in an annual series of reports depicting water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction in the Chicot and Evangeline aquifers in the Houston–Galveston region. The report contains maps showing approximate water-level altitudes for 2012 (calculated from measurements of water levels in wells made during December 2011–February 2012) for the Chicot, Evangeline, and Jasper aquifers; maps showing 1-year (2011–12) water-level-altitude changes for each aquifer; maps showing 5-year (2007–12) water-levelaltitude changes for each aquifer; maps showing long-term (1990–2012 and 1977–2012) water-level-altitude changes for the Chicot and Evangeline aquifers; a map showing long-term (2000–12) water-level-altitude change for the Jasper aquifer; a map showing locations of borehole extensometer sites; and graphs showing measured compaction of subsurface sediments at the extensometers from 1973 (or later) through 2011. Tables listing the data that were used to construct each water-level map for each aquifer and the cumulative compaction graphs are included. In 2012, water-level-altitude contours for the Chicot aquifer ranged from 250 feet (ft) below North American Vertical Datum of 1988 (hereinafter, datum) in a small area in southwestern Harris County to 200 ft above datum in westerncentral Montgomery County. Water-level-altitude changes during 2012 in the Chicot aquifer ranged from a 48-ft decline to an 18-ft rise. Contoured 5-year and long-term changes in water-level altitudes in the Chicot aquifer ranged from a 60-ft decline to a 40-ft rise (2007–12), from a 100-ft decline to an 80-ft rise (1990–2012), and from a 100-ft decline to a 200-ft rise (1977–2012). In 2012, water-level-altitude contours for the Evangeline aquifer ranged from 300 ft below datum in isolated areas located in south-central Montgomery County and north-central Harris County, in southwest Harris County, and in northeastern Fort Bend County to 200 ft above datum near the county boundary intersection of Waller, Montgomery, and Grimes Counties. Water-level-altitude changes for 2012 in the Evangeline aquifer ranged from a 90-ft decline to a 39-ft rise. Contoured 5-year and long-term changes in waterlevel altitudes in the Evangeline aquifer ranged from an 80-ft decline to an 80-ft rise (2007–12), from a 220-ft decline to a 220-ft rise (1990–2012), and from a 360-ft decline to a 260-ft rise (1977–2012). In 2012, water-level-altitude contours for the Jasper aquifer ranged from 250 ft below datum in south-central Montgomery County to 250 ft above datum in northwest Montgomery County. Water-level-altitude changes for 2012 in the Jasper aquifer ranged from a 74-ft decline to a 4-ft rise. Contoured changes in water-level altitudes in the Jasper aquifer ranged from a 120-ft decline to no change (2007–12), and from a 220-ft decline to no change (2000–12). Compaction of subsurface sediments (mostly in the clay layers) composing the Chicot and Evangeline aquifers was recorded continuously at 13 borehole extensometers at 11 sites. For the period of record beginning in 1973 (or later) and ending in December 2011, cumulative compaction data collected from the 13 extensometers ranged from 0.102 ft at the Texas City–Moses Lake site to 3.621 ft at the Addicks site. The rate of compaction varies from site to site because of differences in groundwater withdrawals near each site and differences among sites in the clay-to-sand ratio in the subsurface sediments. Therefore, it is not possible to extrapolate or infer a rate of compaction for adjacent areas on the basis of the rate of compaction measured at a nearby extensometer.

Texas↗

USGS2016

Submarine glacial landforms on the Bay of Fundy–northern Gulf of Maine continental shelf

The Bay of Fundy–northern Gulf of Maine region surrounds the southern part of Nova Scotia, encompassing, from west to east, the Bay of Fundy, Grand Manan Basin, German Bank, Browns Bank, Northeast Channel and northeastern Georges Bank (Fig. 1a, b). During the last glacial maximum ( c. 24–20 14 C ka BP), the SE margin of the Laurentide Ice Sheet (LIS) occupied the study area, the rest of the Gulf of Maine and the continental Scotian Shelf off Atlantic Canada (see Dyke et al. 2002 , fig. 1; Shaw et al. 2006 , fig. 8; Hundert & Piper 2008 , fig. 16). Early mapping of the glaciated region on the Scotian Shelf using side-scan sonar imagery and seismic-reflection profiles revealed topographic features interpreted to be recessional moraines indicative of retreat of the LIS ( King et al. 1972 ; King 1996 ). Subsequently, multibeam sonar seafloor mapping of local-scale glacial landforms on the inner Scotian Shelf off Halifax, Nova Scotia (Fig. 1b) provided further information on the dynamics of the advance and retreat of the ice sheet ( Loncarevic et al. 1994 ). Interpretation of seismic-reflection profiles across Georges Bank revealed that the surficial sediment is a veneer of glacial debris transported to Georges Bank by the LIS during the late Pleistocene from continental areas to the north ( Shepard et al. 1934 ; Knott & Hoskins 1968 ; Schlee 1973 ; Twichell et al. 1987 ; Fader et al. 1988 ). Recent high-resolution multibeam sonar surveys of German Bank and the Bay of Fundy mapped a complex of ice-advance and ice-retreat features attributed to the activity of the LIS ( Todd et al. 2007 ; Todd & Shaw 2012 ).

Geological Society, London, Memoirs↗

USGS1997

Early diagenetic partial oxidation of organic matter and sulfides in the Middle Pennsylvanian (Desmoinesian) Excello Shale Member of the Fort Scott Limestone and equivalents, northern Midcontinent region, USA

A process of early diagenetic partial oxidation of organic matter and sulfides has altered the chemical composition of the Middle Pennsylvanian Excello Shale Member of the Fort Scott Limestone and equivalents in the northern Midcontinent region. This process was identified by comparison of organic carbon contents, Rock-Eval hydrogen indices, organic carbon ??13C and element compositions of core and surface mine samples of the Excello Shale Member with analyses of three other underlying and overlying organic-matter-rich marine shales (offshore shale lithofacies) from southern Iowa, northern Missouri, eastern Kansas and northeastern Oklahoma. The end product of the partial oxidation process is shale with relatively low contents of hydrogen-poor, C13-enriched organic matter, lower contents of sulfur and sulfide-forming elements, and relatively unchanged contents of phosphorus and many trace elements (e.g. Cr, Ni, and V). However, because of lower organic carbon contents, element/organic carbon ratios are greatly increased. The partial oxidation process apparently took place during subaerial exposure of the overlying marine carbonate member (Blackjack Creek Member of the Fort Scott Limestone) following a marine regression when meteoric waters percolated down to the level of the Excello muds allowing oxidation of organic matter and sulfides. This hypothesis is supported by earlier workers, who have identified meteoric carbonate cements within, and soil horizons at the top of the Blackjack Creek Member. The period of oxidation is constrained in that organic matter and sulfides in the Little Osage Shale Member of the Fort Scott Limestone and equivalents (immediately overlying the Blackjack Creek Member) appear unaltered. Similar alteration of other shales in the Middle and Upper Pennsylvanian sections may be local to regional in extent and would depend on the extent and duration of the marine regression and be influenced by local variations in permeability and topography. The partial oxidation process has likely led to a redistribution of sulfur and sulfide-forming elements into other organic-rich lithologies in the section. The altered/oxidized shales are nongenerative with respect to hydrocarbon generation.

Chemical Geology↗

USGS2000

Geochronology and geology of late Oligocene through Miocene volcanism and mineralization in the western San Juan Mountains, Colorado

Twenty-five new 40Ar/39Ar ages from volcanic rocks and veins in the western San Juan Mountains clarify relationships between volcanism and mineralization in this classic area. Five calc-alkaline ash-flow sheets erupted from caldera sources (Ute Ridge, Blue Mesa, Dillon Mesa, Sapinero Mesa, and Crystal Lake Tuffs) from 28.6 to 27.6 Ma. This is a much more restricted time interval than previously thought and indicates that the underlying batholith rose and evolved very rapidly beneath the western San Juan Mountains. The new ages and geologic relations constrain the timing of joint resurgence of the Uncompahgre and San Juan calderas to between 28.2 and 27.6 Ma. The collapse of the Silverton caldera produced a set of strong ring fractures that intersected with graben faults on the earlier resurgent dome to produce the complex set of structures that localized the mid-Miocene epithermal gold veins. Later calc-alkaline monzonitic to quartz monzontic plutons solidified at 26.5-26.0 Ma as the underlying batholith rose through its volcanic cover. A new age from lavas near Uncompahgre Peak supports earlier interpretations that these lavas were fed by nearby 26 Ma monzonite intrusions. Nearly all of these intrusions are associated with subeconomic Mo and Cu mineralization and associated alteration, and new ages of 26.40 and 25.29 Ma from the Ute-Ulay and Lilly veins in the Lake City region show that some of the most important silver and base-metal veins were temporally and possibly genetically connected to these plutons. In addition, the Golden Fleece telluride vein cuts all of the post-Uncompahgre caldera volcanics in the area and is probably temporally related to this cycle, though its age of 27.5 ? 0.3 Ma was determined by less precise U/Pb methods. The 22.9 Ma Lake City caldera collapsed within the older Uncompahgre caldera structure but is petrologically unrelated to the older calc-alkaline activity. The distinctive suite of high-silica rhyolite tuff and alkaline resurgent intrusions indicates that it is closely related to the early stages of bimodal high-silica rhyolite-alkali basalt volcanism that accompanied the onset of extensional tectonism in the region. Both 40Ar/39Ar ages and paleomagnetic data confirm that the entire caldera sequence formed in less than 330,000 years. Only weak quartz vein mineralization is present in the center of the caldera, and it appears to be related to leaching of metals from the intracaldera tuffs above the resurgent intrusion. Massive alunitization and weak Mo and Cu mineralization along the eastern ring fracture are associated with calc-alkaline lavas and stocks related to late stages of the caldera cycle. These calc-alkaline stocks also appear to be genetically and temporally linked to a radial pattern of barite-precious metal veins on the northeastern margin of the Lake City caldera.

Colorado↗

USGS1941

Geologic structure and occurrence of gas in part of southwestern New York. Part 1, Structure and gas possibilities of the Oriskany sandstone in Steuben, Yates, and parts of the adjacent counties

The area covered by this report is in southwestern New York and includes a little more than 3,000 square miles in Steuben and Yates counties and parts of the six adjacent counties. This area has been mapped to determine the structural attitude of the exposed rocks, so as to aid those interested in prospecting for natural gas in the Oriskany sandstone of Lower Devonian age. Because of the gentle regional dip toward the southwest, the youngest beds are exposed in the southwest corner of the area, and progressively older beds crop out northeastward in successive bands that strike generally northwest. All the exposed rocks are of Upper Devonian age except those in a narrow belt at the extreme north edge of the area, where a small thickness of Middle Devonian rocks crops out. The maximum thickness of beds so exposed is nearly 4,000 feet, of which the lower part is predominantly soft dark shale and the upper part predominantly fine-grained sandstone and gray shale. All the beds are marine except a few tongues of continental deposits red shale and sandstone and gray mudstone in the youngest beds. All the beds thicken southeastward, so that there is a northwestward convergence between any two lithologic units in the series. More than 30 key horizons that are persistent and distinctive were mapped, and altitudes on these key horizons served as a basis for constructing the structure contour map. Many of the key horizons are formation or member boundaries, but others are the tops or bottoms of limestone or sandstone beds within formations. All the stratigraphic units mapped are purely lithologic. (See pi. 2.) The Tully limestone, which crops out along the northern border of the area, is an easily recognizable and therefore valuable key bed for subsurface correlations in this part of the State. Below the Tully limestone is a thick body of Middle Devonian shales of the Hamilton group which rests on another valuable key bed, the hard, cherty Onondaga limestone, also of Middle Devonian age. Below the Onondaga limestone is the Lower Devonian Oriskany sandstone, which is the gas-producing bed. Unlike the Onondaga, the Oriskany is locally thin or absent. The structure of most of the area is shown by contour lines at 25-foot intervals, but where key horizons are lacking the structure is indicated by dip symbols. Upon the regional south and southwest dip are superposed numerous gentle folds whose axes trend approximately northeastward in the greater part of the area but more nearly eastward in the eastern part. The folds generally tend to become narrower and steeper, and therefore more closely spaced, southwestward. Many of the anticlines fork southwestward, whereas the synclines. tend to fork northeastward. All the folds have a westward or southwestward plunge. Throughout the area the rocks are jointed in two dominant sets one that trends northwest and the other east or northeast. No evident relation .between these joints, which were measured only in the hard, relatively brittle beds, and the individual folds or domes was discernible. The faults are concentrated in the northeastern and southwestern parts of the area and trend either northeastward or northwestward. Some are nearly vertical normal faults; others are steep reverse faults. Subsurface data show that most of the faults increase in throw downward and also that many subsurface faults do not reach the surface. A group of faults in the northwestern part of the Greenwood quadrangle and the southwestern part of the Hornell quadrangle were active during Upper Devonian time, while the Gowanda shale and overlying beds were being deposited. At this stratigraphic horizon the beds in a zone a few hundred feet thick are highly deformed in a wide belt on both sides of the faults. Sandstone layers are thinned out into long stringers or swollen into thick masses and in places are bent acutely without fracture. Thin layers of shale, coquina, and sand have flowed together into intricately plicated zones that lack cleavage and joints. These features show that the sediments were deformed while wet and plastic and buried only a little way below the sea floor. The beds that were laid down over these disturbed zones were not involved in this deformation. Many of the sharper flexures and most of the faults are not evident in the beds several hundred feet stratigraphically higher. Accordingly, broad, gentle folds in these higher beds in parts of the area south and west of the northwest corner of the Greenwood quadrangle may conceal, at considerable depths below them, narrow folds separated by abrupt flexures or faults. Several of the larger streams and rivers occupy strike valleys, and their j courses swing to follow the changing strike of the rocks where they cross ( successive folds. But, with few exceptions, the small streams are not adjusted to the bedrock structure. Domes likely to serve as traps for natural gas are concentrated in the northeastern and southwestern parts of the area. The Wayne-Dundee gas field is in the northeastern part. All the other potentially valuable domes in this part of the area have been drilled and found valueless except one small structural feature in the southern part of the Ovid quadrangle, which, if the Oriskany is present, may trap a small quantity of gas. In the Greenwood quadrangle in the southwestern part of the area there is one gas field and four well-defined domes, all of which may be productive if the Oriskany sandstone is present. In the northwest corner of the quadrangle the dips indicate at least two domes that can be adequately defined and evaluated only by geophysical prospecting. The State Line gas field is in tbe Wellsville quadrangle. In the southeast corner of this quadrangle there are three other domes of comparable size that may also be productive if underlain by the Oriskany sandstone. At other places in the Wellsville quadrangle the dips suggest several anticlinal axes on which analogous productive domes maybe found. The structural features in this quadrangle, however, are defined by contours only in the southeastern part. In the Woodhull quadrangle a large dome east of Jasper may be productive, and the western top of the large Woodhull dome in the southwestern part of the quadrangle seems to warrant drilling, despite the absence of the Oriskany in a well on the eastern top. Two wells drilled in 1936 and 1937 a little northeast of a broad, nearly flat-topped dome in the Hornell quadrangle, a few miles east of Hornell,, struck small flows of gas, suggesting that wells drilled higher on this dome may be productive. In much of the southwestern part of the area seismograph surveys should be of great value in determining the structure at the Tully and Onondaga horizons. Without abundant subsurface control of this sort, the danger of drilling into subsurface faults can hardly be overemphasized. Three closed or nearly closed synclines in the Greenwood and Wellsville quadrangles appear to be favorable places to drill for oil in the shallow sands presumably parts of the Dunkirk sandstone.

New York↗

USGS2004

Petroleum geology and resources of the Amu-Darya basin, Turkmenistan, Uzbekistan, Afghanistan, and Iran

The Amu-Darya basin is a highly productive petroleum province in Turkmenistan and Uzbekistan (former Soviet Union), extending southwestward into Iran and southeastward into Afghanistan. The basin underlies deserts and semideserts north of the high ridges of the Kopet-Dag and Bande-Turkestan Mountains. On the northwest, the basin boundary crosses the crest of the Karakum regional structural high, and on the north the basin is bounded by the shallow basement of the Kyzylkum high. On the east, the Amu-Darya basin is separated by the buried southeast spur of the Gissar Range from the Afghan-Tajik basin, which is deformed into a series of north-south-trending synclinoria and anticlinoria. The separation of the two basins occurred during the Neogene Alpine orogeny; earlier, they were parts of a single sedimentary province. The basement of the Amu-Darya basin is a Hercynian accreted terrane composed of deformed and commonly metamorphosed Paleozoic rocks. These rocks are overlain by rift grabens filled with Upper Permian-Triassic rocks that are strongly compacted and diagenetically altered. This taphrogenic sequence, also considered to be a part of the economic basement, is overlain by thick Lower to Middle Jurassic, largely continental, coal-bearing rocks. The overlying Callovian-Oxfordian rocks are primarily carbonates. A deep-water basin surrounded by shallow shelves with reefs along their margins was formed during this time and reached its maximum topographic expression in the late Oxfordian. In Kimmeridgian-Tithonian time, the basin was filled with thick evaporites of the Gaurdak Formation. The Cretaceous-Paleogene sequence is composed chiefly of marine clastic rocks with carbonate intervals prominent in the Valanginian, Barremian, Maastrichtian, and Paleocene stratigraphic units. In Neogene time, the Alpine orogeny on the basin periphery resulted in deposition of continental clastics, initiation of new and rejuvenation of old faults, and formation of most structural traps. A single total petroleum system is identified in the Amu-Darya basin. The system is primarily gas prone. Discovered gas reserves are listed by Petroconsultants (1996) at about 230 trillion cubic feet, but recent discoveries and recent reserve estimates in older fields should increase this number by 40 to 50 trillion cubic feet. Reserves of liquid hydrocarbons (oil and condensate) are comparatively small, less than 2 billion barrels. Most of the gas reserves are concentrated in two stratigraphic intervals, Upper Jurassic carbonates and Neocomian clastics, each of which contains about one-half of the reserves. Reserves of other stratigraphic units?from Middle Jurassic to Paleogene in age?are relatively small. Source rocks for the gas are the Lower to Middle Jurassic clastics and coal and Oxfordian basinal black shales in the east-central part of the basin. The latter is probably responsible for the oil legs and much of the condensate in gas pools. Throughout most of the basin both source-rock units are presently in the gas-window zone. Traps are structural, paleogeomorphic, and stratigraphic, as well as a combination of these types. The giant Dauletabad field is in a combination trap with an essential hydrodynamic component. Four assessment units were identified in the total petroleum system. One unit in the northeastern, northern, and northwestern marginal areas of the basin and another in the southern marginal area are characterized by wide vertical distribution of hydrocarbon pools in Middle Jurassic to Paleocene rocks and the absence of the salt of the Gaurdak Formation. The other two assessment units are stratigraphically stacked; they occupy the central area of the basin and are separated by the regional undeformed salt seal of the Gaurdak Formation. The largest part of undiscovered hydrocarbon resources of the Amu-Darya basin is expected in older of these assessment units. The mean value of total assessed resources of the Amu-Darya basin is estimated

Bulletin↗

USGS1982

Submarine valleys in the northeastern Gulf of Alaska: Characteristics and probable origin

The continental shelf of the northeastern Gulf of Alaska Between Prince William Sound and Cross Sound is cut by at least eight major valleys. From west to east, these are Hinchinbrook Seavalley, Egg Island Trough, Kayak Trough, Bering Trough, Pamplona Troughs, Yakutat Valley, Alsek Valley and Yakobi Valley. Evidence common to most of these troughs or valleys indicating that the present morphology is due to glacial processes includes: (1) a pre-Holocene subbottom erosional surface incised into the underlying lithified strata of the shelf; (2) U-shaped cross sections, both at the sea floor and at the pre-Holocene erosional surface; (3) concave longitudinal sections, commonly shoaling at the seaward end; (4) till-like sediments collected from the walls or outer shelf adjacent to the troughs; and (5) seismic stratigraphy that can be correlated with bottom samples indicative of glacially derived strata. Depressions with tens of meters of relief are present on the pre-Holocene subbottom erosional surface beneath most of these valleys. These depressions have been partially filled by a seaward-thinning wedge of Holocene glacial flour (clayey silt) that is filling the valleys and blanketing the inner shelf at rates as high as 15 mm/yr (based on 210 Pb measurements). Although glaciation played a dominant role in the modern morphology of these sea valleys, structural features, including structurally controlled topographic highs on the shelf (e.g. Tarr Bank, Kayak Island, Pamplona Spur and Fairweather Ground) influenced the flow directions of the glacial lobes.

Alaska↗

USGS2004

Stable isotopic constraints on fluid-rock interaction and Cu-PGE-S redistribution in the Sonju Lake intrusion, Minnesota

The Sonju Lake intrusion, part of the 1.1 Ga Midcontinent rift-related Beaver Bay Complex, is a 1,200-m-thick, strongly differentiated, layered sequence of mafic cumulates located in northeastern Minnesota. Basal melatroctolite and dunite layers are overlain by troctolite, gabbro, Fe-Ti oxide-rich gabbro, apatite diorite, and monzodiorite. Stratigraphic intervals rich in Pt + Pd, Cu, and S occur over ???500 m in the Fe-Ti oxide-rich gabbro and apatite diorite units. Peak concentrations show offsets that are similar to those found in other tholeiitic layered intrusions. Concentrations of Pd in excess of 100 ppb are confined to the lowermost 25 m of the interval. Copper shows a sharp increase to 630 ppm above the Pd-rich interval. Sulfur contents are low (<375 ppm) in the Cu-rich interval, but they increase to values as high as 3,150 ppm above in the apatite diorite. Disseminated sulfides in the intrusion have ??34S values that range from -2.2 to 3 per mil Vienna-Canyon Diablo Troilite (V-CDT) and suggest that contamination by country rock sulfur was not an important process in the formation of the metal-rich interval. ??18O values of plagioclase from the intrusion range from 5.6 to 12.0 per mil (V-SMOW) and indicate that a relatively low-18O fluid (??18O ???3-5 ???) interacted with the rocks of the intrusion at temperatures less than ???275??C. Clinopyroxene and Fe-Ti oxides (ilmenite with minor amounts of titanomagnetite) show much more restricted ranges in ??18O values (4.6-5.7 and 5.5-6.7 per mil, respectively) and attest to the kinetic control of the oxygen isotope exchange process. The externally derived fluid that interacted with rocks now enriched in platinum group elements (PGE) + Cu- and Fe-sulfide minerals locally liberated sulfur and replaced chalcopyrite and pyrite with goethite. In the Cu-rich zone, goethite that replaces chalcopyrite may contain up to 8.5 weight percent Cu. It is evident that hydrothermal alteration resulted in a decoupling of copper and sulfur, with sulfur being transferred out of the Cu-rich interval. Interaction between rocks in the PGE-Cu-S interval of the Sonju Lake intrusion and an externally derived fluid at low temperatures modified what appears to have been a primary stratigraphic metal-sulfur zonation. The effects of hydrothermal alteration on PGE and base-metal sulfide mobility and redistribution must be understood before models of primary zonation processes can be meaningfully applied. ?? 2004 by Economic Geology.

Economic Geology↗

USGS2010

Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model area in Michigan and Indiana. Overall water use for the selected areas in Wisconsin, Michigan, Indiana, and Illinois was less during early time intervals than during more recent intervals, with large increases beginning around the 1960s. Total estimated groundwater withdrawals for model input range from 18.01 million gallons per day (Mgal/d) for interval 1 (1864-1900) to 1,280.25 Mgal/d for interval 12 (2001-5). Withdrawals for the public-supply category make up the majority of the withdrawals in each of the four states. In Wisconsin and Michigan, the second largest withdrawals are for the irrigation category; in Indiana and Illinois, industrial withdrawals account for the second largest withdrawal amounts. The smallest withdrawals are for miscellaneous uses in Wisconsin and irrigation uses in Indiana and Illinois. Estimated groundwater withdrawals in the Southern Lower Peninsula of Michigan, Northeastern Illinois, and the farfield model area are generally larger than in the other model subregions. Withdrawals in Michigan and Indiana are predominantly from the Quaternary aquifer system, whereas withdrawals in Illinois are predominantly from the Cambrian-Ordovician aquifer systems. Withdrawals in Wisconsin are about equal from the Quaternary and Cambrian-Ordovician aquifer systems. Estimated groundwater withdrawals in Michigan and Indiana are predominantly from the unconfined unconsolidated aquifer type. Withdrawals in Illinois are largely from the deep confined bedrock aquifer type, although they decreased considerably in more recent time intervals. Wisconsin withdrawals are about equal from unconfined unconsolidated and deep confined bedrock aquifer types. Groundwater-withdrawal estimates in Wisconsin were compiled for the 47 easternmost counties within the boundary of the Lake Michigan Basin model, of which 32 counties, though not entirely contained, are at least partly within the Lake Michigan Basin. Overall, 6,457 withdrawal locations were estima

Scientific Investigations Report↗

USGS2020

Interpretation of hydrogeologic data to support groundwater management, Bazile Groundwater Management Area, northeast Nebraska, 2019—A case demonstration of the Nebraska Geocloud

Nitrate, age tracer, and continuous groundwater-level data were interpreted in conjunction with airborne electromagnetic (AEM) survey data to understand the movement of nitrate within the Bazile Groundwater Management Area (BGMA) in northeastern Nebraska. Previously published age tracer data and nitrate data indicated vertical stratification of groundwater quality. Younger groundwater sampled within shallow parts of the aquifer had higher concentrations of nitrate, with 70 percent exceeding the U.S. Environmental Protection Agency maximum contaminant level of 10 milligrams per liter. In contrast, groundwater sampled from deeper parts of the aquifer indicated that nitrate concentrations were less than 2 milligrams per liter and that groundwater likely recharged prior to widespread use of commercial fertilizer. The hydrostratigraphic interpretation of AEM profiles indicated that shallow and deep monitoring wells were often screened within the same homogenous zone of aquifer material. In contrast, test-hole logs indicated that there often are fine-grained layers within these homogenous zones that separate the shallow and deep monitoring well screens, but these fine-grained layers are not detected by the AEM technique because of decreased resolution of the AEM technique with depth. The stratification of groundwater ages and nitrate concentrations likely was caused by groundwater-flow paths of different length, location and time of recharge, and denitrification. Within paleochannels interpreted from AEM and test-hole data, pesticides detected in groundwater generally coincide with elevated nitrate concentrations. Continuous groundwater-level data from four monitoring well nests indicated that groundwater pumping can impose or increase downward hydraulic gradients and facilitate the downward movement of nitrate into deeper parts of the High Plains aquifer. Given the density of irrigation wells within the BGMA, this effect on the hydraulic gradient is likely prevalent in other areas of the BGMA. Understanding seasonal water-level changes can allow water managers to better predict and assess the hydraulic gradient and the vulnerability of groundwater in deeper parts of the High Plains aquifer. Nitrate, age tracer, and continuous groundwater-level data within the BGMA were interpreted in conjunction with AEM data as a case demonstration of the Nebraska Geocloud. The Nebraska Geocloud was initiated to protect taxpayer investments in AEM data collection and realize maximum benefit of these data by creating a publicly available, online digital database for long-term data storage. The Lower Platte North, Lower Platte South, Papio-Missouri River, Nemaha, Lower Loup, Central Platte, Upper Elkhorn, Lower Elkhorn, Lower Niobrara, and Lewis and Clark Natural Resources Districts; the University of Nebraska-Lincoln Conservation and Survey Division, Nebraska Natural Resources Commission, Nebraska Department of Natural Resources; and the U.S. Geological Survey entered a cooperative agreement to begin a program of data management and research aimed at understanding the best use of AEM for groundwater sustainability and management. Resulting case-study interpretations are provided to guide use of the Nebraska Geocloud to assess water-quality conditions and can be used by water managers and staff to address applicable water resource problems.

Nebraska↗

USGS1978

Petrographic differentiation of depositional environments of sandstones of the Pennsylvanian Breathitt Formation, northeastern Kentucky and southwestern West Virginia

Petrographic properties of sandstone samples from the lower and middle Pennsylvanian Breathitt Formation in northeastern Kentucky and southwestern West Virginia were utilized to classify their environments of deposition. Simultaneous consideration of mineral composition, grain size, and sorting of the sandstones by using stepwise multiple discriminant-function analysis has permitted differentiation between alluvial-plain and delta-plain sandstones. Stepwise multiple discriminant-function analysis of the sandstone samples- was performed in two parts. The first analysis assigned the sandstones of all grain-size classes to one of three environments of deposition: (1) alluvial plain, (2) upper delta plain, and (3) lower delta plain. The second analysis classed the sandstones by grain size and subdivided each grain-size class into these environmental groups. Discriminant analyses of sandstones of all grain-size classes led to correct assignment of depositional environment to 76 percent of the samples. Analyses of the sandstones classed by grain size led to correct classification of 67 percent of the medium-grained sandstones, 83 percent of the fine-grained sandstones, and 91 percent of the very fine grained sandstones. The discriminantfunction method was adopted to classify the depositional environment and stratigraphic position of test sandstone samples. Analyses of fine-grained and very fine grained sandstones resulted in correct classification of 84 percent of the test samples to their depositional setting and stratigraphic position. The best discrimination is offered by analyses of fine-grained sandstones, which resulted in correct classification of 88 percent of the samples. Thus, this statistical method can be adopted as a powerful tool for exploration programs that concern delineation of genetically similar sandstones.

Kentucky, Virginia↗

USGS1987

The geochemistry of water near a surficial organic-rich uranium deposit, northeastern Washington State, U.S.A.

The chemistry of three stream, three spring and six near-surface waters in the vicinity of a Holocene organic-rich uranium deposit is described, with particular emphasis on the chemistry of U. Results characterize the solution behavior of uranium as U-bearing water interacts with relatively undecomposed, surficial organic matter. Of the measured major and trace chemical species, only U is consistently highly enriched (17-318 ppb) relative to reported values for regional waters, or to literature values for waters in largely granitic terrains. R-mode factor analysis of the chemical data suggests that most U is present in a soluble form, but that some U is also associated with fine suspended particulates of clay, organic matter, or hydrous oxides. Calculations that apply thermodynamic data to predict U speciation in solution indicate the relative importance of uranyl carbonate and uranyl phosphate complexes. Analysis of more finely filtered samples (0.05 ??m vs. 0.45 ??m), and direct radiographic observations using fission-track detectors suspended in the waters indicate the presence of some uraniferous particulate matter. Application of existing thermodynamic data for uranous- and uranyl-bearing minerals indicates that all waters are undersaturated with U minerals as long as ambient Eh ??? +0.1 v. If coexisting surface and near-surface waters are sufficiently oxidizing, initial fixation of U in the deposit should be by a mechanism of adsorption. Alternatively, more reducing conditions may prevail in deeper pore waters of the organic-rich host sediments, perhaps leading to direct precipitation or diagenetic formation of U4+ minerals. A 234U 238U alpha activity ratio of 1.08 ?? 0.02 in a spring issuing from a hillslope above the deposit suggests a relatively soluble source of U. In contrast, higher activity ratios of 234U 238U (??? 1.3) in waters in contact with the uraniferous valley-fill sediments suggest differences in the nature of interaction between groundwater and the local, U-rich source rocks. ?? 1987.

Chemical Geology↗

USGS2023

Roles of regional structures and country-rock facies in defining mineral belts in central Idaho mineral province with detail for Yellow Pine and Thunder Mountain mining districts

The central Idaho metallogenic province hosts numerous mineral deposit types. These include Late Cretaceous precious-polymetallic vein deposits, amagmatic Paleocene–Eocene breccia-hosted gold-tungsten-antimony deposits, and Eocene mercury deposits in metasedimentary roof pendants and in Late Cretaceous granitoids. Hot-springs gold deposits in Eocene volcanic rocks are also included in the central Idaho province. New sensitive high mass-resolution ion microprobe (SHRIMP) uranium-lead (U-Pb) ages for igneous rocks and for detrital zircon analyses of metasedimentary rocks along with geologic mapping clarify the geologic framework of the mineral deposits. This framework includes (1) structural controls for regional distribution of mining districts, (2) progressive structural development of individual districts, (3) regional sedimentary facies and their control of metals associations resulting in regional belts, and (4) influences of the several regional magmatic events. In central Idaho, 15 mining districts form two clusters that are grouped about a 200-kilometer (km) long system of normal faults. The northwestern cluster is in the regional hanging wall west of large, west-side-down faults, and the mineral deposits are located along smaller faults and fractures that cut the regional hanging wall. The southeastern cluster is in the regional hanging wall east of a linked large east-side-down fault and along and controlled by related hanging wall faults. At the southern extent of the regional fault system, the Yellow Pine-Thunder Mountain districts span a nearly 24-km-wide, east-tilted crustal block of normal-fault dominoes, exposing original crustal depths from 5 to 10 km deep on the west in the Late Cretaceous to shallow-surface depths on the east in the Eocene. Ore deposition in the northwestern district cluster was primarily Late Cretaceous and related to Idaho batholith plutons with only a single deposit related to a small Eocene intrusion; in the southeastern cluster, most deposits were initiated in the Late Cretaceous but with varying manifestations of overprinted Eocene mineralization activity. In the Yellow Pine-Thunder Mountain districts at the southern extent of the southern cluster, several mineralizing pulses occurred during hanging-wall collapse, such that (1) early deposits were multiply overprinted and (2) deposit depths, ages, and structural characteristics change progressively eastward. Originally deep-seated western Yellow Pine district deposits are Late Cretaceous viscoplastic mesothermal veins overprinted by Paleocene and Eocene breccia-hosted epithermal deposits. Central Yellow Pine district deposits contain early deeper vein systems but are primarily Paleocene and Eocene breccia-hosted epithermal deposits in Late Cretaceous plutonic rocks and Proterozoic–Paleozoic roof pendant rocks. Eastern district deposits are Eocene hot-springs-related deposits in the roof pendant. Thunder Mountain deposits farthest east are near-surface hot-springs deposits in Eocene volcanic and volcaniclastic rocks that overlie buried Cretaceous igneous and older roof pendant rocks. The mining district clusters are sited across several northwest-striking paleostratigraphic belts that are exposed in roof pendants and are offset by the regional normal fault system. A northeastern belt is Mesoproterozoic strata associated with gold-silver-copper±cobalt deposits. A central belt of Neoproterozoic rocks is not associated with mineral deposits in the central Idaho mineral province. A southwestern belt composed of probable Paleozoic deep-water miogeoclinal slope rocks and late Paleozoic epicratonic basinal rocks is thin and narrowly exposed but associated with gold-silver-antimony-tungsten±mercury deposits. These metasedimentary rocks (and their metal associations) are parts of regional mineral belts in which metal endowments are related to particular sedimentary facies belts and their Cretaceous thrust-fault juxtaposition and where these features have proximity to Late Cretaceous or Eocene igneous rocks. Offset and preservation or erosional stripping of these facies belts, thrust plates, igneous settings, and the associated regional mineral belts were controlled by the sense and magnitude of displacements across the regional normal-fault system.

Idaho↗

USGS1986

Pb-, Sr- and Nd-Isotopic systematics and chemical characteristics of cenozoic basalts, Eastern China

Forty-eight Paleogene, Neogene and Quaternary basaltic rocks from northeastern and east-central China have been analyzed for major-element composition, selected trace-element contents, and Pb, Sr and Nd isotopic systematics. The study area lies entirely within the marginal Pacific tectonic domain. Proceeding east to west from the continental margin to the interior, the basalts reveal an isotopic transition in mantle source material and/or degree of crustal interaction. In the east, many of the rocks are found to merge both chemically and isotopically with those previously reported from the Japanese and Taiwan island-arc terrains. In the west, clear evidence exists for component(s) of Late Archean continental lithosphere to be present in some samples. A major crustal structure, the Tan-Lu fault, marks the approximate boundary between continental margin and interior isotopic behaviors. Although the isotopic signature of the western basalts has characteristics of lower-crustal contamination, a subcrustal lithosphere, i.e. an attached mantle keel, is probably more likely to be the major contributor of their continental “flavor”. The transition from continental margin to interior is very pronounced for Pb isotopes, although Sr and Nd isotopes also combine to yield correlated patterns that deviate strikingly from the mid-ocean ridge basalt (MORB) and oceanic-island trends. The most distinctive chemical attribute of this continental lithosphere component is its diminished U/Pb as reflected in the Pb isotopic composition when compared to sources of MORB, oceanic-island and island-arc volcanic rocks. Somewhat diminished Sm/Nd and elevated Rb/Sr, especially in comparison to the depleted asthenospheric mantle, are also apparent from the Nd- and Sr-isotopic ratios.

Chemical Geology↗

USGS1975

Mercury distribution in ancient and modern sediment of northeastern Bering Sea

Reconnaissance sampling of surface and subsurface sediment to a maximum depth of 80 m below the sea floor shows that typical values of 0.03 p.p.m. and anomalies of 0.2-1.3 p.p.m. mercury have been present in northeastern Bering Sea since Early Pliocene time. Values are highest in modern beach (maximum 1.3 and mean 0.22 p.p.m. Hg) and nearshore subsurface gravels (maximum 0.6 and mean 0.06 p.p.m. Hg) along the highly mineralized Seward Peninsula and in clayey silt rich in organic matter (maximum 0.16 and mean 0.10 p.p.m. Hg) throughout the region. Although gold mining may be partly responsible for high mercury levels in the modern beach near Nome, Alaska (maximum 0.45 p.p.m.), equally high or greater concentrations of mercury occur in buried Pleistocene sediments immediately offshore (maximum 0.6 p.p.m.) and in modern unpolluted beach sediments at Bluff (maximum 1.3 p.p.m.); this suggests that the contamination effects of mining may be no greater than natural concentration processes in the Seward Peninsula region. The mercury content of offshore surface sediment, even adjacent to mercury-rich beaches, corresponds to that of unpolluted marine and fresh-water sediment elsewhere. The normal values that prevail offshore may be attributable to entrapment of mercury-bearing heavy minerals on beaches near sources and/or dilution effects of offshore sedimentation. The few minor anomalies offshore occur in glacial drift derived from mercury source regions of Chukotka (Siberia) and Seward Peninsula; Pleistocene shoreline processes have reworked the drift to concentrate the heavy metals. The distribution pattern of mercury indicates that particulate mercury-bearing minerals have not been widely dispersed from onland deposits in quantities sufficient to increase mercury levels above normal in offshore sediments of Bering Sea; however, it shows that natural sedimentary processes can concentrate this mercury in beaches of the coastal zone where there already is concern because of potential pollution from man's activities.

Marine Geology↗

USGS2009

Occurrence and distribution of iron, manganese, and selected trace elements in ground water in the glacial aquifer system of the northern United States

Dissolved trace elements, including iron and manganese, are often an important factor in use of ground water for drinking-water supplies in the glacial aquifer system of the United States. The glacial aquifer system underlies most of New England, extends through the Midwest, and underlies portions of the Pacific Northwest and Alaska. Concentrations of dissolved trace elements in ground water can vary over several orders of magnitude across local well networks as well as across regions of the United States. Characterization of this variability is a step toward a regional screening-level assessment of potential human-health implications. Ground-water sampling, from 1991 through 2003, of the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey determined trace element concentrations in water from 847 wells in the glacial aquifer system. Dissolved iron and manganese concentrations were analyzed in those well samples and in water from an additional 743 NAWQA land-use and major-aquifer survey wells. The samples are from monitoring and water-supply wells. Concentrations of antimony, barium, beryllium, cadmium, chromium, cobalt, copper, iron, lead, manganese, molybdenum, nickel, selenium, strontium, thallium, uranium, and zinc vary as much within NAWQA study units (local scale; ranging in size from a few thousand to tens of thousands of square miles) as over the entire glacial aquifer system. Patterns of trace element concentrations in glacial aquifer system ground water were examined by using techniques suitable for a dataset with zero to 80 percent of analytical results reported as below detection. During the period of sampling, the analytical techniques changed, which generally improved the analytical sensitivity. Multiple reporting limits complicated the comparison of detections and concentrations. Regression on Order Statistics was used to model probability distributions and estimate the medians and other quantiles of the trace element concentrations. Strontium and barium were the most frequently detected and usually were present in the highest concentrations. Iron and manganese were the next most commonly detected and next highest in concentrations. Iron concentrations were the most variable with respect to the range of variations (both within local networks and aquifer-wide) and with respect to the disparity between magnitude of concentrations (detections) and the frequency of samples below reporting limits (nondetections). Antimony, beryllium, cadmium, silver, and thallium were detected too infrequently for substantial interpretation of their occurrence or distributions or potential human-health implications. For those elements that were more frequently detected, there are some geographic patterns in their occurrence that primarily reflect climate effects. The highest concentrations of several elements were found in the West-Central glacial framework area (High Plains and northern Plains areas). There are few important patterns for any element in relation to land use, well type, or network type. Shallow land-use (monitor) wells had iron concentrations generally lower than the glacial aquifer system wells overall and much lower than major-aquifer survey wells, which comprise mostly private- and public-supply wells. Unlike those for iron, concentration patterns for manganese were similar among shallow land-use wells and major-aquifer survey wells. An apparent relation between low pH and relatively low concentrations of many elements, except lead, may be more indicative of the relatively low dissolved-solids content in wells in the Northeastern United States that comprise the majority of low pH wells, than of a pH dependent pattern. Iron and manganese have higher concentrations and larger ranges of concentrations especially under more reducing conditions. Dissolved oxygen and well depth were related to iron and manganese concentrations. Redox conditions also affect several trace elements such as arsenic and copper; however, a comparison of redox categories, based in part on iron and manganese concentrations, indicated that the concentrations of many redox-sensitive elements were not significantly different among redox categories. Some of the redox-related patterns were not what would be expected on the basis of solubility constraints. Furthermore, barium is affected by redox conditions in at least one well network even though it is not a redox-sensitive element. Concentrations of barium in portions of the glacial aquifer system are limited by sulfate, which is strongly affected by redox conditions. Few samples had concentrations of any trace element that exceeded drinking-water standards (Maximum Contaminant Levels), for compounds regulated in drinking water or Health-Based Screening Levels for unregulated trace elements. More unregulated trace elements had concentrations greater than benchmarks than regulated trace elements. More samples had manganese concentrations greater its benchmark than any other element in the glacial aquifer system wells. Of the 1,590 wells sampled for manganese, only 556 are for private or public drinking-water supplies, and of those, 9.9 percent (55) exceeded the manganese Lifetime Health Advisory. Concentrations of arsenic, selenium, and uranium less frequently exceeded Maximum Contaminant Levels. There are 29 wells that had 2 element concentrations that exceeded their respective benchmarks. Most concentrations that exceeded a health-based benchmark were from wells in the West-Central area (Iowa, Minnesota, North and South Dakota, Nebraska, and Kansas); however, there is little geographical pattern to the wells with element concentrations of concern.

Scientific Investigations Report↗

USGS2022

Mapping areas of groundwater susceptible to transient contamination events from rapid infiltration into shallow fractured-rock aquifers in agricultural regions of the conterminous United States

Current time-invariant groundwater vulnerability assessments may not capture intermittent contamination events in landscape areas that experience rapid infiltration following precipitation or snowmelt. Occurrences of rapid infiltration and intermittent degradation of groundwater quality are frequently reported in fractured-rock aquifers. This investigation identifies landscape areas underlain by fractured rock within the conterminous United States (CONUS) that may be susceptible to rapid infiltration and where groundwater is a principal source of water supply to the population. Our analysis shows that approximately 27 percent of the CONUS, corresponding to a population of approximately 150 million people, is both underlain by fractured rock and denoted as an area of significant groundwater use. The results of this survey identified shallow fractured-rock aquifers underlying glacial sediments in the upper Midwest and northeastern United States as areas that may be subject to rapid infiltration events. Additionally, aquifers associated with the early Mesozoic basins located in the northeastern and mid-Atlantic United States and bands of carbonate aquifers in the southeastern United States show high susceptibility to rapid infiltration. Index values used in this investigation indicate isolated areas in the western half of the United States also show high susceptibility to rapid infiltration. The isolated areas in Oklahoma, Texas, Arkansas, and southwestern Missouri correspond to karst regions of carbonate aquifers. The isolated areas showing high susceptibility to rapid infiltration and contamination from agricultural sources are locations where more detailed investigations of transient contamination events are warranted. This survey also addresses the potential for contaminant longevity in fractured-rock aquifers stemming from intermittent contamination events. Contaminants that can dissolve into the groundwater following infiltration may be introduced into fractures, and the dissolved constituents can diffuse from fractures into the porosity of the adjacent rock matrix. These constituents can then diffuse back into permeable fractures and adversely affect groundwater quality at downgradient locations over an extended time frame. Rock types with larger matrix porosities have the capacity to retain and then release larger quantities of dissolved constituents, resulting in longer residence times for dissolved groundwater contaminants. The magnitude of the dissolved contaminant concentration infiltrating to the water table will also dictate whether the contaminant concentration in the groundwater exceeds limits for human consumption over the duration of a contamination event. In general, sedimentary- and carbonate-rock aquifers have larger matrix porosities in comparison to igneous- and metamorphic-rock aquifers, and thus, they are more susceptible to longer contaminant residence times. Aquifers composed of sedimentary or carbonate rock constitute approximately 51 percent of the CONUS, and 19 percent of the CONUS is associated with sedimentary- or carbonate-rock aquifers that are of significance for groundwater use. Depending on the contaminants of concern and the concentration of the contaminants introduced into the groundwater from infiltrating water, it would be beneficial for investigations of susceptibility to rapid infiltration to also consider the potential for contaminant longevity. This investigation identifies areas of rapid infiltration into fractured rock using index values applied to the attributes (1) depth to the water table, (2) depth to bedrock, and (3) percentage of sand in soil, where larger index values indicate a greater susceptibility to rapid infiltration. These attributes are selected as the most likely factors that affect rapid infiltration to the water table. The combination of depth to water table and depth to bedrock highlight those aquifer settings that are characterized as shallow fractured-rock aquifers, where the water table may reside either in the bedrock or in overlying unconsolidated geologic materials. In addition, we consider the percentage of agricultural use as a land-use attribute when formulating an index of susceptibility to rapid infiltration and contamination. Agricultural areas are well recognized as nonpoint sources of contaminants that can affect groundwater quality because of seasonal amendments applied to the land surface. Rural agricultural areas are also characterized by septic tanks and leach fields for onsite treatment of wastewater, which may also be a source of contamination that may be introduced into the groundwater following precipitation or snowmelt events.

Conterminous United States↗

USGS2015

Potential groundwater recharge for the State of Minnesota using the Soil-Water-Balance model, 1996-2010

Groundwater recharge is one of the most difficult components of a water budget to ascertain, yet is an important boundary condition necessary for the quantification of water resources. In Minnesota, improved estimates of recharge are necessary because approximately 75 percent of drinking water and 90 percent of agricultural irrigation water in Minnesota are supplied from groundwater. The water that is withdrawn must be supplied by some combination of (1) increased recharge, (2) decreased discharge to streams, lakes, and other surface-water bodies, and (3) removal of water that was stored in the system. Recent pressure on groundwater resources has highlighted the need to provide more accurate recharge estimates for various tools that can assess the sustainability of long-term water use. As part of this effort, the U.S. Geological Survey, in cooperation with the Minnesota Pollution Control Agency, used the Soil-Water-Balance model to calculate gridded estimates of potential groundwater recharge across Minnesota for 1996‒2010 at a 1-kilometer (0.621-mile) resolution. The potential groundwater recharge estimates calculated for Minnesota from the Soil-Water Balance model included gridded values (1-kilometer resolution) of annual mean estimates (that is, the means for individual years from 1996 through 2010) and mean annual estimates (that is, the mean for the 15-year period 1996−2010). The Soil-Water-Balance model uses a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated on a daily basis. A key advantage of this approach includes the use of commonly available geographic information system data layers that incorporate land cover, soil properties, and daily meteorological data to produce temporally and spatially variable gridded estimates of potential recharge. The Soil-Water-Balance model was calibrated by using a combination of parameter estimation techniques, making manual adjustments of model parameters, and using parameter values from previously published Soil-Water-Balance models. Each calibration simulation compared the potential recharge estimate from the model against base-flow estimates derived from three separate hydrograph separation techniques. A total of 35 Minnesota watersheds were selected for the model calibration. Meteorological data necessary for the model included daily precipitation, minimum daily temperature, and maximum daily temperature. All of the meteorological data were provided by the Daymet dataset, which included daily continuous surfaces of key climatological data. Land-cover data were provided by the 2001 and 2006 National Land Cover Database: the 2001 classification was used from 1994 through 2003, and the 2006 classification was used from 2004 through 2010. Soil data used in the model included hydrologic soils group and the available soil-water capacity. These soil data were obtained from the Natural Resources Conservation Service Soil Survey Geographic (SSURGO) database and the State Soil Geographic (STATSGO) database. The statewide mean annual potential recharge rate from 1996–2010 was 4.9 inches per year. Potential recharge estimates increased from west to east across Minnesota. The mean annual potential recharge estimates across Minnesota at a 1-km resolution for the overall simulation period (1996–2010) ranged from less than 0.1 to 17.8 inches per year. Some of the lowest potential recharge rates for the simulation period were in the Red River of the North Basin of northwestern Minnesota, and generally were between 1.0 and 1.5 inches per year. The highest potential recharge rates were in northeastern Minnesota and the Anoka Sand Plain in central Minnesota. Eighty-eight percent of the potential recharge rates (by grid cell) were between 2 and 8 inches per year from 1996–2010. Only about 3 percent of all the potential recharge estimates (by grid cell) were less than 2 inches per year, and 9 percent of estimates were greater than 8 inches per year. On an annual basis, however, potential recharge rates were as high as 27.2 inches per year. The highest annual mean recharge estimate across the State was for 2010, and the lowest mean recharge estimate was for 2003. Although precipitation variability partially explained the annual differences in potential recharge estimates, precipitation alone did not account for these differences, and other factors such as antecedent moisture conditions likely were important. Also, because precipitation gradients across the State can vary from year to year, the dominant land-cover class and hydrologic soil group combinations for a particular region had a large effect on the resulting potential recharge value. During 1996–2010, April had the greatest monthly mean potential recharge compared to all other months, accounting for a mean of 30 percent of annual potential recharge in this single month.

Minnesota↗