Geology topics

Martha G. Nielsen

Publications and source records attributed to Martha G. Nielsen.

At least 19 recordsLinked to original sources

USGS2026

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↗

USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗

USGS2023

Updated estimates of water budget components for the Mississippi Embayment Region using a soil-water-balance model, 2000–2020

A Soil-Water-Balance (SWB) model for the Mississippi embayment region in Arkansas, Tennessee, Mississippi, and Louisiana was constructed and calibrated to gain insight into potential recharge patterns for the Mississippi River Valley alluvial aquifer, which has had substantial drawdown under intense pumping stress over the last several decades. An analysis of the net infiltration term from the SWB model combined with newly gathered airborne electromagnetic geophysical data on the surficial sediments in a calibrated modular three-dimensional finite-difference (MODFLOW 6) groundwater flow model of one area in the alluvial plain found that the distribution of net infiltration was significantly different from the recharge that gets to the water table through the complicated silt and clay stratigraphy of the unsaturated zone. The net infiltration of water through the rooting zone as simulated by SWB ranges from 5.7 to 12.3 inches per year in the alluvial plain part of the model domain, and is fairly evenly distributed within local areas. Recharge to the underlying aquifer is less and is much more focused in particular zones where the connectivity through the upper layers of the unsaturated zone above the water table is greater, indicating possible horizontal flow and perched water table conditions in the unsaturated zone. Runoff and net infiltration together account for 32 percent of the incoming precipitation overall and somewhat higher percentages in the alluvial plain area on an annual basis. These terms are much higher in the fall and winter than in the summer. Actual evapotranspiration accounts for between 62 and 72 percent on average of the annual precipitation but dominates all other terms in the summer months. Without irrigation, summertime net infiltration and runoff would be near zero in the crop-dominated alluvial plain area. The SWB model reproduced reported irrigation rates for corn, soybeans, rice, and cotton on an annual basis fairly well. The SWB model for the Mississippi embayment region was calibrated using more than 15,000 observations representing four parts of the calculated water budget: actual evapotranspiration, surface runoff, net infiltration, and irrigation. Using a Monte Carlo approach to determine the uncertainty in the model results stemming from the uncertainty in the model parameters used in the calibration, the uncertainty in the annual actual evapotranspiration values was around 5 percent, whereas the uncertainty in the irrigation, net infiltration, and runoff was around 20 percent.

Mississippi Embayment Region↗

USGS2023

Documentation of a pilot workflow for reanalyzing the U.S. Geological Survey principal aquifers datasets and prototype principal aquifer version 2 dataset for three aquifer systems

A pilot workflow to refine the principal aquifers of the United States as defined in the Ground Water Atlas of the United States and create a new version of the principal aquifers (referred to as “version 2”) is documented in this report. The workflow incorporates decision points for creating finer scale spatial data for the principal aquifers and refining the original principal aquifer definitions if warranted. This workflow was applied to four principal aquifers in the upper Midwest region of the United States that were not previously refined as part of a U.S. Geological Survey regional groundwater availability study: the Cambrian-Ordovician aquifer system, the Jacobsville aquifer, the Silurian-Devonian aquifer, and the upper carbonate aquifer. The refinement resulted in the consolidation of two of these aquifers (the Silurian-Devonian and upper carbonate aquifers), an expansion of the Jacobsville aquifer into a larger newly defined Midcontinent Rift sandstone aquifers unit, and a slight refinement of the Cambrian-Ordovician aquifer system to exclude Precambrian units. The U.S. Geological Survey State Geologic Map Compilation geodatabase provided the base data used in the refined version 2 dataset, which are published in an accompanying U.S. Geological Survey data release as a prototype version 2 shapefile and include attributes describing the aquifer, data lineage, and source of the originally defined principal aquifer.

Illinois, Indiana, Iowa, Kentucky, Michigan, Minne↗

USGS2023

Simulation of groundwater flow at the former Badger Army Ammunition Plant, Sauk County, Wisconsin

To help support remedial efforts at the former Badger Army Ammunition Plant the U.S. Geological Survey built and calibrated a transient groundwater flow model using the Newton Raphson formulation (MODFLOW–NWT) of the U.S. Geological Survey’s modular three-dimensional finite-difference code. The model simulates the groundwater flow system at the site from 1984 to 2020. The former Badger Army Ammunition Plant is a 7,275-acre site in Sauk County, Wisconsin. The plant produced smokeless gunpower and solid rocket propellent as munitions components. Peak production periods were during World War II, the Korean War, and the Vietnam War. Subsequent groundwater contamination investigations have found four plumes at the site. A health risk assessment identified at least one contaminant of concern for human health risk present in three of the plumes: the propellant burning ground plume, the deterrent burning ground plume, and the central plume. A cooperative study began between the U.S. Army Environmental Command and U.S. Geological Survey to better understand the groundwater flow system at the former Badger Army Ammunition Plant. Field data, including aquifer tests, streamflow measurements, continuous groundwater elevations, and groundwater gradients with the Wisconsin River were collected and used to inform and calibrate the groundwater flow model. The model was used to assess the variability of the groundwater system over the study period, the components of the groundwater budget, and groundwater flow directions from identified source areas towards the Wisconsin River. Model performance assessment focused on using particle tracking to compare groundwater flowpaths that originate in the contaminant source areas to the observed plume footprints. This focus on plume behavior geometry should help constrain the advective component of a future groundwater transport model of the site.

Wisconsin↗

USGS2021

Initial estimates of net infiltration and irrigation from a soil-water-balance model of the Mississippi Embayment Regional Aquifer Study Area

The Mississippi embayment encompasses about 100,000 square miles and covers parts of eight States. In 2016, the U.S. Geological Survey began updating previous work for a part of the embayment known as the Mississippi Alluvial Plain to support informed water use and agricultural policy in the region. Groundwater, water use, economic, and other related models are being combined with field surveys and observations to create a quantitative framework for evaluating regional groundwater withdrawals and their effects on long-term water availability in the Mississippi Alluvial Plain. As part of this effort, the U.S. Geological Survey’s Soil-Water-Balance code (version 2.0) is being used to model potential groundwater recharge and irrigation water use, as necessary inputs to the long-term groundwater modeling efforts. The Soil-Water-Balance code is designed to estimate the distribution and timing of net infiltration leaving the root zone. Soil-Water-Balance makes use of gridded datasets of elevation, soils, land use (including specific crop types), and daily weather datasets to calculate other components of the root-zone water balance, including soil moisture, reference, actual evapotranspiration, snowfall, snowmelt, and canopy interception. Parameters on plant height and growing-season water needs are used to estimate crop-water demand and potential irrigation water use. This report documents the initial construction, calibration, and application of a Soil-Water-Balance model of the Mississippi Embayment Regional Aquifer Study area for simulations running from 1915 to 2017. Further refinements of the model calibration for an expanded model area are planned.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗

USGS2020

Estimated groundwater withdrawals from principal aquifers in the United States, 2015

In 2015, about 84,600 million gallons per day (Mgal/d) of groundwater were withdrawn in the United States for various uses including public supply, self-supplied domestic, industrial, mining, thermoelectric power, aquaculture, livestock, and irrigation. Of this total, about 94 percent (79,200 Mgal/d) was withdrawn from principal aquifers, which are defined as regionally extensive aquifers or aquifer systems that have the potential to be used as sources of water of suitable quality and quantity to meet various needs. The remaining 6 percent (5,400 Mgal/d) was withdrawn from other, nonprincipal aquifers in the United States. Sixty-six principal aquifers belonging to 5 major lithologic groups have been identified and delineated in the United States, including Puerto Rico and the U.S. Virgin Islands. Of the water withdrawn from principal aquifers in 2015, 81 percent (63,900 Mgal/d) was from the unconsolidated and semiconsolidated sand and gravel lithologic group, 7.1 percent (5,630 Mgal/d) was from the igneous and metamorphic-rock lithologic group, 6.8 percent (5,360 Mgal/d) was from the carbonate-rock lithologic group, 3.4 percent (2,680 Mgal/d) was from the sandstone lithologic group, and 2.2 percent (1,710 Mgal/d) was from the sandstone and carbonate-rock lithologic group. The most heavily pumped of the 24 principal aquifers and aquifer systems within the unconsolidated and semiconsolidated sand and gravel lithologic group were the High Plains aquifer (12,300 Mgal/d), Mississippi River Valley alluvial aquifer (12,100 Mgal/d), Central Valley aquifer system (11,100 Mgal/d), and Basin and Range basin-fill aquifers (7,390 Mgal/d). Withdrawals for irrigation were 48,100 Mgal/d and accounted for 75 percent of the total withdrawals from this lithologic group. Although unconsolidated sand and gravel aquifers are widely distributed and were used as sources of water in all States except Hawaii and the U.S. Virgin Islands, 56 percent of the total withdrawn from unconsolidated and semiconsolidated sand and gravel aquifers was in just four States: California (15,600 Mgal/d), Arkansas (9,560 Mgal/d), Nebraska (5,570 Mgal/d), and Texas (4,830 Mgal/d). The most heavily pumped of the seven principal aquifers within the igneous and metamorphic-rock lithologic group were the Snake River Plain (2,930 Mgal/d) and Columbia Plateau basaltic-rock aquifers (1,080 Mgal/d), which are located in the northwestern United States and together accounted for 71 percent of the water withdrawn from this lithologic group. Withdrawals for irrigation were 4,190 Mgal/d and accounted for more than 74 percent of the total withdrawals from this lithologic group. Seventy-eight percent of the withdrawals from igneous and metamorphic-rock aquifers were in three States: Idaho (3,230 Mgal/d), Washington (614 Mgal/d), and Oregon (528 Mgal/d). The most heavily pumped of the 15 principal aquifers and aquifer systems within the carbonate-rock lithologic group were the Floridan aquifer system (3,180 Mgal/d) and the Biscayne aquifer (679 Mgal/d), which are in the southeastern United States and together accounted for almost 72 percent of the withdrawals from this lithologic group. Withdrawals for public supply (2,440 Mgal/d) and irrigation (1,610 Mgal/d) together accounted for almost 76 percent of the total withdrawals from this lithologic group. Although water was withdrawn from carbonate-rock aquifers in 35 States, 71 percent of the total withdrawn was in Florida (3,020 Mgal/d) and Georgia (785 Mgal/d). The most heavily pumped of the 15 principal aquifers within the sandstone lithologic group was the Cambrian-Ordovician aquifer system (921 Mgal/d), which is in the north-central United States and accounted for 34 percent of the water withdrawn from this lithologic group. Withdrawals for public supply were 1,030 Mgal/d and accounted for 38 percent of the total withdrawals from this lithologic group. Although sandstone aquifers were used as sources of water in 32 States, 45 percent of the total withdrawn from sandstone aquifers was in five States: Minnesota (321 Mgal/d), Wisconsin (319 Mgal/d), Kansas (193 Mgal/d), Illinois (187 Mgal/d), and Pennsylvania (179 Mgal/d). The most heavily pumped of the five principal aquifers and aquifer systems within the sandstone and carbonate-rock lithologic group were the Edwards-Trinity aquifer system (661 Mgal/d) in the south-central United States and the Valley and Ridge aquifers (551 Mgal/d) of the eastern United States, which together accounted for 71 percent of total withdrawals from this lithologic group. Withdrawals from sandstone and carbonate-rock aquifers for public-supply (713 Mgal/d), irrigation (469 Mgal/d), and self-supplied domestic (253 Mgal/d) uses accounted for about 84 percent of the total withdrawals from this lithologic group. Although water was withdrawn from sandstone and carbonate-rock aquifers in 25 States, 65 percent of the total withdrawn was in Texas (651 Mgal/d), Pennsylvania (238 Mgal/d), and Florida (223 Mgal/d).

Circular↗

USGS2019

Groundwater recharge estimates for Maine using a Soil-Water-Balance model—25-year average, range, and uncertainty, 1991 to 2015

To address the lack of information on the spatial and temporal variability of recharge to groundwater systems in Maine, a study was initiated in cooperation with the Maine Geological Survey to use the U.S. Geological Survey Soil-Water-Balance model to evaluate annual average potential recharge across the State over a 25-year period from 1991 to 2015. The Maine Soil-Water-Balance model was calibrated using annual observations of recharge, runoff, and evapotranspiration for 32 calibration watersheds in the State during 2001–12 (902 total observations). Observations of recharge, runoff, and evapotranspiration were developed for each watershed to reduce the possibility of nonunique combinations of model parameters during the calibration. The Maine Soil-Water-Balance model was run using an optional evapotranspiration calculation method that provides more control for calibration than the standard method. The model was calibrated using the Parameter ESTimation software suite. The overall mean model error (average of all annual residuals for recharge, runoff, and precipitation) was 0.39 inch. The mean of the absolute value of the residuals, or the mean absolute error, was 2.32 inches. The root mean squared error for the calibrated model overall was 3.14 inches. Statistical tests indicated that the model residuals are normally distributed. To determine the potential uncertainty in the median annual potential recharge that results from uncertainty in the parameters as they relate to information contained in the observations, 300 alternate model realizations were run, and the standard deviation of the median potential recharge value at every pixel was calculated. Simulated 25-year median potential recharge across the State is widely variable; this variability closely follows patterns of precipitation, with additional variability contributed by the patchwork nature of the combinations of land-use class and hydrologic soil group inputs, and distribution of available water capacity in the soil across the State. Overall, the 25-year median annual potential recharge across the State is 7.5 inches, ranging from a low of about 5 inches to over 30 inches. The statewide range in the 25-year minimum values is from just over 2 inches to just over 20 inches. The statewide range in the 25-year maximum potential recharge is between 15 and 48 inches per year. The model areas with the highest simulated median potential recharge include areas underlain by type A soils (sandy and well drained), particularly those that also have land uses with low or little vegetation (blueberry barrens, developed, open space, scrub/shrub, and cropland, for example). The potential recharge values for these areas are similar to previously published values for comparable soil types. The 25-year average potential recharge grids were compared to recharge evaluated through groundwater-flow models or other methods in four hydrogeologic settings at six study areas in the State. A key factor in the ability of the Soil-Water-Balance model to reproduce the earlier study results was whether the available water-capacity data were an appropriate match for the hydrologic soil groups. The Maine Soil-Water-Balance model does a good job in representing an accurate potential recharge under circumstances where the surficial mapped soils extend below the surface to the water-table aquifer and where the available water-capacity data are in an appropriate range for the hydrologic soil group. One hydrogeologic setting that was challenging for the model was where a silt and clay layer was below a shallow soil unit that did not have available water-capacity data that were appropriate for the hydrologic soil group. In these cases, typically the available water-capacity data were very low, not accounting for the impedance of water flow provided by the underlying soil. The model also does not simulate well areas where bedrock surfaces are above the water table but below the plant rooting zone. The data products accompanying this report are intended to be used to provide first-cut estimates of recharge for geographic areas no smaller than the smallest watersheds used in the calibration of the model—or about 1.5 square miles. It is recommended that the grids are used to calculate an area-wide average potential recharge for any given area of study, and an uncertainty around the mean should be calculated from the standard deviation grid at the same time.

Maine↗

USGS2018

The Connecticut Streamflow and Sustainable Water Use Estimator—A decision-support tool to estimate water availability at ungaged stream locations in Connecticut

Freshwater streams in Connecticut are subject to many competing demands, including public water supply; agricultural, commercial, and industrial water use; and ecosystem and habitat needs. In recent years, drought has further stressed Connecticut’s water resources. To sustainably allocate and manage water resources among these competing uses, Federal, State, and local water-resource managers require data and modeling tools to estimate the water availability at a variety of temporal and spatial scales for planning purposes. The Connecticut Streamflow and Sustainable Water Use Estimator (CT SSWUE), developed by the U.S. Geological Survey in cooperation with the Connecticut Department of Energy and Environmental Protection, is a decision-support tool for estimating daily unaltered streamflow and sustainable water use at ungaged sites in Connecticut. The CT SSWUE estimates unaltered daily mean streamflow and water-use-adjusted streamflow for the period from October 1, 1960, to September 30, 2015, and the monthly sustainable net withdrawal at ungaged sites in Connecticut. Unaltered streamflow is the estimated daily mean streamflow in a drainage basin in the absence of any water withdrawals or wastewater discharges and with minimal human development. Sustainable net withdrawal is the maximum net withdrawal (withdrawal minus wastewater discharges) that can be drawn from a basin without critically depleting the water available through natural streamflow patterns. Sustainable net withdrawal is defined for this study as the difference between the unaltered daily mean streamflow and a user-defined target minimum streamflow. Weighted least squares and Tobit regression techniques were used to develop equations for estimating streamflow at ungaged sites at 19 streamflow quantiles with exceedance probabilities ranging from 0.005 to 99.995 percent. Regressions were based on streamflow quantiles and basin characteristics from 36 reference streamgages in and around Connecticut. Four basin characteristics—drainage area, mean of the soil permeability, mean of the average annual precipitation, and ratio of the length of streams that overlay sand and gravel deposits to the total length of streams in the basin—are used as explanatory variables in the equations. At an ungaged site, interpolation between the streamflow quantiles estimated from the regression equations produces a continuous flow-duration curve. A time series of daily mean streamflow at an ungaged site is then estimated by assuming that for each day, the streamflow quantile occurs on the same date at both a reference streamgage and the ungaged site. In a remove-one cross validation, estimated unaltered daily mean streamflow agreed well with observed values at reference streamgages, with a few exceptions. Nash Sutcliffe efficiency ranged from −0.43 to 0.97 with a median value of 0.88. The normalized root-mean-square error ranged from 16.6 to 120.4 percent with a median value of 34.5 percent. An empirical method for estimating 95-percent prediction intervals for unaltered daily and monthly mean streamflow was developed and tested by using the cross-validation data. Prediction intervals for unaltered daily mean streamflow at the cross-validation reference streamgages performed well in most cases. Gaged streamflow values from the cross-validation data fell within the prediction intervals a median 96.6 percent of the time for daily mean time series and 93.9 percent of the time for monthly mean time series. The CT SSWUE computes water-use-adjusted streamflow using spatially referenced water-use information provided by the Connecticut Department of Energy and Environmental Protection. Available water-use information included permitted and registered water withdrawals and permitted wastewater discharges during 1998 to 2015 for the Thames River Basin and central coastal drainage basins. Water-use information was incorporated into the U.S. Geological Survey StreamStats web application for Connecticut and can be used for computing water-use-adjusted streamflow and sustainable net withdrawal at selected points of interest. Altered daily streamflow is computed by applying average daily withdrawals and wastewater discharges to the water balance equation. Average daily surface water withdrawals and wastewater discharges are applied directly to the daily water balance equation. Time-lagged alterations on streamflow from groundwater withdrawals or wastewater discharges are estimated by using a response-coefficient method developed from results of previously published, calibrated groundwater models.

Connecticut↗

USGS2017

Groundwater-level trends in the U.S. glacial aquifer system, 1964-2013

The glacial aquifer system in the United States is a major source of water supply but previous work on historical groundwater trends across the system is lacking. Trends in annual minimum, mean, and maximum groundwater levels for 205 monitoring wells were analyzed across three regions of the system (East, Central, West Central) for four time periods: 1964-2013, 1974-2013, 1984-2013, and 1994-2013. Trends were computed separately for wells in the glacial aquifer system with low potential for human influence on groundwater levels and ones with high potential influence from activities such as groundwater pumping. Generally there were more wells with significantly increasing groundwater levels (levels closer to ground surface) than wells with significantly decreasing levels. The highest numbers of significant increases for all four time periods were with annual minimum and/or mean levels. There were many more wells with significant increases from 1964 to 2013 than from more recent periods, consistent with low precipitation in the 1960s. Overall there were low numbers of wells with significantly decreasing trends regardless of time period considered; the highest number of these were generally for annual minimum groundwater levels at wells with likely human influence. There were substantial differences in the number of wells with significant groundwater-level trends over time, depending on whether the historical time series are assumed to be independent, have short-term persistence, or have long-term persistence. Mean annual groundwater levels have significant lag-one-year autocorrelation at 26.0% of wells in the East region, 65.4% of wells in the Central region, and 100% of wells in the West Central region. Annual precipitation across the glacial aquifer system, on the other hand, has significant autocorrelation at only 5.5% of stations, about the percentage expected due to chance.

Journal of Hydrology↗

USGS2015

Simulation of groundwater flow and streamflow depletion in the Branch Brook, Merriland River, and parts of the Mousam River watersheds in southern Maine

Watersheds of three streams, the Mousam River, Branch Brook, and Merriland River in southeastern Maine were investigated from 2010 through 2013 under a cooperative project between the U.S. Geological Survey and the Maine Geological Survey. The Branch Brook watershed previously had been deemed “at risk” by the Maine Geological Survey because of the proportionally large water withdrawals compared to estimates of the in-stream flow requirements for habitat protection. The primary groundwater withdrawals in the study area include a water-supply well in the headwaters of the system and three water-supply wells in the coastal plain near the downstream end of the system. A steady-state groundwater flow model was used to understand the movement of water within the system, to evaluate the water budget and the effect of groundwater withdrawals on streamflows, and to understand streamflow depletion in relation to the State of Maine’s requirements to maintain in-stream flows for habitat protection. Delineation of the simulated groundwater divides compared to the surface-water divides suggests that the groundwater divides in the headwater areas do not exactly correspond to the surface-water divides. Under both pumping and non-pumping conditions, groundwater flows from the headwaters of the Branch Brook watershed into the Mousam River watershed. Pumping in the Mousam River watershed captures a small amount of groundwater from the Branch Brook basin. The cumulative effect of groundwater withdrawals on base flows in two rivers in the study area (Branch Brook and the Merriland River) was evaluated using the groundwater flow model. Streamflow depletion in the headwaters of Branch Brook was 0.12 cubic feet per second (ft 3 /s) for the steady-state simulation, or about 10 percent of the average base flow at that location. Downstream on Branch Brook, the total streamflow depletion from all the wells was 0.59 ft 3 /s, or 3 percent of the average base flow at that location. In the Merriland River downstream from the Merriland River well, the total amount of streamflow depletion was 0.6 ft 3 /s, or about 7 percent of the average base flow. The groundwater model was used to evaluate several different scenarios that could affect streamflow and groundwater discharging to the rivers and streams in the study area. The scenarios were (1) no pumping from the water-supply wells; (2) current pumping from the water-supply wells, but simulated drought conditions (25 percent reduction in recharge); (3) current recharge, but with increased pumping from the large water-supply wells; and (4) drought conditions and increased pumping combined. Simulations of increased pumping in the water-supply wells resulted in streamflow depletion in the headwaters of Branch Brook increasing to 16 percent of the headwater base flow. Simulated increases in the pumping in the coastal plain wells increased the amount of streamflow depletion to 6 percent of the flow in Branch Brook and to 8 percent of the flow in the Merriland River. The additional stress of a drought imposed on the model (25 percent less recharge) had a substantial impact on streamflows, as expected. If the simulated drought occurred simultaneously with an increase in pumping, the base flows would be reduced 48 percent in the headwaters of Branch Brook, compared to the no-pumping scenario. Downstream in Branch Brook, the total reduction in flow would be 29 percent of the simulated base flows in the no-pumping scenario, and in the Merriland River, the reduction would be 33 percent of the base flows in the no-pumping scenario. The study evaluated two different methods of calculating in-stream flow requirements for Branch Brook and the Merriland River—a set of statewide equations used to calculate monthly median flows and the MOVE.1 record-extension technique used on site-specific streamflow measurements. The August median in-stream flow requirement in the Merriland River was calculated as 7.18 ft 3 /s using the statewide equations but was 3.07 ft 3 /s using the MOVE.1 analysis. In Branch Brook, the August median in-stream flow requirements were calculated as 20.3 ft 3 /s using the statewide equations and 11.8 ft 3 /s using the MOVE.1 analysis. In each case, using site-specific data yields an estimate of in-stream flow that is much lower than an estimate the statewide equations provide.

Maine↗

USGS2014

Flow and sorption controls of groundwater arsenic in individual boreholes from bedrock aquifers in central Maine, USA

To understand the hydrogeochemical processes regulating well water arsenic (As) evolution in fractured bedrock aquifers, three domestic wells with [As] up to 478 μg/L are investigated in central Maine. Geophysical logging reveals that fractures near the borehole bottom contribute 70-100% of flow. Borehole and fracture water samples from various depths show significant proportions of As (up to 69%) and Fe (93-99%) in particulates (>0.45 μm). These particulates and those settled after a 16-day batch experiment contain 560-13,000 g/kg of As and 14-35% weight/weight of Fe. As/Fe ratios (2.5-20 mmol/mol) and As partitioning ratios (adsorbed/dissolved [As], 20,000-100,000 L/kg) suggest that As is sorbed onto amorphous hydrous ferric oxides. Newly drilled cores also show enrichment of As (up to 1300 mg/kg) sorbed onto secondary iron minerals on the fracture surfaces. Pumping at high flow rates induces large decreases in particulate As and Fe, a moderate increase in dissolved [As] and As(III)/As ratio, while little change in major ion chemistry. The δD and δ 18 O are similar for the borehole and fracture waters, suggesting a same source of recharge from atmospheric precipitation. Results support a conceptual model invoking flow and sorption controls on groundwater [As] in fractured bedrock aquifers whereby oxygen infiltration promotes the oxidation of As-bearing sulfides at shallower depths in the oxic portion of the flow path releasing As and Fe; followed by Fe oxidation to form Fe oxyhydroxide particulates, which are transported in fractures and sorb As along the flow path until intercepted by boreholes. In the anoxic portions of the flow path, reductive dissolution of As-sorbed iron particulates could re-mobilize As. For exposure assessment, we recommend sampling of groundwater without filtration to obtain total As concentration in groundwater.

Maine↗

USGS2013

High-water marks from tropical storm Irene for selected river reaches in northwestern Massachusetts, August 2011

A Presidential Disaster Declaration was issued for Massachusetts, with a focus on the northwestern counties, following flooding from tropical storm Irene on August 28–29, 2011. Three to 10 inches of rain fell during the storm on soils that were susceptible to flash flooding because of wet antecedent conditions. The gage height at one U.S. Geological Survey (USGS) streamgage rose nearly 20 feet in less than 4 hours because of the combination of saturated soils and intense rainfall. Eight of 16 USGS long-term streamgages in western Massachusetts set new peaks of record on August 28 or 29, 2011. To document the historic water levels of the streamflows from tropical storm Irene, the USGS identified, flagged, and surveyed 323 high-water marks in the Deerfield and Hudson- Hoosic River basins in northwestern Massachusetts. Areas targeted for high-water marks were generally upstream and downstream from structures along selected river reaches. Elevations from high-water marks can be used to confirm peak river stages or help compute peak streamflows, to calibrate hydraulic models, or to update flood-inundation and recovery maps. For areas in western Massachusetts that flooded as a result of tropical storm Irene, high-water marks surveyed for this study have helped to confirm or determine instantaneous peak river gage heights at several USGS streamgages.

Massachusetts↗

USGS2013

Changes in nitrogen loading to the Northeast Creek Estuary, Bar Harbor, Maine, 2000 to 2010

Since 1999, the U.S. Geological Survey and the National Park Service have been monitoring land use and nitrogen loading in a 26.3-square-kilometer (10-square-mile) estuarine watershed at Acadia National Park, Mount Desert Island, Maine. The initial study linking land use and nitrogen loads entering the Northeast Creek estuary was completed in 2000, and findings were used to develop simulations of nitrogen loading to the estuary, thereby helping to inform local land-use planning decisions. At that time, the amount of nitrogen entering the Northeast Creek estuary was relatively small, and no evidence of nutrient-related degradation was observed in the Ruppia -dominated estuarine ecosystem. A new round of water-quality monitoring and streamflow measurements was conducted to determine nitrogen loading from 2008 to 2011 as a means to evaluate the effects of increased rural residential housing within the watershed since 2000. On the basis of a 2.6-percent increase in residential-housing land use in the watershed from 2000 to 2010, simulations of nitrogen export predicted a 7-percent increase in nitrogen loading to Northeast Creek. The measurement-based loads estimated for the Northeast Creek tributaries, however, increased much more than predicted, from 1.89 kilograms per hectare per year (kg/ha/yr) in 2000 to 3.12 kg/ha/yr in the time period centered on 2010—a 66-percent increase. This increase is likely primarily a result of the prevalence of much wetter conditions during the 2008–11 sampling period than during the earlier sampling period. In addition to increasing the physical transport of nitrogen in the watershed, wet climatic conditions have been shown in other studies to increase the rates of biotic and abiotic processes that control nitrogen export from northern-latitude forested watersheds. The new loading estimates, however, also support the possibility that some portion of the increase in nitrogen loading results from the observed land-use changes, and that the increase in residential housing has, in fact, contributed to the observed increase in nitrogen loading.

Maine↗

USGS2013

Estimates of future inundation of salt marshes in response to sea-level rise in and around Acadia National Park, Maine

Salt marshes are ecosystems that provide many important ecological functions in the Gulf of Maine. The U.S. Geological Survey investigated salt marshes in and around Acadia National Park from Penobscot Bay to the Schoodic Peninsula to map the potential for landward migration of marshes using a static inundation model of a sea-level rise scenario of 60 centimeters (cm; 2 feet). The resulting inundation contours can be used by resource managers to proactively adapt to sea-level rise by identifying and targeting low-lying coastal areas adjacent to salt marshes for conservation or further investigation, and to identify risks to infrastructure in the coastal zone. For this study, the mapping of static inundation was based on digital elevation models derived from light detection and ranging (LiDAR) topographic data collected in October 2010. Land-surveyed control points were used to evaluate the accuracy of the LiDAR data in the study area, yielding a root mean square error of 11.3 cm. An independent accuracy assessment of the LiDAR data specific to salt-marsh land surfaces indicated a root mean square error of 13.3 cm and 95-percent confidence interval of ± 26.0 cm. LiDAR-derived digital elevation models and digital color aerial photography, taken during low tide conditions in 2008, with a pixel resolution of 0.5 meters, were used to identify the highest elevation of the land surface at each salt marsh in the study area. Inundation contours for 60-cm of sea-level rise were delineated above the highest marsh elevation for each marsh. Confidence interval contours (95-percent,± 26.0 cm) were delineated above and below the 60-cm inundation contours, and artificial structures, such as roads and bridges, that may present barriers to salt-marsh migration were mapped. This study delineated 114 salt marshes totaling 340 hectares (ha), ranging in size from 0.11 ha (marshes less than 0.2 ha were mapped only if they were on Acadia National Park property) to 52 ha, with a median size of 1.0 ha. Inundation contours were mapped at 110 salt marshes. Approximately 350 ha of low-lying upland areas adjacent to these marshes will be inundated with 60 cm of sea-level rise. Many of these areas are currently freshwater wetlands. There are potential barriers to marsh migration at 27 of the 114 marshes. Although only 23 percent of the salt marshes in the study are on ANP property, about half of the upland areas that will be inundated are within ANP; most of the predicted inundated uplands (approximately 170 ha) include freshwater wetlands in the Northeast Creek and Bass Harbor Marsh areas. Most of the salt marshes analyzed do not have a significant amount of upland area available for migration. Seventy-five percent of the salt marshes have 20 meters or less of adjacent upland that would be inundated along most of their edges. All inundation contours, salt marsh locations, potential barriers, and survey data are stored in geospatial files for use in a geographic information system and are a part of this report.

Maine↗

USGS2012

Simulation of groundwater conditions and streamflow depletion to evaluate water availability in a Freeport, Maine, watershed

In order to evaluate water availability in the State of Maine, the U.S. Geological Survey (USGS) and the Maine Geological Survey began a cooperative investigation to provide the first rigorous evaluation of watersheds deemed "at risk" because of the combination of instream flow requirements and proportionally large water withdrawals. The study area for this investigation includes the Harvey and Merrill Brook watersheds and the Freeport aquifer in the towns of Freeport, Pownal, and Yarmouth, Maine. A numerical groundwater- flow model was used to evaluate groundwater withdrawals, groundwater-surface-water interactions, and the effect of water-management practices on streamflow. The water budget illustrates the effect that groundwater withdrawals have on streamflow and the movement of water within the system. Streamflow measurements were made following standard USGS techniques, from May through September 2009 at one site in the Merrill Brook watershed and four sites in the Harvey Brook watershed. A record-extension technique was applied to estimate long-term monthly streamflows at each of the five sites. The conceptual model of the groundwater system consists of a deep, confined aquifer (the Freeport aquifer) in a buried valley that trends through the middle of the study area, covered by a discontinuous confining unit, and topped by a thin upper saturated zone that is a mixture of sandy units, till, and weathered clay. Harvey and Merrill Brooks flow southward through the study area, and receive groundwater discharge from the upper saturated zone and from the deep aquifer through previously unknown discontinuities in the confining unit. The Freeport aquifer gets most of its recharge from local seepage around the edges of the confining unit, the remainder is received as inflow from the north within the buried valley. Groundwater withdrawals from the Freeport aquifer in the study area were obtained from the local water utility and estimated for other categories. Overall, the public-supply withdrawals (105.5 million gallons per year (Mgal/yr)) were much greater than those for any other category, being almost 7 times greater than all domestic well withdrawals (15.3 Mgal/yr). Industrial withdrawals in the study area (2.0 Mgal/yr) are mostly by a company that withdraws from an aquifer at the edge of the Merrill Brook watershed. Commercial withdrawals are very small (1.0 Mgal/yr), and no irrigation or other agricultural withdrawals were identified in this study area. A three-dimensional, steady-state groundwater-flow model was developed to evaluate stream-aquifer interactions and streamflow depletion from pumping, to help refine the conceptual model, and to predict changes in streamflow resulting from changes in pumping and recharge. Groundwater levels and flow in the Freeport aquifer study area were simulated with the three-dimensional, finite-difference groundwater-flow modeling code, MODFLOW-2005. Study area hydrology was simulated with a 3-layer model, under steady-state conditions. The groundwater model was used to evaluate changes that could occur in the water budgets of three parts of the local hydrologic system (the Harvey Brook watershed, the Merrill Brook watershed, and the buried aquifer from which pumping occurs) under several different climatic and pumping scenarios. The scenarios were (1) no pumping well withdrawals; (2) current (2009) pumping, but simulated drought conditions (20-percent reduction in recharge); (3) current (2009) recharge, but a 50-percent increase in pumping well withdrawals for public supply; and (4) drought conditions and increased pumping combined. In simulated drought situations, the overall recharge to the buried valley is about 15 percent less and the total amount of streamflow in the model area is reduced by about 19 percent. Without pumping, infiltration to the buried valley aquifer around the confining unit decreased by a small amount (0.05 million gallons per day (Mgal/d)), and discharge to the streams increased by about 8 percent (0.3 Mgal/d). A 50-percent increase in pumping resulted in a simulated decrease in streamflow discharge of about 4 percent (0.14 Mgal/d). Streamflow depletion in Harvey Brook was evaluated by use of the numerical groundwater-flow model and an analytical model. The analytical model estimated negligible depletion from Harvey Brook under current (2009) pumping conditions, whereas the numerical model estimated that flow to Harvey Brook decreased 0.38 cubic feet per second (ft 3 /s) because of the pumping well withdrawals. A sensitivity analysis of the analytical model method showed that conducting a cursory evaluation using an analytical model of streamflow depletion using available information may result in a very wide range in results, depending on how well the hydraulic conductivity variables and aquifer geometry of the system are known, and how well the aquifer fits the assumptions of the model. Using the analytical model to evaluate the streamflow depletion with an incomplete understanding of the hydrologic system gave results that seem unlikely to reflect actual streamflow depletion in the Freeport aquifer study area. In contrast, the groundwater-flow model was a more robust method of evaluating the amount of streamflow depletion that results from withdrawals in the Freeport aquifer, and could be used to evaluate streamflow depletion in both streams. Simulations of streamflow without pumping for each measurement site were compared to the calibratedmodel streamflow (with pumping), the difference in the total being streamflow depletion. Simulations without pumping resulted in a simulated increase in the steady-state flow rate of 0.38 ft 3 /s in Harvey Brook and 0.01 ft 3 /s in Merrill Brook. This translates into a streamflow-depletion amount equal to about 8.5 percent of the steady-state base flow in Harvey Brook, and an unmeasurable amount of depletion in Merrill Brook. If pumping was increased by 50 percent and recharge reduced by 20 percent, the amount of streamflow depletion in Harvey Brook could reach 1.41 ft 3 /s.

Maine↗

USGS2011

Inventory and protection of salt marshes from risks of sea-level rise at Acadia National Park, Maine

Recent U.S. Geological Survey (USGS) climate studies in the northeastern United States have shown substantial evidence of climate-related changes during the last 100 years, including earlier snowmelt runoff, decreasing occurrence of river ice, and decreasing winter snowpack. These studies related to climate change are being expanded to include investigation of coastal wetlands that might be at risk from sealevel rise. Coastal wetlands, particularly salt marshes, are important ecosystems that provide wildlife nursery and breeding habitat, migratory bird habitat, water quality enhancement, and shoreline erosion control. The USGS is investigating salt marshes in Acadia National Park with the goal of determining which salt marshes may be threatened by sea-level rise and which salt marshes may be able to adapt to sea-level rise by migrating into adjacent low-lying lands.

Fact Sheet↗

USGS2011

Simulation of streamflow in the Pleasant, Narraguagus, Sheepscot, and Royal Rivers, Maine, using watershed models

The U.S. Geological Survey (USGS) began a study in 2008 to investigate anticipated changes in summer streamflows and stream temperatures in four coastal Maine river basins and the potential effects of those changes on populations of endangered Atlantic salmon. To achieve this purpose, it was necessary to characterize the quantity and timing of streamflow in these rivers by developing and evaluating a distributed-parameter watershed model for a part of each river basin by using the USGS Precipitation-Runoff Modeling System (PRMS). The GIS (geographic information system) Weasel, a USGS software application, was used to delineate the four study basins and their many subbasins, and to derive parameters for their geographic features. The models were calibrated using a four-step optimization procedure in which model output was evaluated against four datasets for calibrating solar radiation, potential evapotranspiration, annual and seasonal water balances, and daily streamflows. The calibration procedure involved thousands of model runs that used the USGS software application Luca (Let us calibrate). Luca uses the Shuffled Complex Evolution (SCE) global search algorithm to calibrate the model parameters. The calibrated watershed models performed satisfactorily, in that Nash-Sutcliffe efficiency (NSE) statistic values for the calibration periods ranged from 0.59 to 0.75 (on a scale of negative infinity to 1) and NSE statistic values for the evaluation periods ranged from 0.55 to 0.73. The calibrated watershed models simulate daily streamflow at many locations in each study basin. These models enable natural resources managers to characterize the timing and amount of streamflow in order to support a variety of water-resources efforts including water-quality calculations, assessments of water use, modeling of population dynamics and migration of Atlantic salmon, modeling and assessment of habitat, and simulation of anticipated changes to streamflow and water temperature resulting from changes forecast for air temperature and precipitation.

Scientific Investigations Report↗