Lead in Lake Michigan and Green Bay surficial sediments

Sediment samples were collected in 1987–1990 from Green Bay and in 1994–1996 from Lake Michigan. Surficial sediments (0–1 cm) from both locations were analyzed for lead for the purpose of describing the horizontal variation of lead in 1994–1996 Lake Michigan and 1987–1990 Green Bay sediments, estimating lead fluxes to surficial sediments, and comparing results to earlier studies. With Lake Michigan concentrations ranging from below the method detection limit to 180 μg/g, the surficial sediments had mean and median lead concentrations of 70 μg/g and 64 μg/g, respectively. Lead concentrations in Green Bay surficial sediments were similar to those in Lake Michigan and ranged between the method detection limit and 160 μg/g. For the bay, mean and median concentrations were 58 and 59 μg/g, respectively. Surficial lead concentrations were highest in the Southern, Waukegan, and Grand Haven basins of Lake Michigan and in the central region of Green Bay in the vicinity of Chambers Island. For Lake Michigan and Green Bay, dated sediment cores illustrate the decline in lead concentrations during the last 30 and 10 years, respectively. Lead fluxes ranged between < 0.049 and 7.2 μg/cm²/yr for Green Bay and between 0.47 and 20 μg/cm²/yr for Lake Michigan. Lead fluxes to Lake Michigan were lower than those reported for 1972. These are the most comprehensive fluxes of lead to Lake Michigan and Green Bay surficial sediments reported to date.     

Shell-free Biomass and Population Dynamics of Dreissenids in Offshore Lake Michigan, 2001–2003

The USGS-Great Lakes Science Center has collected dreissenid mussels annually from Lake Michigan since zebra mussels (Dreissena polymorpha) became a significant portion of the bottomtrawl catch in 1999. For this study, we investigated dreissenid distribution, body mass, and recruitment at different depths in Lake Michigan during 2001–2003. The highest densities of dreissenid biomass were observed from depths of 27 to 46 m. The biomass of quagga mussels (Dreissena bugensis) increased exponentially during 2001–2003, while that of zebra mussels did not change significantly. Body mass (standardized for a given shell length) of both species was lowest from depths of 27 to 37m, highest from 55 to 64 m, and declined linearly at deeper depths during 2001–2003. Recruitment in 2003, as characterized by the proportion of mussels < 11 mm in the catch, varied with depth and lake region. For quagga mussels, recruitment declined linearly with depth, and was highest in northern Lake Michigan. For zebra mussels, recruitment generally declined non-linearly with depth, although the pattern was different for north, mid, and southern Lake Michigan. Our analyses suggest that quagga mussels could overtake zebra mussels and become the most abundant mollusk in terms of biomass in Lake Michigan.     

Status and trends of the Asian clam (Corbicula fluminea) in the lower Fox River and Green Bay

In the upper Great Lakes region, survival and population growth of the non-indigenous Asian clam Corbicula fluminea has been limited by cold climates that cause severe overwinter mortality. At these northern latitudes, Asian clam populations are often limited to thermal refugia – particularly warmwater discharges from industrial facilities. Several such facilities exist in the lower Fox River in Green Bay. Asian clams were first documented in the lower Fox River in 1999 and were extensively surveyed near the river mouth in 2011, but the few individuals found were restricted to the warmwater discharge from the Pulliam Power Plant. We performed a follow-up survey during 2017 to re-assess the population status of Asian clams in the lower Fox River at four industrial discharges, including the Pulliam Power Plant. We found more widespread evidence of Asian clams throughout the lower Fox River than previous surveys, but only one live individual was captured. We suspect that the back-to-back severe winters of 2013–2014 and 2014–2015 caused widespread overwinter mortality. Our investigation highlights the significant challenges for establishment of Asian clam populations in the upper Great Lakes region, and provides an example of a potential invasive species struggling to establish a viable population in a hostile climate.     

A comparative examination of recent changes in nutrients and lower food web structure in Lake Michigan and Lake Huron

The lower food webs of Lake Huron and Lake Michigan have experienced similar reductions in the spring phytoplankton bloom and summer populations of Diporeia and cladocerans since the early 2000s. At the same time phosphorus concentrations have decreased and water clarity and silica concentrations have increased. Key periods of change, identified by using a method based on sequential t-tests, were 2003–2005 (Huron) and 2004–2006 (Michigan). Estimated filtration capacity suggests that dreissenid grazing would have been insufficient to directly impact phytoplankton in the deeper waters of either lake by this time (mid 2000s). Despite some evidence of decreased chlorophyll:TP ratios, consistent with grazing limitation of phytoplankton, the main impact of dreissenids on the offshore waters was probably remote, e.g., through interception of nutrients by nearshore populations. A mass balance model indicates that decreased phosphorus loading could not account for observed in-lake phosphorus declines. However, model-inferred internal phosphorus dynamics were strongly correlated between the lakes, with periods of increased internal loading in the 1990s, and increased phosphorus loss starting in 2000 in Lake Michigan and 2003 in Lake Huron, prior to dreissenid expansion into deep water of both lakes. This suggests a limited role for deep populations of dreissenids in the initial phosphorus declines in the lakes, and also suggests a role for meteorological influence on phosphorus dynamics. The high synchrony in lower trophic level changes between Lake Michigan and Lake Huron suggests that both lakes should be considered when investigating underlying causal factors of these changes.     

Evidence of persistent,recurring summertime hypoxia in Green Bay,Lake Michigan

Six years (2009–2015) of temperature and dissolved oxygen profile data show hypoxic conditions are common in the bottom waters of southern Green Bay, Lake Michigan during the summer. Depleted oxygen concentrations (<5?mg?L^?1) affect nearly 70% of the 38 stations sampled representing an area of ~500–600?km². Stratification typically lasts 2+?months, from late June to early September, and some stations exhibit bottom water hypoxia (<2?mg?L^?1) at a frequency of nearly 25% when sampled during this period. A monitoring program initiated in 1986 by the Green Bay Metropolitan Sewerage District has provided a 23?year, recreational season record (May–September) of continuous (15?min interval) in situ bottom water oxygen and temperature measurements at the Entrance Light station of the Green Bay navigational channel. The duration of the hypoxic season ranges from 2?weeks to over 3?months at this shallow 7?m offshore site. This variability likely results from a combination of thermal stratification, oxygen consumption in deeper waters of the bay, and physical forcing mechanisms that drive cool, oxygen depleted, bottom waters on a southerly trajectory across this sensor. These data suggest the duration of hypoxic conditions may have increased during the stratified season in recent years. Hypoxia in the bay would also appear to be sensitive to relatively small changes in these forces, particularly changes in organic carbon loading and the duration of stratification.     

A mesocosm investigation of the effects of quagga mussels (Dreissena rostriformis bugensis) on Lake Michigan zooplankton assemblages

Dreissenid mussels are known to disrupt the base of the food web by filter feeding on phytoplankton; however, they may also directly ingest zooplankton thereby complicating their effects on plankton communities. The objective of this study was to quantify the effects of quagga mussel feeding on the composition and size structure of Lake Michigan zooplankton assemblages. Two mesocosm (six 946 L tanks) experiments were conducted in summer 2013, using quagga mussels and zooplankton collected near Beaver Island, MI, to examine the response of zooplankton communities to the presence and absence of mussels (experiment 1) and varying mussel density (experiment 2). Mesocosms were sampled daily and zooplankton taxa were enumerated and sized using microscopy and FlowCAM® imaging. In experiment 1, the presence of quagga mussels had a rapid negative effect on veliger and copepod nauplii abundance, and a delayed negative effect on rotifer abundance. In experiment 2, mussel density had a negative effect on veliger, nauplii, and copepodite abundance within 24 h. Multivariate analyses revealed a change in zooplankton community composition with increasing mussel density. Ten zooplankton taxa decreased in abundance and frequency as quagga mussel density increased: except for the rotifer Trichocerca sp., treatments with higher mussel densities (i.e., 1327, 3585, and 5389 mussels/m²) had the greatest negative effect on small-bodied zooplankton (≤ 128 μm). This study confirms results from small-scale (≤ 1 L) experiments and demonstrates that quagga mussels can alter zooplankton communities at mesoscales (~ 1000 L), possibly through a combination of direct consumption and resource depletion.     

By-catch in the Saginaw Bay,Lake Huron commercial trap net fishery

This study provides species-specific catch and baseline mortality estimates of non-target species (by-catch) for the Saginaw Bay, Lake Huron commercial trap net fishery. By-catch can represent a significant mortality source that is often unknown. By-catch and by-catch mortality rates in the Saginaw Bay commercial trap net fishery, which primarily targets lake whitefish (Coregonus clupeaformis), yellow perch (Perca flavescens), and channel catfish (Ictalurus punctatus), are currently unknown. From May through August 2010, we observed onboard commercial trap net vessels and took species-specific counts of by-catch and estimated initial by-catch mortality (i.e., morbid or floating fish). The high levels of walleye (Sander vitreus) catch and mortality observed within inner Saginaw Bay have not been previously documented in the Laurentian Great Lakes. Walleye by-catch averaged 127.3 individuals per trap net lift and 42% of those caught were morbid. The levels of lake trout (Salvelinus namaycush) catch observed were within the range observed in previous studies, but mortality (percent) was higher than has been previously observed. Lake trout by-catch averaged 39.4 individuals per lift and 39.2% of those were morbid. Through the use of generalized linear models, this analysis also indicated factors that most influenced catch of non-target species including time of year and soak time (i.e., time interval between trap net lifts). Surface water temperature and trap net depth most influenced mortality. These results may inform fishers and fisheries managers and highlight the need for comprehensive by-catch monitoring throughout the Great Lakes.     

Variations in chemical speciation and reactivity of phosphorus between suspended-particles and surface-sediment in seasonal hypoxia-influenced Green Bay

Water, suspended-particles, and surface-sediment samples were collected from Green Bay, Lake Michigan, for the measurements of phosphorus (P) species, including dissolved/particulate-P, inorganic/organic-P, and five different forms of particulate-P, namely exchangeable- or labile-P (Ex-P), iron-bound-P (Fe-P), biogenic-apatite and/or CaCO₃-associated-P (CFA-P), organic-P (Org-P) and detrital-apatite-P (Detr-P) to elucidate their reactivity and transformation pathways in the water column. Suspended particles contained mainly Ex-P (25?±?15%), Fe-P (28?±?12%) and Org-P (29?±?7%). In contrast, Detr-P (34?±?10%) and Org-P (36?±?12%) were the predominant P species in surface sediment. Contents of Ex-P, Fe-P, Org-P and CFA-P decreased consistently from suspended-particles to surface-sediment, but an increase was observed for the Detr-P, indicating a net loss of Ex-P, Fe-P, Org-P and CFA-P from particulate into dissolved phase. Such active regeneration of P in the water column between particulate and dissolved phases may serve as an internal phosphate source in Green Bay, especially under hypoxic conditions. Degradation of organic matter in south central bay areas seemed to promote hypoxia and enhance the reductive-dissolution of Fe-P and preservation of Org-P under low-oxygen conditions in the central bay. Overall, Ex-P, Fe-P, CFA-P and Org-P species, which comprised up to 50–90% of total particulate-P, can be collectively considered as potentially-bioavailable-P (BAP). Under low-phosphate (0.022?±?0.014?μM in Green Bay) and summer low-oxygen/hypoxic conditions, suspended-particles may release up to 71% of their BAP before deposited in sediment although the BAP regeneration decreased along the south-north transect in Green Bay.     

Stock structure,dynamics, demographics,and movements of walleyes spawning in four tributaries to Green Bay

To test assumptions related to the current conceptual model for walleye Sander vitreus management in Green Bay, we evaluated whether: 1) spawning aggregations in the Fox, Menominee, Oconto, and Peshtigo rivers represent genetically distinct stocks; 2) population dynamics and demographics vary among walleye spawning at these locations; 3) walleye spawning in these rivers contribute to the fishery in northern Green Bay, and 4) walleye spawning in these rivers exhibit spawning site fidelity or if they stray among rivers. Genetic differentiation among the four tributaries was low and sex-specific total length (TL), mean TL at age 5, and age-class diversity were generally similar among rivers and observed differences were not consistent. Movements of walleye inferred from angler tag returns suggest that walleye spawning (and tagged) in the four tributaries typically remain within southern Green Bay; however, this assertion may be confounded by the distribution of angling effort that provides tag recoveries. Straying rates among rivers ranged from 0 to 23% and were likely sufficient to preclude genetic differentiation among stocks. Collectively, results suggest that walleye spawning in the Fox, Menominee, Oconto, and Peshtigo rivers do not function as separate stocks and do not significantly contribute to the fishery outside of southern Green Bay. The primary assumption of the current conceptual model that remains to be tested is whether the walleye fishery in southern Green Bay is supported primarily by fish spawning in these four rivers, or if there are substantial contributions from fish spawning at other unknown locations.     

Spatial genetic structure and recruitment dynamics of burbot (Lota lota) in Eastern Lake Michigan and Michigan tributaries

Burbot (Lota lota) are the only freshwater member of the Cod like (Lotidae) family that have a circumpolar distribution and occupy the widest geographic distribution of all Laurentian Great Lakes fish species. Information regarding burbot spatial genetic structure and recruitment dynamics is critical for the development of effective management strategies. Although burbot are a species of conservation concern throughout their range, little demographic or behavioral information exists. We estimated levels of genetic diversity within, and the degree of spatial population structure between samples collected from Lake Michigan and tributaries of the Manistee River, MI. Measures of genetic diversity across 10 microsatellite loci were moderately high. Disparities between adult groups sampled in Lake Michigan and the Manistee River were notable for observed heterozygosity (0.662 vs 0.488) and allelic richness (11.7 vs 6.6). Significant levels of inter-population variance in microsatellite allele frequencies (F_ST 0.154 to 0.208) were detected between Lake Michigan and the Manistee River samples. Results indicate reproductive isolation between what plausibly may be riverine and lacustrine spawning life history types. Pedigree analyses for three cohorts sampled in the Manistee River revealed that a sizeable number of adults contributed reproductively to multiple cohorts, indicating spawning philopatry. While data were collected from restricted areas in lacustrine and river habitats, analyses revealing microgeographic genetic structuring, potentially attributed to life history polymorphisms, have significant implications for burbot management in the Great Lakes.     

A short-term look at potential changes in Lake Michigan slimy sculpin diets

Diporeia hoyi and Mysis relicta are the most important prey items of slimy sculpins (Cottus cognatus) in the Great Lakes. Slimy sculpins were collected from dreissenid-infested bottoms off seven Lake Michigan ports at depths of 27–73 m in fall 2003 to study their lake-wide diets. Relatively large dreissenid biomass occurred at depths of 37- and 46-m. Quagga mussels (Dreissena bugnesis) composed at least 50% of dreissenid biomass at Manistique, Saugatuck, and Sturgeon Bay. Mysis accounted for 82% of the sculpin diet by dry weight at eastern Lake Michigan while Diporeia composed 54–69% of the diet at western Lake Michigan and dominated the diets of slimy sculpins at all sites deeper than 46 m. In northern Lake Michigan, this diet study in new sites showed that slimy sculpin consumed more prey with low energy contents, especially chironomids, than Mysis and Diporeia in shallow sites (depth <55 m). We recommend diet studies on sedentary benthic fishes to be conducted along perimeters of the Great Lakes to observe changes in their diets that may be impacted by changing benthic macroinvertebrate communities.     

Dreissenid veliger dynamics along a nearshore to offshore transect in Lake Michigan

Dreissenid mussel veligers compose a substantial component of pelagic biomass in the Great Lakes, yet their dynamics are poorly understood. To evaluate seasonal, spatial, and inter-annual variation in veliger density, we used a 64-μm mesh plankton net (2008, 2013–2016) and a 153-μm mesh plankton net (2007–2016) to collect dreissenid veligers at nearshore (15–25?m depth), transitional (45?m) and offshore (93–110?m) sites in southeast Lake Michigan during March–December. We also evaluated trends in density of recently settled mussels relative to veliger abundance and the density of the standing stock of adult mussels. Veliger density peaked during both summer and fall at all sites, but peak densities in summer were generally higher nearshore, whereas peak densities in the fall were generally higher offshore. The density of veligers in the 153-μm net was overall 28% of that in the 64-μm net, but there was high variability in this comparison among months. Smaller veligers were much more abundant in the 64-μm net, but there was little difference in the size distribution and abundance between nets for the 210–300?μm size classes. Thus, the 153-μm net could still be a useful tool for assessing density trends of larger veligers just prior to their settlement. Newly settled mussels (≤2?mm) were most abundant in summer or fall at the nearshore and offshore sites but were nearly absent at the transitional site despite the high density of veligers there. Factors other than veliger density must play an important role in mussel recruitment.     

Distribution of seston and nutrient concentrations in the eastern basin of Lake Erie pre- and post-dreissenid mussel invasion

Increased human population growth, reduction of phosphorus (P) loading, and the invasion of dreissenid mussels may have changed the spatial pattern and relationships between the nearshore and the offshore seston and nutrient concentrations in the eastern basin of Lake Erie over the past 30 years. We compared seston characteristics, nutrient concentrations, and phytoplankton nutrient status between nearshore and offshore zones in years before (1973–1985) and after (1990–2003) the dreissenid invasion. In 1973 (the only pre-dreissenid year nearshore data was collected), chlorophyll a (chla) and nutrient concentrations were higher nearshore than offshore. In post-dreissenid years, nearshore chla concentrations became significantly lower than the offshore, while carbon (C):chla ratios became higher, which was related to mussel grazing and possibly photoacclimation. Phosphorus deficiency in the phytoplankton increased over the 30-year period, and in the post-dreissenid years was less acute in the nearshore than offshore. Mean water column irradiance became higher in the nearshore relative to the offshore in the post-dreissenid years. The nutrient changes and phytoplankton physiology were consistent with the expected effects of nutrient cycling by mussels and diminished demand by phytoplankton despite increased demand from benthic algae in the nearshore. This basin-scale study suggests that dreissenid mussel invasion can be associated with alterations in the spatial pattern of water column properties in large lakes even on open coasts with vigorous circulation and exchange.     

Trends in Mysis diluviana abundance in the Great Lakes, 2006–2016

With the large Diporeia declines in lakes Michigan, Huron, and Ontario, there is concern that a similar decline of Mysis diluviana related to oligotrophication and increased fish predation may occur. Mysis density and biomass were assessed from 2006 to 2016 using samples collected by the Great Lakes National Program Office's biomonitoring program in April and August in all five Great Lakes. Summer densities and biomasses were generally greater than spring values and both increased with bottom depth. There were no significant time trends during these 10–11 years in lakes Ontario, Michigan, or Huron, but there was a significant increase in Lake Superior. Density and biomass were highest in lakes Ontario and Superior, somewhat lower in Lake Michigan, and substantially lower in Lake Huron. A few Mysis were collected in eastern Lake Erie, indicating a small population in the deep basin of that lake. On average, mysids contributed 12–18% (spring-summer, Michigan), 18–14% (spring-summer, Superior), 30–13% (spring-summer, Ontario), and 3% (Huron) of the total open-water crustacean biomass. Size distributions consisted of two peaks, indicating a 2-year life cycle in all four of the deep lakes. Mysis were larger in Lake Ontario than in lakes Michigan, Superior, and Huron. Comparisons with available historic data indicated that mysid densities were higher in the 1960s–1990s (5 times higher in Huron, 2 times higher in Ontario, and around 40% higher in Michigan and Superior) than in 2006–2016.