Open Access
Translator Disclaimer
1 June 2014 Benthic and planktonic primary production along a nutrient gradient in Green Bay, Lake Michigan, USA
Bryan Althouse, Scott Higgins, M. Jake Vander Zanden
Author Affiliations +

Primary production in lakes occurs in both planktonic (water column) and benthic (bottom) habitats. How whole-lake primary production is distributed between these 2 habitats—referred to as autotrophic structure—is a key ecosystem property. Empirical research examining the balance between benthic and planktonic primary production in lakes is scarce, and how autotrophic structure changes across depth, nutrient, water clarity, and biological invasion gradients is unclear. Therefore, we are ill equipped to anticipate ecosystem-level responses to environmental change. We assessed the magnitude of offshore planktonic, nearshore planktonic, and nearshore benthic gross primary production (GPP) along a gradient of nutrients, water clarity, and Cladophora biomass in Green Bay, Lake Michigan, USA, during summer 2010, 2011, and 2012. Benthic and planktonic GPP varied strongly along the trophic gradient. Planktonic GPP increased with nutrients status, whereas benthic GPP decreased. From shore to 10 m depth, autotrophic structure shifted from planktonic dominance near the mouth of the Fox River (95% planktonic) to a mix of benthic and planktonic GPP 35 km from the mouth of the Fox (∼40% benthic). The steep bathymetry at more-distant sites reduced the relative importance of benthic GPP at the whole-ecosystem level. Our work highlights the dual-pathway (i.e., benthic and planktonic) nature of lentic food webs from the perspective of GPP, and shows how both trophic status and bathymetry affect autotrophic structure.

Researchers studying lakes have emphasized measurement of planktonic primary production and processes (Vadeboncoeur et al. 2002) with particular emphasis on planktonic primary production in studies of eutrophication (Carpenter et al. 1998, Correll 1998). However, recognition is growing that substantial primary production occurs in both planktonic and benthic habitats. How whole-lake primary production (planktonic + benthic) is distributed between these 2 habitats is a poorly understood, but fundamental ecosystem property referred to as the autotrophic structure (Higgins et al, in press). Little understanding exists of how primary production is partitioned between these habitats, how partitioning varies along environmental gradients, and the effects of environmental change on autotrophic structure.

Benthic primary production can be a significant contributor to littoral zone, and sometimes to whole-lake, primary production (Vadeboncoeur et al. 2001, 2008, DeNicola et al. 2003, van de Bogert et al. 2007, Sadro et al. 2011). Even in lakes where benthic habitats are a minor contributor to whole-lake primary production, top predators generally assimilate C fixed by both benthic and planktonic autotrophs (Vander Zanden and Vadeboncoeur 2002, Vander Zanden et al. 2011), often leading to a disproportionately high importance of benthic primary production as a source of energy for higher trophic levels relative to its contribution to wholelake primary production (Vander Zanden et al. 2006).

Nearshore benthic habitats in Lake Erie and Lake Ontario are primary production hotspots, with areal rates of benthic primary production equal to or exceeding planktonic primary production down to depths of 12 m (DeNicola et al. 2003, Davies and Hecky 2005, Malkin et al. 2010). These high rates of primary production are associated with mats of the filamentous green macroalga Cladophora glomerata, which are a widespread problem throughout coastal areas of the lower Laurentian Great Lakes (DeNicola et al. 2003, Davies and Hecky 2005, Malkin et al. 2010). In these ecosystems, the growth and biomass accrual of Cladophora is controlled primarily by the availability of hard substratum (e.g., rocky lake bottom) for attachment of filaments, water clarity, and the supply of bioavailable P (Higgins et al. 2008). Cladophora blooms were a common problem in the Great Lakes from the 1950s through the early 1980s, but successful implementation of strict P-abatement programs as part of the Canada-USA Great Lakes Water Quality agreement reduced total P (TP) loading and concentrations to levels that significantly reduced growth rates and biomass (Higgins et al. 2008). The resurgence of widespread Cladophora blooms in the Great Lakes has been attributed to the arrival of the nonnative zebra mussel (Dreissena polymorpha) and quagga mussel (Dreissena rostriformis bugensis) (Hecky et al. 2004, Auer et al. 2010). Dreissenids have a well described set of effects on lakes: reduced phytoplankton biomass and primary production, increased light penetration, increased bioavailability of P at the sediment-water interface, and increased biomass and primary production of benthic primary producers (Lowe and Pillsbury 1995, Hecky et al. 2004, Higgins et al. 2008, Higgins and Vander Zanden 2010).

Some investigators have measured the balance between benthic and planktonic primary production in lakes (Malkin et al. 2010, Sadro et al. 2011), but most investigators have not considered how autotrophic structure varies along key ecological gradients, such as nutrient status and water clarity. In the coming decades, nutrient-abatement efforts paired with global environmental change will alter nutrient regimes, shift temperature and precipitation patterns, and bring new species invasions (Carpenter et al. 2011). Our limited understanding of lake autotrophic structure leaves researchers and resource managers ill equipped to anticipate how lake ecosystems will respond to such changes.

Green Bay, Lake Michigan, USA, presents a unique opportunity to study how autotrophic structure varies across key gradients (e.g., nutrient levels and water clarity) within a single ecosystem. Large inputs of nutrients and sediment from the Fox River lead to nutrient-enriched and turbid waters within inner Green Bay. However, trophic status decreases and water clarity increases with distance from the Fox River. Differences in basin morphology, which can influence the available habitat for benthic autotrophs, also occur with increasing distance from the Fox River. Water depths are shallow near the mouth of the Fox River but give way to deeper water and steep slopes in middle and outer Green Bay. We assessed how autotrophic structure varies along these gradients in a large lake basin during the summer stratification period.


Study sites

We measured benthic and planktonic gross primary production (GPP) at 4 transects along Green Bays nutrient gradient. GPP transects were 1.2, 12.8, 22.4, and 34.8 km from the mouth of the Fox River (transects 1, 2, 3, and 4, respectively; Fig. 1) at locations selected to span the full trophic gradient of Green Bay. At each transect, we used chamber incubations to estimate nearshore benthic GPP (GPPB) and free-water gas dynamics to estimate nearshore planktonic GPP (GPPnp) and offshore planktonic GPP (GPPOP). At transect 1, benthic substrate consisted almost entirely of mud and sand. Transect 2 had a mixture of sand, gravel, and small rocks. Substrate at transects 3 and 4 was composed almost entirely of medium to large rocks.

Water chemistry

We measured total P (TP), soluble reactive P (SRP), total N (TN), chlorophyll a (chi a), total suspended solids (TSS), and the light attenuation coefficient (Kd) at 5 offshore locations along the trophic gradient every 2 wk during June–August 2010 and 2011 (14–25 sampling events/transect). We transported water samples to the laboratory on ice where samples were frozen (SRP, TSS, chl a) or treated with a weak solution of HCL (TP, TN). Samples were processed by the US Geological Survey (USGS)-certified water-chemistry laboratory at University of Wisconsin-Madison within 3 wk of field collection. We estimated Kd from vertical light profiles, taken at 1-m depth intervals with a LI-COR 250A underwater light meter with a flat-plate (LI-192) sensor (LI-COR, Lincoln, Nebraska).

Figure 1.

Map of Green Bay, Lake Michigan, USA. Dashed lines represent transects 1 to 4.


Cladophora biomass

We estimated Cladophora % cover at 1-m depth once a month from late May to late August (with the exception of transect 1, which lacked Cladophora) with a 0.33-m quadrat at 10 locations at each transect. We estimated Cladophora biomass from an empirical relationship between Cladophora % cover and biomass (sensu Higgins et al. 2005) developed by direct comparison of these variables measured at a subset of locations and dates. We estimated Cladophora % cover visually within each quadrat, harvested Cladophora filaments by hand, and transferred them to a mesh bag. We blotted Cladophora filaments with a paper towel and weighed them to estimate wet mass (WM; g). We dried filaments overnight at 60°C and reweighed them to estimate dry mass (DM; g). The relationship between % cover and WM was linear: biomass = 0.88 (% cover) (CV = 0.22), with no significant difference in the multiplier among categories (5, 10, and 100%) of % cover. We converted estimated WM to DM with our empirically derived equation: ln(DM) = 0.861n(WM) - 1.37 (r 2 = 0.88, p < 0.0001, F = 653.9, n = 85, df = 86), and expressed values as g/m. We averaged estimates of DM/m2 for each transect and used them in subsequent calculations (see CPPB estimates below).

Planktonic primary production

We estimated GPPOP and GPPNP from June to August 2010, 2011, and 2012 from diurnal fluctuations of O2 and estimates of air-water gas exchange. At each transect, we deployed a YSI multiparameter sonde (model 6600 V2-4; Yellow Springs Instruments, Yellow Springs, Ohio) at 1.5-m depth in the middle of the Bay (OP) and ∼ 1 km from shore (NP). We measured dissolved O2 (DO) and temperature every 15 min. We used the bookkeeping method to estimate GPP from free-water DO (Odum 1956)


where NEP = net ecosystem production, R = ecosystem respiration, and D = rate of diffusion between lake and atmosphere (see Table 1 for symbol definitions) (Odum 1956). GPP estimates can be converted from O2 to C units by multiplying by 0.375 and can be converted from C to O2 units by multiplying by 2.667.

We calculated instantaneous rates of volumetric NEP (NEP,) over 15-min intervals with the equations outlined by Staehr et al. (2010) so that


Table 1.

Definition of symbols used in the model.


where t = time interval, k is the coefficient of gas exchange between the water surface and the atmosphere, z = depth of integration, Y t = DO concentration at time t, and O2sat was derived using the equation of Weiss (1970) and corrected for barometric pressure using the methods in USGS memoranda 81.11 and 81.15 (USGS 1981a, b).




(Wanninkhof 1992) and


(Cole and Caraco 1998) where T = temperature and U10 = wind speed 10 m above the ground. We obtained wind data from the National Oceanic and Atmospheric Administration’s (NOAA) National Data Buoy Center for a location (lat 45°5′45″N, long 87°35′24″W) near Marinette, Wisconsin ( =mnmm4). We summed NEP, obtained from each 15-min time step to obtain daily estimates of NEP (g C m-1 d-1). We estimated daytime R from the average nocturnal R for the evenings before and after each daytime period. We calculated daily GPP (g C m-3 d-1) by summing daily NEP and daily R. We converted volumetric rates to areal rates (g Cm-2 d-1) by multiplying volumetric rates by the depth of the mixed layer or the water depth (if mixed layer > water depth). Where water depth was shallower than the mixed layer depth, we corrected planktonic GPP for the benthic signal by adding the depth-specific benthic net primary production (NPP; see following section). For each transect, we calculated the average GPPOP and GPPNP. As noted elsewhere (e.g., Staehr et al. 2010), our approach is based on the assumption that daytime and nocturnal R are equivalent. If daytime R exceeds nocturnal R, our estimates of NEP will be unaffected, but estimates of GPP will be elevated.


We did benthic chamber incubations 3 to 4 times at each transect. At transects 2, 3, and 4, we collected rocks with attached Cladophora from 1-m depth and placed them in 2-L clear Plexiglas® chambers. We incubated a total of 6 chambers across a range of light levels (created by screening material, and including >1 completely opaque chamber) for 45 to 90 min under natural sunlight at ambient lake temperature. We measured ambient and withinchamber photosynthetically active radiation (PAR) with a Li-Cor 250A light meter every 30 min and temperature and DO inside each chamber with a YSI ProODO probe inserted through a sampling port every 15 min. Before the first DO and temperature readings, we placed all chambers in the dark with the top gasket removed for 5 min. Water in the chambers was mixed gently via the natural movement of the surrounding water and by a mixing paddle built into the chambers.

Benthic substrata at transect 1 consisted of sand and mud and lacked Cladophora. We collected sediment cores with ∼1-L (5.5-cm diameter) clear Plexiglas chambers and measured temperature and DO in situ. We sealed the cores on both ends and incubated them with the method described above. We agitated cores gently by hand every 15 min rather than with a paddle. At the end of the incubation, we measured temperature and DO in each chamber and compared them to initial conditions measured in situ.

We estimated NPP in each chamber from the net change in DO during the incubation, corrected for chamber volume. We converted values from O2 to C units by multiplying by 0.375 (photosynthetic quotient of 1; Davies and Hecky 2005). At transects 2, 3, and 4, we normalized C fixation to Cladophora DM in each chamber, after collecting, drying (37.7°C, 24 h), and weighing collected material.

We calculated GPPB with data obtained from the incubations. We set light intensity at the onset of photosaturation (IK ) to 300 µmol photons m-2 s-1 , which is near the middle of IK values reported in a literature review of Cladophora (Higgins et al. 2008). We calculated the maximum rate of GPPB (BPmax) as the average GPPB in chambers receiving >300 pmol photons m-2 s-2. We used this approach because it yielded more robust estimates of BPmax than did dynamically fitting incubation data to saturationtangent functions used to fit photosynthetic-irradiation (PI) curves to incubation data. We modeled daily GPPB at each transect by inserting average BPmax into the equation


(Vadeboncoeur et al. 2008) where Kd is the average light attenuation coefficient at that transect, depth is 1 m, and surface light (I 0) was modeled as a sine function and estimated at 15-min intervals with the equations by Fee (1990). We summed rates of GPPB corresponding to each modeled solar-irradiance value to estimate daily GPPB (g C m-2 d-1). To scale biomass-specific GPPB up to areal GPPB at transects 2, 3, and 4, we multiplied biomass-specific GPPB by average Cladophora DM/m2 at 1 m depth for each transect.

Data manipulation

We used estimates of benthic and planktonic GPP to model areal (g C m-2 d-1) GPP-depth curves from the 0 to 10-m depths at each transect. We modeled GPPB every 0.5 m with average BPmax from each transect and the equations used to estimate GPPB (see previous section; Fee 1990, Vadeboncoeur et al. 2008). We estimated areal planktonic GPP at each 0.5-m depth by integrating the average daily volumetric GPPOP (g C m-2 d-1) at that transect to depth. At each depth, we summed benthic and planktonic GPP to estimate depth-specific total GPP (benthic + planktonic). We calculated the benthic fraction (Bfi of total GPP as Bf = GPPB/total GPP.

To evaluate how GPPB, planktonic GPP, and their relative importance (autotrophic structure) varied along the trophic gradient, we integrated the area under benthic and planktonic GPP-depth curves. This produced an overall picture of the contributions of benthic and planktonic production to total primary production at each transect, but did not incorporate among-site differences in bathymetry. We call this approach the GPP-depth integration method, which is based on the approach of Higgins et al. (2012).

We expanded upon the GPP-depth integration method by incorporating transect-specific bathymetric data. We obtained bathymetric data from a NOAA nautical map of Green Bay (OceanGrafix chart 14910; We determined depth every 100 m from shore to the middle of Green Bay along each transect. We estimated total transect GPPB, planktonic GPP (g C/d), and Bf by summing along the transect with the same approach as the GPP-depth integration method, but with the real bathymetric data. We did all modeling and statistical analyses in R (version 2.15.0; R Project for Statistical Computing, Vienna, Austria).


Nutrients, water clarity, and Cladophora

Nutrients, chl a, TSS, and Kd were highest at transect 1 and decreased sharply with distance from the Fox River (Fig. 2A–F). Mean Cladophora biomass varied significantly among sites (ANOVA, F = 24.49, df = 1, p < 0.0001). Cladophora was absent at transect 1 and averaged 52.1 to 73.6 g DM/m2 at transects 2 to 4 (Table 2).

Planktonic and benthic GPP measurements and depth trends

Volumetric GPPqp decreased from 4.7 to 0.4 g C m-3 d-1 from transect 1 to 4 (ANOVA, F = 165.84, df = 1, p < 0.0001; Table 2). Areal GPPOP decreased from 11.6 g C m-2 d-1 to 3.3 g C m-2 d-1 from transect 1 to 4 (ANOVA, F = 96.327, p < 0.0001, df = 1; Table 2). Volumetric GPPNP ranged from 0.5 and 1 g C m-3 d-1 and did not vary with location (ANOVA, F = 0.021, df = 1 ,p = 0.8857). Areal GPPNP varied among transects (ANOVA, F = 5.927, df = 1, p = 0.0153) and had the highest rates at transects 2 and 3 (Table 2). BPmax was ∼4× higher at transect 1 (soft sediment) than transect 2 (hard substrate). BPmax was highest at transects 3 and 4 (Table 2).

For each transect, we used our empirical GPP measurements to model how areal benthic and planktonic GPP varied with depth (Fig. 3A–C). GPPB was highest in shallow water and decreased with depth (Fig. 3A). GPPB at depths <1 m were almost 4× higher at transect 1 than transect 2, and declined to 0 g C m-2 d-1 by 2 m at both transects. Transects 3 and 4 had higher GPPB at shallow depths and declined to 0 g C m-2 d-1 at 6 m for transect 3, and 8 m for transect 4. In contrast, planktonic GPP increased linearly with depth, with the highest rates at transect 1 and lowest rates at transects 3 and 4 (Fig. 3B). Bf declined rapidly as a function depth at transects 1 and 2, and extended to greater depths at transects 3 and 4 (Fig. 3C).

Benthic and planktonic contributions to total transect GPP

We used 2 methods to estimate benthic and planktonic contribution to total GPP (g C/d) at the 4 transects. In the 1st method, GPP-depth integration, the areas under the benthic and planktonic GPP-depth curves (Fig. 3A, B) were compared. We plotted these estimates of benthic and planktonic GPP (g C/d) against potential trophic drivers (Fig. 4A–H). Planktonic GPP increased linearly with trophic variables, whereas GPPB values suggested a threshold response between transects 2 and 3. Bf was low at transects 1 and 2 (0.04 and 0.02, respectively) but increased substantially at transects 3 and 4 (0.3 and 0.38 respectively). Planktonic GPP was highest at transect 1 and lowest at transect 4 (Table 2, Fig. 5A). In the 2nd method, local bathymetry was incorporated to estimate benthic and planktonic contribution to total transect GPP for each transect. With this method, the benthic contribution to total transect GPP was highest at transect 1 (0.10) and lower at transects 2, 3, and 4 (0.02, 0.05, and 0.04, respectively; Table 2, Fig. 5B).

The 2 methods produced different estimates of benthic and planktonic contribution to total transect GPP (Fig. 5A, B). The GPP-depth integration method showed decreasing total GPP and increasing Bf (to a maximum near 0.4) with distance from the Fox River. However, the method based on transect-specific bathymetry showed little change in total GPP among transects, with a maximum Bf of 0.1 (transect 1) and a minimum of 0.02 (transect 2).


Studies of lake primary production, and particularly research on eutrophication, have traditionally been focused on the dynamics and consequences of primary production in the pelagic zone of lakes (Carpenter et al. 1998, Correll 1998). However, a broader view of lake ecology and eutrophication would consider both benthic and planktonic autotrophs. Such an approach is of particular relevance for the Laurentian Great Lakes and other systems impacted by nutrient enrichment and dreissenid mussels. Our study highlights the broader implications of eutrophication for autotrophic structure in the Laurentian Great Lakes. Though we consider autotrophic structure along a spatial trophic gradient, our approach and findings are applicable to other situations where changes in nutrient status occur either across space or through time.

Benthic and planktonic GPP as a function of depth

GPPB changed with depth at all transects but the rate of change depended on nutrient status. GPPB declined rapidly with depth at the 2 most eutrophic transects (1 and 2), and continued to much greater depths (6–8 m) at the least eutrophic transects (3 and 4). We attribute the large variation in maximum depth of GPPB to differences in light attenuation between transects in inner Green Bay (1 and 2) and those in middle/outer Green Bay (3 and 4).

GPPB at depths <2 m at transects with Cladophora were among the highest reported for Cladophora and were similar to primary production rates in highly productive tropical coral reefs (Carpenter 1985, Gattuso et al. 1998). Below 2 m at transects 3 and 4, our estimates of GPPB were comparable to previous estimates of GPPB in Cladophora mats in the Great Lakes. Malkin et al. (2010) reported rates of GPPB between 0.8 and 1.6 g C m-2 d-1 in Lake Ontario, and Davies and Hecky (2005) reported GPPB rates of 0.9 g C m-2 d-1 in Lake Erie Cladophora mats. Average Cladophora biomass (∼ 60–80 g DM/m2) at transects 3 and 4 were similar to, but slightly lower than reports for Lakes Erie and Ontario, whereas our maximum biomass estimates were similar to previous reports (Higgins et al. 2005, 2008). Our estimates of BPmax at transects with Cladophora (transects 2, 3, and 4 had mean BPmax of 2.1 mg C m-2 d-1) agree well with reported values for other locations in the Great Lakes (BPmax ranging from 1-4 mg C m-2 d-1 reported by Malkin et al. (2010).

Figure 2.

Mean total P (TP) (A), soluble reactive P (SRP) (B), total N (TN) (C), chlorophyll a (chi a) (D), total suspended solids (TSS) (E), and light attenuation coefficient (Kd) (F) at 5 locations along the trophic gradient in Green Bay, Lake Michigan, USA. All data were collected during 2010 and 2011 at offshore planktonic (OP) locations. Water-quality sampling sites do not correspond with gross primary production (GPP) transects except at distances 1.2 km (transect 1) and 34.8 km (transect 4).


Autotrophic structure and environmental variables

Trophic status played a role in structuring the relative contribution of benthic and planktonic primary production along the trophic gradient in Green Bay. As nutrient concentration decreased, total planktonic primary production decreased whereas total GPPB increased (GPP-depth integration method). Planktonic GPP was more responsive to nutrient status than GPPB. This result is not surprising given that phytoplankton are generally nutrient limited (Hecky et al. 1993). In contrast, benthic algae tend to be light rather than nutrient limited (Vadeboncoeur et al. 2001), and this limitation explains the more muted response of benthic algae to water-column nutrient status.

Nutrients were correlated with changes in benthic and planktonic GPP, but bathymetry strongly influenced patterns of GPP. We anticipated that oligotrophic sites (i.e., transect 4) would have the largest Bf to total GPP because of increased water clarity, but we found the opposite to be true. The most eutrophic location (transect 1), had the greatest benthic contribution (Bf = 0.1). Thus, despite high areal rates of GPPB at shallow depths at transect 3 and 4, the steep bathymetry restricted GPPB to a narrow strip of the littoral zone. Meanwhile, the high water clarity at these sites (Fig. 2F) allowed planktonic primary production to extend to greater depths. Our results indicate that trophic status may be an important factor in determining the relative importance of benthic and planktonic primary production at fine spatial scales, but bathymetry can be a driving factor at the ecosystem level in large lakes. Considering that much of Lake Michigan has even higher ratios of planktonic volume to littoral surface area than transects 3 and 4, benthic habitats are likely to be important contributors only in nearshore areas, whereas planktonic primary production is the dominant contributor of fixed energy and nutrients at the ecosystem scale.

Table 2.

Planktonic and benthic primary production characteristics along the trophic gradient. Planktonic (plank) gross primary production (GPP) rates are mean (±1 SD), whereas benthic (benth) GPP is reported as the mean production rate for each individual transect. Cladophora biomass at each transect is shown as mean (±1 SD). Maximum benthic production (BPmax) is the average BPmax rate (mg C m-2 h-1) at each transect. All estimates of planktonic and benthic production using M1 and M2 are in units of g C/d. B f = the benthic fraction of total GPP using either M1 or M2.


Figure 3.

Benthic gross primary production (GPPB) (A), planktonic GPP (B), and the benthic fraction of total GPP (Bf ) (C) as a function of depth at transects 1 to 4. GPPB was derived from chamber estimates, whereas planktonic GPP was derived from free water O2 dynamics and adjusted for the benthic signal.


In addition to nutrient status and bathymetry, biological invasion also has the potential to shape autotrophic structure. Establishment of dreissenid mussels in the Great Lakes has sharply reduced phytoplankton biomass and light attenuation (Brezonik et al. 2005, Binding et al. 2007) and redirected nutrients to nearshore benthic habitats (Hecky et al. 2004) to yield increased GPPB (Higgins et al. 2008, Auer et al. 2010). These changes have led to increased reliance on benthic resources in Great Lake food webs (Sierszen et al. 2006, Rennie et al. 2009) despite the dominance of planktonic primary production at the whole-lake level. GPPB may not be a significant contributor to total ecosystem primary production, but on an areal basis, rates of GPPB are equivalent to those occurring in coral reefs (Carpenter 1985, Gattuso et al. 1998). This narrow band of extremely high productivity may explain why these shallow littoral areas are also hotspots of biodiversity in large lakes across the globe (Vadeboncoeur et al. 2011). Furthermore, Great Lake food webs may rely disproportionately on benthicderived C and could be sensitive to changes in rates of GPPB in nearshore and littoral habitats. For these reasons, changes in GPPB consequent to eutrophication, dreissenid mussels, and other drivers could have a large influence on lake food webs and biodiversity without substantially affecting total lake primary production.

Figure 4.

Total daily transect planktonic (A, C, E, G) and benthic (B, D, F, H) gross primary production (GPP) as a function of soluble reactive P (SRP) (A, B), total P (TP) (C, D), total N (TN) (E, F), and the light attenuation coefficient (Kd) (G, H). Both benthic and planktonic total transect GPP were calculated as the area under the GPP vs depth curves from the depth-integration method (Fig. 3A–C). Data points are indicated by transect numbers (1–4).


Figure 5.

Total gross primary production (GPP) and benthic contribution to total GPP calculated by the depth-integration method (A) and using transect-specific bathymetry (B) as a function of distance from the Fox River. Points represent total daily GPP in each transect, vertical grey lines indicate planktonic contribution to total GPP, and black diagonal lines indicate benthic contribution to total GPP.


Benthic and planktonic primary production in Green Bay: past trends and future scenarios

The cumulative effects of eutrophication and biological invasions have undoubtedly influenced rates of both planktonic and benthic primary production over the past 25 y. TP entering Green Bay via the Fox River increased from a mean of ∼100 pg/L in 1986 to a peak of nearly 210 µg/L between 2001–2004 before dropping to current levels near 150 µg/L (Green Bay Metropolitan Sewage District; This increase in P load corresponded to increases in areal planktonic primary production from 2.5 g C m-2 d-1 in 1986 (Auer and Canale 1986) to 3.7 g C m-2 d-1 in 1990 (Millard and Sager 1994). Volumetric rates of planktonic production also increased during that time period from 3.0 g C m-3 d-1 in 1986 to near 6 g C m-3 d-1 in 1990.

Current TP concentrations in inner Green Bay are similar to levels in 1990, and our estimates of volumetric planktonic production are comparable to those reported in the early 1990s (Millard and Sager 1994). At more distant locations from the Fox River (e.g., our transect 4) TP declined from ∼80 µg/L in 1986 to ∼5 µg/L in 2010–2011, presumably as a result of filtration by dreissenid mussels. Though TP in middle Green Bay has declined, planktonic primary production appears to have increased since 1990. Part of this discrepancy is a result of differences in methods used to calculate planktonic GPP. However, although dreissenids reduced phytoplankton biomass in the Great Lakes through grazing, they also may have increased phytoplankton-specific growth rates (Heath et al. 1995). Dreissenid-induced increases in water clarity in the Great Lakes (Budd et al. 2001, Binding et al. 2007) may have offset reductions in phytoplankton biomass, resulting in increased biomass-specific rates of planktonic primary production and higher rates of volumetric and areal GPP (Millard and Sager 1994).

At present, nutrient abatement offers the best opportunity to control nuisance growth of both benthic and planktonic algae in Green Bay. However, reductions in suspended sediment and nutrient loading from the Fox River could increase phytoplankton primary production near transect 1 because of increased light. In contrast, P reduction may do little to reduce nuisance Cladophora at transects 3 and 4, where dreissenid density is high and watercolumn TP and light attenuation are currently at levels indicative of an oligotrophic state.


On an areal basis, the benthic contribution to total primary Bf in Green Bay was high in shallow water and declined with depth, but how the benthic contribution changed with depth varied strongly along the trophic gradient. GPPB made relatively minor contributions to total primary production at the transect or whole-lake level in Green Bay (and presumably other Great Lakes) during summer stratification because of steep bathymetry and large ecosystem size. However, on an areal basis, the narrow band of high GPPB is probably important for maintenance of biodiversity hotspots in the littoral zones of large lakes (Vadeboncoeur et al. 2011). Furthermore, GPPB contributes a disproportionately large amount of fixed energy and nutrients to higher trophic levels in food webs of large lakes invaded by dreissenid mussels (Rennie et al. 2013). Thus, food webs in large lakes may be sensitive to changes in GPPB. The relative contributions of benthic and planktonic habitats to total primary production (autotrophic structure) at an individual location or the ecosystem level is a fundamental ecosystem attribute that is poorly explored but demands future research.


This work was funded by the University of Wisconsin Sea Grant Institute under grants from the National Sea Grant College Program, NO A A, US Department of Commerce, and the State of Wisconsin (Federal grant number NA100AR4170070, project number R/HCE-7). Special thanks to the following for their contributions to this research: Bart DeStasio, Paul Hanson, Tony Ives, Jim Hurley, Jake Walsh, Kaity Taylor, Aaron Heimann, Zach Ruff, and Margaret Porco.



M. T. Auer , and R. P. Canale . 1986. Mathematical modelling of primary production in Green Bay (Lake Michigan, USA), a phosphorus- and light-limited system. Hydrobiological Bulletin 20:195–211. Google Scholar


M. T. Auer , L. M. Tomlinson , S. N. Higgins , S. Y. Malkin , E. T. Howell , and H. A. Bootsma . 2010. Great Lakes Cladophora in the 21st century: same algae—different ecosystem. Journal of Great Lakes Research 36:248–255. Google Scholar


C. E. Binding , J. H. Jerome , R. P. Bukata , and W. G. Booty . 2007. Trends in water clarity of the lower Great Lakes from remotely sensed aquatic color. Journal of Great Lakes Research 33:828–841. Google Scholar


P. Brezonik , K. D. Menken , and M. Bauer . 2005. Landsat-based remote sensing of lake water quality characteristics, including chlorophyll and colored dissolved organic matter (CDOM). Lake and Reservoir Management 21:373–382. Google Scholar


J. W. Budd , T. D. Drummer , T. F. Nalepa , and G. L. Fahnenstiel . 2001. Remote sensing of biotic effects: zebra mussels (Dreissena polymorpha) influence on water clarity in Saginaw Bay, Lake Huron. Limnology and Oceanography 46:213–223. Google Scholar


R. C. Carpenter 1985. Relationships between primary production and irradiance in coral-reef algal communities. Limnology and Oceanography 30:784–793. Google Scholar


S. R. Carpenter , N. F. Caraco , D. L. Correll , R. W. Howarth , A. N. Sharpley , and V. H. Smith . 1998. Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecological Applications 8:559–568. Google Scholar


S. R. Carpenter , E. H. Stanley , and M. J. Vander Zanden . 2011. State of the worlds freshwater ecosystems: physical, chemical, and biological changes. Pages 75–99 in A. Gadgil and D. M. Liverman (editors). Volume 36. Annual review of environment and resources. Annual Reviews, Palo Alto, California. Google Scholar


J. Cole , and N. F. Caraco . 1998. Atmospheric exchange of carbon dioxide in a low-wind oligotrophic lake measured by the addition of SF6. Limnology and Oceanography 43:647–656. Google Scholar


D. L. Correll 1998. The role of phosphorus in the eutrophication of receiving waters: a review. Journal of Environmental Quality 27:261–266. Google Scholar


J. M. Davies , and R. E. Hecky . 2005. Initial measurements of benthic photosynthesis and respiration in Lake Erie. Journal of Great Lakes Research 31:195–207. Google Scholar


D. M. DeNicola , E. De Eyto , A. Wemaere , and K. Irvine . 2003. Production and respiration of epilithic algal communities in Irish lakes of different trophic status. Archiv für Hydrobiologie 157:67–87. Google Scholar


E. J. Fee 1990. Computer programs for calculating in-situ phytoplankton photosynthesis. Canadian Technical Report of Fisheries and Aquatic Sciences No. 1740. Department of Fisheries and Oceans, Central and Arctic Region, Winnipeg, Manitoba. Google Scholar


J. P. Gattuso , M. Frankignoulle , and R. Wollast . 1998. Carbon and carbonate metabolism in coastal aquatic ecosystems. Annual Review of Ecology and Systematics 29:405–434. Google Scholar


R. T. Heath , G. L. Fahnenstiel , W. S. Gardner , J. F. Cavaletto , and S. J. Hwang . 1995. Ecosystem-level effects of zebra mussels (Dreissena polymorpha): an enclosure experiment in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 21:501–516. Google Scholar


R. E. Hecky , P. Campbell , and L. L. Hendzel . 1993. The stoichiometry of carbon, nitrogen, and phosphorus in particulate matter in lakes and oceans. Limnology and Oceanography 38:709–724. Google Scholar


R. E. Hecky , R. E. H. Smith , D. R. Barton , S. J. Guildford , W. D. Taylor , M. N. Charlton , and T. Howell . 2004. The nearshore phosphorus shunt: a consequence of ecosystem engineering by dreissenids in the Laurentian Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 61:1285–1293. Google Scholar


S. N. Higgins , B. Althouse , S. P. Devlin , Y. Vadeboncoeur , and M. J. Vander Zanden . Potential for large-bodied zooplankton and dreissenids to alter the productivity and autotrophic structure of lakes. Ecology (in press). Google Scholar


S. N. Higgins , R. E. Hecky , and S. J. Guildford . 2005. Modeling the growth, biomass, and tissue phosphorus concentration of Cladophora glomerata in eastern Lake Erie: model description and field testing. Journal of Great Lakes Research 31:439–455. Google Scholar


S. N. Higgins , S. Y. Malkin , E. T. Howell , S. J. Guildford , L. Campbell , V. Hiriart-Baer , and R. E. Hecky . 2008. An ecological review of Cladophora glomerata (Chlorophyta) in the Laurentian Great Lakes. Journal of Phycology 44:839–854. Google Scholar


S. N. Higgins , C. M. Pennuto , E. T. Howell , T. W. Lewis , and J. C. Makarewicz . 2012. Urban influences on Cladophora blooms in Lake Ontario. Journal of Great Lakes Research 38: 116–123. Google Scholar


S. N. Higgins , and M. J. Vander Zanden . 2010. What a difference a species makes: a meta-analysis of dreissenid mussel impacts on freshwater ecosystems. Ecological Monographs 80(2):179–196. Google Scholar


R. L. Lowe , and R. W. Pillsbury . 1995. Shifts in benthic algal community structure and function following the appearance of zebra mussels (Dreissena polymorpha) in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 21:558–566. Google Scholar


S. Y. Malkin , S. A. Bocaniov , R. E. Smith , S. J. Guildford , and R. E. Hecky . 2010. In situ measurements confirm the seasonal dominance of benthic algae over phytoplankton in nearshore primary production of a large lake. Freshwater Biology 55:2468–2483. Google Scholar


E. S. Millard , and P. E. Sager . 1994. Comparison of phosphorus, light climate, and photosynthesis between two culturally eutrophied bays—Green Bay, Lake Michigan, and the Bay of Quinte, Lake Ontario. Canadian Journal of Fisheries and Aquatic Sciences 51:2579–2590. Google Scholar


H. T. Odum 1956. Primary production in flowing waters. Limnology and Oceanography 1:102–117. Google Scholar


M. D. Rennie , D. O. Evans , and J. D. Young . 2013. Increased dependence on nearshore benthic resources in the Lake Simcoe ecosystem after dreissenid invasion. Inland Waters 3: 297–310. Google Scholar


M. D. Rennie , W. G. Sprules , and T. B. Johnson . 2009. Resource switching in fish following a major food web disruption. Oecologia (Berlin) 159:789–802. Google Scholar


S. Sadro , J. M. Melack , and S. MacIntyre . 2011. Spatial and temporal variability in the ecosystem metabolism of a high-elevation lake: integrating benthic and pelagic habitats. Ecosystems 14:1123–1140. Google Scholar


M. E. Sierszen , G. S. Peterson , and J. V. Scharold . 2006. Depthspecific patterns in benthic-planktonic food web relationships in Lake Superior. Canadian Journal of Fisheries and Aquatic Sciences 63:1496–1503. Google Scholar


P. A. Staehr , D. Bade , M. C. van de Bogert , G. R. Koch , C. Williamson , P. Hanson , J. J. Cole , and T. Kratz . 2010. Lake metabolism and the diel oxygen technique: state of the science. Limnology and Oceanography: Methods 8:628–644. Google Scholar


USGS (US Geological Survey). 1981a. New tables of dissolved oxygen saturation values. Water Quality Technical Memorandum 81.11. US Geological Survey, Reston, Virginia. Google Scholar


USGS (US Geological Survey). 1981b. New tables of dissolved oxygen saturation values; amendment of quality of water technical memorandum no. 81.11. Water Quality Technical Memorandum 81.15. US Geological Survey, Reston, Virginia. Google Scholar


Y. Vadeboncoeur , D. M. Lodge , and S. R. Carpenter . 2001. Whole-lake fertilization effects on distribution of primary production between benthic and pelagic habitats. Ecology 82: 1065–1077. Google Scholar


Y. Vadeboncoeur , P. B. McIntyre , and M. J. Vander Zanden . 2011. Borders of biodiversity: life at the edge of the world's large lakes. BioScience 61:526–537. Google Scholar


Y. Vadeboncoeur , G. Peterson , M. J. Vander Zanden , and J. Kalff . 2008. Benthic algal production across lake size gradients: interactions among morphometry, nutrients, and light. Ecology 89:2542–2552. Google Scholar


Y. Vadeboncoeur , M. J. Vander Zanden , and D. M. Lodge . 2002. Putting the lake back together: reintegrating benthic pathways into lake food web models. BioScience 52:44–54. Google Scholar


M. C. van de Bogert , S. R. Carpenter , J. J. Cole , and M. L. Pace . 2007. Assessing pelagic and benthic metabolism using free water measurements. Limnology and Oceanography: Methods 5:145–155. Google Scholar


M. J. Vander Zanden , S. Chandra , S. K. Park , Y. Vadeboncoeur , and C. R. Goldman . 2006. Efficiencies of benthic and pelagic trophic pathways in a subalpine lake. Canadian Journal of Fisheries and Aquatic Sciences 63:2608–2620. Google Scholar


M. J. Vander Zanden , and Y. Vadeboncoeur . 2002. Fishes as integrators of benthic and pelagic food webs in lakes. Ecology 83:2152–2161. Google Scholar


M. J. Vander Zanden , Y. Vadeboncoeur , and S. Chandra . 2011. Fish reliance on littoral-benthic resources and the distribution of primary production in lakes. Ecosystems 14:894–903. Google Scholar


R. Wanninkhof 1992. Relationship between wind-speed and gas exchange over the ocean. Journal of Geophysical Research: Oceans 97:7373–7382. Google Scholar


R. F. Weiss 1970. Solubility of nitrogen, oxygen and argon in water and seawater. Deep Sea Research 17:721–735. Google Scholar
© 2014 by The Society for Freshwater Science.
Bryan Althouse, Scott Higgins, and M. Jake Vander Zanden "Benthic and planktonic primary production along a nutrient gradient in Green Bay, Lake Michigan, USA," Freshwater Science 33(2), 487-498, (1 June 2014).
Received: 13 August 2013; Accepted: 6 November 2013; Published: 1 June 2014

autotrophic structure
benthic primary production
Green Bay
Lake Michigan
Get copyright permission
Back to Top