Estuarine, Coastal and Shelf Science 117 (2013) 107e124
Contents lists available at SciVerse ScienceDirect
Estuarine, Coastal and Shelf Science
journal homepage: www.elsevier.com/locate/ecss
Sediment flux modeling: Calibration and application for coastal systems
Damian C. Brady a,*, Jeremy M. Testa b, Dominic M. Di Toro c, Walter R. Boynton d, W. Michael Kemp b a School of Marine Sciences, University of Maine, 193 Clark Cove Road, Walpole, ME 04573, USA b Horn Point Laboratory, University of Maryland Center for Environmental Science, 2020 Horns Point Rd., Cambridge, MD 21613, USA c Department of Civil and Environmental Engineering, University of Delaware, 356 DuPont Hall, Newark, DE 19716, USA d Chesapeake Biological Laboratory, University of Maryland Center for Environmental Science, P.O. Box 38, Solomons, MD 20688, USA article info abstract
Article history: Benthicepelagic coupling in shallow estuarine and coastal environments is an important mode of Received 19 June 2012 particle and solute exchange and can influence lag times in the recovery of eutrophic ecosystems. Links Accepted 7 November 2012 between the water column and sediments are mediated by particulate organic matter (POM) deposition Available online 21 November 2012 to the sediment and its subsequent decomposition. Some fraction of the regenerated nutrients are returned to the water column. A critical component in modeling sediment fluxes is the organic matter Keywords: flux to the sediment. A method is presented for estimating POM deposition, which is especially difficult sediments to measure, by using a combination of long term time series of measured inorganic nutrient fluxes and modeling fl e Chesapeake Bay a mechanistic sediment ux model (SFM). A Hooke Jeeves pattern search algorithm is used to adjust fi þ fl ammonium organic matter deposition to t ammonium (NH4 ) ux at 12 stations in Chesapeake Bay for up to 17 sediment oxygen demand years. diagenesis POM deposition estimates matched reasonably well with sediment trap estimates on average (within 10%) and were strongly correlated with previously published estimates based upon sediment chloro- fi þ fl phyll-a (chl-a) distribution and degradation. Model versus eld data comparisons of NH4 ux, sediment 2 þ oxygen demand (SOD), sulfate (SO4 ) reduction rate, porewater NH4 , and sedimentary labile carbon concentrations further validated model results and demonstrated that SFM is a powerful tool to analyze solute fluxes. Monthly model-field data comparisons clearly revealed that POM decomposition in the original SFM calibration was too rapid, which has implications for lagged ecosystem responses to nutrient management efforts. Finally, parameter adjustments have been made that significantly improve model-field data comparisons and underscore the importance of revisiting and recalibrating models as long time series (>5 years) become available. Ó 2012 Elsevier Ltd. All rights reserved.
1. Introduction demand (Kemp and Boynton, 1992). Despite the importance of particulate organic matter (POM) deposition to coastal sediment Strong coupling between pelagic and benthic habitats is processes, few direct measurements have been made because of a fundamental feature of most shallow coastal aquatic ecosystems the inherent difficulty in deriving depositional rates from sediment (Nixon, 1981; Kamp-Nielsen, 1992; Kemp and Boynton, 1992; traps (e.g., Buesseler et al., 2007), especially in shallow coastal Soetaert and Middelburg, 2009). A primary mechanism that waters where resuspended sediment material can accumulate in connects benthic and pelagic habitats is the downward transport of traps (Kemp and Boynton, 1984; Hakanson et al., 1989; Kozerski, organic material, which may be composed of senescent algal cells, 1994; Gust and Kozerski, 2000). zooplankton fecal pellets, and relatively refractory organic Eutrophication is a common phenomenon in coastal aquatic compounds imported from marshes or rivers. Organic matter that ecosystems where increased nutrient loading from the adjacent reaches the sediment drives sediment biogeochemical processes watershed promotes increased phytoplankton growth that and resultant nutrient fluxes (Jensen et al., 1990), feeds benthic degrades water quality (e.g. hypoxia and harmful algal blooms, organisms (Heip et al., 1995), and can control sediment oxygen Cloern, 2001; Diaz and Rosenberg, 2008; Nixon, 1995). Typically, nutrient stimulated phytoplankton growth leads to higher rates of POM deposition, and thus higher rates of sediment oxygen demand e fl * Corresponding author. (SOD) and sediment water nutrient uxes (e.g., Kelly and Nixon, E-mail address: [email protected] (D.C. Brady). 1984; Kelly et al., 1985). Eutrophication mitigation via reductions
0272-7714/$ e see front matter Ó 2012 Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.ecss.2012.11.003 108 D.C. Brady et al. / Estuarine, Coastal and Shelf Science 117 (2013) 107e124 in nutrient and organic carbon loads to estuaries has become a goal description of the development, calibration, theory, and applica- of most water management agencies (Boesch et al., 2001; Conley tion of SFM. et al., 2002). To evaluate and predict the consequences of such The purpose of this study is to: (1) utilize a stand-alone version eutrophication mitigation strategies, observational data are linked of SFM to support and analyze sediment process observations at to models of varying complexity that include linked hydrodynamic a variety of Chesapeake Bay stations; (2) use SFM to estimate and biogeochemical processes, including sediment process models organic matter deposition to sediments in Chesapeake Bay; and (3) (Cerco and Cole, 1993; Lin et al., 2007). calibrate and validate SFM predictions with a multi-decadal time- In shallow coastal ecosystems, the need for sediment process series of water column concentrations and sedimentewater models is especially important. In very shallow ecosystems (<5 m), exchanges of nutrients and O2 at 12 stations in Chesapeake Bay. þ fl sediments may be populated by submerged vascular plants and/or We focus our process studies on NH4 ,NO3 and O2 uxes. SFM also benthic algal communities, both of which modify sediment computes phosphorus and silica fluxes and the calibration and biogeochemical reactions via nutrient uptake and sediment validation of these fluxes will be the focus of a companion paper oxygenation (Miller et al., 1996; McGlathery et al., 2007). In (Testa et al., in press). It is suggested that SFM can serve as a widely moderately shallow systems (5e50 m), pelagicebenthic coupling applicable tool for a variety of coastal ecosystems to aid in sediment involves sediments as sites for high rates of nutrient recycling processes and water quality investigations. (Cowan et al., 1996) and oxygen consumption (Kemp et al., 1992). Depleted bottom-water O , a common feature of eutrophication, 2 2. Methods complicates predictions of sediment nutrient exchange because low O concentrations exert control over sediment biogeochemical 2 SFM was previously calibrated and validated in Chesapeake Bay reactions (Middelburg and Levin, 2009; Testa and Kemp, 2012). for the years 1985e1988, using sedimentewater flux and Models of sediment diagenetic processes can be used to simulate overlying-water nutrient and O data that were available at that exchanges of particulate and dissolved substances between the 2 time for 8 sites (Chapter 14 in Di Toro, 2001). More than two water column and sediments. Such models are particularly useful decades later, data availability has expanded. New calibration and to simulate nutrient storage in sediments and resulting delays in validation can improve SFM simulation skill for Chesapeake Bay ecosystem recovery following eutrophication abatement (Conley and can demonstrate its utility as a stand-alone tool available for et al., 2009; Soetaert and Middelburg, 2009). use in other systems. In this paper we describe and analyze SFM In general, sediment process model structures range from performance in the upper portion of Chesapeake Bay using field relatively simple empirical relationships to more complex process data from 12 stations (Fig. 1) over 4e17 year time series. SFM can be simulations that include time-varying state variables (Boudreau, run on a personal computer, executing a 25-year run on the time- 1991b; Cerco and Noel, 2005). Simple model representations include assigning a constant O2 or nutrient sedimentewater flux (e.g., Scully, 2010) or using basic parameterizations of sedimente water flux as a function of overlying water conditions (Fennel et al., 2006; Hetland and DiMarco, 2008). More complex process models may simulate one or two layers, each of which represent a particular chemical environment, and such models tend to demand limited computational effort (Emerson et al., 1984; Di Toro, 2001). Process models may also resolve depth into several layers, allowing for simulations of pore-water constituent vertical profiles (Boudreau, 1991a; Dhakar and Burdige, 1996; Soetaert et al., 1996; Cai et al., 2010). Depth resolution in such complex models is usually associated with high computational demand, therefore less vertically-resolved models are commonly used when sediment biogeochemical models are coupled to water column models (Cerco and Noel, 2005; Sohma et al., 2008). A sediment biogeochemical model was developed to link to spatially articulated water column models describing biogeo- chemical processes at a limited computational cost (Di Toro, 2001). This sediment flux model (SFM) separates sediment reac- tions into two layers and generates fluxes of key solutes (e.g., nitrogen, phosphorus, silica, dissolved oxygen (O2), sulfide, and methane). SFM has been successfully integrated into water quality models in many coastal systems, including Massachusetts Bay (Jiang and Zhou, 2008), Chesapeake Bay (Cerco and Cole, 1993; Cerco, 1995; Cerco and Noel, 2005), Delaware’s Inland Bays (Cerco and Seitzinger, 1997), and Long Island Sound (Di Toro, 2001). SFM can also be used as a stand-alone diagnostic tool in sediment process studies, especially over seasonal to decadal time scales. For example, it has been used to simulate sediment dynamics in the Rhode Island Marine Ecosystem Research Laboratory (MERL) mesocosms (Di Toro and Fitzpatrick, 1993; Di Toro, 2001). Such a stand-alone sediment modeling tool, when combined with high-quality time series observations, can be used for parameter Fig. 1. Map of northern Chesapeake Bay, on the east coast of the United States (inset), optimization, scenario analysis, and process investigations. showing the locations where sediment flux model (SFM) simulations were compared We direct the reader to Di Toro (2001) for a comprehensive to Sediment Oxygen and Nutrient Exchange (SONE) observations. D.C. Brady et al. / Estuarine, Coastal and Shelf Science 117 (2013) 107e124 109 scale of seconds, and a MATLAB interface is available for input approach is used to compute temporal changes in concentrations generation, model execution, post-processing, and plotting. of particulate and dissolved constituents, which in turn produce fluxes across the volume boundaries. 2.1. General model description The general forms of the equations are presented in Table 1 (Eqs. (1) and (2)). For example, one can replace CT1 with NH4,1 to compute þ What follows is a brief summary of the structure and initial the change in NH4 concentration in the aerobic layer (Eq. (1) in calibration of SFM (Di Toro and Fitzpatrick, 1993; Di Toro, 2001). Table 1). The governing equations are mass balance equations that k2= The modeling framework encompasses three processes: (1) the include biogeochemical reactions ðð 1 KL01ÞCT1Þ and dissolved and u sediment receives depositional fluxes of particulate organic carbon particulate mixing (i.e., burial; 2CT1, diffusion of dissolved mate- and nitrogen from the overlying water; (2) the mineralization of rial from the overlying water column; KL01(fd0CT0 fd1CT1) and particulate organic matter produces soluble intermediates that are between layers; KL12(fd2CT2 fd1CT1), and mixing of particulate quantified as diagenesis fluxes; and (3) solutes react in the aerobic material between layers;u12(fp2CT2 fp1CT1): see Eqs. (1) and (2) in and anaerobic layers of the sediment and portions are returned to Table 1). the overlying water (Fig. 2). The model assumes that organic matter Dissolved O2 concentrations are not modeled as a state variable mineralization is achieved by denitrification, sulfate reduction, and in SFM, but concentrations in each layer are determined as follows. methanogenesis. Carbon is first utilized by denitrification (which is In layer 2, the anaerobic layer, O2(2) is assumed to be zero. Because calculated independently) and any remaining carbon is processed O2(0), the overlying water concentration, is a boundary condition, 2 by sulfate reduction until SO4 becomes limiting, after which a linear decline is assumed from O2(0) to zero (O2(2)) in the aerobic methanogenesis occurs. To model these processes, SFM numeri- layer, thus O2(1) ¼ (O2(0) þ O2(2))/2 ¼ O2(0)/2. cally integrates mass-balance equations for chemical constituents A key to this model is an algorithm that continually updates the in two functional layers: an aerobic layer near the sediment water thickness of the aerobic layer, H1, by computing the product of the fi interface of variable depth (H1) and an anaerobic layer below that is diffusion coef cient ðDO2 Þ and the ratio of overlying-water O2 equal to the total sediment depth (10 cm) minus the depth of H1 concentration ([O2(0)]) to sediment oxygen demand (SOD). (Figs. 2e4). The SFM convention is to use “0” when referring to the overlying water, “1” when referring to the aerobic layer, and “2” ½O ð0Þ H ¼ D 2 (1) when referring to the anaerobic layer. A strict mass-balance 1 O2 SOD
This relationship was perhaps first suggested by Grote (1934) e quoted by Hutchinson (1957) e and verified by measurements (Jorgensen and Revsbech,1985; Cai and Sayles, 1996). The inverse of the second term on the right hand side of Equation (1) is the surface fi = = mass transfer coef cient (KL01 ¼ðSOD ½O2ð0Þ Þ ¼ ðDO2 H1Þ; Eq. (5) in Table 1). The surface mass transfer coefficient controls solute exchange between the aerobic layer and the overlying water column. The model uses the same mass transfer coefficient for all solutes since differences in the diffusion coefficients between solutes are subsumed in the kinetic parameters that are fit to data (Chapter 3 in Di Toro, 2001; Fennel et al., 2009). The SOD in Equation (1), which is computed as described below, also depends on the thickness of the aerobic layer. Therefore, the equation for SOD is nonlinear. It is solved numerically at each time step of the simulation. Dissolved and particle mixing between layers 1 and 2 (KL12 and u12, Eqs. (6) and (7) in Table 1) are modeled as a function of passive transport and proxies for the activities of benthic organisms. KL12 enhancement due to benthic activity is parameterized directly, that is, the dissolved mixing coefficient (Dd) was increased such that KL12 is 2e3 times molecular diffusion (Matisoff and Wang, 1998). The rate of mixing of sediment particles (u12; Eq. (7) in Table 1)by macrobenthos is quantified by estimating the apparent particle diffusion coefficient (Dp; Chapter 13 in Di Toro, 2001). In the model, particle mixing is controlled by temperature (1st term in Eq. (7) in Table 1; Balzer, 1996), carbon input (2nd term in Eq. (7) in Table 1; Robbins et al., 1989), and oxygen (3rd term in Eq. (7) in Table 1; Diaz and Rosenberg, 1995). In order to make the model self-consistent, that is to use only internally-computed variables in the parame- terizations, the model assumes that benthic biomass and therefore, particle mixing is correlated with the amount of labile carbon (e.g., POC1) present in the sediment. However, if excess carbon loading creates unfavorable oxygen conditions, the model accumulates this as benthic stress (S in Eq. (7) in Table 1). As O2 decreases, 1 kSS approaches zero. After the stress has passed, the minimum is carried forward for the rest of the year in order to simulate the Fig. 2. Generic schematic diagram of the sediment flux model (SFM), including state variables, transport and biogeochemical processes, and boundary conditions. Note that observation that benthic communities do not recover until the the depths of the aerobic (H1) and anaerobic (H2) layers vary over time. following year (Diaz and Rosenberg, 1995). 110 D.C. Brady et al. / Estuarine, Coastal and Shelf Science 117 (2013) 107e124
+ - MUINOMMA NH (0) NO (0) ETARTIN N (g) JPON 4 3 2 NMULOC RETAW NMULOC RETAW NMULOC
KL01 KL01 NH +(1) - 4 NO3 (1) κ NH ,1 + 4 - DIAGENESIS: NONE SOURCE: NH4 (1) NO3 (1) NITRIFICATION κ κ LAYER 1 LAYER AEROBIC NH ,1 NO ,1 + 4 - - 3 N (g) REACTION: NH4 (1) NO3 (1) REACTION: NO3 (1) 2 NOITACIFIRTIN NOITACIFIRTINED NH +(1) - 4 NO3 (1)
H1 KL12 KL12 SEDIMENT NH +(2) - 4 NO3 (2)
JN2 + DIAGENESIS: PON NH4 (2) SOURCE: NONE AMMONIFICATION
LAYER 2 LAYER κ
ANAEROBIC NO3,2 - N (g) REACTION: NONE REACTION: NO3 (2) 2 DENITRIFICATION a b H2
fl þ Fig. 3. Schematic representation of nitrogen transport and kinetics in the sediment ux model (SFM). Panels a & b represent the dynamics of the NH4 and NO3 models, respectively. Note: (1) there is no diagenesis (ammonification) in layer 1 (panel a). (2) there is no source of nitrate in the anaerobic layer, as no O2 is present (Panel b).
2.2. Ammonium flux kinetics (Eq. (8) in Table 1). Mass-transfer coefficients are þ employed to model diffusion of NH4 between the anaerobic and þ NH4 concentrations are computed for both aerobic and aerobic layers (KL12) and between the aerobic layer and the anaerobic layers via mass balances of biogeochemical and phys- overlying water (KL01). þ fi þ ical processes. Fig. 3a shows the sources and sinks of NH4 in the Calculation of nitri cation (i.e., the oxidation of NH4 to NO3 )is þ model. NH4 is produced by organic matter diagenesis (JN2 in a good example of a process that depends on the depth of H1. The fi þ Fig. 3a, JT2 in Eq. (10) in Table 1) in the anaerobic layer, while in nitri cation rate expression in the mass balance equations for NH4 þ the aerobic layer, NH is converted to nitrate ðNO Þ via nitrifi- is a product of the aerobic layer nitrification rate ðk þ; Þ and the 4 3 NH4 1 cation via a first order reaction rate ðk þ; Þ (Fig. 3a) with Monod depth of the aerobic layer (Fig. 3a): NH4 1
CH4(g) 2- RUFLUS SO (0) H S(0) ENAHTEM CH (aq)(0) 4 2 J 4 JPOC POC NMULOCRETAW NMULOCRETAW
KL01 KL01 2- SO4 (1) H2S(1) CH4(aq)(1) π 1 PARTITIONING: H2S(1) FeS(1) SOURCE: NONE κ H2Sd1 H S(1) + 2 O SO 2-(2) κ REACTIONS: 2 2 4 CH4,1 κ REACTION: CH (aq)(1) CO
LAYER 1 LAYER 4 2 AEROBIC H2Sp1 2- OXIDATION FeS(1) + 9/4O2 1/2Fe2O3 + SO4 (1) + 2H2O
FeS(1) SO 2-(1) H S(1) 4 2 CH4(aq)(1) ω 12 H1 KL12 KL12 SEDIMENT 2- FeS(2) SO4 (2) H2S(2) CH4(aq)(2)
DIAGENESIS: POC CH4(aq)(2) DIAGENESIS: METHANOGENESIS 2- 2CH2O + SO4 (2) 2CO2 + H2S(2) + 2H2O SULFATE REDUCTION CH (aq)(2) > CH (sat)
LAYER 2 LAYER 4 4
ANAEROBIC π 2 PARTITIONING: H2S(2) FeS(2) REACTION: CH4(aq)(2) CH4(g) ω BUBBLE FORMATION 2 a b H2
Fig. 4. Schematic representation of the sulfide model, with both solid and particulate phases and associated transport (Panel a), as well as the methane model which includes oxidation, diffusion, and bubble formation (Panel b). The oxygen uptake associated with these two models is added to nitrification-associated oxygen uptake to calculate sediment oxygen demand. D.C. Brady et al. / Estuarine, Coastal and Shelf Science 117 (2013) 107e124 111
D þ k þ; 2.4. Sulfate NH4 NH4 1 kNHþ;1H1 ¼ (2) 4 KL01 Sulfate reduction is an important diagenetic process in coastal The product D þ k þ; is made up of two coefficients, neither fi NH4 NH4 1 marine sediments and is the source of sul de, the oxidation of of which is readily measured. The diffusion coefficient in a milli- which is a large portion of SOD (Fig. 4a; Howarth and Jorgensen, meter thick layer of sediment at the sediment water interface may 2 2 1984). Thus, SO4 is explicitly modeled in SFM (Fig. 4a). SO4 in be larger than the diffusion coefficient in the bulk of the sediment the overlying water-column is calculated via a linear regression due to the effects of overlying water shear. It is therefore conve- 2 with overlying-water salinity. SO4 is produced in the aerobic layer nient to subsume two relatively unknown parameters into one as a result of particulate and dissolved sulfide oxidation and parameter that is calibrated to data called k þ; (Fig. 3a): NH4 1 removed in the anaerobic layer via sulfate reduction (Fig. 4a, Eq. qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 (16) in Table 1). Because anaerobic layer SO4 may become limiting k ; ¼ D þ k þ; (3) NH4 1 NH4 NH4 1 to sulfate reduction in low-salinity sediments, the model computes 2 a depth of SO4 penetration (HSO4 , Eq. (17) in Table 1). Since 2 which is termed the reaction velocity for NH4 oxidation, since its SO4 (2) ¼ 0atHSO4 one can solve for HSO4 using a one dimensional fi 2 dimensions are length per unit time. Nitri cation is then computed SO4 mass balance equation in the anaerobic layer (Eq. (17) in k þ Table 1; Chapter 11 in Di Toro (2001)). Because SO2 is often less as a function of NH4;1, temperature, and NH4 and O2 limitation (Eq. 4 2 2 (8) in Table 1). than the total sediment depth, KL12 is scaled to SO4 to yield a SO4 þ fi fi The mass balance NH4 equations for each layer are: speci c mass-transfer coef cient, KL12,S (Eq. (17) in Table 1). n h io þ 2 h i h 2.5. Sulfide d H1 NH4 ð1Þ k ; ¼ NH4 1 NHþð1Þ K NHþð1Þ dt K 4 L01 4 L01 i h i The product of sulfate reduction in the anaerobic layer is dis- þ þ þ fi NH4 ð0Þ þ KL12 NH4 ð2Þ NH4 ð1Þ solved sul de (H2S). Note that the carbon available for sulfate reduction is computed as the depositional flux of carbon minus the (4) n h i organic carbon utilized during denitrification (Eq. (16) in Table 1). þ h i 2 d H2 NH 2 fi 4 þ þ SO4 limitation on sul de production (via sulfate reduction) is ¼ KL12 NH4 ð2Þ NH4 ð1Þ þ JN2 (5) e 2 dt modeled via a Michaelis Menten dependency on SO4 (2) (Eq. (20) in Table 1). A portion of sulfide reacts with oxidized iron in the þ þ þ þ fi where NH4 (0), NH4 (1), and NH4 (2) are the NH4 concentrations in sediment to form particulate iron sul de (FeS in Fig. 4a; Morse the overlying water, layer 1 and layer 2, respectively. Note that for et al., 1987). SFM distinguishes between the solid and dissolved þ k fi fi fi NH4 , 2 ¼ 0 and JN1 ¼ 0. It should be noted that in the time varying phases of sul de via a layer-speci c partition coef cient (Eqs. (3) solution, H1 and H2 are within the derivative since the depth of the and (4) in Table 1, Table 2). The oxidation of both particulate and layers are variable and can entrain and lose mass when they are dissolved sulfide is the product of a first-order kinetic constant 2 changing (see Chapter 13 in Di Toro, 2001 for the method to (k ) and a linear dissolved O2 dependency (Eq. (18) in Table 1) H2S;d;p;1 account for these mass transfers). The depth of the anaerobic layer consistent with experimental data (Di Toro, 2001; Chapter 9). (H2) is then calculated as the difference between the total sediment e þ fl depth (10 cm) and H1. Finally, the sediment water NH4 ux 2.6. Methane þ þ þ (denoted as J½NH4 ) is equal to KL01ð½NH4 ð1Þ NH4 ð0Þ Þ. One important consequence of this equation is that the gradient in Methanogenesis is computed as a fraction of total diagenesis, 2 solute concentration between the sediment and the overlying which is controlled by anaerobic-layer SO4 (Eqs. (20) and (23) in water-column is not the sole control on solute flux during certain Table 1). For example, methanogenesis is equal to total diagenesis times of the year, as fluxes are primarily controlled by mixing via 2 when SO4 (2) is extremely low relative to the constant KM;SO4 (Eq. KL01. Mixing is primarily controlled by temperature, overlying- (23) in Table 1). The dissolved methane ([CH4] (aq)) produced by water dissolved [O2(0)], and the reactions that control the depth carbon diagenesis can either be oxidized in the aerobic layer via the of the aerobic layer. product of a first-order reaction constant and MichaeliseMenten O2 dependency (Fig. 4b, Eq. (21) in Table 1), released as a flux to 2.3. Nitrate flux the overlying water column, or becomes gaseous (J[CH4] (g)) when its concentration exceeds saturation ([CH4](sat); Fig. 4b and Eq. (24) There are two sources of NO3 in SFM: (1) NO3 enters from the in Table 1). J[CH4](g) is computed via a mass-balance of CH4(aq) overlying water column as controlled by surface mass transfer and sources and sinks when CH4(aq) CH4(sat) (Eq. (24) in Table 1). þ concentration gradient and (2) NH4 is oxidized in the aerobic layer (i.e., nitrification) (Eq. (13) in Table 1, Fig. 3b; see Chapter 4 of Di 2.7. Sediment oxygen demand Toro, 2001). In turn, NO3 can be returned to the overlying water fl column as NO3 ux (J[NO3]) or converted to nitrogen gas (i.e., Sediment oxygen demand (SOD) is computed by summing all denitrification, Eqs. (11) and (12) in Table 1). There is no biogeo- model processes that consume molecular O2, including the oxida- þ chemical nitrate source in the anaerobic layer (Fig. 3). Although tion of NH4 (nitrogenous SOD ¼ NSOD; Eqs. (28) and (29) in unconventional (see Gypens et al., 2008 for a discussion of models Table 1) and the oxidation of dissolved sulfide þ particulate sulfide fi fi fi of denitri cation con ned to the anaerobic layer), denitri cation (SODH2S) and dissolved methane (SODCH4 ). Carbonaceous SOD e occurs in both the aerobic and anaerobic layers (Fig. 3b). The close (CSOD) is equivalent to SODH2S þ SODCH4 (Eqs. (25) (27) in coupling between nitrification and denitrification has been sug- Table 1). In the case of the sulfur and methane models, all solutes gested by Blackburn et al. (1994) and there is considerable evidence are modeled in units of stoichiometric O2 equivalents (O2*, Tables 1 for denitrification in the oxic layer, for example, within anoxic and 2). This model formulation allows the model to track reduced microsites (Jenkins and Kemp, 1984; Berg et al., 1998). The end-products (e.g., sulfide) in the summer for re-oxidization during þ fl remainder of the formulation parallels the NH4 ux model dis- the winter, a potentially important process when determining the cussed in Section 3.2 and is diagrammed in Fig. 3b. lag between loading and sediment recovery. 112 D.C. Brady et al. / Estuarine, Coastal and Shelf Science 117 (2013) 107e124
Table 1 The model equations are listed below. The solutions are found by numerically integrating the equations (Di Toro, 2001). Parameter definitions are located at the bottom of this table.a
Mass balance equations (these are the general equations for CT1 and CT2):
k2 dfH1CT1g 1 u u 1 ¼ CT1 þ KL01ðfd0CT0 fd1CT1Þþ 12ðfp2CT2 fp1CT1ÞþKL12ðfd2CT2 fd1CT1Þ 2CT1 þ JT1 dt KL01 dfH C g 2 2 T2 ¼ k C u ðf C f C Þ K ðf C f C Þþu ðC C ÞþJ dt 2 T2 12 p2 T2 p1 T1 L12 d2 T2 d1 T1 2 T1 T2 T2 1 ; 3 where: fdi ¼ p i ¼ 1 2 1 þ mi i
4 fpi ¼ 1 fdi
DO2 SOD 5 KL01 ¼ ¼ H1 ½O2ð0Þ
D qðT 20Þ d Dd 6 KL12 ¼ ðH1 þ H2Þ=2
ðT 20Þ Dpq Dp POC1 u ¼ ð1 ksSÞ min 12 H þ H POC ðeach yearÞ 7 1 2 R
dS KM;Dp where ¼ ksS þ = dt KM;Dp þ½O2ð0Þ 2
The kinetic and source terms for the ammonium, nitrate, sulfide, methane, and oxygen systems are listed next.
Ammonium ðNHþÞ 4 0 1 ðT 20Þ K ; q M NH4 KM;NH ½O ð0Þ =2 k2 k2 qðT 20Þ@ 4 A 2 8 1 ¼ NH ;1 NH 4 4 ðT 20Þ K ; ; O 0 =2 K ; q NH 1 M NH4 O2 þ½ 2ð Þ M NH4 K ; þ½ 4ð Þ M NH4 k 9 2 ¼ 0
10 JT1 ¼ 0 X2 qðT 20Þ Jc ¼ kPOC;i POC;i POCiH2 i ¼ 1 JT2 ¼ aN;C Jc