Ecology of larval fishes and large zooplankton in the Keweenaw Current region of Lake Superior, with special focus on lake herring, Coregonus artedi

I assessed the influence of the Keweenaw Current and spring thermal bar on the distribution of larval fishes and large zooplankton in Lake Superior. In 1998 and 1999, samples were collected from inshore (0.2 – 3.0 km from shore) and offshore (5.0 – 9.0 km from shore) locations on three transects off the western coast of the Keweenaw Peninsula, Michigan. For larval fishes, density and size distribution patterns of lake herring (Coregonus artedi), rainbow smelt (Osmerus mordax), burbot (Lota lota), deepwater sculpin (Myoxocephalus thompsoni), and spoonhead sculpin (Cottus ricei) suggest a seasonal inshore to offshore movement. For zooplankton, seasonal warming appeared to be the major factor that limited planktonic catches of the primarily benthic Mysisrelicta and Diporeia spp., while simultaneously stimulated growth and reproduction of the cladocerans Daphnia spp., Holopedium gibberum, and Bythotrephes cederstroemi. In contrast, calanoid copepods as a group were abundant throughout the entire sampling season. The greatest abundances of zooplankton were generally encountered offshore, even for the cladocerans, which apparently expanded from inshore to offshore locations with seasonal warming. In 2000, sampling efforts focused on lake herring. Samples were collected from surface waters at 0.1 – 17.0 km from shore on two transects. Lake herring larvae were also reared in the laboratory from eggs in order to validate the use of otolith microstructure for aging. Increment deposition was not statistically different from a daily rate starting from 28 days after hatching, near the time of yolk-sac absorption, but larvae with lower growth rates could not be aged as accurately. In Lake Superior, lake herring tended to be slightly more abundant, larger, and older at inshore locations, but a dense patch of younger larvae was also encountered 7 – 13 km from shore. The distribution iiipatterns suggest that larvae were transported by prevailing currents into the study region, possibly from the more productive spawning regions in western Lake Superior. Growth rates were suppressed at offshore locations where temperatures were less than 8°C. These results indicate that lake herring larvae may be transported far distances from spawning concentrations by longshore currents, and water temperatures may largely control their growth.

Open top chambers and infrared lamps : a comparison of heating efficacy and CO₂/CH₄ dynamics in a Lake Superior coastal peatland

Experimental warming provides a method to determine how an ecosystem will respond to increased temperatures. Northern peatland ecosystems, sensitive to changing climates, provide an excellent setting for experimental warming. Storing great quantities of carbon, northern peatlands play a critical role in regulating global temperatures. Two of the most common methods of experimental warming include open top chambers (OTCs) and infrared (IR) lamps. These warming systems have been used in many ecosystems throughout the world, yet their efficacy to create a warmer environment is variable and has not been widely studied. To date, there has not been a direct, experimentally controlled comparison of OTCs and IR lamps. As a result, a factorial study was implemented to compare the warming efficacy of OTCs and IR lamps and to examine the resulting carbon dioxide (CO2) and methane (CH4) flux rates in a Lake Superior peatland. IR lamps warmed the ecosystem on average by 1-2 #°C, with the majority of warming occurring during nighttime hours. OTC's did not provide any long-term warming above control plots, which is contrary to similar OTC studies at high latitudes. By investigating diurnal heating patterns and micrometeorological variables, we were able to conclude that OTCs were not achieving strong daytime heating peaks and were often cooler than control plots during nighttime hours. Temperate day-length, cloudy and humid conditions, and latent heat loss were factors that inhibited OTC warming. There were no changes in CO2 flux between warming treatments in lawn plots. Gross ecosystem production was significantly greater in IR lamp-hummock plots, while ecosystem respiration was not affected. CH4 flux was not significantly affected by warming treatment. Minimal daytime heating differences, high ambient temperatures, decay resistant substrate, as well as other factors suppressed significant gas flux responses from warming treatments.

Modeling contaminant behavior in Lake Superior : a comparison of PCBs, PBDEs, and mercury

A mass‐balance model for Lake Superior was applied to polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and mercury to determine the major routes of entry and the major mechanisms of loss from this ecosystem as well as the time required for each contaminant class to approach steady state. A two‐box model (water column, surface sediments) incorporating seasonally adjusted environmental parameters was used. Both numerical (forward Euler) and analytical solutions were employed and compared. For validation, the model was compared with current and historical concentrations and fluxes in the lake and sediments. Results for PCBs were similar to prior work showing that air‐water exchange is the most rapid input and loss process. The model indicates that mercury behaves similarly to a moderately‐chlorinated PCB, with air‐water exchange being a relatively rapid input and loss process. Modeled accumulation fluxes of PBDEs in sediments agreed with measured values reported in the literature. Wet deposition rates were about three times greater than dry particulate deposition rates for PBDEs. Gas deposition was an important process for tri‐ and tetra‐BDEs (BDEs 28 and 47), but not for higher‐brominated BDEs. Sediment burial was the dominant loss mechanism for most of the PBDE congeners while volatilization was still significant for tri‐ and tetra‐BDEs. Because volatilization is a relatively rapid loss process for both mercury and the most abundant PCBs (tri‐ through penta‐), the model predicts that similar times (from 2 ‐ 10 yr) are required for the compounds to approach steady state in the lake. The model predicts that if inputs of Hg(II) to the lake decrease in the future then concentrations of mercury in the lake will decrease at a rate similar to the historical decline in PCB concentrations following the ban on production and most uses in the U.S. In contrast, PBDEs are likely to respond more slowly if atmospheric concentrations are reduced in the future because loss by volatilization is a much slower process for PBDEs, leading to lesser overall loss rates for PBDEs in comparison to PCBs and mercury. Uncertainties in the chemical degradation rates and partitioning constants of PBDEs are the largest source of uncertainty in the modeled times to steady‐state for this class of chemicals. The modeled organic PBT loading rates are sensitive to uncertainties in scavenging efficiencies by rain and snow, dry deposition velocity, watershed runoff concentrations, and uncertainties in air‐water exchange such as the effect of atmospheric stability.

Tracing the source of groundwater for three different coastal peatlands along Lake Superior

The goal of this project was to investigate the influence of a large inland lake on adjacent coastal freshwater peatlands. The specific aim was to determine the source of groundwater for three differently formed peatlands located on the southern shore of Lake Superior. The groundwater study was conducted at Bete Grise, a peatland complex in a dune-swale system; Pequaming, a peatland developed in the swale of a tombolo; and Lightfoot Bay, a peatland developed in a barrier beach wetland complex. To determine the source of groundwater in the peatlands, transects of six groundwater monitoring wells were established at each study site, covering distinctly different vegetation zones. At Pequaming and Lightfoot Bay the transects monitored two vegetation zones: transition zone from upland and open fen. At Bete Grise, the transects monitored dunes and swales. Additionally, at all three sites, upland groundwater was monitored using three wells that were installed into the adjacent upland forest. Biweekly measurements of well water pH and specific conductance were carried out from May to October of 2010. At each site, vegetation cover, peat depths and surface elevations were determined and compared to Lake Superior water levels. From June 14 – 17, July 20 – 21 and September 10 – 12, stable isotopes of oxygen (18O/16O) ratios were measured in all the wells and for Lake Superior water. A mixing model was used to estimate the percentage of lake water influencing each site based on the oxygen isotope ratios. During the sampling period, groundwater at all three sites was supported primarily by upland groundwater. Pequaming was approximately 80 % upland groundwater supported and up to 20 % Lake water supported in the uppermost 1 m layer of peat column of the transition zone and open fen. Bete Grise and Lightfoot Bay were 100 % upland groundwater supported throughout the season. The height of Lake Superior was near typical levels in 2010. In years when the lake level is higher, Lake water could intrude into the adjacent peatlands. However, under typical hydrologic conditions, these coastal peatlands are primarily supported by upland groundwater.

Mechanisms contributing to the deep chlorophyll maximum in Lake Superior

The seasonal appearance of a deep chlorophyll maximum (DCM) in Lake Superior is a striking phenomenon that is widely observed; however its mechanisms of formation and maintenance are not well understood. As this phenomenon may be the reflection of an ecological driver, or a driver itself, a lack of understanding its driving forces limits the ability to accurately predict and manage changes in this ecosystem. Key mechanisms generally associated with DCM dynamics (i.e. ecological, physiological and physical phenomena) are examined individually and in concert to establish their role. First the prevailing paradigm, “the DCM is a great place to live”, is analyzed through an integration of the results of laboratory experiments and field measurements. The analysis indicates that growth at this depth is severely restricted and thus not able to explain the full magnitude of this phenomenon. Additional contributing mechanisms like photoadaptation, settling and grazing are reviewed with a one-dimensional mathematical model of chlorophyll and particulate organic carbon. Settling has the strongest impact on the formation and maintenance of the DCM, transporting biomass to the metalimnion and resulting in the accumulation of algae, i.e. a peak in the particulate organic carbon profile. Subsequently, shade adaptation becomes manifest as a chlorophyll maximum deeper in the water column where light conditions particularly favor the process. Shade adaptation mediates the magnitude, shape and vertical position of the chlorophyll peak. Growth at DCM depth shows only a marginal contribution, while grazing has an adverse effect on the extent of the DCM. The observed separation of the carbon biomass and chlorophyll maximum should caution scientists to equate the DCM with a large nutrient pool that is available to higher trophic levels. The ecological significance of the DCM should not be separated from the underlying carbon dynamics. When evaluated in its entirety, the DCM becomes the projected image of a structure that remains elusive to measure but represents the foundation of all higher trophic levels. These results also offer guidance in examine ecosystem perturbations such as climate change. For example, warming would be expected to prolong the period of thermal stratification, extending the late summer period of suboptimal (phosphorus-limited) growth and attendant transport of phytoplankton to the metalimnion. This reduction in epilimnetic algal production would decrease the supply of algae to the metalimnion, possibly reducing the supply of prey to the grazer community. This work demonstrates the value of modeling to challenge and advance our understanding of ecosystem dynamics, steps vital to reliable testing of management alternatives.

Study of particle movement in the nearshore region of Lake Superior with radioisotope tracers

Biogeochemical processes in the coastal region, including the coastal area of the Great Lakes, are of great importance due to the complex physical, chemical and biological characteristics that differ from those on either the adjoining land or open water systems. Particle-reactive radioisotopes, both naturally occurring (210Pb, 210Po and 7Be) and man-made (137Cs), have proven to be useful tracers for these processes in many systems. However, a systematic isotope study on the northwest coast of the Keweenaw Peninsula in Lake Superior has not yet been performed. In this dissertation research, field sampling, laboratory measurements and numerical modeling were conducted to understand the biogeochemistry of the radioisotope tracers and some particulate-related coastal processes. In the first part of the dissertation, radioisotope activities of 210Po and 210Pb in a variability of samples (dissolved, suspended particle, sediment trap materials, surficial sediment) were measured. A completed picture of the distribution and disequilibrium of this pair of isotopes was drawn. The application of a simple box model utilizing these field observations reveals short isotope residence times in the water column and a significant contribution of sediment resuspension (for both particles and isotopes). The results imply a highly dynamic coastal region. In the second part of this dissertation, this conclusion is examined further. Based on intensive sediment coring, the spatial distribution of isotope inventories (mainly 210Pb, 137Cs and 7Be) in the nearshore region was determined. Isotope-based focusing factors categorized most of the sampling sites as non- or temporary depositional zones. A twodimensional steady-state box-in-series model was developed and applied to individual transects with the 210Pb inventories as model input. The modeling framework included both water column and upper sediments down to the depth of unsupported 210Pb penetration. The model was used to predict isotope residence times and cross-margin fluxes of sediments and isotopes at different locations along each transect. The time scale for sediment focusing from the nearshore to offshore regions of the transect was on the order of 10 years. The possibility of sediment longshore movement was indicated by high inventory ratios of 137Cs: 210Pb. Local deposition of fine particles, including fresh organic carbon, may explain the observed distribution of benthic organisms such as Diporeia. In the last part of this dissertation, isotope tracers, 210Pb and 210Po, were coupled into a hydrodynamic model for Lake Superior. The model was modified from an existing 2-D finite difference physical-biological model which has previously been successfully applied on Lake Superior. Using the field results from part one of this dissertation as initial conditions, the model was used to predict the isotope distribution in the water column; reasonable results were achieved. The modeling experiments demonstrated the potential for using a hydrodynamic model to study radioisotope biogeochemistry in the lake, although further refinements are necessary.

FACTORS AFFECTING DIPOREIA GROWTH RATES IN LAKE SUPERIOR

An ability to predict population dynamics of the amphipod Diporeia is important in understanding how energy pathways in the Lake Superior food web might be altered by disturbances to the ecosystem. Estimating growth rates for this prominent prey item for fish requires information on the physiological effects of changes to its environment. These effects have been investigated for Diporeia in other Great Lakes, but little is known about Lake Superior populations. The primary objective of this study is to obtain quantitative data for rates of Diporeia respiration and consumption that can be incorporated into a bioenergetics model for Lake Superior. Benthic communities in Lake Superior were sampled bimonthly from April through September during 2011 and 2012 to investigate spatial and temporal trends of Diporeia abundances as well as size class structures of the population. Additional samples of Diporeia were collected and kept alive in natural sediment for laboratory experiments. Respiration rates for Diporeia were measured by monitoring dissolved oxygen concentrations in microcosoms using microelectrodes. Additionally, a series of experiments to estimate consumption rates based on food availability were conducted using 14C-labeled algae (Selenastrum capricornutum). Amphipod population densities are highest between 30-110 m (slope) compared to 0-30 m (shelf) or >110 m (profundal) regions in Lake Superior. This heterogeneous distribution of Diporeia in Lake Superior is an important component to quantifying lake-wide biomass. Rates of oxygen consumption by Diporeia range from 32.0 to 44.7 mgO2*gDW-1*d-1, and do not vary significantly with body size per individual. The predicted consumption rate corresponding to average Lake Superior algal carbon fluxes was 0.08 ± SE mgC*gDW-1*d-1. Data on Lake Superior Diporeia biomass and bioenergetics found in this study can be incorporated in a model used to estimate the viability of this population under potential future environmental stressors.