Geology Reports⌕ Search

Geology topics

John A. Colman

Publications and source records attributed to John A. Colman.

At least 19 recordsLinked to original sources
USGS2024

Evidence of nitrate attenuation in intertidal and subtidal groundwater in a subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2015–16

Nitrogen dynamics in intertidal and nearshore subtidal groundwater (subterranean estuary) adjacent to the Seacoast Shores peninsula, Falmouth, Massachusetts, were investigated during 2015–16 by the U.S. Geological Survey. The peninsula is a densely populated residential area with septic systems and cesspools that are substantial sources of nitrogen to groundwater. The study area is in the Eel River, an estuarine saltwater embayment connected to the ocean adjacent to the western shore of the peninsula, that was the subject of an earlier study by Colman and others (2018, https://doi.org/10.3133/sir20185095 ) on nitrogen transport and transformations in groundwater between onshore and offshore locations. The previous study documented the distribution of nitrate concentrations and nitrate attenuation reactions in fresh groundwater beneath the peninsula and the estuary. The current study extended those observations with more detailed sampling and analysis of shallow groundwater from wells near discharge sites beneath the estuary. The current field investigation included sampling of existing wells and installation and sampling of clusters of wells and temporary sampling points in the subterranean estuary, including (1) shallow transects 0.3 to 1.2 meters (m) deep extending from 1 to 13.5 m offshore and (2) deeper wells (from 1.83 to 4.88 m deep) extending from 4.3 to 14.3 m offshore. Measurements of hydraulic-head gradients 2–5 m below the sediment/water interface in the intertidal and nearshore subtidal zones indicated that groundwater flow generally was upwards (towards the estuary) under all tide conditions in October 2016. The magnitude of the gradient was greatest during low tide conditions, indicating that groundwater discharge likely decreased during high tides. Measurements of specific conductance in shallow groundwater in the subterranean estuary in three transects perpendicular to shore were consistent with the existence of saltwater flow cells (infiltration of overlying saline water, mixing with fresh groundwater, and discharge to the overlying saline water) in the intertidal and nearshore subtidal regions. The size of these flow cells was variable in space and time and dependent on the elevation of the tide (spring or neap). At this location in the Eel River subterranean estuary, and offshore to at least 13.5 m, offshore flow of fresh groundwater apparently prevented a deeper saltwater wedge from discharging to the surface. Nitrate concentrations in shallow groundwater (30 to 122 centimeters [cm] depth) were variable in space and time, ranging from not detectable to 600 micromoles per liter (μmol/L) (8.4 milligrams per liter as N), and were highest in June 2016 at depths from 61 to 122 cm below the sediment/water interface and from 4 to 9 m offshore. Nitrate generally was not detectable in saline shallow groundwater at 30-cm depth or at any depth from 30 to 122 cm from 10 to 13.5 m offshore. Dissolved oxygen concentrations were suboxic (less than 16 μmol/L) in 60 percent of the sampled subterranean groundwater beneath the intertidal and subtidal zones. In the remaining sites, the range of dissolved oxygen concentrations was from 18 to 272 μmol/L and the median concentration was 43 μmol/L. Evidence for microbial nitrate reduction (denitrification and possibly anammox) was provided by the distribution of the reaction product nitrogen gas (excess N 2 , or N 2MIC ), as determined from analysis of the dissolved nitrogen gas and argon gas (Ar) concentrations in groundwater samples. Excess nitrogen gas provided evidence for nitrate reduction in shallow groundwater below the subtidal and, to a lesser extent, intertidal zones adjacent to the Seacoast Shores peninsula. These zones, where evidence for nitrate reduction was detected, were in fresh and brackish groundwater near subtidal or intertidal saltwater cells where discharging fresh groundwater mixed with infiltrating saline water. Infiltrating seawater may have supplied organic carbon, one of several potential electron donors that are required for denitrification. Other potential electron donors, such as organic carbon, iron, manganese, hydrogen, methane, ammonium, elemental sulfur, or sulfide phases, may have been supplied by the estuarine sediments. Drainage from surface runoff near the shore also may have supplied organic carbon to fresh groundwater near the intertidal saltwater cell. The highest amounts of nitrate converted to excess nitrogen gas were estimated to be in the range of 230 to 430 μmol/L in nearly fresh groundwater near the subtidal saltwater cell at depths of 61 to 122 cm below the sediment/water interface and from 10 to 13.5 m offshore. Evidence of denitrification within 10 m of the shore was sparse (generally limited to less than 50 μmol/L of N 2 -N) despite the presence of high nitrate concentrations. The spatial distribution of estimated nitrate reduction in the intertidal and nearshore subtidal fresh and brackish groundwater may be related to local variability in the distribution of reactive electron donors in those zones. Variations in the amount of nitrate reduction to nitrogen gas were not clearly related to potential aqueous electron donors such as dissolved organic carbon, nor to potential reaction products such as alkalinity, but may have been controlled by combinations of aqueous and solid-phase reactants. The distribution of relatively shallow fresh groundwater containing nitrate could indicate potential nitrate discharge areas in the lower intertidal zone and uncertain locations farther offshore; however, the data did not extend all the way to the sediment/water interface or to the offshore freshwater limit. This study confirmed substantial loss of nitrate from some of the fresh and brackish groundwater in shallow subestuarine sediments prior to discharge but did not quantify how much nitrate eventually discharged to the estuary.

Massachusetts↗
USGS2021

Assessment of water quality and discharge in the Herring River, Wellfleet, Massachusetts, November 2015 to September 2017

The U.S. Geological Survey, Cape Cod National Seashore of the National Park Service, and Friends of Herring River cooperated from 2015 to 2017 to assess nutrient concentrations and fluxes across the ocean-estuary boundary at a dike on the Herring River in Wellfleet, Massachusetts. The purpose of this assessment was to characterize environmental conditions prior to a future removal of the dike, which has restricted saltwater inputs into the Herring River watershed for more than 100 years. Water temperature, dissolved oxygen, pH, and specific conductance were monitored continuously, and flow-weighted composite samples were collected approximately twice per month at the ocean-estuary boundary. Bidirectional discharge was computed for the U.S. Geological Survey Herring River at Chequessett Neck Road at Wellfleet, Massachusetts, streamgage (011058798) by using a stage-area rating and index-velocity ratings developed with acoustic Doppler current profile measurements made upstream and downstream from the dike. LOADEST regression modeling software was used to estimate nutrient fluxes (loads) from composite, paired nutrient concentration and discharge data in conjunction with continuous discharge data. Temperature, dissolved oxygen, pH, and specific conductance were also monitored continuously on two tributaries to the Herring River, Pole Dike Creek and Bound Brook, from late-May 2016 to mid-June 2017. Composite or discrete water samples were collected from the tributaries approximately twice per month in most months from late-May 2016 to mid-June 2017 and analyzed for total nitrogen, total phosphorus, and dissolved organic carbon. Flow-weighted concentrations of ammonium, nitrate, and total nitrogen on the Herring River at the dike on the ebb tide generally varied between 0.01 and 0.1, 0.003 and 0.03, and 0.3 and 0.7 milligram per liter as nitrogen, respectively. Flow-weighted concentrations of orthophosphate, total dissolved phosphorus, and total phosphorus generally varied between 0.002 and 0.02, 0.003 and 0.06, and 0.03 and 0.1 milligram per liter as phosphorus, respectively, on the ebb tide. Flow-weighted concentrations of silicate and dissolved organic carbon on the ebb tide generally varied between 0.08 and 3.0 milligrams per liter of silica (silicon dioxide), and 1.7 and 5.6 milligrams per liter of carbon, respectively. Ebb tide concentrations of nitrate were highest in winter and lowest in summer. By contrast, ebb tide concentrations of phosphorus species were highest in late summer and early fall and lowest in winter. Silica and dissolved organic carbon did not exhibit systematic variation in seasonal concentrations. There was uncertainty in estimates of nutrient fluxes, but the LOADEST-estimated fluxes indicated that annual (and in almost all cases seasonal) exports (ebb tides) exceeded inputs (flood tides). Ebb tide concentrations of ammonium, nitrate, total nitrogen, and silica were positively correlated with antecedent cumulative 7-day precipitation.

Massachusetts↗
USGS2018

Geochemical conditions and nitrogen transport in nearshore groundwater and the subterranean estuary at a Cape Cod embayment, East Falmouth, Massachusetts, 2013–14

Nitrogen transport and transformation were studied during 2013 to 2014 by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, in a subterranean estuary beneath onshore locations on the Seacoast Shores peninsula, a residential area in Falmouth, Massachusetts, served by septic systems and cesspools, and adjacent offshore locations in the Eel River, a saltwater embayment connected to the ocean. The field investigation included installation and sampling of clusters of wells and temporary sampling points near a transect extending from about 35 meters (m) onshore to 18 m offshore. The fresh groundwater at the study site formed a lens about 11 m thick at the shoreline that was underlain by saline groundwater. Groundwater flow in the water-table aquifer was oriented northwestward toward the embayment. Nitrate concentrations in the fresh groundwater at a site about 35 m onshore increased in the downward direction from less than 500 micromoles per liter near the water table to about 1,700 micromoles per liter just above the freshwater/saltwater transition zone. Dissolved oxygen was largely absent in the onshore fresh groundwater. Distributions of salinity, dissolved oxygen, and nitrate at the shoreline and offshore generally were similar to those onshore; at some locations, however, shallow saline water was present above the freshwater, and there were scattered occurrences of elevated dissolved oxygen concentrations. Geochemical indicators of nitrate reduction, including concentrations of the reaction product nitrogen gas, stable isotope ratios of nitrate and nitrogen gas, and changes in alkalinity, provided evidence for nitrate reduction in two zones separated vertically by a zone 7–8 m thick with no evidence of nitrate reduction. The shallow nitrate-reduction zone was near the water table in fresh groundwater onshore, where nitrate reduction may be related to particular recharge conditions at nearby sources. The shallow nitrate-reduction zone also may be related to an interval of fine-grained sediments at about the same altitude (−1 to −6 m relative to the National Geodetic Vertical Datum of 1929), where flow is slower and reactive electron donors such as solid organic carbon, iron, or sulfide phases may be present to drive the reduction. The deep nitrate-reduction zone was near the freshwater/saltwater transition zone, where nitrate reduction may be related to mixing of freshwater containing nitrate and saltwater containing dissolved organic carbon and ammonium, or to fine-grained sediments near the transition zone. The maximum amount of nitrate converted to nitrogen gas was estimated to be less than or equal to 300 micromoles per liter in both nitrate-reduction zones. The presence of nitrate and low dissolved oxygen concentrations in the 7–8-meter-thick zone between the shallow and deep nitrate-reduction zones are conditions that could permit nitrate reduction. The absence of evidence of nitrate reduction in the high-nitrate zone may have resulted from the lack of reactive electron donors in that depth interval. The high-nitrate zone dissipated somewhat in the offshore direction, but the current study did not extend far enough to encompass the fresh groundwater discharge area or determine how much of the nitrate was removed prior to discharge. A shallow intertidal saltwater cell was formed during a spring tide by saltwater infiltration during tidal run-up on the beach. Nitrate reduction might have occurred if nitrate-containing fresh groundwater discharging to the estuary mixed with the saltwater containing dissolved organic carbon in this zone, but samples collected from the intertidal saltwater cell during this study were not analyzed for indicators of nitrate reduction. Elevated dissolved oxygen concentrations in fresh groundwater 9 m offshore may indicate that groundwater flow was partly oblique to the sampling transect or that groundwater from a regional flow system was converging under the river near the study area. Flow directions also may have been affected by aquifer heterogeneity such as the shallow fine-grained sediments onshore and at the bottom of the Eel River. Improved understanding of the fate of nitrate in this type of complex setting might be gained by including additional characterization of aquifer heterogeneity and groundwater flow and extending investigations of nitrate reduction to the shallow sediments in the intertidal saltwater cell and adjacent subtidal zone and to locations farther offshore beneath the estuary.

Massachusetts↗
USGS2015

Simulation of nitrogen attenuation in a subterranean estuary, representative of the southern coast of Cape Cod, Massachusetts

A two-dimensional model was developed by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, to assess flow and chemical reaction associated with groundwater discharge through the subterranean estuary representative of coastal salt ponds of southern Cape Cod. The model simulated both the freshwater and saltwater flow systems and accounted for density-dependent flow, tidal fluctuation, and chemical reactivity among oxygen, dissolved organic carbon, nitrate, and ammonia. Not previously incorporated into one model, the interaction of these effects can now be simulated in the subterranean estuary context. An analysis of the flow system under mean-tide conditions was conducted first to provide the initial conditions for a subsequent analysis that included the effects of tidal fluctuations. Tidal fluctuations were simulated with a repeated couplet that represented a high tide-low tide sequence and alternating locations of head-dependent flux boundaries placed along the simulated seabed, above and below the levels of the respective high and low tides. Boundary conditions for chemical species included nitrate in recharge, and oxygen and organic matter (including organic nitrogen) in infiltrating solutions of head-dependent boundaries. Reaction chemistry was limited to oxidative degradation of organic matter (including remineralization of ammonia) with oxygen or nitrate as electron acceptors and nitrification of ammonia in the presence of oxygen. Simulations using the SEAWAT-2000 computer program resulted in two mixing zones-between freshwater and saltwater in a deep saltwater wedge and in an intertidal salt zone, which results from tidal fluctuation. The mixing zones are the principal locations where nitrogen attenuation reactions occurred-between organic matter in the saltwater zones of the aquifer and nitrate in the freshwater zone. In mean-tide PHT3D model simulations, 15 percent of nitrogen that is recharged was attenuated because of reaction with dissolved organic matter, a denitrification reaction that reduces nitrate to nitrogen gas. When a fluctuating tide was simulated, the amount of recharged nitrogen that was denitrified increased to 20 percent. Chemical reaction was controlled by the rate of mixing of freshwater and saltwater, which contained the reactants nitrate and dissolved organic matter, respectively, necessary for nitrogen attenuation reactions to take place. Reaction occurred in both the deep saltwater wedge and in an increased denitrification. However, mixing may also have been enhanced partly by numerical dispersion.

Massachusetts↗
USGS2015

Mercury in Pacific bluefin tuna (Thunnus orientalis):bioaccumulation and trans-Pacific Ocean migration

Pacific bluefin tuna (Thunnus orientalis) have the largest home range of any tuna species and are well known for the capacity to make transoceanic migrations. We report the measurement of mercury (Hg) concentrations in wild Pacific bluefin tuna (PBFT), the first reported with known size-of-fish and capture location. The results indicate juvenile PBFT that are recently arrived in the California Current from the western Pacific Ocean have significantly higher Hg concentrations in white muscle (0.51 ug/g wet mass, wm) than PBFT of longer California Current residency (0.41 ug/g wm). These new arrivals are also higher in Hg concentration than PBFT in farm pens (0.43 ug/g wm) that were captured on arrival in the California Current and raised in pens on locally derived feed. Analysis by direct Hg analyzer and attention to Hg by tissue type and location on the fish allowed precise comparisons of mercury among wild and captive fish populations. Analysis of migration and nearshore residency, determined through extensive archival tagging, bioaccumulation models, trophic investigations, and potential coastal sources of methylmercury, indicates Hg bioaccumulation is likely greater for PBFT juvenile habitats in the western Pacific Ocean (East China Sea, Yellow Sea) than in the eastern Pacific Ocean (California Current). Differential bioaccumulation may be a trophic effect or reflect methylmercury availability, with potential sources for coastal China (large hypoxic continental shelf receiving discharge of three large rivers, and island-arc volcanism) different from those for coastal Baja California (small continental shelf, no large rivers, spreading-center volcanism).

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2011

Determination of dilution factors for discharge of aluminum-containing wastes by public water-supply treatment facilities into lakes and reservoirs in Massachusetts

Dilution of aluminum discharged to reservoirs in filter-backwash effluents at water-treatment facilities in Massachusetts was investigated by a field study and computer simulation. Determination of dilution is needed so that permits for discharge ensure compliance with water-quality standards for aquatic life. The U.S. Environmental Protection Agency chronic standard for aluminum, 87 micrograms per liter (μg/L), rather than the acute standard, 750 μg/L, was used in this investigation because the time scales of chronic exposure (days) more nearly match rates of change in reservoir concentrations than do the time scales of acute exposure (hours). Whereas dilution factors are routinely computed for effluents discharged to streams solely on the basis of flow of the effluent and flow of the receiving stream, dilution determination for effluents discharged to reservoirs is more complex because (1), compared to streams, additional water is available for dilution in reservoirs during low flows as a result of reservoir flushing and storage during higher flows, and (2) aluminum removal in reservoirs occurs by aluminum sedimentation during the residence time of water in the reservoir. Possible resuspension of settled aluminum was not considered in this investigation. An additional concern for setting discharge standards is the substantial concentration of aluminum that can be naturally present in ambient surface waters, usually in association with dissolved organic carbon (DOC), which can bind aluminum and keep it in solution. A method for dilution determination was developed using a mass-balance equation for aluminum and considering sources of aluminum from groundwater, surface water, and filter-backwash effluents and losses caused by sedimentation, water withdrawal, and spill discharge from the reservoir. The method was applied to 13 reservoirs. Data on aluminum and DOC concentrations in reservoirs and influent water were collected during the fall of 2009. Complete reservoir volume was determined to be available for mixing on the basis of vertical and horizontal aluminum-concentration profiling. Losses caused by settling of aluminum were assumed to be proportional to aluminum concentration and reservoir area. The constant of proportionality, as a function of DOC concentration, was established by simulations in each of five reservoirs that differed in DOC concentration. In addition to computing dilution factors, the project determined dilution factors that would be protective with the same statistical basis (frequency of exceedance of the chronic standard) as dilutions computed for streams at the 7-day-average 10-year-recurrence annual low flow (the 7Q10). Low-flow dilutions are used for permitting so that receiving waters are protected even at the worst-case flow levels. The low-flow dilution factors that give the same statistical protection are the lowest annual 7-day-average dilution factors with a recurrence of 10 years, termed 7DF10s. Determination of 7DF10 values for reservoirs required that long periods of record be simulated so that dilution statistics could be determined. Dilution statistics were simulated for 13 reservoirs from 1960 to 2004 using U.S. Geological Survey Firm-Yield Estimator software to model reservoir inputs and outputs and present-day values of filter-effluent discharge and aluminum concentration. Computed settling velocities ranged from 0 centimeters per day (cm/d) at DOC concentrations of 15.5 milligrams per liter (mg/L) to 21.5 cm/d at DOC concentrations of 2.7 mg/L. The 7DF10 values were a function of aluminum effluent discharged. At current (2009) effluent discharge rates, the 7DF10 values varied from 1.8 to 115 among the 13 reservoirs. In most cases, the present-day (2009) discharge resulted in receiving water concentrations that did not exceed the standard at the 7DF10. Exceptions were one reservoir with a very small area and three reservoirs with high concentrations of DOC. Maximum permissible discharges were determined for water-treatment plants by adjusting discharges upward in simulations until the 7DF10 resulted in reservoir concentrations that just met the standard. In terms of aluminum flux, these discharges ranged from 0 to 28 kilograms of aluminum per day.

Massachusetts↗
USGS2009

Mercury sources to Lake Ozette and Lake Dickey: Highly contaminated remote coastal lakes, Washington State, USA

Mercury concentrations in largemouth bass and mercury accumulation rates in age-dated sediment cores were examined at Lake Ozette and Lake Dickey in Washington State. Goals of the study were to compare concentrations in fish tissues at the two lakes with a larger statewide dataset and examine mercury pathways to the lakes. After accounting for fish length, tissue concentrations at the lakes were significantly higher than other Washington State lakes. Wet deposition and historical atmospheric monitoring from the area show no indication of enhanced local or regional deposition. Sediment core records from the lakes indicate rising sedimentation rates coinciding with logging in the lakes’ drainages has greatly increased the net flux of mercury to the waterbodies.

Washington↗
USGS2008

Submarine Hydrogeological Data from Cape Cod National Seashore

In order to test hypotheses about ground water flow under and into estuaries and the Atlantic Ocean, geophysical surveys, geophysical probing, submarine ground-water sampling, and sediment coring were conducted by U.S. Geological Survey (USGS) scientists at Cape Cod National Seashore (CCNS) in Massachusetts from 2004 through 2006. This USGS Open-File Report presents the data collected as a result of these field efforts.

Open-File Report↗
USGS2007

Monitoring Ground-Water Quality in Coastal Ecosystems

INTRODUCTION The Cape Cod National Seashore (CACO) extends along more than 70 km of Atlantic Ocean open-beach coastline and includes three large saltwater bays - Wellfleet Harbor, Nauset Marsh, and Pleasant Bay (fig. 1). CACO encompasses about 18,000 ha of uplands, lakes, wetlands, and tidal lands (Godfrey and others, 1999) including most habitats typical of the sandy coast in National seashores and parks extending southward from Massachusetts to Florida. In 1995, CACO was selected by the National Park Service (NPS) as a prototype park typifying the Atlantic and Gulf Coast biogeographic region for long-term coastal ecosystem monitoring. The U.S. Geological Survey (USGS) is currently (2007) assisting the NPS in the development of protocols for a Long-Term Coastal Ecosystem Monitoring Program at the CACO in Massachusetts. The overall purpose of the monitoring program is to characterize both natural and human-induced change in the biological resources of the CACO, over a time scale of decades, in the context of a changing global ecosystem.

Open-File Report↗
USGS2006

Effects of aquifer travel time on nitrogen transport to a coastal embayment

Effects of aquifer travel time on nitrogen reaction and loading to Popponesset Bay, a eutrophic coastal embayment on western Cape Cod, Massachusetts, are evaluated through hydrologic analysis of flow and transport. Approximately 10% of the total nitrogen load to the embayment is intercepted by fresh water ponds and delivered to the coast by connecting streams. For the nitrogen load not intercepted by ponds, we compare two steady-state methods of analyzing nitrogen loss in the aquifer, one using a constant-loss factor and the other time-dependent loss rates. The constant-loss method, which assumes that all similar land uses have the same per unit area loading rate to surface water regardless of location within the watershed, predicts that 42% of the nonpond watershed nitrogen load originated within the zero to 2 yr time-of-travel zone, which is 40% of the contributing area. The time-of-travel loss method calculates loss rates based on aquifer travel times and denitrification reaction kinetics, evaluated separately for carbon-unlimited and carbon-limited cases. Time-of-travel loss calculations for percent of nonpond load that originated within the area of < 2 yr aquifer residence time are 64% when carbon is not limiting, but only 49% when carbon limitation is included, not greatly different from the constant-loss method. A feature of the kinetics used is that carbon (and the denitrified nitrogen) is lost rather quickly in the aquifer travel path, after which carbon limitation stops denitrification altogether. Carbon limitation causes the time-of-travel loss model to approximate the constant-loss model such that in most of the watershed, a nearly constant fraction of the nitrogen input is lost in both models.

Groundwater↗
USGS2005

Submarine groundwater discharge to a small estuary estimated from radon and salinity measurements and a box model

Submarine groundwater discharge was quantified by a variety of methods for a 4-day period during the early summer of 2004, in Salt Pond, adjacent to Nauset Marsh, on Cape Cod, USA. Discharge estimates based on radon and salinity took advantage of the presence of the narrow channel connecting Salt Pond to Nauset Marsh, which allowed constructing whole-pond mass balances as water flowed in and out due to tidal fluctuations. The data suggest that less than one quarter of the discharge in the vicinity of Salt Pond happened within the pond itself, while three quarters or more of the discharge occurred immediately seaward of the pond, either in the channel or in adjacent regions of Nauset Marsh. Much of this discharge, which maintains high radon activities and low salinity, is carried into the pond during each incoming tide. A box model was used as an aid to understand both the rates and the locations of discharge in the vicinity of Salt Pond. The model achieves a reasonable fit to both the salinity and radon data assuming submarine groundwater discharge is fresh and that most of it occurs either in the channel or in adjacent regions of Nauset Marsh. Salinity and radon data, together with seepage meter results, do not rule out discharge of saline groundwater, but suggest either that the saline discharge is at most comparable in volume to the fresh discharge or that it is depleted in radon. The estimated rate of fresh groundwater discharge in the vicinity of Salt Pond is 3000-7000 m 3 d -1 . This groundwater flux estimated from the radon and salinity data is comparable to a value of 3200-4500 m 3 d -1 predicted by a recent hydrologic model (Masterson, 2004; Colman and Masterson, 2004), although the model predicts this rate of discharge to the pond whereas our data suggest most of the groundwater bypasses the pond prior to discharge. Additional work is needed to determine if the measured rate of discharge is representative of the long-term average, and to better constrain the rate of groundwater discharge seaward of Salt Pond.

Biogeosciences Discussions↗
USGS2005

Response curves for phosphorus plume lengths from reactive-solute-transport simulations of onland disposal of wastewater in noncarbonate sand and gravel aquifers

Surface-water resources in Massachusetts often are affected by eutrophication, excessive plant growth, which has resulted in impaired use for a majority of the freshwater ponds and lakes and a substantial number of river-miles in the State. Because supply of phosphorus usually is limiting to plant growth in freshwater systems, control of phosphorus input to surface waters is critical to solving the impairment problem. Wastewater is a substantial source of phosphorus for surface water, and removal of phosphorus before disposal may be necessary. Wastewater disposed onland by infiltration loses phosphorus from the dissolved phase during transport through the subsurface and may be an effective disposal method; quantification of the phosphorus loss can be simulated to determine disposal feasibility. In 2003, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, initiated a project to simulate distance of phosphorus transport in the subsurface for plausible conditions of onland wastewater disposal and subsurface properties. A coupled one-dimensional unsaturated-zone and three-dimensional saturated-zone reactive-solute-transport model (PHAST) was used to simulate lengths of phosphorus plumes. Knowledge of phosphorus plume length could facilitate estimates of setback distances for wastewater-infiltration sites from surface water that would be sufficient to protect the surface water from eutrophication caused by phosphorus transport through the subsurface and ultimate discharge to surface water. The reactive-solute-transport model PHAST was used to simulate ground-water flow, solute transport, equilibrium chemistry for dissolved and sorbed species, and kinetic regulation of organic carbon decomposition and phosphate mineral formation. The phosphorus plume length was defined for the simulations as the maximum extent of the contour for the 0.015 milligram-per-liter concentration of dissolved phosphorus downgradient from the infiltration bed after disposal cessation. Duration of disposal before cessation was assumed to be 50 years into an infiltration bed of 20,000 square feet at the rate of 3 gallons per square foot per day. Time for the maximum extent of the phosphorus plume to develop is on the order of 100 years after disposal cessation. Simulations indicated that phosphorus transport beyond the extent of the 0.015 milligram-per-liter concentration contour was never more than 0.18 kilogram per year, an amount that would likely not alter the ecology of most surface water. Simulations of phosphorus plume lengths were summarized in a series of response curves. Simulated plume lengths ranged from 200 feet for low phosphorus-concentration effluents (0.25 milligram per liter) and thick (50 feet) unsaturated zones to 3,400 feet for high phosphorus-concentration effluents (14 milligrams per liter) discharged directly into the aquifer (unsaturated-zone thickness of 0 feet). Plume length was nearly independent of unsaturated-zone thickness at phosphorus concentrations in the wastewater that were less than 2 milligrams per liter because little or no phosphorus mineral formed at low phosphorus concentrations. For effluents of high phosphorus concentration, plume length varied from 3,400 feet for unsaturated-zone thickness of 0 to 2,550 feet for unsaturated-zone thickness of 50 feet. Model treatments of flow and equilibrium-controlled chemistry likely were more accurate than rates of kinetically controlled reactions, notably precipitation of iron-phosphate minerals; the kinetics of such reactions are less well known and thus less well defined in the model. Sensitivity analysis indicated that many chemical and physical aquifer properties, such as hydraulic gradient and model width, did not affect the simulated plume length appreciably, but duration of discharge, size of infiltration bed, amount of dispersion, and number of sorption sites on the aquifer sediments did affect plume length appreciably. Because simulation of plume length in carbonate-mineral sediments indicated that the plume would be substantially longer than in noncarbonate-mineral sediments, the application of the response curves in locations with carbonate-mineral sediments would be inappropriate. The effect of carbonate minerals in sediments is to increase pH, which causes decreased sorption of phosphorus on aquifer sediments. Phosphorus removal from solution by precipitation onto aquifer sediments is more efficient at high concentrations of disposed phosphorus than at low concentrations. At very low phosphorus concentrations, the solubility product of phosphorus minerals is not exceeded and no phosphorus mineral forms. An important consequence is that removal of dissolved phosphorus from the plume by processes in the subsurface is decreased the more that removal efforts are applied in treatment before wastewater is disposed. Model simulations indicate that removal of phosphorus from wastewater disposed through septic systems would have the advantage of efficient phosphorus removal in the subsurface because phosphorus concentrations are high in septic-system effluent. Short plume lengths result from wastewater disposal through septic systems because of the efficient phosphorus removal and because of the low volume of wastewater involved. The simulation results for small-volume systems are not quantitative, however, because wastewater-infiltration rates are much lower than those of the higher-volume system that was used to calibrate the model and to create the plume-length response curves. The response curves for phosphorus plume lengths, as defined by the maximum extent of the 0.015 milligram-per-liter concentration contour, is clearly defined in the model simulations, although the relation between simulated plume length and protective setback distance is subject to interpretation. Phosphorus does move beyond the point at which the simulated 0.015 milligram-per-liter concentration contour has stopped, so that a determination of protective setback distance must include a consideration of whether that continued flux, or some other flux amount, is appropriate. Also, simulations indicate that phosphorus plumes do not reach their full extent until 50 to 200 years after disposal cessation, depending on concentration of phosphorus disposed. No phosphorus plume has been monitored for that long after cessation, so there is no way to verify the long-term simulation results.

Massachusetts↗
USGS2003

Reactive-transport simulation of phosphorus in the sewage plume at the Massachusetts Military Reservation, Cape Cod, Massachusetts

The subsurface transport of phosphorus introduced by the disposal of treated sewage effluent to ground-infiltration disposal beds at the Massachusetts Military Reservation on western Cape Cod was simulated with a three-dimensional reactive-transport model. The simulations were used to estimate the load of phosphorus transported to Ashumet Pond during operation of the sewage-treatment plant from 1936 to 1995 and for 60 years following cessation of sewage disposal. The model accounted for spatial and temporal changes in water discharge from the sewage-treatment plant, ground-water flow, transport of associated chemical constituents, and a set of chemical reactions, including phosphorus sorption on aquifer materials, dissolution and precipitation of iron- and manganese-oxyhydroxide and iron phosphate minerals, organic carbon sorption and decomposition, cation sorption, and irreversible denitrification. The flow and transport in the aquifer were simulated by using parameters consistent with those used in previous flow models of this area of Cape Cod, except that numerical dispersion was much larger than the physical dispersion estimated in previous studies. Sorption parameters were fit to data derived from phosphorus sorption and desorption laboratory column experiments. Rates of organic carbon decomposition were adjusted to match the location of iron concentrations in an anoxic iron zone within the sewage plume. The sensitivity of the simulated load of phosphorus transported to Ashumet Pond was calculated for a variety of processes and input parameters. Model limitations included large uncertainties associated with the loading of the sewage beds, the flow system, and the chemistry and sorption characteristics in the aquifer. The results of current model simulations indicate a small load of phosphorus transported to Ashumet Pond during 1965-85, but this small load was particularly sensitive to model parameters that specify flow conditions and the chemical process by which non-desorbable phosphorus is incorporated in the sediments. The uncertainties were large enough to make it difficult to determine whether loads of phosphorus transported to Ashumet Pond in the 1990s were greater or less than loads during the previous two decades. The model simulations indicate substantial discharge of phosphorus to Ashumet Pond after about 1965. After the period 2000-10 the simulations indicate that the load of phosphorus transported to Ashumet Pond decreases continuously, but the load of phosphorus remains substantial for many decades. The current simulations indicate a peak in phosphorus discharge to Ashumet Pond of about 1,000 kilograms per year during the 1990s; however, comparisons of simulated phosphorus concentrations with measured concentrations in 1993 indicate that the peak in phosphorus load transported to Ashumet Pond may be larger and moving more quickly in the model simulations than in the aquifer. The results of the three-dimensional reactive-transport simulations are consistent with the loading history, experimental laboratory data, and field measurements. The results of the simulations adequately reproduce the spatial distribution of phosphorus concentrations measured in 1993, the magnitude of changes in phosphorus concentration with time in a profile near the disposal beds following cessation of sewage disposal, the observed iron zone in the sewage plume, the approximate flow of treated sewage effluent into Ashumet Valley, and laboratory-column data for phosphorus sorption and desorption.

Massachusetts↗
USGS2003

Natural remediation potential of arsenic-contaminated ground water

Migration of leachate from a municipal landfill in Saco, Maine has resulted in arsenic concentrations in ground water as high as 647 μg/L . Laboratory experimental data indicate the primary source of arsenic to be reductive dissolution of arsenic-enriched iron oxyhydroxides in the aquifer by organic carbon in landfill leachate. A core from an uncontaminated part of the aquifer yielded no dissolved iron or arsenic when leached with oxic ground water. Eluent ground water spiked with organic carbon in order to create reducing conditions mobilized both ferrous iron and arsenite from this core. The landfill was capped in early 1998 to eliminate the source of leachate. Cores from the contaminated portion of the aquifer were collected and leached with uncontaminated ground water in the laboratory to simulate natural remediation conditions. Data from these experiments show that significant concentrations of labile organic carbon have accumulated on aquifer solids, causing significant biological oxygen demand. In laboratory leaching experiments of the most contaminated core, the organic carbon caused complete consumption of the influent dissolved oxygen (6 mg/L) for 220 pore volumes. Arsenic leaching from contaminated cores rapidly decreased in concentration initially in response to flushing with uncontaminated ground water. Subsequent leaching produced more gradual decreases in dissolved arsenic concentrations, controlled by a combination of reductive dissolution of arsenic-enriched iron oxyhydroxides and adsorption/desorption. In leachate from the most contaminatedcore, arsenic concentrations exceeded the new United States Environmental Protection Agency drinking-water standard of 10 μg/L for more than 200 pore volumes. A geochemical model simulated the concentration of selected constituents as uncontaminated ground water eluted through contaminated aquifer solids. Concentrations of dissolved oxygen, arsenic, and iron, in leachate from one core were used to calibrate the model. This model was validated by successfully simulating constituent concentrations in leachate from cores collected from other contaminated areas of this aquifer.

Saco↗
USGS2001

Methodology and significance of studies of atmospheric deposition in highway runoff

Atmospheric deposition and the processes that are involved in causing and altering atmospheric deposition in relation to highway surfaces and runoff were evaluated nationwide. Wet deposition is more easily monitored than dry deposition, and data on wet deposition are available for major elements and water properties (constituents affecting acid deposition) from the inter-agency National Atmospheric Deposition Program/ National Trends Network (NADP/NTN). Many trace constituents (metals and organic compounds) of interest in highway runoff loads, however, are not included in the NADP/NTN. Dry deposition, which constitutes a large part of total atmospheric deposition for many constituents in highway runoff loads, is difficult to monitor accurately. Dry-deposition rates are not widely available. Many of the highway-runoff investigations that have addressed atmospheric-deposition sources have had flawed investigative designs or problems with methodology. Some results may be incorrect because of reliance on time-aggregated data collected during a period of changing atmospheric emissions. None of the investigations used methods that could accurately quantify the part of highway runoff load that can be attributed to ambient atmospheric deposition. Lack of information about accurate ambient deposition rates and runoff loads was part of the problem. Samples collected to compute the rates and loads were collected without clean-sampling methods or sampler protocols, and without quality-assurance procedures that could validate the data. Massbudget calculations comparing deposition and runoff did not consider loss of deposited material during on-highway processing. Loss of deposited particles from highway travel lanes could be large, as has been determined in labeled particle studies, because of resuspension caused by turbulence from passing traffic. Although a cause of resuspension of large particles, traffic turbulence may increase the rate of deposition for small particles and gases by impaction, especially during precipitation periods. Ultimately, traffic and road maintenance may be determined to be the source of many constituents measured in highway runoff previously attributed to ambient atmospheric deposition. An investigative design using tracers of ambient deposition that are not present in highway traffic sources could determine conclusively what fraction of highway runoff load is contributed by ambient atmospheric deposition.

Open-File Report↗