Geology Reports⌕ Search

Geology topics

Jared J. Trost

Publications and source records attributed to Jared J. Trost.

At least 19 recordsLinked to original sources

Simulations of Minnesota water budget components using the Soil-Water-Balance model, past (1981-2022) and future (2040-2059 and 2080-2099)

This study presents historical estimates and future projections of net infiltration (potential recharge), actual evapotranspiration (ET), surface runoff, and related water budget components across Minnesota using the USGS Soil Water Balance model, version 2 (SWB). The model was calibrated to streamflow and actual ET observations from 2000–2022, achieving an overall r² of 0.973 with a standard error of 0.002 inches per year (in/yr). Two primary sets of simulations were produced. First, historical water budget components for 1981–2022 were estimated using observed daily PRISM weather data. During this period, Minnesota’s mean annual precipitation was 28.5 in/yr. The calibrated SWB model simulated mean annual actual ET of 21.6 in/yr, mean annual surface runoff of 3.0 in/yr, and mean annual net infiltration of 4.0 in/yr for Minnesota. Second, water budget component projections were simulated for 1995–2014, 2040–59, and 2080–99 using 4 kilometer, dynamically downscaled climate models representing intermediate, high, and very high emissions scenarios. All future climate models showed warming relative to the historical period, leading to consistent responses in temperature‑sensitive components on an annual basis: less snowfall and snowmelt and higher reference ET. Seasonal shifts were also apparent, with most models projecting increases in winter precipitation, actual ET, runoff, and net infiltration, and decreases in summer and fall precipitation, actual ET, and net infiltration relative to the historical period. The accompanying data products are intended for estimating multi-year mean water budget components at spatial scales no smaller than the smallest watershed included in model calibration (about 20 square miles). These updated historical estimates and future projections provide decision relevant information for water resource managers and offer a framework for developing long range projections of water budget variables.

EarthArXiv↗

Treatability study to evaluate bioremediation of trichloroethene at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, 2020–22

Executive Summary Chlorinated solvents, including trichloroethene (TCE) and other chlorinated volatile organic compounds (cVOCs), are widespread contaminants that can be treated by bioremediation approaches that enhance anaerobic reductive dechlorination. Reductive dechlorination can be enhanced either through the addition of an electron donor (biostimulation) or the addition of a known dechlorinating culture (bioaugmentation) along with an electron donor. Although bioremediation has been applied at many TCE-contaminated groundwater sites, application in source zones at sites where residual dense nonaqueous phase liquid (DNAPL) is present is more limited. In this study, laboratory and field treatability tests were completed to evaluate the potential application of anaerobic bioremediation for a shallow groundwater plume containing TCE in a perched alluvial aquifer at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, which was on the National Priorities List as the New Brighton/Arden Hills Superfund site until 2019. In addition to the presence of residual DNAPL at the site, temporal variability in groundwater flow directions and input of oxygenated recharge were possible complicating factors for the application of enhanced anaerobic biodegradation in the shallow plume. The Site K plume extends beneath the footprint of Building 103, which was demolished in 2006, and soil excavations to a maximum depth of 6 feet (ft) below ground surface in 2014 were known to leave some deeper contaminated soil in place in the TCE source area. Groundwater treatment at the site, formalized as part of the 1997 Record of Decision, has been in operation since 1986 and consists of an extraction trench at the downgradient edge of the plume to collect groundwater, which is then pumped to an on-site air stripper. Groundwater concentrations in the plume have been relatively stable since treatment began, indicating a continued source of TCE in the aquifer. The desire for a destructive remedy that would enhance the removal of cVOCs in the aquifer at Site K and shorten the remediation timeframe led the U.S. Army to request that the U.S. Geological Survey conduct a groundwater treatability study to assess bioremediation. This report describes the U.S. Geological Survey bioremediation treatability study conducted during 2020–22, including pre-design site characterization to assist in formulating the bioremediation approach, laboratory experiments to support the design of the field pilot test, and implementation and 1-year performance monitoring results for the pilot test. Pre-design site characterization included the collection of soil cores for cVOC analysis and lithologic descriptions and the re-installment of three wells to obtain hydrologic measurements and initial groundwater chemistry. Relatively flat head gradients were measured at the site, and substantial decreases in water-level elevations occurred from spring to summer (May–July 2021). Continuous water-level monitoring indicated a rapid response to precipitation. Groundwater flow velocities were consistently less than 0.5 foot per day, and the pilot bioremediation test was therefore designed with short lateral distances (about 5 ft) between injection and individual monitoring points. Soil analyses confirmed that high volatile organic compound contamination was left in place in the source area. The highest concentrations were near or in clay at the base of the perched aquifer. Concentrations of cVOCs measured in the replaced wells were consistent with historical data and had a maximum TCE concentration of 57,700 micrograms per liter (μg/L), indicative of nearby residual DNAPL based on the general rule of observed concentrations exceeding 1 percent of solubility. The primary TCE daughter product detected was 1,2-cis-dichloroethene (cisDCE), which indicated limited reductive dechlorination in the plume. Groundwater in both the source and downgradient areas was relatively reducing during the pre-design characterization, particularly in the source area where methane concentrations greater than 400 μg/L were measured. Initial laboratory tests conducted using native aquifer microorganisms from the three replacement wells showed that anaerobic TCE biodegradation rates were low when biostimulated with the addition of sodium lactate as an electron donor, also known as a carbon donor, and resulted in the production of only cisDCE. Addition of a known dechlorinating culture, WBC-2, however, resulted in rapid biodegradation and production of ethene, verifying complete reductive dechlorination of TCE. Microcosms constructed with aquifer soil collected from the site were used to evaluate other electron donors besides lactate to support reductive dechlorination by WBC-2, including corn syrup as an alternative fast-release compound and whey, soy-based vegetable oil, and 3-D Microemulsion (Regenesis, San Clemente, California) as slow-release compounds. First-order rate constants for total organic chlorine removal in these WBC-2 amended microcosms were greatest with either lactate or vegetable oil as the donor, ranging between 0.061 and 0.047 per day or corresponding half-lives of 11–15 days. Testing of commercial products in other WBC-2-bioaugmented microcosms led to selection for the field pilot test of an emulsified vegetable oil product that also contained some sodium lactate as a fast-release donor. Delaying the addition of WBC-2 relative to the donor in the microcosms resulted in the most rapid overall biodegradation rates. The selected design for the pilot test utilized three separate test plots, each about 30-ft wide and 60-ft long: plots GS1 and GS2 in the source area of the plume and plot GS3 in the downgradient area of the plume near the excavation trench. Each test plot had one injection well, one monitoring well upgradient from the injection point, and 12 surrounding monitoring wells in a grid to capture variable groundwater flow directions. Donor injections, which included a bromide tracer, were completed in October 2021, immediately following baseline sampling, and the WBC-2 culture was injected about 40 days later, between November 30 and December 2, 2021. Performance monitoring conducted until December 2022 included hydrologic measurements and analyses of cVOCs, redox-sensitive constituents, dissolved organic carbon, bromide, volatile fatty acids, compound-specific carbon isotopes, and microbial communities. The biogeochemical data collected during the pilot tests in the three treatment plots showed that enhanced, complete reductive dechlorination of cVOCs in the groundwater was achieved in the GS1 and GS3 plots. In contrast, evidence of distribution of the injected amendments and subsequent biodegradation was limited in GS2, which was in an area of more heterogeneous soil lithology and low water table elevations. The molar composition of volatile organic compounds in the GS1 and GS3 plots was dominated by ethene in wells that were reached by the injected amendments by the end of the monitoring period. In the GS1 and GS3 plots, similar patterns were observed of cVOC concentrations decreasing to near detection levels, or below, at some wells sampled in July and October 2022, whereas ethene became dominant and indicated sustained complete reductive dechlorination. Baseline cVOC concentrations were more than a factor of 10 higher in the groundwater in the GS1 plot than in GS3, but no apparent inhibition of complete dechlorination occurred. As expected from the initial pre-design site data and the laboratory experiments, enhanced dissolution of residual DNAPL coupled to biodegradation was evident in the GS1 plot, where a marked increase in dichloroethene (DCE) above the initial baseline and upgradient TCE and DCE concentrations occurred. DCE concentrations subsequently declined where DNAPL dissolution was evident, concurrent with production of vinyl chloride and then predominantly ethene. Thus, overall biodegradation rates outpaced the DNAPL dissolution and desorption and DCE production in the source area. This success in complete degradation to predominantly ethene was achieved even in areas where the DCE concentrations reached a maximum of about 30,000 μg/L. Compound specific isotope analysis of carbon in TCE, cisDCE, trans-1,2-dichloroethene, and vinyl chloride was conducted to provide another line of evidence of the occurrence and extent of anaerobic biodegradation. Along a flow path in each plot that was affected by the injected amendments, carbon isotopes in the TCE and daughter cVOCs in the groundwater became isotopically heavier, indicating biodegradation. Enhanced biodegradation rates calculated from the field tests in GS1 and GS3 showed half-lives of 36.9–75.3 days for DCE degradation and 9.48–38.5 days for ethene production. Notably, these ethene production rates calculated from the field tests are consistent with the results of WBC-2-bioaugmented microcosms amended with either lactate or vegetable oil, which had half-lives for total organic chlorine removal that ranged from 11 to 15 days. These rates indicated rapid enhanced biodegradation, which is promising for application of a full-scale bioremediation remedy. Ultimately, however, the mass of residual or sorbed TCE in the aquifer that remains accessible for dissolution and biodegradation would likely control the time required for a full-scale bioremediation effort to achieve performance goals for TCE and cisDCE specified in the Record of Decision for Site K. The field pilot tests showed that the relatively low hydraulic head gradients and temporal changes in groundwater flow directions in the shallow aquifer would add complexity to a full-scale bioremediation effort. The radius of influence (ROI) at GS1 and GS3 (16.3 ft and 12.7 ft, respectively) were close to the design ROI of 15 ft. The estimated ROI at GS2 was about four times the design ROI, but may be less reliable at this location owing to groundwater flow direction. In addition, the low temperatures following WBC-2 injection in late November to early December 2021, in combination with the low hydraulic head gradients, were probably major factors in the delay observed before the onset of enhanced biodegradation following injection of the culture. Additional test injections could be beneficial to optimize the timing of donor and culture injections with the variable temperatures and hydraulic head in the shallow aquifer.

Minnesota↗

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii↗

Photochemical mobilization of dissolved hydrocarbon oxidation products from petroleum contaminated soil into a shallow aquifer activate human nuclear receptors

Elevated non-volatile dissolved organic carbon (NVDOC) concentrations in groundwater monitoring wells under oil-contaminated hydrophobic soils originating from a pipeline rupture at the National Crude Oil Spill & Natural Attenuation Research Site near Bemidji, MN are documented.. We hypothesized the elevated NVDOC is comprised of water-soluble photooxidation products transported from the surface to the aquifer. We use field and laboratory samples in combination with complementary analytical methods to test this hypothesis and determine the biological response to these products. Observations from optical spectroscopy and ultrahigh-resolution mass spectrometry reveal a significant correlation between the chemical composition of NVDOC leached from photochemically weathered soils and groundwater monitoring wells with high NVDOC concentrations measured in the aquifer beneath the contaminated soil. Conversely, the chemical composition from the uncontaminated soil photoleachate, matches the NVDOC observed in the uncontaminated wells. Contaminated groundwater and photodissolution leachates from contaminated soil activated biological targets indicative of xenobiotic metabolism and exhibited potential for adverse effects. Newly formed hydrocarbon oxidation products (HOPs) from fresh oil could be distinguished from those downgradient. This study illustrates another pathway for dissolved HOPs to infiltrate groundwater and potentially affect human health and the environment.

Minnesota↗

Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18

Confined (or buried) aquifers of glacial origin overlain by till confining units provide drinking water to hundreds of thousands of Minnesota residents. The sustainability of these groundwater resources is not well understood because hydraulic properties of till that control vertical groundwater fluxes (leakage) to underlying aquifers are largely unknown. The U.S. Geological Survey, Iowa State University, Minnesota Geological Survey, and Minnesota Department of Health investigated hydraulic properties and groundwater flow through till confining units using field studies and heuristic MODFLOW simulations. Till confining units in the following late-Wisconsinan stratigraphic units (with locations in parentheses) were characterized: Des Moines lobe till of the New Ulm Formation (Litchfield, Minnesota), Superior lobe till of the Cromwell and Aitkin Formations (Cromwell, Minn.), and Wadena lobe till of the Hewitt Formation (hydrogeology field camp [HFC] near Akeley, Minn.). Pre-Illinoian till of the Good Thunder formation (Olivia, Minn.) was also characterized. Hydraulic and geochemical field data were collected from sediment cores and a series of five piezometer nests. Each nest consisted of five to eight piezometers screened at short vertical intervals in hydrostratigraphic units including (if present) surficial aquifers, till confining units, confined/buried aquifers, and underlying bedrock. Till hydraulic conductivity was estimated from slug tests (horizontal [ K h ]) and constant-rate aquifer tests in the confined aquifer (vertical [ K v ]). Travel times through the till were evaluated with Darcy’s law and stable isotope concentrations. A series of heuristic MODFLOW simulations were used to evaluate groundwater fluxes through till across the range of till hydraulic properties and pumping rates observed at the field sites. The field data demonstrated variability in hydraulic properties between and within till stratigraphic units horizontally and vertically. The variability in hydraulic properties within and between sites resulted in substantial differences in groundwater flux through till. A conceptual understanding that emerges from the vertical till profiles is that they are not homogeneous hydrostratigraphic units with uniform properties; rather, each vertical sequence is a heterogeneous mixture of glacial sediment with differing abilities to transmit water. Till thicknesses varied from 60 to 166 feet, and till textures ranged from a sandy loam (Hewitt Formation, HFC site) to a silt loam/clay loam (Good Thunder formation, Olivia site). Till Kh varied by one to three orders of magnitude within each piezometer nest. Four piezometer nests had downward hydraulic gradients ranging from 0.04 to 0.56, and one nest had a slight upward hydraulic gradient of 0.02. The Cromwell, HFC, and Litchfield 1 sites were examples of “leaky” tills with high Kv (0.001 to 1.1 feet per day [ft/d]) and geometric mean Kh (0.03 to 0.07 ft/d) and extensive vertical hydraulic connectivity between the confined aquifer and the overlying till. Estimated groundwater travel times through these sites ranged from 1 to 81 years, and two of these sites had tritium throughout their till profiles. The tills at the other two sites, Olivia and Litchfield 2, were effective confining units that had low Kv (0.001 to 0.0005 ft/d) and geometric mean Kh (0.0002 to 0.004 ft/d). The till piezometers at these sites had no drawdown response to short-term (up to 10 hours for Olivia and up to 5 days for Litchfield) high-capacity pumping from the confined aquifer. Estimated groundwater travel times through the tills at these sites ranged from 165 to nearly 1,800 years, and tritium was only detected in the upper one-third of these till profiles. Across all sites, the till vertical anisotropy (ratio of K h to K v ) ranged by four orders of magnitude from 0.05 at the Cromwell nest to 70 at the Litchfield 1 nest. Stable isotopes of oxygen and hydrogen indicate that groundwater throughout all five till profiles is younger than the last glacial advance into Minnesota at about 11,000 years ago. The heuristic modeling demonstrated that, for understanding sustainability of groundwater pumping from confined aquifers, knowledge of till hydraulic properties is just as important as knowledge of aquifer hydraulic properties. Substantial differences in groundwater fluxes into and through till were observed across hydrogeologic settings representative of the field sites. Over long periods of time (hundreds of years), pumping-induced hydraulic gradients are established in confined aquifer systems and, even in low hydraulic conductivity tills, these pumping-induced hydraulic gradients increase leakage into and through till compared to ambient conditions. In conclusion, groundwater flowing vertically downward through till confining units (leakage) replenishes water pumped from confined aquifers. Till hydraulic properties, such as those presented in this report, provide important information that can be used to quantify leakage rates through till. Till hydraulic properties are variable over short distances and profoundly affect leakage rates, demonstrating the importance of site-specific till hydraulic data for evaluating the sustainability of groundwater withdrawals from confined aquifers.

Minnesota↗

Effects of a crude-oil recovery remediation system operated 1999–2003 on groundwater plumes and unsaturated-zone vapor concentrations at a crude-oil spill site near Bemidji, Minnesota

A crude-oil spill occurred in 1979 when a pipeline burst near Bemidji, Minnesota. More than 70 percent of the 1.7 million liters of spilled crude oil was removed shortly thereafter. In response to a requirement by the State regulatory agency to remove the remaining crude to a sheen in all wells, in 1998, the pipeline company installed a dual-pump recovery system at the site. This additional remediation from 1999 to 2003 resulted in removal of about 115,000 liters of crude oil, representing between 36 and 41 percent of the volume of oil (281,000–317,000 liters) estimated to be present in 1998. Effects of the 1999–2003 remediation on groundwater plumes and unsaturated-zone vapor concentrations were evaluated by the U.S. Geological Survey using several methods including measurements of oil thicknesses in wells; field water-quality properties of dissolved oxygen, specific conductance, temperature, and pH in groundwater; and vapor concentrations of methane, carbon dioxide, nitrogen, and oxygen in the unsaturated zone. Although the recovery system decreased oil thicknesses near the remediation wells, average oil thicknesses measured in all wells at the site were not reduced substantially. Dissolved oxygen and specific conductance measurements indicate that a secondary plume was created during the remediation, caused by the disposal of pumped water from the remediation wells in an upgradient infiltration gallery. This plume expanded rapidly immediately after the start of the remediation in 1999, resulting in expansion of the anoxic zone of groundwater upgradient and beneath the existing natural attenuation plume. Beginning in 2000–1, for example, specific conductance concentrations noticeably increased in many wells at the north oil pool from about 400 to more than 700 microsiemens per centimeter. The rapid expansion of the anoxic and elevated specific conductance plume indicates that the remediation contributed substantial amounts of biodegradable dissolved organic carbon to groundwater through the infiltration gallery. The trends in vapor data collected before, during, and after the remediation generally support the research hypothesis that crude-oil removal would have an insignificant effect on vapor concentrations in the unsaturated zone. Although there were some small changes in the concentration of methane, carbon dioxide, nitrogen, and oxygen in the unsaturated zone, these changes were not coincident with the beginning or cessation of the remediation and are therefore thought to be the result of other factors affecting biodegradation rates. A decrease in methane concentrations in one representative well, for example, is thought to be the result of reduced rates of biodegradation and methane production from the increasingly more weathered crude oil. Oil-phase recovery at this site was determined to be challenging and resulted in considerable volumes of mobile and entrapped oil remaining in the subsurface despite remediation efforts.

Minnesota↗

Geophysical mapping of plume discharge to surface water at a crude oil spill site: Inversion versus machine learning

The interpretation of geophysical survey results to answer hydrologic, engineering, and geologic questions is critical to diverse problems for management of water, energy, and mineral resources. Although geophysical images provide valuable qualitative insight into subsurface architecture and conditions, translating geophysical images into quantitative information (e.g., saturation, concentration, and hydraulic properties) often involves substantial nonuniqueness and uncertainty owing to the limited resolution of geophysical imaging and uncertainty in petrophysical relations. We have developed a machine-learning approach to address these challenges in the context of a field-based investigation to map zones where a hydrocarbon plume was discharging to surface water at the National Crude Oil Spill Fate and Natural Attenuation Research Site in Bemidji, Minnesota, USA. The two-step approach combines multiple types of geophysical and direct information and effectively bypasses inversion and its associated assumptions. Integrating multifrequency electromagnetic induction, ground-penetrating radar, and fluid-sampling data, we first identify discharge zones and second estimate specific conductance versus depth. Compared with conventional inversion results, the machine-learning results (1) directly address the study objectives (delineating the discharge zones); (2) better extract depth-dependent information from the data, for which sensitivity diminishes rapidly with depth; and (3) quantify the uncertainty of the predictions (i.e., discharge versus nondischarge zones), rather than the uncertainty of the geophysical estimates (i.e., the standard error of estimation for the logarithm of electrical conductivity).

Minnesota↗

UZIG research: Measurement and characterization of unsaturated zone processes under wide-ranging climates and changing conditions

Unsaturated zone properties and processes are central to understanding the interacting effects of land-use change, contamination, and hydroclimate on our ability to grow food, sustain clean water supplies, and minimize loss of life and property. Advances in unsaturated zone science are being achieved through collaborations across traditional boundaries where information from biological, physical, and chemical disciplines is combined for new insights. The Unsaturated Zone Interest Group (UZIG) is an organization that exists principally to promote multidisciplinary collaborations and the sharing of ideas, expertise, and technical assets. Here we summarize key findings from 14 papers, several of which originated from a meeting convened by UZIG in 2017 at the University of Florida in Gainesville titled “Land-Use Change, Climate Change, and Hydrologic Extremes: Unsaturated Zone Responses and Feedbacks.” This special section of Vadose Zone Journal contains multidisciplinary research in three general categories relevant to measuring and understanding unsaturated zone responses to changing land uses and climate: (i) unsaturated zone properties and processes; (ii) soil–plant–atmosphere interactions; and (iii) novel field sampling devices. A strong cross-cutting theme in these papers is the value of continuous monitoring data and ways of utilizing them to discover novel hydrologic, biologic, and pedologic information. As climatic and land-use conditions change and demands for resources and stresses on ecosystems continue to intensify, it is vital to improve our fundamental understanding of the processes at work in the unsaturated zone. Toward that goal, we discuss the need for improved ground-based unsaturated zone monitoring networks.

Vadose Zone Journal↗

A direct-push freezing core barrel for sampling unconsolidated subsurface sediments and adjacent pore fluids

Contaminants passing through the unsaturated zone can undergo changes in narrow reaction zones upon reaching saturated sediments. Understanding these reactions requires sampling of sediment together with adjacent water and microbes in a manner that preserves in situ redox conditions. Use of a basket-type core catcher for saturated, noncohesive sediments results in redistribution or loss of fluids during sample retrieval. Previously developed sample-freezing drive shoes for hollow-stem auger drilling rigs lessened fluid redistribution and retained all material that entered the core barrel in noncohesive sediment cores by freezing the base of the core with liquid CO 2 . This technology has not previously been compatible with direct-push rigs that are commonly used for contaminated site assessments. Here, we describe a freezing core barrel designed for direct-push rigs that is compatible with commercially available tool strings. The device can be used interchangeably with unsaturated-zone direct-push tool strings, enabling core collection for studies of contaminant transport and transformation spanning unsaturated to saturated profiles. In all 10 attempts during testing near Bemidji, MN, the device froze a 10- to 15-cm (4–6-in) plug that retained fluids and sediments in a 1.2-m (4-ft)-long, 5.0-cm (2.0-in)-diameter polyvinyl chloride (PVC) sleeve. Cores were collected from variably saturated sediments spanning the capillary fringe through the upper 2 m of the saturated zone in sandy glacial outwash sediments. The median recovery was 81% of the drive length, similar to a sample-freezing drive shoe developed for a wire-line piston core sampler operated with a hollow-stem auger drill rig.

Vadose Zone Journal↗

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗

Generalized hydrogeologic framework and groundwater budget for a groundwater availability study for the glacial aquifer system of the United States

The glacial aquifer system groundwater availability study seeks to quantify (1) the status of groundwater resources in the glacial aquifer system, (2) how these resources have changed over time, and (3) likely system response to future changes in anthropogenic and environmental conditions. The glacial aquifer system extends from Maine to Alaska, although the focus of this report is the part of the system in the conterminous United States east of the Rocky Mountains. The glacial sand and gravel principal aquifer is the largest source of public and self-supplied industrial supply for any principal aquifer and also is an important source for irrigation supply. Despite its importance for water supply, water levels in the glacial aquifer system are generally stable varying with climate and only locally from pumping. The hydrogeologic framework developed for this study includes the information from waterwell records and classification of material types from surficial geologic maps into likely aquifers dominated by sand and gravel deposits. Generalized groundwater budgets across the study area highlight the variation in recharge and discharge primarily driven by climate.

Glacial Aquifer System↗

Weathering of oil in a surficial aquifer

The composition of crude oil in a surficial aquifer was determined in two locations at the Bemidji, MN, spill site. The abundances of 71 individual hydrocarbons varied within 16 locations sampled. Little depletion of these hydrocarbons (relative to the pipeline oil) occurred in the first 10 years after the spill, whereas losses of 25% to 85% of the total measured hydrocarbons occurred after 30 years. The C 6‐30 n ‐alkanes, toluene, and o ‐xylene were the most depleted hydrocarbons. Some hydrocarbons, such as the n‐ C 10–24 cyclohexanes, tri‐ and tetra‐ methylbenzenes, acyclic isoprenoids, and naphthalenes were the least depleted. Benzene was detected at every sampling location 30 years after the spill. Degradation of the oil led to increases in the percent organic carbon and in the δ 13 C of the oil. Another method of determining hydrocarbon loss was by normalizing the total measured hydrocarbon concentrations to that of the most conservative analytes. This method indicated that the total measured hydrocarbons were depleted by 47% to 77% and loss of the oil mass over 30 years was 18% to 31%. Differences in hydrocarbon depletion were related to the depth of the oil in the aquifer, local topography, amount of recharge reaching the oil, availability of electron acceptors, and the presence of less permeable soils above the oil. The results from this study indicate that once crude oil has been in the subsurface for a number of years there is no longer a “starting oil concentration” that can be used to understand processes that affect its fate and the transport of hydrocarbons in groundwater.

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↗

Seasonal Variability in Vadose zone biodegradation at a crude oil pipeline rupture site

Understanding seasonal changes in natural attenuation processes is critical for evaluating source-zone longevity and informing management decisions. The seasonal variations of natural attenuation were investigated through measurements of surficial CO2 effluxes, shallow soil CO2 radiocarbon contents, subsurface gas concentrations, soil temperature, and volumetric water contents during a 2-yr period. Surficial CO2 effluxes varied seasonally, with peak values of total soil respiration (TSR) occurring in the late spring and summer. Efflux and radiocarbon data indicated that the fractional contributions of natural soil respiration (NSR) and contaminant soil respiration (CSR) to TSR varied seasonally. The NSR dominated in the spring and summer, and CSR dominated in the fall and winter. Subsurface gas concentrations also varied seasonally, with peak values of CO2 and CH4 occurring in the fall and winter. Vadose zone temperatures and subsurface CO2 concentrations revealed a correlation between contaminant respiration and temperature. A time lag of 5 to 7 mo between peak subsurface CO2 concentrations and peak surface efflux is consistent with travel-time estimates for subsurface gas migration. Periods of frozen soils coincided with depressed surface CO2 effluxes and elevated CO2 concentrations, pointing to the temporary presence of an ice layer that inhibited gas transport. Quantitative reactive transport simulations demonstrated aspects of the conceptual model developed from field measurements. Overall, results indicated that source-zone natural attenuation (SZNA) rates and gas transport processes varied seasonally and that the average annual SZNA rate estimated from periodic surface efflux measurements is 60% lower than rates determined from measurements during the summer.

Vadose Zone Journal↗

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↗

Land-cover effects on the fate and transport of surface-applied antibiotics and 17-beta-estradiol on a sandy outwash plain, Anoka County, Minnesota, 2008–09

A plot-scale field experiment on a sandy outwash plain in Anoka County in east-central Minnesota was used to investigate the fate and transport of two antibiotics, sulfamethazine (SMZ) and sulfamethoxazole (SMX), and a hormone, 17-beta-estradiol (17BE), in four land-cover types: bare soil, corn, hay, and prairie. The SMZ, SMX, and 17BE were applied to the surface of five plots of each land-cover type in May 2008 and again in April 2009. The cumulative application rate was 16.8 milligrams per square meter (mg/m 2 ) for each antibiotic and 0.6 mg/m2 for 17BE. Concentrations of each chemical in plant-tissue, soil, soil-water, and groundwater samples were determined by using enzyme-linked immunosorbent assay (ELISA) kits. Soil-water and groundwater sampling events were scheduled to capture the transport of SMZ, SMX, and 17BE during two growing seasons. Soil and plant-tissue sampling events were scheduled to identify the fate of the parent chemicals of SMZ, SMX, and 17BE in these matrices after two chemical applications. Areal concentrations (mg/m 2 ) of SMZ and SMX in soil tended to decrease in prairie plots in the 8 weeks after the second chemical application, from April 2009 to June 2009, but not in other land-cover types. During these same 8 weeks, prairie plots produced more aboveground biomass and had extracted more water from the upper 125 centimeters of the soil profile compared to all other land-cover types. Areal concentrations of SMZ and SMX in prairie plant tissue did not explain the temporal changes in areal concentrations of these chemicals in soil. The areal concentrations of SMZ and SMX in the aboveground plant tissues in June 2009 and August 2009 were much lower, generally two to three orders of magnitude, than the areal concentrations of these chemicals in soil. Pooling all treatment plot data, the median areal concentration of SMZ and SMX in plant tissues was 0.01 and 0.10 percent of the applied chemical mass compared to 22 and 12 percent in soil, respectively. Furthermore, areal concentrations of SMZ and SMX in plant-tissue samples were variable, and did not differ significantly between control and treatment plots within each land-cover type. SMZ was detected in 23 percent of soil-water samples and in 16 percent of groundwater samples collected between October 2008 and October 2009 in treatment plots, indicating that surface-applied SMZ leached below the rooting zone and reached groundwater. SMX was detected in only 1 percent of soil-water and groundwater samples during this same time period. In contrast to the antibiotics, 17BE was not reliably detected in soil samples. Additionally, ELISA-determined 17BE concentrations in plant-tissue, soil-water, and groundwater samples indicated the presence of chemicals that were not applied as part of this experiment [17BE from an external source or other chemical(s) that interfered with the 17BE ELISA kits].

Minnesota↗

Energy and conservation benefits from managed prairie biomass

Marginally productive land, such as that enrolled in the Conservation Reserve Program (CRP), may provide acreage and economic incentives for cellulosic energy production. Improving the yields from these lands will help establish a biomass producer's position in the marketplace. The effects of water and nitrogen on biomass yields were investigated in both a plot-scale experiment and a broad-scale survey of CRP lands. The plot-scale experiment demonstrated that irrigation improved mixed-species prairie biomass yields more than nitrogen fertilizer on coarse-textured, marginally productive soils. Experimental plots amended with both irrigation and moderate (but not high) nitrogen produced more biomass than other treatment combinations, but this trend was not statistically significant. The survey of biomass yields on CRP lands across four Midwestern States indicates that yields are better correlated with June rainfall than any other individual month. Applying nutrient-enriched water such as agricultural runoff could benefit prairie yields if applied at appropriate times.

Conference Paper↗