No description available.

No description available.

No description available.

The Great Lakes has been impaired by elevated concentrations of contaminants, and different regulatory actions are taken to decrease theses concentrations. Recent research of Dr. Azim suggests a decline in mercury levels in the Lake Erie, and an increase in the Hg concentration in the fish communities. This research is a follow up study to examine whether the contaminant levels are indeed increasing. In the first chapter we evaluate the temporal trends of polychlorinated biphenyls (PCBs) in Lake Erie fish over the last 30 years by conducting the Exponential and Dynamic Linear Modeling analysis. DLM analysis in the chapter two detects the PCB and Hg contaminant trends in a broader range. There is a high variability in two contaminants trends. These patterns may be because of factors such as fish physiological parameters, and the impact of invasive species. The results of my research will contribute to the Toxics Reduction Strategy.

The Great Lakes has been impaired by elevated concentrations of contaminants, and different regulatory actions are taken to decrease theses concentrations. Recent research of Dr. Azim suggests a decline in mercury levels in the Lake Erie, and an increase in the Hg concentration in the fish communities. This research is a follow up study to examine whether the contaminant levels are indeed increasing. In the first chapter we evaluate the temporal trends of polychlorinated biphenyls (PCBs) in Lake Erie fish over the last 30 years by conducting the Exponential and Dynamic Linear Modeling analysis. DLM analysis in the chapter two detects the PCB and Hg contaminant trends in a broader range. There is a high variability in two contaminants trends. These patterns may be because of factors such as fish physiological parameters, and the impact of invasive species. The results of my research will contribute to the Toxics Reduction Strategy.

Variable fluorescence of chlorophyll a (Fv/Fm), measured by pulse amplitude modulated (PAM) fluorometers, is an attractive target for phytoplankton-related water quality management. Lowered Fv/Fm is believed to reflect the magnitude of nutrient sufficiency or deficiency in phytoplankton. This rapid and specific metric is relevant to Lake Erie, which often experiences problematic Cyanobacteria blooms. It is unknown whether PAMs reliably measure phytoplankton nutrient status or if different PAMs provide comparable results. Water samples collected from Lake Erie and two Lake Ontario sites in July and September 2011 were analysed using alkaline phosphatase assay (APA), P-debt, and N-debt to quantify phytoplankton nutrient status and with three different PAM models (PhytoPAM, WaterPAM and DivingPAM) to determine Fv/Fm. The Lake Ontario, Lake Erie East and Central Basin sites were all N- and P-deficient in July, but only the East and Central Basin and one Lake Ontario site were P-deficient in September. The West Basin sites were P-deficient in July and one West Basin site and a river site were N-deficient in September. Between-instrument Fv/Fm comparisons did not show the expected 1:1 relationship. Fv/Fm from the PhytoPAM and WaterPAM were well-correlated with each other but not with nutrient deficiency. DivingPAM Fv/Fm did not correlate with the other PAM models, but correlated with P-deficiency. Spectral PAM fluorometers (PhytoPAM) can potentially resolve Fv/Fm down to phytoplankton group by additionally measuring accessory pigment fluorescence. The nutrient-induced fluorescent transient (NIFT) is the observation that Fv/Fm drops immediately and recovers when the limiting nutrient is reintroduced to nutrient-starved phytoplankton. A controlled laboratory experiment was conducted on a 2x2 factorial mixture design of P-deficient and P-sufficient Asterionella formosa and Microcystis aeruginosa cultures. Patterns consistent with published reports of NIFT were observed for P-deficient M. aeruginosa in mixtures; the pattern for A. formosa was less clear. This thesis showed that Fv/Fm by itself was not a reliable metric of N or P deficiency and care must be taken when interpreting results obtained by different PAM fluorometers. NIFT analysis using spectral PAM fluorometers may be able to discriminate P-deficiency in M. aeruginosa, and possibly other Cyanobacteria, in mixed communities.

Lake Erie has undergone a substantial amount of ecosystem changes over the past century; including cultural eutrophication and several invasions by industrious exotic species. Simple mass balance models for phosphorus have been useful in guiding policy decisions that led to reduced eutrophication, but new, confounding threats to the ecological health of Lake Erie continue to appear and lake managers continue to need useful tools to better understand the lake. As more complex ecological questions are asked to guide future management decisions, more complex ecological models are developed in an effort to provide some clues. The walleye fishery in Lake Erie is economically very important. Walleye recruitment has been highly variable from year to year since the 1990s. Modelling zooplankton is desired as a diagnostic tool for elucidating the quality of habitat – spatially and temporally – that is available to walleye in their vulnerable larval state. ELCOM-CAEDYM (or ELCD) is a 3-dimensional, coupled hydrodynamic and ecological model, which has been successfully applied to Lake Erie to model the nutrients and phytoplankton. The objectives of this study were to better understand the ELCOM-CAEDYM model of Lake Erie through a sensitivity analysis (SA), which has not been done before, and to explicitly simulate zooplankton in this model. An SA is important for determining which of the uncertain parameters have the greatest impact on the output variables. Due to the complexity of the CAEDYM model and the highly interdependent functions and variables modelled, a local SA (comparing changes in output by perturbing parameters one-at-a-time from some baseline configuration) was not desirable. Local SA’s ignore the possibility of a parameter’s effects being correlated to the status of other parameters. However, quantitative global methods are enormously computationally expensive for a complex model. The Lake Erie ELCD model simulates temperature, mixing, nutrient cycles, and phytoplankton dynamics. Phytoplankton are represented by 5 functional groups. With the explicit inclusion of 2 functional groups of zooplankton (copepods and cladocerans), the model uses over 300 function parameters in addition to requiring meteorological data and river inflow characteristics throughout the simulation. The model is set up with a 2-km grid over 40 layers with a 5-minute timestep from April 11 to September 1. This full simulation takes 6 days to complete. A quantitative global method to evaluate all parameters potentially significant to zooplankton would be impossible. The Morris method was selected for its streamlined global sampling procedure combined with the manageable computational demands of a one-at-a-time analysis. This method provides the relative sensitivity of diagnostic outputs to perturbed parameters. Ninety-one parameters were selected to be evaluated in 3680 simulations for the Morris SA. The selection of which parameters to evaluate and their assigned ranges are critical components in any SA. The ranges for parameters that represent a measurable quantity were assessed based on observed values in Lake Erie and other relevant studies. For some parameters, a measured realistic range was unknown. In these cases, values from relevant published models or judgements based on experience with the ELCD-Erie model were used to choose a suitable range. To assess the sensitivity of CAEDYM variables to parameters, DYRESM was substituted for ELCOM to vastly decrease the computation time of a single run. DYRESM is not suitable to model the entire lake due to the large size and irregular shape of the entire lake. Therefore, only the West Basin was modelled and analysed using DYCD. The West Basin was of special interest for a sensitivity analysis of CAEDYM parameters with respect to zooplankton because it is an important area for walleye larval development. DYCD output profiles for temperature, total chlorophyll a (TChla as a surrogate for total phytoplankton concentration) were similar in magnitudes and seasonal dynamics compared to ELCD outputs in deep West Basin stations. The sensitivity of zooplankton, TChla and TP to each parameter was assessed using two single value diagnostics: the simulated seasonal maximum and the simulated day on which peak maximum was reached. Zooplankton were sensitive to almost all of the zooplankton parameters perturbed in the analysis. This may indicate that modelling zooplankton is extremely complex, relying on many dynamic processes, or that evaluated ranges were not constrained well enough. An example of sensitivity to a poorly known parameter is the messy feeding coefficient. Reducing the uncertainty of this parameter would improve the confidence in the zooplankton assimilation submodel. Other parameters that stood out for being especially significant to zooplankton were: the respiration rate, mortality rate, internal phosphorus to carbon ratio, the temperature multiplier and standard temperature for feeding dynamics, and the half saturation constant. Most of these are easily explained as they directly aid or impede growth or they directly affect zooplankton losses. The most significant phytoplankton parameters on TChla and zooplankton outputs were, not surprisingly, the maximum growth rate and the mortality and respiration coefficient. Some particulate matter parameters proved to be important to outputs as well. More than 2500 of the 3680 parameter configurations resulted in unrealistic zooplankton simulations: peak values that did not much exceed initial conditions on the first day of the simulation. The SA exercise pinpointed a few configurations that resulted in reasonable peak zooplankton values and timing; these runs were used as a starting point for calibrating the ELCD model. Parameters were further manually adjusted by quickly checking their impacts on DYCD before applying them to ELCD. Post SA and minor calibration, the modelled zooplankton results were dramatically better than initial modelling attempts prior to the SA. Zooplankton concentrations throughout the lake were close to measured ranges and in some parts of the lake seasonal patterns were also similar to measured patterns. Modelled zooplankton results were least consistent with observations in the south west area of the lake: zooplankton were overestimated in late June-early July and they subsequently crashed and were underestimated in late July-August. It is supposed that this is due to higher grazing pressure from fish larvae in that area of the lake, which is not explicitly modelled. Although it is not anticipated that the south west seasonal zooplankton patterns will improve through parameter calibration (since predator effects are uniformly characterized throughout the lake by the same mortality factor) , further calibration is needed to improve results in the rest of the lake since copepods are generally overestimated and cladocerans generally underestimated. Phytoplankton groups must also be calibrated simultaneously to ensure that they are still operating within reasonable concentrations given more successful zooplankton simulations.

A lake-wide benthic survey of Lake Erie during summer 2002 indicated that Dreissena bugensis is the dominant dreissenid in Lake Erie, especially in the east basin where this species was found at every station but no Dreissena polymorpha were collected. Mean (±SD) densities of dreissenid mussels were comparable between the west (601±2,110/m2,n=49) and central (635±1,293/m2; n=41) basins, but were much greater in the east basin (9,480±11,173/m2;n=17). The greater variability in mussel density among stations and replicate samples in the central and west basins than in the east basin is attributable to the preponderance of fine-grained substrata in the nearshore, higher episodic rates of sediment deposition and periodic hypoxia in bottom waters. Although there was little change in lake-wide mean dreissenid densities between 1992 and 2002 (declining from ca. 2,636 individuals/m2 to 2,025 individuals/m2), basin-averaged shell-free dry tissue mass increased by almost four-fold from ca. 6.8±15.6 g /m2 to 24.7±71.3 g/m2 in the same interval. Up to 90% of this biomass is in the eastern basin. Other changes in 2002 include the virtual absence of mussels in the 3 to 12 mm size range, probably because of predation by round gobies, and an increase in the average size of mature mussels. The substantial changes observed between 1992 and 2002 suggest that dreissenid populations in Lake Erie were still changing rapidly in abundance and biomass, as well as species composition. The results of this survey suggest that a direct link between Dreissena spp. and hypolimnetic hypoxia in the central basin is unlikely. The dominant organisms of Lake Erie in 2002 were D. rostriformis bugensis (38%), Oligochaeta (33%), Chironomidae (18%), Sphaeriidae (2.7%), Amphipoda (2.3%) and Hydrozoa (2.2%). Mean invertebrate density was greater in the east basin, especially on hard substrates, than either the west or central basin. In the central basin, sites ≥5m supported greater numbers of organisms, than shallow (≤2m) sites in the nearshore wave zone. The greatest number of taxa were observed in the central basin, likely a resutlt of greater sampling effort there. Gammarus fasciatus comprised 80% of all amphipods, being most abundant on Dreissena-dominated hard substrates in the east basin. The introduced species, Echinogammarus ischnus occured at only 11 of 69 sites, and was the only amphipod found at 4 east basin sites, but at relatively low densities. Diporeia were not found in our survey. Hexagenia was collected at only 4 stations, all in the west basin. Chironomids were dominated by Tanytarsus and Chironomus, with Procladius, Dicrotendipes and Polypedilum also being relatively common. Oligochaete Trophic Index for 2002 indicates a similarly mesotrophic condition throughout the lake and marginal nutrient enrichment of sediments between years 1979-2002. Multivariate ordination of community data indicates clear separation of sites by year and basin as expected given the extirpation of Diporeia and the introduction of Dreissena and E. ischnus, but also reveals subtle changes in benthic structure over the last 2 decades. The benthic community of Lake Erie in 2002 does not likely represent an equilibrium condition.

Lake Erie has undergone a substantial amount of ecosystem changes over the past century; including cultural eutrophication and several invasions by industrious exotic species. Simple mass balance models for phosphorus have been useful in guiding policy decisions that led to reduced eutrophication, but new, confounding threats to the ecological health of Lake Erie continue to appear and lake managers continue to need useful tools to better understand the lake. As more complex ecological questions are asked to guide future management decisions, more complex ecological models are developed in an effort to provide some clues. The walleye fishery in Lake Erie is economically very important. Walleye recruitment has been highly variable from year to year since the 1990s. Modelling zooplankton is desired as a diagnostic tool for elucidating the quality of habitat – spatially and temporally – that is available to walleye in their vulnerable larval state. ELCOM-CAEDYM (or ELCD) is a 3-dimensional, coupled hydrodynamic and ecological model, which has been successfully applied to Lake Erie to model the nutrients and phytoplankton. The objectives of this study were to better understand the ELCOM-CAEDYM model of Lake Erie through a sensitivity analysis (SA), which has not been done before, and to explicitly simulate zooplankton in this model. An SA is important for determining which of the uncertain parameters have the greatest impact on the output variables. Due to the complexity of the CAEDYM model and the highly interdependent functions and variables modelled, a local SA (comparing changes in output by perturbing parameters one-at-a-time from some baseline configuration) was not desirable. Local SA’s ignore the possibility of a parameter’s effects being correlated to the status of other parameters. However, quantitative global methods are enormously computationally expensive for a complex model. The Lake Erie ELCD model simulates temperature, mixing, nutrient cycles, and phytoplankton dynamics. Phytoplankton are represented by 5 functional groups. With the explicit inclusion of 2 functional groups of zooplankton (copepods and cladocerans), the model uses over 300 function parameters in addition to requiring meteorological data and river inflow characteristics throughout the simulation. The model is set up with a 2-km grid over 40 layers with a 5-minute timestep from April 11 to September 1. This full simulation takes 6 days to complete. A quantitative global method to evaluate all parameters potentially significant to zooplankton would be impossible. The Morris method was selected for its streamlined global sampling procedure combined with the manageable computational demands of a one-at-a-time analysis. This method provides the relative sensitivity of diagnostic outputs to perturbed parameters. Ninety-one parameters were selected to be evaluated in 3680 simulations for the Morris SA. The selection of which parameters to evaluate and their assigned ranges are critical components in any SA. The ranges for parameters that represent a measurable quantity were assessed based on observed values in Lake Erie and other relevant studies. For some parameters, a measured realistic range was unknown. In these cases, values from relevant published models or judgements based on experience with the ELCD-Erie model were used to choose a suitable range. To assess the sensitivity of CAEDYM variables to parameters, DYRESM was substituted for ELCOM to vastly decrease the computation time of a single run. DYRESM is not suitable to model the entire lake due to the large size and irregular shape of the entire lake. Therefore, only the West Basin was modelled and analysed using DYCD. The West Basin was of special interest for a sensitivity analysis of CAEDYM parameters with respect to zooplankton because it is an important area for walleye larval development. DYCD output profiles for temperature, total chlorophyll a (TChla as a surrogate for total phytoplankton concentration) were similar in magnitudes and seasonal dynamics compared to ELCD outputs in deep West Basin stations. The sensitivity of zooplankton, TChla and TP to each parameter was assessed using two single value diagnostics: the simulated seasonal maximum and the simulated day on which peak maximum was reached. Zooplankton were sensitive to almost all of the zooplankton parameters perturbed in the analysis. This may indicate that modelling zooplankton is extremely complex, relying on many dynamic processes, or that evaluated ranges were not constrained well enough. An example of sensitivity to a poorly known parameter is the messy feeding coefficient. Reducing the uncertainty of this parameter would improve the confidence in the zooplankton assimilation submodel. Other parameters that stood out for being especially significant to zooplankton were: the respiration rate, mortality rate, internal phosphorus to carbon ratio, the temperature multiplier and standard temperature for feeding dynamics, and the half saturation constant. Most of these are easily explained as they directly aid or impede growth or they directly affect zooplankton losses. The most significant phytoplankton parameters on TChla and zooplankton outputs were, not surprisingly, the maximum growth rate and the mortality and respiration coefficient. Some particulate matter parameters proved to be important to outputs as well. More than 2500 of the 3680 parameter configurations resulted in unrealistic zooplankton simulations: peak values that did not much exceed initial conditions on the first day of the simulation. The SA exercise pinpointed a few configurations that resulted in reasonable peak zooplankton values and timing; these runs were used as a starting point for calibrating the ELCD model. Parameters were further manually adjusted by quickly checking their impacts on DYCD before applying them to ELCD. Post SA and minor calibration, the modelled zooplankton results were dramatically better than initial modelling attempts prior to the SA. Zooplankton concentrations throughout the lake were close to measured ranges and in some parts of the lake seasonal patterns were also similar to measured patterns. Modelled zooplankton results were least consistent with observations in the south west area of the lake: zooplankton were overestimated in late June-early July and they subsequently crashed and were underestimated in late July-August. It is supposed that this is due to higher grazing pressure from fish larvae in that area of the lake, which is not explicitly modelled. Although it is not anticipated that the south west seasonal zooplankton patterns will improve through parameter calibration (since predator effects are uniformly characterized throughout the lake by the same mortality factor) , further calibration is needed to improve results in the rest of the lake since copepods are generally overestimated and cladocerans generally underestimated. Phytoplankton groups must also be calibrated simultaneously to ensure that they are still operating within reasonable concentrations given more successful zooplankton simulations.

A lake-wide benthic survey of Lake Erie during summer 2002 indicated that Dreissena bugensis is the dominant dreissenid in Lake Erie, especially in the east basin where this species was found at every station but no Dreissena polymorpha were collected. Mean (±SD) densities of dreissenid mussels were comparable between the west (601±2,110/m2,n=49) and central (635±1,293/m2; n=41) basins, but were much greater in the east basin (9,480±11,173/m2;n=17). The greater variability in mussel density among stations and replicate samples in the central and west basins than in the east basin is attributable to the preponderance of fine-grained substrata in the nearshore, higher episodic rates of sediment deposition and periodic hypoxia in bottom waters. Although there was little change in lake-wide mean dreissenid densities between 1992 and 2002 (declining from ca. 2,636 individuals/m2 to 2,025 individuals/m2), basin-averaged shell-free dry tissue mass increased by almost four-fold from ca. 6.8±15.6 g /m2 to 24.7±71.3 g/m2 in the same interval. Up to 90% of this biomass is in the eastern basin. Other changes in 2002 include the virtual absence of mussels in the 3 to 12 mm size range, probably because of predation by round gobies, and an increase in the average size of mature mussels. The substantial changes observed between 1992 and 2002 suggest that dreissenid populations in Lake Erie were still changing rapidly in abundance and biomass, as well as species composition. The results of this survey suggest that a direct link between Dreissena spp. and hypolimnetic hypoxia in the central basin is unlikely. The dominant organisms of Lake Erie in 2002 were D. rostriformis bugensis (38%), Oligochaeta (33%), Chironomidae (18%), Sphaeriidae (2.7%), Amphipoda (2.3%) and Hydrozoa (2.2%). Mean invertebrate density was greater in the east basin, especially on hard substrates, than either the west or central basin. In the central basin, sites ≥5m supported greater numbers of organisms, than shallow (≤2m) sites in the nearshore wave zone. The greatest number of taxa were observed in the central basin, likely a resutlt of greater sampling effort there. Gammarus fasciatus comprised 80% of all amphipods, being most abundant on Dreissena-dominated hard substrates in the east basin. The introduced species, Echinogammarus ischnus occured at only 11 of 69 sites, and was the only amphipod found at 4 east basin sites, but at relatively low densities. Diporeia were not found in our survey. Hexagenia was collected at only 4 stations, all in the west basin. Chironomids were dominated by Tanytarsus and Chironomus, with Procladius, Dicrotendipes and Polypedilum also being relatively common. Oligochaete Trophic Index for 2002 indicates a similarly mesotrophic condition throughout the lake and marginal nutrient enrichment of sediments between years 1979-2002. Multivariate ordination of community data indicates clear separation of sites by year and basin as expected given the extirpation of Diporeia and the introduction of Dreissena and E. ischnus, but also reveals subtle changes in benthic structure over the last 2 decades. The benthic community of Lake Erie in 2002 does not likely represent an equilibrium condition.