Geology ReportsSearch

Geology topics

Perry M. Jones

Publications and source records attributed to Perry M. Jones.

17 recordsLinked to original sources

Assessing the influence of natural copper-nickel-bearing bedrocks of the Duluth Complex on water quality in Minnesota, 2013–15

The U.S. Geological Survey, in cooperation with the University of Minnesota-Duluth Natural Resources Research Institute, completed an assessment of regional water quality in areas of potential base-metal mining in Minnesota. Bedrock, soil, streambed sediment, and surface-water samples were collected in three watersheds that cross the basal part of the Duluth Complex with different mineral-deposit settings: (1) copper-nickel-platinum group element mineralization (Filson Creek), (2) iron-titanium-oxide mineralization (headwaters of the St. Louis River), and (3) no identified mineralization (Keeley Creek). At least 10 bedrock, 30 soil (2 each from 15 sites), and as many as 13 streambed sediment samples were collected in each watershed and analyzed for 44 major and trace elements, total and inorganic carbon, and 10 loosely bound metals (when possible). Surface-water samples were collected at four to nine locations in each watershed three to four times per year for 2 years (total of 141 environmental samples). Surface-water samples were analyzed for 10 trace metals (total and dissolved concentrations), 8 trace elements, 8 major ions (dissolved concentrations), alkalinity, and total and dissolved organic carbon. Metal and element concentrations in solid media varied by watershed, representing local geology. Copper-nickel sulfide mineralization in the Filson Creek watershed was evidenced in bedrock, soil, and streambed sediments. In the Keeley Creek watershed, silicate mineralogy of underlying bedrock contributed metals to streambed sediments. Thick glacial cover masked potential bedrock contributions to solid media in the St. Louis River watershed. Water-quality data indicate that waters in all three watersheds are dilute. Water quality is more similar between the Filson and Keeley Creek watersheds, compared to the St. Louis River watershed, because of the difference in glacial cover. Metal concentrations (copper and nickel, in particular) in surface-water samples follow similar patterns of concentrations in solid media, indicating the influence of bedrock on water quality in Filson and Keeley Creeks. Data from this study provide a baseline of metal concentrations and general water quality within an area of active mineral exploration.

Minnesota

Multi-year data from satellite- and ground-based sensors show details and scale matter in assessing climate’s effects on wetland surface water, amphibians, and landscape conditions

Long-term, interdisciplinary studies of relations between climate and ecological conditions on wetland-upland landscapes have been lacking, especially studies integrated across scales meaningful for adaptive resource management. We collected data in situ at individual wetlands, and via satellite for surrounding 4-km 2 landscape blocks, to assess relations between annual weather dynamics, snow duration, phenology, wetland surface-water availability, amphibian presence and calling activity, greenness, and evapotranspiration in four U.S. conservation areas from 2008 to 2012. Amid recent decades of relatively warm growing seasons, 2012 and 2010 were the first and second warmest seasons, respectively, dating back to 1895. Accordingly, we observed the earliest starts of springtime biological activity during those two years. In all years, early-season amphibians first called soon after daily mean air temperatures were ≥ 0°C and snow had mostly melted. Similarly, satellite-based indicators suggested seasonal leaf-out happened soon after snowmelt and temperature thresholds for plant growth had occurred. Daily fluctuations in weather and water levels were related to amphibian calling activity, including decoupling the timing of the onset of calling at the start of season from the onset of calling events later in the season. Within-season variation in temperature and precipitation also was related to vegetation greenness and evapotranspiration, but more at monthly and seasonal scales. Wetland water levels were moderately to strongly associated with precipitation and early or intermittent wetland drying likely reduced amphibian reproduction success in some years, even though Pseudacris crucifer occupied sites at consistently high levels. Notably, satellite-based indicators of landscape water availability did not suggest such consequential, intra-seasonal variability in wetland surface-water availability. Our cross-disciplinary data show how temperature and precipitation interacted to affect key ecological relations and outcomes on our study landscapes. These results demonstrate the value of multi-year studies and the importance of scale for understanding actual climate-related effects in these areas.

Minnesota, Wisconsin

A multidisciplinary-based conceptual model of a fractured sedimentary bedrock aquitard: improved prediction of aquitard integrity

A hydrogeologic conceptual model that improves understanding of variability in aquitard integrity is presented for a fractured sedimentary bedrock unit in the Cambrian-Ordovician aquifer system of midcontinent North America. The model is derived from multiple studies on the siliciclastic St. Lawrence Formation and adjacent strata across a range of scales and geologic conditions. These studies employed multidisciplinary techniques including borehole flowmeter logging, high-resolution depth-discrete multilevel well monitoring, fracture stratigraphy, fluorescent dye tracing, and three-dimensional (3D) distribution of anthropogenic tracers regionally. The paper documents a bulk aquitard that is highly anisotropic because of poor connectivity of vertical fractures across matrix with low permeability, but with ubiquitous bed parallel partings. The partings provide high bulk horizontal hydraulic conductivity, analogous to aquifers in the system, while multiple preferential termination horizons of vertical fractures serve as discrete low vertical hydraulic conductivity intervals inhibiting vertical flow. The aquitard has substantial variability in its ability to protect underlying groundwater from contamination. Across widespread areas where the aquitard is deeply buried by younger bedrock, preferential termination horizons provide for high aquitard integrity (i.e. protection). Protection is diminished close to incised valleys where stress release and weathering has enhanced secondary pore development, including better connection of fractures across these horizons. These conditions, along with higher hydraulic head gradients in the same areas and more complex 3D flow where the aquitard is variably incised, allow for more substantial transport to deeper aquifers. The conceptual model likely applies to other fractured sedimentary bedrock aquitards within and outside of this region.

Minnesota

Simulation and assessment of groundwater flow and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2003 through 2013: Chapter B of Water levels and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015

Water levels during 2003 through 2013 were less than mean water levels for the period 1925–2013 for several lakes in the northeast Twin Cities Metropolitan Area in Minnesota. Previous periods of low lake-water levels generally were correlated with periods with less than mean precipitation. Increases in groundwater withdrawals and land-use changes have brought into question whether or not recent (2003–13) lake-water-level declines are solely caused by decreases in precipitation. A thorough understanding of groundwater and surface-water exchanges was needed to assess the effect of water-management decisions on lake-water levels. To address this need, the U.S. Geological Survey, in cooperation with the Metropolitan Council and the Minnesota Department of Health, developed and calibrated a three-dimensional, steady-state groundwater-flow model representing 2003–13 mean hydrologic conditions to assess groundwater and lake-water exchanges, and the effects of groundwater withdrawals and precipitation on water levels of 96 lakes in the northeast Twin Cities Metropolitan Area. Lake-water budgets for the calibrated groundwater-flow model indicated that groundwater is flowing into lakes in the northeast Twin Cities Metropolitan Area and lakes are providing water to underlying aquifers. Lake-water outflow to the simulated groundwater system was a major outflow component for Big Marine Lake, Lake Elmo, Snail Lake, and White Bear Lake, accounting for 45 to 64 percent of the total outflows from the lakes. Evaporation and transpiration from the lake surface ranged from 19 to 52 percent of the total outflow from the four lakes. Groundwater withdrawals and precipitation were varied from the 2003‒13 mean values used in the calibrated model (30-percent changes in groundwater withdrawals and 5-percent changes in precipitation) for hypothetical scenarios to assess the effects of groundwater withdrawals and precipitation on water budgets and levels in Big Marine Lake, Snail Lake, and White Bear Lake. Simulated lake-water levels and budgets for Snail Lake and White Bear Lake were affected by 30-percent changes in groundwater withdrawals and 5-percent changes in precipitation in the area, whereas the water level in Big Marine Lake was mainly affected by 5-percent precipitation changes. The effects of groundwater withdrawals on the lake-water levels depend on the number of wells and amount of withdrawals from wells near the lakes. Although lake-water levels are sensitive to precipitation changes, increases in groundwater withdrawals during dry periods exacerbate lake-water level declines. The calibrated, groundwater-flow model is a tool that water-resources managers can use to address future water management issues in the northeast Twin Cities Metropolitan Area.

Minnesota

Statistical analysis of lake levels and field study of groundwater and surface-water exchanges in the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015: Chapter A of Water levels and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015

Water levels declined from 2003 to 2011 in many lakes in Ramsey and Washington Counties in the northeast Twin Cities Metropolitan Area, Minnesota; however, water levels in other northeast Twin Cities Metropolitan Area lakes increased during the same period. Groundwater and surface-water exchanges can be important in determining lake levels where these exchanges are an important component of the water budget of a lake. An understanding of groundwater and surface-water exchanges in the northeast Twin Cities Metropolitan Area has been limited by the lack of hydrologic data. The U.S. Geological Survey, in cooperation with the Metropolitan Council and Minnesota Department of Health, completed a field and statistical study assessing lake-water levels and regional and local groundwater and surface-water exchanges near northeast Twin Cities Metropolitan Area lakes. This report documents the analysis of collected hydrologic, water-quality, and geophysical data; and existing hydrologic and geologic data to (1) assess the effect of physical setting and climate on lake-level fluctuations of selected lakes, (2) estimate potential percentages of surface-water contributions to well water across the northeast Twin Cities Metropolitan Area, (3) estimate general ages for waters extracted from the wells, and (4) assess groundwater inflow to lakes and lake-water outflow to aquifers downgradient from White Bear Lake. Statistical analyses of lake levels during short-term (2002–10) and long-term (1925–2014) periods were completed to help understand lake-level changes across the northeast Twin Cities Metropolitan Area. Comparison of 2002–10 lake levels to several landscape and geologic characteristics explained variability in lake-level changes for 96 northeast Twin Cities Metropolitan Area lakes. Application of several statistical methods determined that (1) closed-basin lakes (without an active outlet) had larger lake-level declines than flow-through lakes with an outlet; (2) closed-basin lake-level changes reflected groundwater-level changes in the Quaternary, Prairie du Chien, and Jordan aquifers; (3) the installation of outlet-control structures, such as culverts and weirs, resulted in smaller multiyear lake-level changes than lakes without outlet-control structures; (4) water levels in lakes primarily overlying Superior Lobe deposits were significantly more variable than lakes primarily overlying Des Moines Lobe deposits; (5) lake-level declines were larger with increasing mean lake-level elevation; and (6) the frequency of some of these characteristics varies by landscape position. Flow-through lakes and lakes with outlet-control structures were more common in watersheds with more than 50 percent urban development compared to watersheds with less than 50 percent urban development. A comparison of two 35-year periods during 1925–2014 revealed that variability of annual mean lake levels in flow-through lakes increased when annual precipitation totals were more variable, whereas variability of annual mean lake levels in closed-basin lakes had the opposite pattern, being more variable when annual precipitation totals were less variable. Oxygen-18/oxygen-16 and hydrogen-2/hydrogen-1 ratios for water samples from 40 wells indicated the well water was a mixture of surface water and groundwater in 31 wells, whereas ratios from water sampled from 9 other wells indicated that water from these wells receive no surface-water contribution. Of the 31 wells with a mixture of surface water and groundwater, 11 were downgradient from White Bear Lake, likely receiving water from deeper parts of the lake. Age dating of water samples from wells indicated that the age of water in the Prairie du Chien and Jordan aquifers can vary widely across the northeast Twin Cities Metropolitan Area. Estimated ages of recharge for 9 of the 40 wells sampled for chlorofluorocarbon concentrations ranged widely from the early 1940s to mid-1970s. The wide range in estimated ages of recharge may have resulted from the wide range in the open-interval lengths and depths for the wells. Results from stable isotope analyses of water samples, lake-sediment coring, continuous seismic-reflection profiling, and water-level and flow monitoring indicated that there is groundwater inflow from nearshore sites and lake-water outflow from deep-water sites in White Bear Lake. Continuous seismic-reflection profiling indicated that deep sections of White Bear, Pleasant, Turtle, and Big Marine Lakes have few trapped gases and little organic material, which indicates where groundwater and lake-water exchanges are more likely. Water-level differences between White Bear Lake and piezometer and seepage measurements in deep waters of the lake indicate that groundwater and lake-water exchange is happening in deep waters, predominantly downgradient from the lake and into the lake sediment. Seepage fluxes measured in the nearshore sites of White Bear Lake generally were higher than seepage fluxes measured in the deep-water sites, which indicates that groundwater-inflow rates at most of the nearshore sites are higher than lake-water outflow from the deep-water sites.

Minnesota

Assessment of aquifer properties, evapotranspiration, and the effects of ditching in the Stoney Brook watershed, Fond du Lac Reservation, Minnesota, 2006-9

The U.S. Geological Survey, in cooperation with the Fond du Lac Band of Lake Superior Chippewa, assessed hydraulic properties of geologic material, recharge, and evapotranspiration, and the effects of ditching on the groundwater resources in the Stoney Brook watershed in the Fond du Lac Reservation. Geologic, groundwater, and surface-water data were collected during 2006–9 to estimate hydrologic properties in the watershed. Streamflow and groundwater levels in the shallow glacial deposits in the Stoney Brook watershed were analyzed to estimate groundwater-flow directions, groundwater recharge, and evapotranspiration within the watershed and to assess the effect of ditches on surrounding groundwater resources. Groundwater, streamflow, and precipitation data collected during the study (2006–9) can be used to update the U.S. Department of Agriculture’s Natural Resource Conservation Service and Fond du Lac Resource Management Division surface-water models, which are used to evaluate the effect of proposed adjustments to the ditching system on streamflow on wild rice production and aquatic habitats. Specific yields calculated from the well water levels ranged from 0.11 to 0.40, and hydraulic conductivities determined from water levels measured during well slug tests ranged from 1 to 7 feet per day. The values for specific yields were similar to values obtained in other studies done in glacial materials of similar composition in Minnesota. The higher hydraulic conductivity estimate (7 feet per day) was similar to lower hydraulic conductivities estimated in another hydrologic study conducted in Carlton County, Minnesota. The installation of drainage ditches in the Stoney Brook watershed has reduced water levels in lakes connected to the ditch system, and has locally reduced groundwater levels in shallow groundwater adjacent to the ditches and lakes. Differences in near-ditch groundwater hydrographs relative to far-ditch groundwater hydrographs indicate that the effect of the ditches on groundwater is only localized to near-ditch areas. These hydrograph differences resulted in large differences between recharge estimated at wells near and far from ditches. In this study, recharge estimated at wells within 50 feet of a ditch was influenced by ditch-water levels. Annual groundwater recharge estimates from water levels and streamflows during 2006–9 ranged from 0.36 to 34.8 inches, and varied with climate, geology, and well location relative to ditches. The higher recharge estimates were determined from analysis of groundwater levels in wells near the ditches because the shallow groundwater in these wells received both infiltration from ditches and areal groundwater recharge from precipitation. The water-table fluctuation method using a manual groundwater recession approach for wells far from ditches provided the best estimates of areal groundwater recharge to the shallow glacial aquifer because water levels in these wells were not affected by water infiltrating from ditches (bank storage). For wells more than 400 feet from ditches, mean annual areal groundwater recharge estimates using the manual groundwater recession approach for wells screened mostly in outwash sands during 2007, 2008, and 2009 ranged from 4.47 to 18.6 inches (wells 5, 7, 13, 14 and 15), and ranged from 0.43 to 2.85 inches for wells screened mostly in clayey sand or sandy clay (wells 9 and 16). Recharge estimates at wells far from ditches were similar to basinwide recharge estimates from streamflow. Daily fluctuations in water levels in two wells indicated that the evapotranspiration extinction depth in the Stoney Brook watershed is approximately 4.6 to 6 feet below the land surface. A polynomial regression fit of the daily evapotranspiration rates during 2006–9 for well 1 produced a total evapotranspiration estimate of 16.1 inches from June 26 to October 6 for every year. Evapotranspiration estimated from daily water-level fluctuations in wells near ditches is relatively high. The ditch-water surface allowed for relatively high evaporation compared to the land surface, which, with a good hydraulic connection to surrounding groundwater, resulted in relatively high fluctuations in daily groundwater levels near ditches, resulting in high evapotranspiration estimates.

Minnesota

Copper toxicity and organic matter: Resiliency of watersheds in the Duluth Complex, Minnesota, USA

We estimated copper (Cu) toxicity in surface water with high dissolved organic matter (DOM) for unmined mineralized watersheds of the Duluth Complex using the Biotic Ligand Model (BLM), which evaluates the effect of DOM, cation competition for biologic binding sites, and metal speciation. A sediment-based BLM was used to estimate stream-sediment toxicity; this approach factors in the cumulative effects of multiple metals, incorporation of metals into less bioavailable sulfides, and complexation of metals with organic carbon. For surface water, the formation of Cu-DOM complexes significantly reduces the amount of Cu available to aquatic organisms. The protective effects of cations, such as calcium (Ca) and magnesium (Mg), competing with Cu to complex with the biotic ligand is likely not as important as DOM in water with high DOM and low hardness. Standard hardness-based water quality criteria (WQC) are probably inadequate for describing Cu toxicity in such waters and a BLM approach may yield more accurate results. Nevertheless, assumptions about relative proportions of humic acid (HA) and fulvic acid (FA) in DOM significantly influence BLM results; the higher the HA fraction, the higher calculated resiliency of the water to Cu toxicity. Another important factor is seasonal variation in water chemistry, with greater resiliency to Cu toxicity during low flow compared to high flow. Based on generally low total organic carbon and sulfur content, and equivalent metal ratios from total and weak partial extractions, much of the total metal concentration in clastic streambedsediments may be in bioavailable forms, sorbed on clays or hydroxide phases. However, organicrich fine-grained sediment in the numerous wetlands may sequester significant amount of metals, limiting their bioavailability. A high proportion of organic matter in waters and some sediments will play a key role in the resiliency of these watersheds to potential additional metal loads associated with future mining operations.

Minnesota

Groundwater–surface-water exchange and the geologic setting of northern Minnesota's lakes, wetlands, and streams—Modern-day relevance of Tom Winter's legacy

Tom Winter spent nearly 50 years conducting research in earth science, and he specialized in the exchange between groundwater and surface water. Tom's highly productive career began in Minnesota. This fi eld trip revisits many of the places where Tom conducted his early research and demonstrates the continuing relevance of that research. Stops and topics include the groundwater infl uence on the record low stage of White Bear Lake, the contribution of groundwater to continually rising water levels in an abandoned open-pit iron mine, hydrogeology of the Shingobee headwaters aquatic ecosystem research site, hydrogeology of Lake Sallie, geology associated with the Pillager water gap, and the hydrogeology of Little Rock Lake.

GSA Field Guides

Groundwater and surface-water interactions near White Bear Lake, Minnesota, through 2011

The U.S. Geological Survey, in cooperation with the White Bear Lake Conservation District, the Minnesota Pollution Control Agency, the Minnesota Department of Natural Resources, and other State, county, municipal, and regional planning agencies, watershed organizations, and private organizations, conducted a study to characterize groundwater and surface-water interactions near White Bear Lake through 2011. During 2010 and 2011, White Bear Lake and other lakes in the northeastern part of the Twin Cities Metropolitan Area were at historically low levels. Previous periods of lower water levels in White Bear Lake correlate with periods of lower precipitation; however, recent urban expansion and increased pumping from the Prairie du Chien-Jordan aquifer have raised the question of whether a decline in precipitation is the primary cause for the recent water-level decline in White Bear Lake. Understanding and quantifying the amount of groundwater inflow to a lake and water discharge from a lake to aquifers is commonly difficult but is important in the management of lake levels. Three methods were used in the study to assess groundwater and surface-water interactions on White Bear Lake: (1) a historical assessment (1978-2011) of levels in White Bear Lake, local groundwater levels, and their relation to historical precipitation and groundwater withdrawals in the White Bear Lake area; (2) recent (2010-11) hydrologic and water-quality data collected from White Bear Lake, other lakes, and wells; and (3) water-balance assessments for White Bear Lake in March and August 2011. An analysis of covariance between average annual lake-level change and annual precipitation indicated the relation between the two variables was significantly different from 2003 through 2011 compared with 1978 through 2002, requiring an average of 4 more inches of precipitation per year to maintain the lake level. This shift in the linear relation between annual lake-level change and annual precipitation indicated the net effect of the non-precipitation terms on the water balance has changed relative to precipitation. The average amount of precipitation required each year to maintain the lake level has increased from 33 inches per year during 1978-2002 to 37 inches per year during 2003-11. The combination of lower precipitation and an increase in groundwater withdrawals can explain the change in the lake-level response to precipitation. Annual and summer groundwater withdrawals from the Prairie du Chien-Jordan aquifer have more than doubled from 1980 through 2010. Results from a regression model constructed with annual lake-level change, annual precipitation minus evaporation, and annual volume of groundwater withdrawn from the Prairie du Chien-Jordan aquifer indicated groundwater withdrawals had a greater effect than precipitation minus evaporation on water levels in the White Bear Lake area for all years since 2003. The recent (2003-11) decline in White Bear Lake reflects the declining water levels in the Prairie du Chien-Jordan aquifer; increases in groundwater withdrawals from this aquifer are a likely cause for declines in groundwater levels and lake levels. Synoptic, static groundwater-level and lake-level measurements in March/April and August 2011 indicated groundwater was potentially flowing into White Bear Lake from glacial aquifers to the northeast and south, and lake water was potentially discharging from White Bear Lake to the underlying glacial and Prairie du Chien-Jordan aquifers and glacial aquifers to the northwest. Groundwater levels in the Prairie du Chien-Jordan aquifer below White Bear Lake are approximately 0 to 19 feet lower than surface-water levels in the lake, indicating groundwater from the aquifer likely does not flow into White Bear Lake, but lake water may discharge into the aquifer. Groundwater levels from March/April to August 2011 declined more than 10 feet in the Prairie du Chien-Jordan aquifer south of White Bear Lake and to the north in Hugo, Minnesota. Water-quality analyses of pore water from nearshore lake-sediment and well-water samples, seepage-meter measurements, and hydraulic-head differences measured in White Bear Lake also indicated groundwater was potentially flowing into White Bear Lake from shallow glacial aquifers to the east and south. Negative temperature anomalies determined in shallow waters in the water-quality survey conducted in White Bear Lake indicated several shallow-water areas where groundwater may be flowing into the lake from glacial aquifers below the lake. Cool lake-sediment temperatures (less than 18 degrees Celsius) were measured in eight areas along the northeast, east, south, and southwest shores of White Bear Lake, indicating potential areas where groundwater may flow into the lake. Stable isotope analyses of well-water, precipitation, and lake-water samples indicated wells downgradient from White Bear Lake screened in the glacial buried aquifer or open to the Prairie du Chien-Jordan aquifer receive a mixture of surface water and groundwater; the largest surface-water contributions are in wells closer to White Bear Lake. A wide range in oxygen-18/oxygen-16 and deuterium/protium ratios was measured in well-water samples, indicating different sources of water are supplying water to the wells. Well water with oxygen-18/oxygen-16 and deuterium/protium ratios that plot close to the meteoric water line consisted mostly of groundwater because deuterium/protium ratios for most groundwater usually are similar to ratios for rainwater and snow, plotting close to meteoric water lines. Well water with oxygen-18/oxygen-16 and deuterium/protium ratios that plot between the meteoric water line and ratios for the surface-water samples from White Bear Lake consists of a mixture of surface water and groundwater; the percentage of each source varies relative to its ratios. White Bear Lake is the likely source of the surface water to the wells that have a mixture of surface water and groundwater because (1) it is the only large, deep lake near these wells; (2) these wells are near and downgradient from White Bear Lake; and (3) these wells obtain their water from relatively deep depths, and White Bear Lake is the deepest lake in that area. The percentages of surface-water contribution to the three wells screened in the glacial buried aquifer receiving surface water were 16, 48, and 83 percent. The percentages of surface-water contribution ranged from 5 to 79 percent for the five wells open to the Prairie du Chien-Jordan aquifer receiving surface water; wells closest to White Bear Lake had the largest percentages of surface-water contribution. Water-balance analysis of White Bear Lake in March and August 2011 indicated a potential discharge of 2.8 and 4.5 inches per month, respectively, over the area of the lake from the lake to local aquifers. Most of the sediments from a 12.4-foot lake core collected at the deepest part of White Bear Lake consisted of silts, sands, and gravels likely slumped from shallower waters, with a very low amount of low-permeability, organic material.

Minnesota

Evaluation of Methods for Delineating Zones of Transport for Production Wells in Karst and Fractured-Rock Aquifers of Minnesota

Assessment of groundwater-flow conditions in the vicinity of production wells in karst and fractured-rock settings commonly is difficult due in part to the lack of detailed hydrogeologic information and the resources needed to collect it. To address this concern and to better understand the hydrogeology and aquifer properties of karst and fractured-rock aquifers in Minnesota, the U.S. Geological Survey, in cooperation with the Minnesota Department of Health, conducted a study to evaluate methods for delineating zones of transport for 24 production wells in karst and fractured-rock aquifers in Minnesota. Two empirical methods for delineating zones of transport around wells were applied to the 24 production wells that extract groundwater from karst and fractured-rock aquifers in nine Minnesota communities. These methods were the truncated-parabola and modified-ellipse methods, and both methods assume porous-media flow conditions. The 24 wells extracted water from a karst aquifer (Prairie du Chien-Jordan aquifer), porous aquifers interspersed with solution-enhanced fractures (Jordan and Hinckley aquifers), or fractured-bedrock aquifers (Biwabik Iron-Formation and Sioux Quartzite aquifers). Zones of transport delineated using these two empirical methods were compared with zones of transport previously delineated by Minnesota Department of Health hydrologists for the wells using the calculated-fixed-radius method and groundwater-flow models. Large differences were seen in the size and shapes of most zones of transport delineated using the truncated-parabola and modified-ellipse methods compared with the zones of transport delineated by the Minnesota Department of Health. In general, the zones of transport delineated using the truncated-parabola and modified-ellipse methods were smaller in area than those delineated by the Minnesota Department of Health and included only small parts of the Minnesota Department of Health zones of transport. About two-thirds(67 percent) of the individual or composite truncated parabolas and modified ellipses covered less than 50 percent of the area included in zones of transport delineated by the Minnesota Department of Health. The shapes of some of the truncated parabola and modified ellipses did not closely match the zones of transport delineated by the Minnesota Department of Health using the calculated-fixed-radius method and groundwater-flow models. Differences between the zones of transport delineated by the truncated-parabola and modified-ellipse methods and those delineated by the Minnesota Department of Health can be explained by variations inherent to the methods and by the amount of complexity taken into account by different groundwater-flow models. Additional field hydrogeologic studies would be needed at specific sites to support the use of these zone-of-transport delineation methods. Application of the truncated-parabola and modified-ellipse methods to sites for which existing hydrogeologic information is limited can produce questionable results in karst and fractured-rock settings, particularly in areas where many high-capacity wells or active mining operations are nearby.

Scientific Investigations Report

Water quality (2000-08) and historical phosphorus concentrations from paleolimnological studies of Swamp and Speckled Trout Lakes, Grand Portage Reservation, northeastern Minnesota

A paleolimnological approach was taken to aid the Grand Portage Reservation, in northeastern Minnesota, in determining reference conditions for lakes on the reservation. The U.S. Geological Survey, in cooperation with the Grand Portage Band of Chippewa Indians and the Science Museum of Minnesota, conducted a study to describe water quality (2000-08) and historical total phosphorus concentrations (approximately 1781-2006) for Swamp and Speckled Trout Lakes. Results from this study may be used as a guide in establishing nutrient criteria in these and other lakes on the Grand Portage Reservation. Historical phosphorus concentrations were inferred through paleolimnological reconstruction methods involving diatom analysis and lead-210 dating of lake-sediment cores. Historical diatom-inferred total phosphorus concentrations in Swamp Lake ranged from 0.017 to 0.025 milligrams per liter (mg/L) based on diatom assemblages in sediment samples dated 1781-2005. Historical diatom-inferred total phosphorus concentrations in Speckled Trout Lake ranged from 0.008 to 0.014 mg/L based on diatom assemblages in sediment samples dated 1825-2006. In both lakes, historical changes in diatom-inferred total phosphorus concentrations did not exceed model error estimates, indicating that there has been minimal change in total phosphorus concentrations in the two lakes over about two centuries. Nutrient concentrations in monthly water samples collected May through October during 2000, 2002, 2004, 2006, and 2008 were compared to the diatom-inferred total phosphorus concentrations. Total phosphorus concentrations from water samples collected from Swamp Lake during 2000-08 ranged from less than 0.002 to 0.160 mg/L (median= 0.023 mg/L) compared to diatom-inferred total phosphorus concentrations of 0.018 to 0.020 mg/L for 2002 to 2005. Total phosphorus concentrations in water samples collected from Speckled Trout Lake during 2000-08 were similar to those of Swamp Lake, ranging from less than 0.002 to 0.147 mg/L (median=0.012 mg/L), whereas the diatom-inferred total phosphorus concentrations were smaller, ranging from 0.009 to 0.010 mg/L for 2003 to 2006. Differences in total phosphorus concentrations between the two lakes may be because of differences in watershed characteristics, particularly the number of wetlands in the two watersheds. Similarities between recent total phosphorus concentrations in water-quality samples and diatom-inferred total phosphorus indicate that diatom-inferred phosphorus reconstructions might be used to help establish reference conditions. Nutrient criteria for Grand Portage Reservation lakes may be established when a sampling program is designed to ensure representative phosphorus concentrations in water samples are comparable to diatom-inferred concentrations.

Scientific Investigations Report

Groundwater-flow model of the Ozark Plateaus aquifer system, northwestern Arkansas, southeastern Kansas, southwestern Missouri, and northeastern Oklahoma

To assess the effect that increased water use is having on the long-term availability of groundwater within the Ozark Plateaus aquifer system, a groundwater-flow model was developed using MODFLOW 2000 for a model area covering 7,340 square miles for parts of Arkansas, Kansas, Missouri, and Oklahoma. Vertically the model is divided into five units. From top to bottom these units of variable thickness are: the Western Interior Plains confining unit, the Springfield Plateau aquifer, the Ozark confining unit, the Ozark aquifer, and the St. Francois confining unit. Large mined zones contained within the Springfield Plateau aquifer are represented in the model as extensive voids with orders-of-magnitude larger hydraulic conductivity than the adjacent nonmined zones. Water-use data were compiled for the period 1960 to 2006, with the most complete data sets available for the period 1985 to 2006. In 2006, total water use from the Ozark aquifer for Missouri was 87 percent (8,531,520 cubic feet per day) of the total pumped from the Ozark aquifer, with Kansas at 7 percent (727,452 cubic feet per day), and Oklahoma at 6 percent (551,408 cubic feet per day); water use for Arkansas within the model area was minor. Water use in the model from the Springfield Plateau aquifer in 2005 was specified from reported and estimated values as 569,047 cubic feet per day. Calibration of the model was made against average water-level altitudes in the Ozark aquifer for the period 1980 to 1989 and against waterlevel altitudes obtained in 2006 for the Springfield Plateau and Ozark aquifers. Error in simulating water-level altitudes was largest where water-level altitude gradients were largest, particularly near large cones of depression. Groundwater flow within the model area occurs generally from the highlands of the Springfield Plateau in southwestern Missouri toward the west, with localized flow occurring towards rivers and pumping centers including the five largest pumping centers near Joplin, Missouri; Carthage, Missouri; Noel, Missouri; Pittsburg, Kansas; and Miami, Oklahoma. Hypothetical scenarios involving various increases in groundwater-pumping rates were analyzed with the calibrated groundwater-flow model to assess changes in the flow system from 2007 to the year 2057. Pumping rates were increased between 0 and 4 percent per year starting with the 2006 rates for all wells in the model. Sustained pumping at 2006 rates was feasible at the five pumping centers until 2057; however, increases in pumping resulted in dewatering the aquifer and thus pumpage increases were not sustainable in Carthage and Noel for the 1 percent per year pumpage increase and greater hypothetical scenarios, and in Joplin and Miami for the 4 percent per year pumpage increase hypothetical scenarios. Zone-budget analyses were performed to assess the groundwater flow into and out of three zones specified within the Ozark-aquifer layer of the model. The three zones represented the model part of the Ozark aquifer in Kansas (zone 1), Oklahoma (zone 2), and Missouri and Arkansas (zone 3). Groundwater pumping causes substantial reductions in water in storage and induces flow through the Ozark confining unit for all hypothetical scenarios evaluated. Net simulated flow in 2057 from Kansas (zone 1) to Missouri (zone 3) ranges from 74,044 cubic feet per day for 2006 pumping rates (hypothetical scenario 1) to 625,319 cubic feet per day for a 4 percent increase in pumping per year (hypothetical scenario 5). Pumping from wells completed in the Ozark aquifer is the largest component of flow out of zone 3 in Missouri and Arkansas, and varies between 88 to 91 percent of the total flow out of zone 3 for all of the hypothetical scenarios. The largest component of flow into Oklahoma (zone 2) comes from the overlying Ozark confining unit, which is consistently about 45 percent of the total. Flow from the release of water in storage, from general-head boundaries, and from zones 1 and 3 is considerably smaller values that range from 3 to 22 percent of the total flow into zone 2. The largest flow out of the Oklahoma part of the model occurs from pumping from wells and ranges from 52 to 69 percent of the total.

Arkansas, Kansas, Missouri, Oklahoma

Ground-water/surface-water interaction in nearshore areas of Three Lakes on the Grand Portage Reservation, northeastern Minnesota, 2003-04

The U.S. Geological Survey (USGS), in cooperation with the Grand Portage Band of Chippewa Indians, applied three techniques to assess ground-water/surface-water interaction in nearshore areas of three lakes (North, Teal, and Taylor) on the Grand Portage Reservation in northeastern Minnesota. At each lake, analyses of existing aerial photographs, in-situ temperature measurements of shoreline lake sediment, and chemical analyses of surface water and pore water were conducted. Surface-water and pore-water samples were analyzed for major constituents, nutrients, and stable isotopes of oxygen and hydrogen. Bulk precipitation samples were collected and analyzed (1) for nutrient concentrations to determine nutrient input to the lakes through atmospheric deposition and (2) for stable isotope ratios of oxygen and hydrogen to determine a meteoric waterline that was needed for the stable isotope analyses of surface-water and pore-water samples. Total nitrogen concentrations in the precipitation samples ranged from 0.51 to 8.4 mg/L (milligrams per liter) as nitrogen at the North Lake precipitation station and from 0.42 to 2.3 mg/L as nitrogen at the Grand Portage precipitation station. Oxygen-18/oxygen-16 and deuterium/protium isotope ratios for the bulk precipitation samples lie relatively close to a meteoric waterline for northern Wisconsin, except for the ratios for samples collected on May 20, 2004. Analyses of existing aerial photographs, nearshore lake-sediment temperatures, and seasonal isotope ratios of surface-water and pore-water samples were the most valuable data for identifying locations of ground-water inflow and surface-water outseepage. Analyses of existing aerial photographs of the three lakes indicated the location of potential inflow channels and lineaments identifying potential ground-water inflow locations for pore-water sampling. Lake-sediment temperatures at potential ground-water inflow locations ranged from 4 to 16 ºC, varying between lakes, seasons, and climatic conditions. Major constituent chemistry was valuable at Taylor Lake, and to a limited extent at North and Teal Lakes, in confirming results from the isotope and lake-sediment temperature data. Ground-water inflow to North Lake likely occurs along the southwest and south shores, and along portions of the west, southeast, north, and northeast shores. Relatively cool lake-sediment temperatures along the southwest, south, west, and southeast shores, and in isolated beaver channels along the north and northeast shores of North Lake indicate potential ground-water inflow at these locations. Both localized ground-water inflow and surface-water outseepage occurs along portions of the north, northeast, southeast, and south shores, varying seasonally. Conflicting evidence for ground-water flow conditions exist for the northwest and north-northwest pore-water samples. Only minor differences in the major constituent concentrations were seen between the surface-water and pore-water samples from the North Lake area with the exception of iron and manganese concentrations. Ground-water inflow likely takes place along the south-southwest and north shores of Teal Lake, with a mixture of ground-water inflow and surface-water outseepage occurring in other areas of the lake. Cooler lake-sediment temperatures occurred along the south-southwest, west, and northwest shores, portions of the north shore, and in channels identified in aerial photographs throughout the lake, indicating potential ground-water inflow at those locations. Warmer lake-sediment temperatures along the northeast and portions of the southwest and northwest shores of Teal Lake indicate potential locations where surface-water outseepage or little ground- and surface-water interaction occurs. The major constituent concentrations were higher in the pore-water samples collected from the south-southwest and northeast shores of Teal Lake, indicating ground-water inflow. Cation adsorption, cation exchanges with hydrogen ions, and chelation with organic materials occurring in the fen surrounding the lake likely resulted in the low dissolved calcium, magnesium, and sodium concentrations in north, northwest, and west pore-water samples from the Teal Lake area. Pore-water samples from the south-southwest, north, and southwest shores of Teal Lake had isotopic compositions that plotted closest to the meteoric waterline, indicating that little evaporation or transpiration occurred in these samples and that ground-water inflow may be occurring at these locations. Surface-water outseepage from Teal Lake likely occurs along the northeast shore even though major constituent concentrations were high. Major constituent concentrations may be high because of a nearby beaver dam. Ground-water inflow to Taylor Lake likely occurs at the north and south pore-water sampling sites. Higher major constituent concentrations and the least evaporative isotope ratios were found in pore-water samples along the south, north, and west shores of Taylor Lake, indicating potential locations of ground-water inflow. However, a combination of warmer and cooler lake-sediment temperatures along the west lowland indicated that ground-water inflow and surface-water outseepage may occur at that location. Surface-water outseepage likely occurs from Taylor Lake along the south shore through a surface-water drainage channel to a downgradient bog. Warmer lake-sediment temperatures along portions of the south and southeast shores indicate that surface-water outseepage may occur at those locations. Both ground-water inflow and surface-water outseepage may occur along the west, southeast, and east shores of Taylor Lake, varying seasonally and with local precipitation. Knowledge of general water-flow directions in lake watersheds and how they may change seasonally can help water-quality specialists and lake managers address a variety of water-quality and aquatic habitat protection issues for lakes. Results from this study indicate that ground-water and surface-water interactions at the study lakes are complex, and the ability of the applied techniques to identify ground-water inflow and surface-water outseepage locations varied among the lakes. Measurement of lake-sediment temperatures proved to be a reliable and relatively inexpensive reconnaissance technique that lake managers may apply in complex settings to identify general areas of ground-water inflow and surface-water outseepage.

Minnesota

Simulated effects of water-level changes in the Mississippi River and Pokegama Reservoir on ground-water levels, Grand Rapids area, Minnesota

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, used an existing, three-dimensional, numerical ground-water flow model (referred to as the calibrated model) to assess the effects of water-level changes in the Mississippi River and Pokegama Reservoir on ground-water levels in adjacent glaciofluvial aquifers in the Grand Rapids area of north-central Minnesota. Pokegama Reservoir consists of Pokegama Lake, Little Jay Gould Lake, Jay Gould Lake, Cut-off Lake, and Blackwater Lake. Water levels in the Pokegama Reservoir are regulated at Pokegama Dam on the Mississippi River west of Grand Rapids. A steady-state model was used, and simulations represent “worse-case” scenarios for the effects of lowering or raising the river and lake water levels. The simulated ground-water levels represent levels that would result if the river and lake stages permanently declined or rose by the specified amounts. Eight hypothetical scenarios were simulated by varying water levels in the Mississippi River and Pokegama Reservoir from values used in the calibrated model. In the simulations, water levels for the Mississippi River, riverine wetlands of the Mississippi River, and lakes of the Pokegama Reservoir were raised and lowered uniformly by 0.50, 1.00, 2.00, and 3.00 feet from calibrated water levels. The extent of aquifer water-level changes resulting from these river, wetland, and lake water-level changes varied because of the complex hydrogeology of the study area. A 1.00-foot decline in reservoir/river water levels caused a maximum simulated ground-water-level decline in the middle aquifer near Jay Gould and Little Jay Gould Lakes of 1.09 feet and a maximum simulated ground-water-level decline of 1.00 foot in the lower aquifer near Cut-off and Blackwater Lakes. The amount and extent of ground-water-level changes in the middle and lower aquifers can be explained by the thickness, extent, and connectivity of the aquifers. Surface-water/ground-water interactions near wetlands and lakes with water levels unchanged from the calibrated model resulted in small water-table altitude differences among the simulations. Results of the ground-water modeling indicate that lowering of the reservoir and river water levels by 1.00 foot likely will not substantially affect water levels in the middle and lower aquifers.

Minnesota

Simulation of ground-water flow in glaciofluvial aquifers in the Grand Rapids area, Minnesota

A calibrated steady-state, finite-difference, ground-waterflow model was constructed to simulate ground-water flow in three glaciofluvial aquifers, defined in this report as the upper, middle, and lower aquifers, in an area of about 114 mi2 surrounding the city of Grand Rapids in north-central Minnesota. The calibrated model will be used by Minnesota Department of Health and communities in the Grand Rapids area in the development of wellhead protection plans for their water supplies. The model was calibrated through comparison of simulated ground-water levels to measured static water levels in 351 wells, and comparison of simulated base-flow rates to estimated base-flow rates for reaches of the Mississippi and Prairie Rivers. Model statistics indicate that the model tends to overestimate ground-water levels. The root mean square errors ranged from +12.83 ft in wells completed in the upper aquifer to +19.10 ft in wells completed in the middle aquifer. Mean absolute differences between simulated and measured water levels ranged from +4.43 ft for wells completed in the upper aquifer to +9.25 ft for wells completed in the middle aquifer. Mean algebraic differences ranged from +9.35 ft for wells completed in the upper aquifer to +14.44 ft for wells completed in the middle aquifer, with the positive differences indicating that the simulated water levels were higher than the measured water levels. Percentage errors between simulated and estimated base-flow rates for the three monitored reaches all were less than 10 percent, indicating good agreement. Simulated ground-water levels were most sensitive to changes in general-head boundary conductance, indicating that this characteristic is the predominant model input variable controlling steady-state water-level conditions. Simulated groundwater flow to stream reaches was most sensitive to changes in horizontal hydraulic conductivity, indicating that this characteristic is the predominant model input variable controlling steady-state flow conditions.

Minnesota

Characterization of rainfall-runoff response and estimation of the effect of wetland restoration on runoff, Heron Lake Basin, southwestern Minnesota, 1991-97

The U.S. Geological Survey (USGS), in cooperation with the Minnesota Department of Natural Resources and the Heron Lake Watershed District, conducted a study to characterize the rainfall-runoff response and to examine the effects of wetland restoration on the rainfall-runoff response within the Heron Lake Basin in southwestern Minnesota. About 93 percent of the land cover in the Heron Lake Basin consists of agricultural lands, consisting almost entirely of row crops, with less than one percent consisting of wetlands. The Hydrological Simulation Program – Fortran (HSPF), Version 10, was calibrated to continuous discharge data and used to characterize rainfall-runoff responses in the Heron Lake Basin between May 1991 and August 1997. Simulation of the Heron Lake Basin was done as a two-step process: (1) simulations of five small subbasins using data from August 1995 through August 1997, and (2) simulations of the two large basins, Jack and Okabena Creek Basins, using data from May 1991 through September 1996. Simulations of the five small subbasins was done to determine basin parameters for the land segments and assess rainfall-runoff response variability in the basin. Simulations of the two larger basins were done to verify the basin parameters and assess rainfall-runoff responses over a larger area and for a longer time period. Best-fit calibrations of the five subbasin simulations indicate that the rainfall-runoff response is uniform throughout the Heron Lake Basin, and 48 percent of the total rainfall for storms becomes direct (surface and interflow) runoff. Rainfall-runoff response variations result from variations in the distribution, intensity, timing, and duration of rainfall; soil moisture; evapotranspiration rates; and the presence of lakes in the basin. In the spring, the amount and distribution of rainfall tends to govern the runoff response. High evapotranspiration rates in the summer result in a depletion of moisture from the soils, substantially affecting the rainfall-runoff relation. Five wetland restoration simulations were run for each of five subbasins using data from August 1995 through August 1997, and for the two larger basins, Jack and Okabena Creek Basins, using data from May 1991 through September 1996. Results from linear regression analysis of total simulated direct runoff and total rainfall data for simulated storms in the wetland-restoration simulations indicate that the portion of total rainfall that becomes runoff will be reduced by 46 percent if 45 percent of current cropland is converted to wetland. The addition of wetlands reduced peak runoff in most of the simulations, but the reduction varied with antecedent soil moisture, the magnitude of the peak flow, and the presence of current wetlands and lakes. Reductions in the simulated total and peak runoff from the Jack Creek Basin for most of the simulated storms were greatest when additional wetlands were simulated in the North Branch Jack Creek or the Upper Jack Creek Subbasins. In the Okabena Creek Basin, reductions in simulated peak runoff for most of the storms were greatest when additional wetlands were simulated in the Lower Okabena Creek Subbasin.

Minnesota