Skip Links

USGS - science for a changing world

Open-File Report 2015–1085

Prepared in cooperation with the U.S. Environmental Protection Agency

Simulation of Nitrogen Attenuation in a Subterranean Estuary, Representative of the Southern Coast of Cape Cod, Massachusetts

By John A. Colman, Carl S. Carlson, and Clare Robinson

Thumbnail of and link to report PDF (3.45 MB)Abstract

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.

First posted June 9, 2015

For additional information, contact:
Director, New England Water Science Center
U.S. Geological Survey
10 Bearfoot Road
Northborough, MA 01532
(508) 490–5000

Part or all of this report is presented in Portable Document Format (PDF). For best results viewing and printing PDF documents, it is recommended that you download the documents to your computer and open them with Adobe Reader. PDF documents opened from your browser may not display or print as intended. Download the latest version of Adobe Reader, free of charge. More information about viewing, downloading, and printing report files can be found here.

Suggested citation:

Colman, J.A., Carlson, C.S, and Robinson, Clare, 2015, Simulation of nitrogen attenuation in a subterranean estuary, representative of the southern coast of Cape Cod, Massachusetts: U.S. Geological Survey Open-File Report 2015–1085, 35 p.,

ISSN 2331-1258 (online)




Sources and Quality of Model Data

Model Development and Approach

Simulation of Nitrogen Attenuation with a Model of Reactive-Solute Transport

Limitations of Analysis

Summary and Conclusions

References Cited

Accessibility FOIA Privacy Policies and Notices logo U.S. Department of the Interior | U.S. Geological Survey
Page Contact Information: Contact USGS
Page Last Modified: Wednesday, December 07, 2016, 09:42:52 PM