Effect of benthic boundary layer transport on the productivity of Mono Lake, California
© Bruce et al. 2008
Received: 13 December 2005
Accepted: 19 August 2008
Published: 19 August 2008
Skip to main content
© Bruce et al. 2008
Received: 13 December 2005
Accepted: 19 August 2008
Published: 19 August 2008
The significance of the transport of nutrient-rich hypolimnetic water via the benthic boundary layer (BBL) to the productivity of Mono Lake was studied using a coupled hydrodynamic and ecological model validated against field data. The coupled model enabled us to differentiate between the role of biotic components and hydrodynamic forcing on the internal recycling of nutrients necessary to sustain primary productivity. A 4-year period (1991–1994) was simulated in which recycled nutrients from zooplankton excretion and bacterially-mediated mineralization exceeded sediment fluxes as the dominant source for primary productivity. Model outputs indicated that BBL transport was responsible for a 53% increase in the flux of hypolimnetic ammonium to the photic zone during stratification with an increase in primary production of 6% and secondary production of 5%. Although the estimated impact of BBL transport on the productivity of Mono Lake was not large, significant nutrient fluxes were simulated during periods when BBL transport was most active.
The transport of nutrient-rich water from benthic to pelagic regions has been linked to increased levels of primary productivity in stratified lakes [1–3]. Ostrovsky et al.. (1996) suggest that seiche activity in the boundary layer of Lake Kinneret sustained a vertical flux between the hypolimnetic and epilimnetic waters enhancing biological productivity in the lake. MacIntyre et al. (1999) calculated the upward fluxes of ammonium across the nutricline in Mono Lake and suggested nearshore boundary fluxes could be the dominant pathway supplying ammonium to the deep chlorophyll maxima. Eckert et al. (2002) used microstructure measurements of temperature, oxygen and hydrogen sulphide in Lake Kinneret to conclude that following the onset of stratification, the flux of benthic nutrients to the water column controls primary productivity. In this study we have defined BBL transport as that which occurs in the layer bordering the sediments of a lake [4, 5] alternatively referred to as the bottom boundary layer .
The development of basin-scale internal waves arising from wind-induced energy are responsible for large scale water motions and most of the turbulence caused by these large-scale motions occurs in the BBL [7, 8]. In order to differentiate between boundary and internal modes of vertical transport Yeates and Imberger (2004) parameterized the split between mixing in the internal and benthic boundary layer (BBL) using values of Lake number, LN  and Burger number, BN . The LN is a measure of the amplitude of basin-scale internal waves in response to surface wind forcing, and BN describes waves that evolve from simple seiches . Simulations performed on a number of monomictic lakes indicated that fluxes through the BBL were dominant during strong wind events occurring during period of stratification .
A number of studies, aimed at identifying sources and sinks of nutrients in the photic zone have focused on bacterial mineralization , regeneration through planktonic organisms [11–13], nitrogen fixation , hypolimnetic flux and inflows and outflows . Although the occurrence of BBL transport and its potential impact on primary productivity has been examined, the upward mixing of nutrient-rich hypolimnetic waters via the BBL and the consequent effect on lake-wide ecological processes deserves further analysis.
Mono Lake is a nitrogen-limited saline lake with a relatively simple food web  and is subjected to wind-driven boundary-layer mixing events . Yeates and Imberger (2004) simulated a BBL thickness in Mono Lake of 10–15 m during a sequence of strong wind event suggesting an active role in the development of the thermal structure of the lake. These features make it well-suited for examining the role of BBL-supplied nutrients and the influence of these nutrients on the seasonal plankton dynamics and overall productivity of the lake.
The objective of the present study is to investigate the role of BBL transport in the supply of nutrients to the photic zone and its consequent impact on the lake's ecology. A coupled hydrodynamic and ecological model was used to quantify nitrogen biogeochemistry during a 4-yr period from 1991–1994 when the lake mixed to the bottom during the winter. Initially, we calibrated the model parameters and processes to ensure an acceptable representation of the field data. The simulated output was then used to calculate the sources and sinks of nitrogen to the photic zone. A comparison could then be made between the roles of recycled and external sources on the primary and secondary productivity in the lake. To enable quantification of the significance of BBL transport for ecological processes, a series of simulations were run in which this mechanism was switched off allowing a comparison between lake behavior with and without BBL transport.
The planktonic community of Mono Lake has few species as is typical of hypersaline waters. The phytoplankton is dominated by a newly described picoplanktonic (2–3 μm) green alga, Picocystis salinarum Lewin (Lewin et al., 2000), and several bacillarophytes, mainly Nitzschia spp. (20–30 μm) (Lovejoy & Dana, 1977; Mason, 1967). A brine shrimp, Artemia monica Verill, is the only macrozooplankter (Lenz, 1980; Lenz, 1984). While pelagic ciliates and rotifers may also be present at times (Mason, 1967; Jellison et al. 2001), they contribute a negligible amount to the total zooplankton biomass.
There is a strong seasonal pattern in the nutrient and plankton dynamics of Mono Lake . The seasonal patterns are driven by biotic and abiotic forces affecting productivity via bottom-up and top-down controls. Water temperatures of the surface mixed-layer ranged from 2–5°C in winter to 12–22°C in summer. Seasonal stratification and high productivity result in anoxic conditions in the hypolimnion where ammonium accumulates. The flux of this ammonium to the photic zone is limited until winter overturn mixes the whole lake providing nutrients for a pronounced spring algal bloom. Daily primary productivity rates are relatively high (Jellison and Melack 1993a).
The lake's only macrozooplankter, A. monica, produces over-wintering cysts that lie dormant on the bottom during the winter and hatch during early spring (February-April) . A. monica biomass usually peaks in the late spring, remains high during the summer and gradually declines during the autumn as food is scarce and temperatures decline. The spring growth of A. monica biomass is associated with a simultaneous decline in phytoplankton biomass due to grazing and rise in surface concentrations of ammonium from zooplankton excretion. Phytoplankton biomass remains low during the summer and only increases toward the end of the year when grazing pressure is reduced .
As phosphorus concentrations are always high (>400 μM; Jellison et al. 1993), nitrogen limits primary production in the photic zone (Jellison & Melack 1993a, 2001). Nitrogen inputs from inflowing streams and planktonic nitrogen fixation are very low relative to internal fluxes where the main sources are from sediment release in the hypolimnion, phytoplankton and zooplankton excretion, and bacterial mineralization of particulate detrital organic nitrogen. Peak concentrations in the photic zone are observed at the breakdown of stratification as nutrient-rich hypolimnetic waters become entrained into the epilimnion. Towards the end of the mixed period and onset of stratification ammonium levels are generally low. When the zooplankton become abundant in late spring, grazing reduces phytoplankton biomass and internal phytoplankton nitrogen is converted to ammonium via the zooplankton grazing and excretion. Zooplankton excretion and reduced ammonium uptake due to low phytoplankton biomass results in an increase in epilimnetic ammonium concentrations.
The simulated values of phytoplankton biomass follow the low summer concentrations, timing and slope of the autumn recovery and spring decline observed in the field data (Fig. 3). However, the elevated values of phytoplankton biomass observed in the field at the end of the mixing periods (early 1992, 1993 and 1994) are not captured by the model (Fig. 3).
The particulate organic nitrogen (PON) and particulate organic carbon (POC) data observed in the field closely followed that of the phytoplankton, with elevated values during the winter in the absence of grazing and low values during the summer months. These patterns were captured by the model, although elevated levels of PON were underestimated by the model during periods of full circulation (Fig. 3). However, elevated levels of POC were generally captured by the model which suggests that the model overestimated the detrital component of the particulate carbon pool (Fig. 3).
The timing and slope of the early spring peak in A. monica biomass observed in the field was matched in the simulated results across the four year simulation period (Fig. 3).
Simulated concentrations of dissolved oxygen are similar to those measured in the field during stratified periods (Fig. 3). The model, however, under predicted concentrations in late spring for both 1991 and 1992.
Primary productivity in Mono Lake has been estimated using a numerical interpolative model incorporating photosynthetic uptake rates and measured vertical attenuation of PAR . During the non-meromictic conditions of 1989 and 1990, Jellison and Melack (1993a) estimated an average daily productivity of 1.6 g C m-2 d-1. This matches the value simulated by DYRESM-CAEDYM for the 1991–1994 monomictic period. During periods of stratification an average daily productivity of 1.7 g C m-2 d-1 was simulated and 1.3 g C m-2 d-1 during periods of full circulation.
Average rates of lake-wide nitrogen deposition measured in 1986 and 1987 ranged from approximately 5.9 Mg N d-1 (ca. 2.5 mmol m-2 d-1) during the summer to 2.7 Mg N d-1 (ca. 1.2 mmol m-2 d-1) during the winter (Jellison et al. 1993). These rates are similar to those simulated by the model, i.e., 3.7 Mg N d-1 and 2.1 Mg N d-1 averaged during periods of stratification and full circulation, respectively. Areal average lake-wide nitrogen fluxes from the sediments were calculated by the model as 12.5 Mg N m-2 d-1 and 6.2 Mg N m-2 d-1 averaged during period of stratification and full circulation, respectively. Jellison et al. (1993) estimated the rate of ammonia release from the sediments based on sediment cores collected in 1988 as 58–162 Mg N m-2 d-1 (ca. 3.6–10.1 mmol m-2 d-1). Although greater than those predicted by the model these estimates were derived under anoxic conditions so should be taken as an upper estimate.
The calculated values of normalized mean absolute error, correlation coefficient and slope are presented in Table 3 for each of the main state variables over the full simulation period from 1991 to 1994 and compared to the calibration period from 1991 to 1992. Calculations of correlation coefficients are all equal to or greater than 0.8 with the exception of ammonium and dissolved oxygen.
The five parameters that displayed the greatest sensitivity to annual estimates of lake-wide nitrogen fluxes were: (1) release rate of NH4 from sediments (SdNH4); (2) the fraction of zooplankton grazing excreted (fex); (3) background attenuation coefficient (Kd); (4) internal nitrogen to carbon ratio of the phytoplankton (INcon) and (5) the fraction of zooplankton grazing egested (feg). The optimal parameter value (determined by the model best fit), and the upper and lower bounds used to determine alternative parameter sets are listed in Table 4.
Calculations based on model output indicate that for 1991 to 1994 BBL transport was responsible for a 53% increase in upwards flux of ammonium across the thermocline during periods of stratification. For the corresponding periods, the simulated increase in primary production was calculated as 6% and secondary production as 5%. The model results, averaged over periods of autumn and winter mixing for the 4 years of simulation, indicated a reduction in upward ammonium flux of 28% when the BBL transport was active. This corresponded with a simulated decrease in primary production of 7% and negligible increase in secondary production of 1% for the same periods. The estimated net increase for 1991–1992 in ammonium flux across the thermocline due to BBL transport was 9%, primary productivity was 2% and secondary productivity was 3%.
Several aspects of the modeling require further examination. The step temperature function used to represent the process responsible for the hatching and initial growth of over-wintering A. monica cysts simulated well the timing and slope of the early spring peak in A. monica biomass. However, experiments have demonstrated that increases in salinity can influence the hatching process . It is anticipated, therefore, that an additional salinity factor would be required before the model could be used to predict A. monica dynamics under alternative salinities. Although simulated mid-summer concentrations of A. monica compare favorably with those observed in the field, the autumn decline was difficult to simulate well (Fig. 3). The model included three processes responsible for decreases in biomass during this period; limited grazing at low temperatures, end of life mortality and grebe predation. Improved understanding of the combination of triggers responsible for the autumn decline in A. monica will aid in the model representation of these processes. Alternatively a cohort model such as that proposed by  may be required to accurately represent the autumn decline.
Differences between measures of fit comparing the calibration and validation periods are small. Although the ecological dynamics of the model during the validation period are similar to that of the calibration period, this result is an indication of model stability. However, this stability only relates to the representation of the interactions between the main processes responsible for determining the ecological patterns observed in the lake over the four years studied. Comparison to measures of fit for other lake ecosystem models is difficult as quantitative measures are rarely given. However, our overall measurement of NMAE compare favorably to Ross et al. (1994) (0.65) and Bruce et al. (2006) (0.52).
Since we defined sensitivity in relation to estimates of lake-wide nitrogen fluxes, it follows that the parameters showing the most sensitivity are related to the nitrogen cycle. Since the inflow of nitrogen into the lake is negligible, it follows that for Mono Lake, sediment release is a critical source of nitrogen to the water column. Similarly both the fractions of zooplankton grazing that goes into either egestion (the bulk of which is deposited into sediments and thus lost from the photic zone) or excretion (providing nitrogen in a form for primary production) have a direct effect on the proportion of phytoplankton nitrogen that is recycled. The ratio of phytoplankton internal nitrogen to carbon controls both the uptake of inorganic nitrogen by phytoplankton and the flux of nitrogen recycled via the zooplankton grazing and excretion pathway. Background attenuation influences nitrogen fluxes indirectly by controlling the amount of light available for primary productivity.
Model results indicated that during the summer stratified periods the N demand by phytoplankton in the surface to 9-m of Mono Lake is predominantly met by zooplankton excretion, phytoplankton leakage of dissolved organics, and bacterially mediated mineralization. Of these, the model predicted that the dominant source was zooplankton excretion. Since zooplankton biomass was well represented by the model including timing and magnitude of the initial peak, it follows that during these peaks the model has the closest fit to the ammonium data. Midwinter spikes in ammonium during period of reduced phytoplankton biomass were not reproduced in the model output. The model simulated almost constant phytoplankton biomass during the winter months that is inconsistent with the field data. From this we would conclude that the processes of phytoplankton ammonium uptake and release are not well represented by the model during winter conditions of high algal biomass and light-limitation. One of the limitations of this study was the assumption (for simplicity) of a constant internal phytoplankton C:N ratio. It is anticipated that modeling the internal nitrogen as a dynamic variable would improve the simulation of the phytoplankton-ammonium interactions particularly during periods of full circulation.
This study employed optimization techniques to determine a series of parameter sets to best represent field data as described by the processes included in the current model formulation. Some field data were better represented than others and misrepresentation of field data by the simulation will serve to direct improvements in future model generations. Although the modeled fluxes sometimes over or underestimated the measured concentrations, the general seasonal patterns were captured by the simulations and thus used to provide insight into the processes that determine the ecosystem dynamics of Mono Lake.
Bruce et al. (2006) in their study of the role of zooplankton in the nutrient cycles of Lake Kinneret, Israel, found zooplankton excretion to be the dominant source of dissolved nitrogen during winter overturn and sediment release the dominant source during summer stratification. For Mono Lake we also found that zooplankton excretion was most influential in the summer stratified period. Although the simulated rate of ammonium flux from the sediments was higher in Lake Kinneret , the main reason for finding sediment-released nutrients relatively less important in Mono Lake is due to two-fold higher rates of primary productivity and greater recycling due to zooplankton excretion in Mono Lake.
MacIntyre and Jellison (2001) suggested that transport of nutrient-rich hypolimnetic water via the BBL layer is responsible for increased ammonium flux across the thermocline and consequential increase in productivity. By comparing the simulation results from the two scenarios we found that, although the increase in upward ammonium flux across the thermocline during the stratified periods of 1991–1994 due to BBL transport was 53% (± 4%), primary productivity for the same period increased only 6% (± 4%). Since the model suggested that 87% of the N demand by phytoplankton is met by regenerated sources, it is not unexpected that an increase in external supply has a limited impact. MacIntyre et al. (1999) highlighted the importance of the flux of BBL transported ammonium across the thermocline in sustaining primary productivity to the deep chlorophyll maxima. As a percentage of phytoplankton demand during the stratified periods the upward flux of ammonium across the thermocline was calculated as 12% with BBL on and 5% with BBL off. MacIntyre et al. (1999) reached a similar conclusion and, assuming that 5–10% of primary productivity occurs in the deep chlorophyll maximum during the summer, suggested that BBL may be the dominant mechanism supplying ammonium to the deep chlorophyll maximum.
As anticipated, during stratified periods simulation results indicate that BBL transport leads to an increase in ammonium transport across the thermocline and concomitant increase in primary productivity. However, this pattern was reversed during periods of mixing. A greater build up of ammonium in the hypolimnion occurred during stratification in the case where BBL transport is absent (Fig. 6). Although in the absence of BBL transport, less flux was available in the photic zone during stratification, at overturn a greater mass of ammonium led to greater upwards flux of ammonium and concomitant increase in primary productivity. As a result, on an annual average, primary productivity was similar under both scenarios.
Our model results have illustrated the importance of timing of BBL transport and its subsequent effect on primary and secondary productivity. The model used in this study did not include algorithms to represent differences in generations using a stage-structured zooplankton model. Inclusion of a stage structured model might enable us to determine whether the timing of BBL transport events and concomitant increases in primary productivity effect the timing and magnitude of successive generations of A. monica in Mono Lake.
For Lake Kinneret, estimates of monthly primary productivity fall between 0.5 and 1.7 g C m-2 d-1 (Bruce et al. 2006). Mean daily values of LN estimated for Lake Kinneret range from 10-2 to 100 with a period of low LN associated with strong wind events (Yeates and Imberger 2004). Given these ranges, it is predicted that the importance of BBL transport in Lake Kinneret may be greater than the 6% predicted for Mono Lake. In Lake Geneva, a study investigating the effect of internal waves on basin exchange indicated up to 40% of the hypolimnetic volume was exchanged following episodes of strong winds . Primary productivity in Lake Geneva is relatively high . For Lake Constance, estimates of LN during stratification are relatively high (Yeates and Imberger 2004) and productivity is less than 1 g C m-2 d-1  suggesting that for Lake Constance the transport of nutrients through the BBL may be less important to overall lake productivity but potentially significant during episodic events associated with low LN.
The results of this study have indicated that the relative importance of BBL transport as a source of nutrients sustaining productivity in the photic zone is determined by productivity and morphology. Future studies will be focused on comparing the effect of BBL transport on the ecology of other lakes. By differentiating between physical and ecological process we will be able to determine what limnological features alter the importance of BBL transport.
The model used in this study is a modified version of the Computational Aquatic Ecosystem Dynamics Model (CAEDYM) [28, 29] coupled to the Dynamic Reservoir Model (DYRESM) . In DYRESM the lake is represented as a series of homogeneous horizontal layers of variable thickness ; as inflows and outflows enter or leave the lake, the affected layers expand or contract, respectively, and those above move up or down to accommodate the volume change. Mass, including that of the ecological state variables, is adjusted conservatively each time layers expand, contract, merge or are affected by inflows and outflows. The main processes modelled in DYRESM are surface heat, mass and momentum transfers, mixed layer dynamics, hypolimnetic mixing, benthic boundary layer mixing, inflows and outflows.
where BN is the Burger number .
Model process equations. Equations used to describe the processes included in the ecological model CAEDYM
∂Zi/∂t = [GiAi f(Z)i f 1(T)(1-fex-feg) - (Ri+Mi)f 2(T) - Predi]Zi
= (assimilation - excretion - egestion) - (respiration + mortality) - predation
∂P/∂t = [Pmax,j f 1(T)min(f(I),f(P),f(N)) - (Rj)f 2(T) - Predj]Pj ± Sj
= photosynthetic uptake - (respiration + excretion + mortality) - predation ± settling
∂POC/∂t = Σ[Gi f(Z)i f 1(T)i((1-Ai) + Aifeg) + Mi f 2(T)i]Zi + Σ[Rj(1-fres)(1-fDOM)f 2(T)]Pj - PredPOCPOC - RPOC f(DO)f 1(T)POC ± SPOM
= (unassimilated zooplankton food + zooplankton egestion + zooplankton mortality) + phytoplankton mortality - zooplankton predation - POC decomposition ± settling
∂DOC/∂t = Σ[Rj(1-fres)fDOM f 2(T)]Pj + RPOC f POC(POC)f (DO)f 2(T)POC - RDOC f (DO)f 2(T)DOC
= phytoplankton excretion + POC decomposition - DOC mineralisation
∂POP/∂t = Σ[Gi f(Z)i f 1(T)i((1-Ai) + Aifeg) + Mi f 2(T)i]IPZiZi + Σ[Rj(1-fres)(1-fDOM)f 2(T)]IPj - PredPOCPOP - RPOP f(DO)f 1(T)POP ± SPOM
= (unassimilated zooplankton food + zooplankton egestion + zooplankton mortality) + phytoplankton mortality - zooplankton predation - POP decomposition ± settling
∂DOP/∂t = Σ[Rj(1-fres)fDOM f 2(T)]IPj + Σ[AifexGi f(Z)i f 1(T)i]IPZiZi + RPOP f (DO)f 1(T)POP - RDOP f (DO)f 1(T)DOP
= phytoplankton release + zooplankton excretion + POP decomposition - DOP mineralisation
∂PO4/∂t = RDOP f (DO)f 2(T)DOP - Σ[UNmax,j f 1(T)j f(IP)j f(P)j]Pj + SdPO4 f(DO)f 2(T)LA/LV
= DOP mineralisation - phytoplankton uptake + PO4 sediment flux
∂PON/∂t = Σ[Gi f(Z)i f 1(T)i((1-Ai) + Aifeg) + Mi f 2(T)i]INZiZi + Σ[Rj(1-fres)(1-fDOM)f 2(T)]INj - PredPONPON - RPON f(DO)f 1(T)PON ± SPOM
= (unassimilated zooplankton food + zooplankton egestion + zooplankton mortality) + phytoplankton mortality - zooplankton predation - PON decomposition ± settling
∂DON/∂t = Σ[Rj(1-fres)fDOM f 2(T)]INj + Σ[AifexGi f(Z)i f 1(T)i]INZiZi + RPON f(DO)f 2(T)PON - RDON f(DO)f 2(T)DON
= phytoplankton release + zooplankton excretion + PON decomposition - DON mineralisation
∂NH4/∂t = RDON f(DO)f 1(T)DON - Σ[UNmax,jPN f 1(T)j f(IN)j f(N)j]Pj - RNO f(DO)f 2(T)NH4 + SdNH4 f(DO)f 2(T)LA/LV
= PON mineralisation - phytoplankton uptake - nitrification + NH4 sediment flux
∂NO3/∂t = RNO f(DO)f 2(T)NH4 - RN2 f(DO)f 2(T)NO3 - Σ[UNmax,j(1-PN)f 1(T)j f(IN)j f(N)j]Pj
= nitrification - denitrification - phytoplankton uptake
∂DO/∂t = kO2(DO_atm - DO) + Σ[Pmax,j f 1(T)jmin(f(I),f(P),f(N)) - Rj f 2(T)j]PjYO2:C Σ[Ri f 2(T)i]ZiYO2:C - RDOC f (DO)f 1(T)DOCYO2:C - RNO f (DO)f 2(T)NH4 - SdO2 f(DO)f 2(T)LA/LV
= atmospheric flux + (phytoplankton oxygen production - phytoplankton respiratory consumption) - zooplankton respiratory consumption - utilisation of oxygen in mineralisation of DOM - utilisation of oxygen in nitrification - sediment oxygen demand.
f 1(T) = θT-20 - θk(T-a) + b
where k, a and b are constants solved numerically to satisfy the following conditions:
f 1(T) = 1; at T = Tsta
∂f 1(T)/∂T = 0; at T = Topt
f 1(T) = 0; at T = Tmax
f 2(T) = θT-20
f(I)j = I/Is exp(1-I/Is)
f(IP)j = [IPmax/(IPmax-IPmin)] [1-IPmin/IP]
f(IN)j = [INmax/(INmax-INmin)] [1-INmin/IN]
f(DO) = DO/(KDO+DO)
f(P) = PO4/(KPO4+PO4)
f(N) = (NH4+NO3)/(KN2+NH4+NO3)
PN = (NH4 NO3)/[(NH4+KN)(NO3+KN)] + (NH4 KN)/[(NH4+KN)(NO3+KN)]
Sj = (ws/Δz)Pj
SPOM = (g(ρPOM - ρw)(DPOM)2/18μ)/Δz)POM
Predi = Σ(Gkf(Z)kf1(T)kZkPzZOOk,i)
Predj = Σ(Gif(Z)if1(T)iZiPzPHYi,j)
Abbreviations: Z, zooplankton; P, phytoplankton; POC, particulate organic carbon; DOC, dissolved organic carbon; POP, particulate organic phosphorus; PO4, phosphate; PON, particulate organic nitrogen; NH4, ammonium; NO3, nitrate; POM, particulate organic matter (C, N or P); IPzi, zooplankton internal phosphorus; INzi, zooplankton internal nitrogen; IPj, phytoplankton internal phosphorus; INj, phytoplankton internal nitrogen; DO, dissolved oxygen; DOatm, concentration of oxygen in the atmosphere; LA, layer area; LV, layer volume; Δz, layer thickness; ρw, density of water; μ, viscosity of water; kO2, oxygen transfer coefficient. Subscripts: i, zooplankton group; j, phytoplankton group; k, zooplankton predator group.
The major nutrient fluxes represented in CAEDYM are uptake of dissolved inorganic nutrients by phytoplankton, release of dissolved nutrients from phytoplankton excretion, grazing, egestion and excretion of nutrients by zooplankton, nitrification and denitrification of inorganic nitrogen, sedimentation of nutrients in particulate form, mineralization of organic nutrients and release of dissolved nutrients from sediments (Table 1).
Net change in carbon concentration of the phytoplankton at each model time step is calculated as the difference between the increment due to gross primary production and losses due to sedimentation, grazing by zooplankton, respiration, excretion and mortality. These terms are calculated using equations parameterized to represent the physiology of the main phytoplankton species. Losses due to grazing by zooplankton are calculated by multiplying the food assimilation rate for the zooplankton by a preference factor for phytoplankton over detrital POC.
Net zooplankton growth is calculated as a balance between food assimilation and losses from respiration, excretion, egestion, predation and mortality. Food assimilation is calculated as the product of the maximum potential rate of grazing, assimilation efficiency, and temperature and food concentration functions. A constant internal nutrient ratio is assumed and excretion of nutrients calculated to maintain this ratio at each time step. Advective movement of zooplankton is carried out in DYRESM.
Bacteria have not been directly simulated as they were not measured during the study period. However the nutrient pathways mediated by bacteria were included as mineralization of the particulate organic pools (POC, POP and PON). The POC, POP and PON pools available for zooplankton grazing include bacteria. Predation of zooplankton by grebes was included by an additional predation term for the months of August to November estimated from predation studies (Cooper et al. 1984).
The advantage of using a depth resolved model DYRESM linked to the ecological model CAEDYM is that we could explore the effect of transport and mixing between the epilimnion, metalimnion and hypolimnion on the ecological processes in the lake. Of most relevance to this study is the exchange of nutrient-rich hypolimnetic waters to the photic zone via the BBL and its consequential effect on the primary productivity. The separation of internal and BBL cells in the layered structure of the current DYRESM allowed us to differentiate between the transport of nutrients in the internal and BBL and to determine the relative importance of each process on the mixed layer ecological dynamics. The ecological model also returns the attenuation coefficient (as a function of the concentration of both phytoplankton and particulate organic matter) to the hydrodynamic model at each one-hour time step. This variable is used to determine the extent of light and heat penetration that in turn governs the deepening of the surface mixed layer and the timing of winter turnover. In this way the feedback on a sub-daily time scale between the ecological and physical models is instrumental in the application of the model to aid in understanding of the interaction of various lake processes.
Seasonal and year-to-year variations in the physical, chemical, and biotic environments were monitored fortnightly from March through October and monthly during November through January. Water temperature and conductivity were measured at nine buoyed, pelagic stations (2, 3, 4, 5, 6, 7, 8, 10 and 12) (Fig. 1). Profiles were taken with a high-precision, conductivity-temperature-depth profiler (CTD) (Seabird Electronics model Seacat 19) equipped with a submersible photosynthetically available radiation (PAR) (LiCor 191S), fluorescence (695 nm) (WETLabs WETStar miniature fluorometer), and transmissivity (660 nm) (WETlabs C-Star Transmissometer). Specific conductivity, salinity, and density were all calculated based on equations derived from measurements on Mono Lake brine . Dissolved oxygen was measured at one centrally located station (Station 6) with a Yellow Springs Instruments temperature-oxygen meter (YSI, model 58) and probe (YSI, model 5739). The oxygen electrode was calibrated at least once each year against Miller titrations of Mono Lake water (Walker et al. 1970).
Ammonium and chlorophyll profiles were determined by sampling 7–10 discrete depths at two pelagic stations (2 and 7; Fig. 1), while A. monica abundance was determined via vertical net tows collected at 10 (1991–1992) or 20 (1993–1994) pelagic stations (Fig. 1). Chlorophyll a was also determined in the upper water column from samples collected with a 9-m integrating tube sampler at 5 pelagic stations (2, 6, 7, 10, 11; Fig. 1).
Nutrient and phytoplankton samples were immediately passed through a 120-μm net to remove all stages of A. monica and a sub-sample filtered through Gelman A/E glass fiber filters for analysis of nutrients (Jellison and Melack 1993a). Ammonium concentrations were measured with the indophenol blue method as described by Jellison et al. (1993). Nitrate and nitrite concentrations were measured but were always low (< 1 μM) [21, 31] and thus not considered in this study. Phosphate concentrations are orders of magnitude greater than the half saturations constants for phytoplankton so were not considered in this study . Phytoplankton chlorophyll a was determined by spectrophotometric analysis as described by Jellison and Melack (1993a). Conversion of chlorophyll to carbon units were made by assuming a C:Chl a ration of 50 (see Jellison & Melack 2001). Subsamples were filtered onto precombusted Gelman A/E filters for the determination of particulate organic carbon (POC) and nitrogen (PON). Duplicate carbon and nitrogen filters were acid fumed for 12 hours over concentrated HCl, and then dried at 40–50°C before determination by combustion in a Perkin-Elmer 240B elemental analyzer standardized with acetanilide. A. monica were collected using vertical net tows (120-μm mesh) to within 1-m of the bottom or well below the oxycline depending on stratification. A. monica biomass (dry weight) was estimated from stage-specific abundance, adult female length data, and weight-length relationship determined in the laboratory simulating in situ conditions of food and temperature . Conversion from dry weight to carbon assumed 0.4 g C/g dry weight [33, 34].
Model input files included data for initialization, meteorology, inflows and outflows. The initialization file was prepared from field data collected on 13 January 1991. On this day the temperature of the lake ranged from 2.5 to 3°C and the salinity from 88 to 88.5 g kg-1. Inflow data included the daily volume, temperature and salinity for two inflows, one representing total surface inflows (streams and direct runoff) and the other, hydrothermal springs. The volume of the hydrothermal springs was set at 3888 m3 day-1, based on a 3He mass balance of Mono Lake . The surface inflows were calculated based on a water mass balance using measured values of water depth and evaporation calculated by DYRESM. As ammonium, phytoplankton and zooplankton concentrations are negligible in the inflows, they were set to zero (Jellison and Melack 2001). Meteorological input data included hourly short- and long-wave radiation, air temperature, vapor pressure, wind speed and precipitation . Air temperature, vapor pressure (converted from relative humidity), wind speed and precipitation were collected at a meteorological station located on Paoha, a central island (Fig. 1). Radiation data were collected from a meteorological station located approximately 7 km from the southwest shore of the lake (Fig. 1).
The physical parameters used to simulate the hydrodynamics of Mono Lake were either physical constants or ones fixed according to the dimensions of the lake .
Model parameters. Parameters used in CAEDYM to simulate ecological variables in Mono Lake.
Values from field/lit
Background extinction coefficient
a Calculated from unpub data on in-situ light measurements
Values from field/literature
Maximum potential growth rate
Parameter for initial slope of PI curve
Specific attenuation coefficient
m2 g C-1
Half saturation constant for phosphorus uptake
Low value as not P limited
Half saturation constant for nitrogen uptake
Constant internal N ratio
mg N (mg C)-1
Constant internal P ratio
mg P (mg C)-1
Temperature multiplier for growth
Metabolic loss rate coefficient
Temperature multiplier for metabolic loss
Fraction of respiration relative to total metabolic loss
Fraction of metabolic loss rate that goes to DOM
aJellison and Melack 1993a, based on maximum value of carbon uptake measured from lake samples 1983–1990 assuming 50 g C g Chl a -1
b Jellison and Melack 1993a, based on minimum value of IK measured from lake samples 1983–1990.
cJellison and Melack 1993a.
dJellison and Melack 2001, estimated from seston ratios during the summer period from monomictic years 1991–1995 1984
eJellison and Melack 1993a, based on Q10 of 1.95.
fJellison et al.. 1993.
Values from field/literature
g C m-3 (g C m-3)-1 d-1
Close to 1 as filter feeders
Respiration rate coefficient
Mortality rate coefficient
Fecal pellet fraction of grazing
Kfz+kez = 0.36–0.68e
Excretion fraction of grazing
Minimum DO tolerance
Temperature multiplier for growth
Respiration temperature dependence
Half saturation constant for grazing
g C m-3
Internal ratio of nitrogen to carbon.
g N g C-1
Internal ratio of phosphorus to carbon
g P g C-1
Preference of zooplankton for phytoplankton
Preference of zooplankton for POC
cJellison et al.. 1993 (based on survival rate over 30 days)
dDana and Lenz 1986 (based on survival rate over 26 days)
eEvjemo et al.. 2000, Artemia fransiscana.
fDO concentration at depth of deep Chl a maxima (unpub data).
gJellison et al.. 1993 (best fit to temperature function used in model)
h Jellison unpub data 1991–1994. (based on temperature at which total biomass < 0.01 before Spring growth)
i Evjemo and Olsen 1999 (Artemia fransiscana 11 days old, 26–28°C, Holling Type II)
j Jellison unpub data (Females/Males)
Dissolved Oxygen and Nutrients
Values from field/literature
DO sediment exchange rate
Half saturation constant for DO sediment flux
mg O L-1
Half saturation constant for dependence of POM/DOM decomposition on DO
mg O L-1
Mineralisation rate for POC to DOC
Mineralisation rate for POP to DOP
Mineralisation rate for PON to DON
Diameter of POM particles
Density of POM particles
Specific light attenuation coefficient for POC
Mineralisation rate for DOC
Set to 1 to eliminate DOP pool for simplicity
Mineralisation rate for DOP to PO4
Set to 1 to eliminate DOP pool for simplicity
Mineralisation rate for DOP to PO4
Set to 1 to eliminate DOP pool for simplicity
Mineralisation rate for DON to NH4
onset to 1 to eliminate DON pool for simplicity.
Specific light attenuation coefficient of DOC
Denitrification rate coefficient
Temperature multiplier for denitrification
Half saturation constant for denitrification dependence on oxygen
mg N L-1
Nitrification rate coefficient
Temperature multiplier for nitrification
Half saturation constant for nitrification dependence on oxygen
mg O L-1
Temperature multiplier for sediment nutrient fluxes
Release rate of NH4 from sediments
g m-2 d-1
Controls sediment release of NH4 via oxygen – Half saturation constant for sediment NH4 release dependence on DO
aJorgensen and Bendoricchio 2001
b Jellison et al.. 1993
Normalised mean absolute error.
NMAE (lower bound)
NMAE (upper bound)
where s t is the simulated value at time t, o t is the observed value at time t, ō is the mean of the observed values over the simulation period and n is the number of observed values. NMAE is a measure of the absolute deviation of simulated values from observations, normalized to the mean; a value of zero indicates perfect agreement and greater than zero an average fraction of the discrepancy normalized to the mean. To compare the extent of variability within the observed data, we also calculated for each state variable, the standard deviation of observed data normalized to the mean over the simulation period (Table 3). In addition, the correlation coefficients and associated slope for the direct comparison of observed against simulated values for each state variable were calculated.
The period from 1991–1992 was used for initial parameter calibration. Comparisons of field and model data were made for six major model state variables: (1) ammonium (NH4); (2) particulate organic nitrogen (PON); (3) particulate organic carbon (POC); (4) dissolved oxygen (DO); (5) phytoplankton carbon and (6) zooplankton carbon. PON and POC refer to the sum of phytoplankton and detrital particulate nitrogen and carbon, respectively. Lake-wide averages of the surface to 9-m integrated concentrations of NH4, PON, POC, DO and phytoplankton carbon were compared. For zooplankton, lake-wide averaged biomass (g C m-2) as determined by vertical net tows were compared to vertically-integrated model output.
A manual calibration procedure was initially applied whereby individual parameters were adjusted and the model response observed. The particular features of the observed data that were used to adjust individual parameters were dependent on the parameter adjusted. For example, the minimum temperature for A. monica growth was adjusted to gain best representation of the timing of the spring zooplankton peak and the grazing rate adjusted to gain best representation of the magnitude of this peak. Individual parameters were adjusted in this way until an overall model average NMAE (calculated using the field data from the five variables listed above) of less than 0.5 was achieved.
Once a reasonable fit was achieved through trial-and-error, the local parameter space optima was determined by applying a Levenberg-Marquant (L-M) method of optimization  using a predefined Matlab® function (The MathWorks Inc., Natick, MA). In this function, parameters were adjusted to optimize the sum of the NMAE values for the same five variables listed above. Additional bounds were placed on simulated values of primary productivity and nitrogen sedimentation to fall within the ranges of those estimated by Jellison et al. (1993). The years 1993–1994 were used for model validation.
Uncertainly in model predictions arises from different sources including those associated with process representation, parameter estimation, uncertainty in inputs and observed data [39–42]. While a full analysis of model uncertainty is beyond the scope of this paper, we made an estimate of the uncertainty associated with parameter estimation by comparing output from simulations using ten different parameter sets. Initially, the five most sensitive parameters to model output were established by the sensitivity analysis described below. We then determined alternative parameter sets by fixing the lower and upper bound of each parameter and then optimizing the remaining parameters via the L-M method described above until appropriate calibration was achieved. A benchmark NMAE value of 0.5 was selected so that calibration was deemed successful if the NMAE was less than 0.5 (Table 4). The lower and upper bounds for each sensitive parameter were determined by experimental or literature ranges. Model output from this suite of parameter sets was then used as an estimate of the relative uncertainty in model output. As the upper and lower bound of the five most sensitive parameters were used, this should provide a conservative estimate of the model uncertainty associated with parameter estimation.
where Δc j is the change in output variable i from the reference value c i and Δβ j is the change in parameter j from the reference value β j . Because this study is concerned with the role of physical transport mechanisms on lake-wide nitrogen fluxes, we focused on the response of the five major nitrogen fluxes (phytoplankton uptake, sediment flux, zooplankton regeneration, settling and upward flux into surface mixed layer) to parameter manipulation. Each parameter was adjusted by ± 10% or by ± 0.01 in the case of the temperature multipliers. A time variant array of sensitivity parameters was calculated for each flux and then the average taken and used to rank the parameters according to sensitivity.
To enable us to quantify the significance of BBL transport on the ecological processes of the lake, a series of simulations were run in which this mechanism was switched off allowing a comparison between lake behavior with and without BBL transport. In DYRESM the lake-wide vertical fluxes are partitioned into internal and BBL contributions (see eq.36 in Yeates and Imberger (2004)). The BBL contribution was set to zero with all other parameters (physical and ecological) remaining the same. The simulated output of nitrogen fluxes and primary and secondary production were then analyzed and compared against base line output.
Benthic boundary layer
Computational aquatic ecosystem dynamic model
Dynamic reservoir simulation model.
The first author was funded by an Australian Postgraduate (Industry) scholarship, sponsored by the Wheatbelt Development Commission. Data were collected as part of a long-term limnological monitoring program supported by grants from the Los Angeles Department of Water and Power. Facilities at the Sierra Nevada Aquatic Research Laboratory (Natural Reserve System, University of California) were utilized for laboratory analyses and field sampling. We thank Sally MacIntyre from USCB who has provided valuable advice and direction to the research. The Contract Research Group at the Centre for Water Research (CWR) and in particular Jose Romero, Peter Yeates and Matt Hipsey provided technical support for the project.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.