Diet complexity of Lake Michigan salmonines: 2015–2016

In Lake Michigan, the unintended introduction of invasive species (e.g., zebra mussel, Dreissena polymorpha; quagga mussel, D. rostriformis bugensis; round goby, Neogobius melanostomus) and reduced nutrient loading has altered nutrient dynamics, system productivity, and community composition over the past two decades. These factors, together with sustained predation pressure, have contributed to declines of several forage fish species, including alewife (Alosa pseudoharengus), which has dominated diets of the five primary salmonine species of Lake Michigan for the last 50 years. Salmonines that have inflexible, less complex diets may struggle if alewife declines continue. We analyzed stomach contents of salmonines collected throughout the main basin of Lake Michigan in 2015 and 2016 to investigate diet composition, diet diversity, and individual variation of alewife lengths consumed. Chinook salmon (Oncorhynchus tshawytscha) almost exclusively consumed alewife and had lower diet diversities compared to the other four species, which consumed relatively high frequencies of round goby (brown trout, Salmo trutta; lake trout, Salvelinus namaycush), aquatic invertebrates (coho salmon, Oncorhynchus kisutch) and terrestrial invertebrates (rainbow trout, Oncorhynchus mykiss) along with alewife. Although clear spatio-temporal feeding patterns existed, much of the variation in diet composition and diet diversity was expressed at the individual level. Salmonine populations consumed the entire size range of alewife that were available, whereas individual stomachs tended to contain a narrow range of alewife sizes. Due to their reliance on alewife, it is likely that Chinook salmon will be more negatively impacted than other salmonine species if alewife abundance continues to decline in Lake Michigan.     

Total and methyl mercury accumulation in 1994–1995 Lake Michigan lake trout and forage fish

As part of the Lake Michigan Mass Balance Project, total and methyl mercury were determined for lake trout (Salvelinus namaycush) and five forage fish species collected from Lake Michigan near Saugatuck, Michigan, and Port Washington, Sheboygan Reef, and Sturgeon Bay, Wisconsin, between 1994 and 1995. With a mean concentration of 179 ng/g wet wt., whole lake trout total mercury (Hg_T) concentrations ranged between 27.6 and 348 ng/g wet wt. For combined sites, 1–4 yrs, 5–6 yrs, 7–11 yrs, and 12–15 yrs lake trout mean Hg_T concentrations were 73.7, 130, 212, and 280 ng/g, respectively. Forage fish species alewife (Alosa pseudoharengus), bloater (Coregonus hoyi), slimy sculpin (Cottus cognatus), deepwater sculpin (Myoxocephalus thompsoni), and rainbow smelt (Osmerus mordax) had mean Hg_T concentrations of 63.8, 55.3, 36.7, 51.4, and 35.2 ng/g wet wt., respectively. With the exception of alewife, bloater, and slimy sculpin, all fish species contained approximately 100% methyl mercury (MeHg). Field bioaccumulation factors (BAF) were consistent with a Lake Michigan food chain that is more efficient at transferring MeHg to higher trophic levels than some inland lakes. This and other studies of lake trout from Lake Michigan document decreasing Hg_T concentrations in lake trout from 1971 to 1985 and constant or increasing concentrations between 1985 and 2000. These observations were supported by a similar trend in Lake Michigan Hg sediment fluxes. To our knowledge, this is the most intense two year study of mercury in fish for any Great Lake or other large fresh water system and is one of the most complete studies of mercury cycling in the Lake Michigan food chain.     

Seasonal variation in light penetration and subsurface chlorophyll-α in southern Lake Michigan observed by a glider

The Cooperative Institute for Great Lakes Research (CIGLR) in collaboration with the Great Lakes Observing System and National Oceanic and Atmospheric Administration, Great Lakes Environmental Research Laboratory (NOAA GLERL) deployed an autonomous underwater glider in southern Lake Michigan several times per year between 2012 and 2019 to collect offshore (>30 m depth) limnological measurements, including temperature, photosynthetically active radiation (beginning during 2015), and chlorophyll-α fluorescence. From these data, we calculated mixed layer depth, several measures of light penetration (diffuse attenuation coefficient, first optical depth, euphotic zone depth), and depth of the subsurface chlorophyll-α maxima. During summer, mean offshore mixed layer depth was typically 10–15 m, K_d for PAR was 0.1–0.17 m^?1, first optical depth was 6–9 m, euphotic zone depth was 35–40 m, and depth of subsurface chlorophyll-α maxima was 30–35 m. We also observed substantial spatial and temporal variation in these values across the basin and within and among seasons. Glider-based observations provide a wider horizontal and vertical perspective than other methods (e.g., ship- and satellite-based observations, buoys, and fixed moorings), and are therefore a valuable, complementary tool for Great Lakes limnology. The set of observations reported here provide seasonal and basin-scale information that may help to identify anomalies useful for future glider-assisted investigation into the role of biophysical processes in Great Lakes limnology and ecology.     

Mathematical modeling of bacteria–virus interactions in Lake Michigan incorporating phosphorus content

Viruses are ubiquitous in aquatic ecosystems and play a critical role in nutrient cycling because viral lysis of cells releases phosphorus, an essential and often limiting nutrient. Previous models have examined dynamics of bacterial and viral communities, but with limited analysis and without explicit consideration of nutrients. A recent model (Fuhrman et al., 2011; Math. Comp. Model. 53, 716-730) incorporated internal nutrient content of bacteria and viruses. In the present study, we modified and tested the model with data from natural planktonic communities from Lake Michigan. Replicate 20 L water samples (135 μm screened), were either untreated or enriched with 8 μM phosphate, and bacterial and viral abundance, chl a fluorescence, and phosphorus (total and dissolved inorganic) were monitored for two weeks. Fuhrman et al.'s model (modified to include phytoplankton) was applied to the data, fitting burst size, lytic latent period and virus decay rate. For enriched samples, model fits were good and parameters were consistent with measurements in other freshwater ecosystems. However, for unenriched samples, where nutrient concentrations approached detection limits, model fits were relatively poorer. The model predicted similar viral decay rates but higher burst sizes and longer latent period in phosphorus-limited versus enriched conditions, underlining the potential importance of nutrients in host–virus interactions. The model is likely to be most useful in meso- to eutrophic systems; requirements for future model development and parameter estimation for application to oligotrophic lakes are discussed.     

Fleet Dynamics of the Commercial Lake Trout Fishery in Michigan Waters of Lake Superior during 1929–1961

Understanding fishing fleet dynamics is important when using fishery dependent data to infer the status of fish stocks. We analyzed data from mandatory catch reports from the commercial lake trout (Salvelinus namaycush) fishery in Michigan waters of Lake Superior during 1929–1961, a period when lake trout populations collapsed through the combined effects of overfishing and sea lamprey (Petromyzon marinus) predation. The number of full-time fishermen increased during 1933–1943 and then decreased during 1943–1957. Addition of new fishermen was related to past yield, market prices, World War II draft exemptions, and lost fishing opportunities in Lake Huron and Lake Michigan. Loss of existing fishermen was related to declining lake trout density. Large mesh (≥ 114-mm stretch-measure) gill net effort increased during 1929–1951 because fishermen fished more net inshore as lake trout density declined, even though catch per effort (CPE) was often higher in deeper waters. The most common gill net mesh size increased from 114-mm to 120-mm stretch-measure during 1929–1957, as lake trout growth increased. More effort was fished inshore than offshore and the amount of inshore effort was less variable over time than offshore effort. Relatively stable yield was maintained by increasing gill net effort and by moving some effort to better grounds. Because fishing-up caused yield and CPE to remain high despite declining lake trout abundance, caution must be used when basing goals for lake trout restoration on historical fishery indices.     

Fatty acids reveal salmonine – prey relationships in Lake Michigan

Lake Michigan salmon and trout populations are important species for recreational fisheries and food web management, and are largely supported through stocking efforts, with varying degrees of natural recruitment. Ongoing fisheries management of these salmonine populations is dictated by relationships between predator and prey abundance as well as community structure within the lake. However, while prey fish biomass has declined, and species composition has changed in recent decades, knowledge of prey consumption by the salmonine community has lagged. Herein, we explore trophic relationships using fatty acids profiles, which offer insights into the foraging habits and energy pathways relied on over weeks to months prior to collection. Fatty acids of the prey base for salmonines in Lake Michigan indicate a gradient of foraging habits that range from pelagic (typified by alewife and rainbow smelt) versus benthic (i.e., slimy sculpin and round goby) resource use. Fatty acids implied that there was more variation in foraging habits among individual lake trout and brown trout compared to Chinook salmon, coho salmon and rainbow trout, which appeared to all rely almost exclusively on pelagic prey. Fatty acid profiles also indicated size-based shifts in foraging habits; for example, larger lake trout consuming a greater proportion of benthic prey than smaller individuals. Data herein suggest that Chinook and coho salmon, as well as rainbow trout, are more likely to experience competitive interactions during times of low pelagic prey-fish abundance in Lake Michigan, whereas brown and lake trout are able to utilize benthic resources to a greater degree.     

Trends of PCB concentrations in Lake Michigan coho and chinook salmon, 1975–2010

The manufacture and use of polychlorinated biphenyls (PCBs) was banned in the United States in 1977 after it was determined that these compounds adversely affect animals and humans. The Wisconsin Department of Natural Resources has quantified total PCB concentrations in Lake Michigan chinook (n = 765) and coho (n = 393) salmon (Oncorhynchus tshawytscha and Oncorhynchus kisutch, respectively) filets since 1975. We analyzed these data to estimate trends in PCB concentrations in these fish (1975–2010). We used generalized linear models with a gamma error distribution and log link fit to the untransformed concentrations. Trend patterns were examined using graphical smoothing and generalized additive models. We identified a candidate set of models that included time trend and other predictor variables. Using the Akaike Information Criterion to select among models we found the best models for both species included piecewise linear time trends, total body length, % lipid, and collection season as predictor variables. The intersection of the two trends was 1985 for chinook salmon and 1984 for coho salmon. PCB concentrations in both species increased with body length and % lipid, and were higher for individuals caught in the fall. Our data reveals a dramatic decline in PCB concentrations of − 16.7% and − 23.9% per year for chinook and coho, respectively, up until the intersection year likely reflecting implementation of restrictions on Aroclor-based PCBs. After the intersection year to 2010, PCB concentrations declined at an annual rate of − 4.0% (95% CI: − 4.4% to − 3.6%) and − 2.6% (95% CI: − 3.3% to − 1.9%) for chinook and coho, respectively.     

Landscape-scale modeling of water quality in Lake Superior and Lake Michigan watersheds: How useful are forest-based indicators?

The Great Lakes watersheds have an important influence on the water quality of the nearshore environment, therefore, watershed characteristics can be used to predict what will be observed in the streams. We used novel landscape information describing the forest cover change, along with forest census data and established land cover data to predict total phosphorus and turbidity in Great Lakes streams. In Lake Superior, we modeled increased phosphorus as a function of the increase in the proportion of persisting forest, forest disturbed during 2000–2009, and agricultural land, and we modeled increased turbidity as a function of the increase in the proportion of persisting forest, forest disturbed during 2000–2009, agricultural land, and urban land. In Lake Michigan, we modeled increased phosphorus as a function of ecoregion, decrease in the proportion of forest disturbed during 1984–1999 and watershed storage, and increase in the proportion of urban land, and we modeled increased turbidity as a function of ecoregion, increase in the proportion of forest disturbed during 2000–2009, and decrease in the proportion softwood forest. We used these relationships to identify priority areas for restoration in the Lake Superior basin in the southwestern watersheds, and in west central and southwest watersheds of the Lake Michigan basin. We then used the models to estimate water quality in watersheds without observed instream data to prioritize those areas for management. Prioritizing watersheds will aid effective management of the Great Lakes watershed and result in efficient use of restoration funds, which will lead to improved nearshore water quality.     

Status of the amphipod Diporeia spp. in Lake Superior, 2006–2016

The amphipod Diporeia spp. has historically been an important component of the benthic food web of the Laurentian Great Lakes. The Great Lakes Water Quality Agreement included its population density as an indicator of ecological condition for Lake Superior, with target values of 220–320 m^?2 in nearshore areas (≤100 m depth) and 30–160 m^?2 in offshore areas (>100 m). To assess the status of Diporeia in Lake Superior, we used a probability-based lake-wide survey design to obtain estimates of Diporeia density and biomass in 2006, 2011 and 2016. A PONAR grab sampler was used to collect Diporeia at 50–53 sites each year, with approximately half in the nearshore (<100 m depth) region of the lake and half in the offshore. The mean area-weighted lake-wide density was 395 ± 56 (SE) m^?2 in 2006, 756 ± 129 m^?2 in 2011, and 502 ± 60 m^?2 in 2016. For all years, both density and biomass were greater in the nearshore than in the offshore stratum. The densities for 2006–2016 were 3–5 times higher than those reported from a lake-wide survey conducted in 1973 by the Canada Centre for Inland Waters. The severe declines in Diporeia populations observed in the other Great Lakes during recent decades have apparently not occurred in Lake Superior. Further research is needed to understand spatial and temporal variability of Diporeia populations in Lake Superior to enhance the utility of Diporeia density as an indicator of benthic condition.     

Temporal and spatial differences in deposition of organic matter and black carbon in Lake Michigan sediments over the period 1850–2010

Inorganic carbon (IC), total organic carbon (OC), and black carbon (BC) were analyzed in eight sediment cores obtained from deep water (>30 m) sediments in the Chippewa and south Chippewa basins, as well as Green Bay in Lake Michigan. These cores were segmented at high resolution and radio-dated to reconstruct a detailed history of deposition to the lake both spatially and temporally since ca. 1850 CE. To help interpret the depositional record, cores were also characterized for stable isotopes (¹³C and ¹⁵N), as well as particle size distribution, density, organic matter (OM), and other parameters. Fine (silt and clay) sediment particles contained OM of primarily lacustrine algal biomass origin. Sedimentation fluxes showed large increases in OM and OC fluxes through much of the lake during the onset of industrialization and the period of rapid industrialization to onset of Great Lakes environmental legislation. In contrast, fluxes and loading of BC increased dramatically in the southern basin until the 1930's, then decreased substantially after the 1940's. This observation was due largely to results from site M009 nearest the steel mills and industrial zones of Chicago and northern Indiana. Together, whole lake loadings of OM and BC provide evidence that changing industrial activity and legislation intended to curb air pollution in the Great Lakes region have had a fairly rapid and dramatic impact. In contrast, legislation intended to decrease eutrophication through reductions in nutrient loading to the lake have not had a similar impact on sedimentation of OM in the lake.     

Pink Salmon Populations in the U.S. Waters of Lake Superior, 1981–1984

Pink salmon were introduced to Lake Superior in very small numbers at Thunder Bay, Ontario, Canada in 1956. Since that time, they have established themselves as permanent members of the Lake Superior ecosystem. Peak spawner densities were observed in the U.S. tributaries of Lake Superior in 1979 when many streams experienced runs of 10,000 fish or more. To better understand the population dynamics of pink salmon, population characteristics, fecundity and egg deposition of spawners, egg and larval survival, and fry outmigration were studied in several Michigan and Minnesota tributaries in 1980. Population estimates indicate a significant decline of pink salmon in the U. S. waters of Lake Superior since 1979. This decline appears to be related to instream fry survival, which is in turn related to the hydrological conditions of the streams in fall and winter. Lake Superior pink salmon are phenotypically similar to their Pacific Coast counterparts, the major difference being age structure and smaller size of adults. Approximately 90% of the adults mature as 2-year-olds while the remaining 10% mature as 3-year-olds. Pink salmon maturing as 3-year-olds have different growth patterns than those maturing at 2 years of age. Two-year-old pink salmon average 390 mm in length, which is about 30% smaller than 2-year-olds from the Pacific Ocean. Three-year-old females have a lower fecundity and a poorer egg quality than 2-year-old females.     

Physical–chemical measurements in the water column along a transect through a tilapia cage fish farm in Lake Malawi,Africa

The initiation of cage aquaculture in the shallow southeast arm of Lake Malawi has raised concerns about its possible impact on the surrounding environment and the highly diverse fish community. To evaluate the impact of the cage operation on the surrounding environment, observations were made over an annual cycle in 2007 at a production capacity of ~ 200 tonnes fish/year. Impacts of the cage wastes in the water column in the vicinity of the cages were minimal despite the substantial discharges from the cages. No significant differences were observed in concentrations of dissolved and particulate nutrients (ammonia, nitrate, phosphate, particulate C,N and P), chlorophyll, chlorophyll fluorescence, dissolved oxygen, total suspended solids, Secchi depths and extinction coefficient of photosynthetically active radiation between the cage site and the control stations upstream or downstream of the fish farm. Although sedimentation rates measured in traps were higher under the cages than at control sites, the sediment flux was a small percentage of total feeds added to the cages. Apparently cage wastes were efficiently dispersed by water currents which averaged 9.3 cm/s below the fish cages. Consumption of the wastes by wild fish species which aggregated around the fish cages and their movement in the vicinity of the cages also contributed to the dispersion of the cage wastes and served to dilute impact of the cages. In combination, these physical and biological processes reduced the immediate impact of the cage farming operation and must be considered in the siting of future cage farms.     

Water Quality Changes in Lake Erie, 1968–1980

Data on chlorophyll a, total phosphorus, and hypolimnetic dissolved oxygen concentration, which were collected by Canada Centre for Inland Waters during the period 1968 to 1980, are statistically analyzed to evaluate the changes in the water quality of Lake Erie. There are strong evidences of a decreasing trend in the value of chlorophyll a and total phosphorus in the western, central, and eastern basins of the lake between 1970 and 1980. Furthermore, a statistical model is developed for the hypolimnetic dissolved oxygen concentration in the central basin. The model shows that the increase in depletion is related to the increase in the level of total phosphorus. Hence, it is possible to improve the anoxic conditions in the lake by controlling total phosphorus loadings.     

City storm-flood events in China, 1984–2015

This study uses web-based information to explore the spatial and temporal trends of reported city flood events for all mainland China cities from 1984 to 2015. Panel data were compiled on flooding (or the lack thereof) within China’s cities for every year in this period, and the relationship between the annual number of cities with reported flood events and possible influencing variables was analyzed. Few cities experienced flood events before 2001, but this situation then increased dramatically after 2010. In 2015, approximately 60% of China’s cities experienced a storm-flood event, which is three times as many as the government estimated in 2011.     

Ice–water heat exchange during ice growth in Lake Baikal

Using a custom-made thermistor chain frozen into the ice cover we obtained the first detailed information on distribution of temperature within ice and structure of the ice–water boundary layer during ice growth in Lake Baikal. A mathematical model of the heat transport in a multilayer ice–water system (Stefan problem) was developed and verified on results of in situ measurements. Effective coefficients of thermal diffusivity and ice–water heat fluxes were estimated from the inverse solution of the model and compared with direct flux estimates from the flux-gradient method. Both estimations agreed on flux values of 1–10 W m^− 2 and demonstrated strong synoptic variability in ice–water heat exchange. We estimated the thickness of viscous laminar sublayer under ice, as well as the thickness of the transitional layer on top of the turbulent water column. The thickness of the viscous sublayer of 1–1.5 cm in Lake Baikal was several times smaller than values reported previously from small lakes, suggesting high magnitudes of convective velocities and/or of the under-ice currents in Lake Baikal. Significant growth of the thermal diffusivity coefficient with increasing distance from the ice bottom was detected: its value at the top of the transition layer of under ice water was 10–40 times higher compared with its value in viscous laminar sublayer. This is also significantly higher than previous estimations in smaller freshwater lakes.     

Studies of Large-Scale Currents in Lake Erie, 1979–80

Currents and water temperatures were recorded at a large-scale grid of fixed moorings in Lake Erie from May 1979 through June 1980. Currents measured in the lower half of the central basin water column were mostly return flows (beneath the surface wind drift) driven by the surface pressure gradient. Often observed was a complex system of Lake Erie circulation gyres as predicted by models. Another common occurrence was for one of the central basin gyres to become dominant and envelop the whole basin in either uniform clockwise or anticlockwise flow. It is not fully certain why one of the circulation cells grows as opposed to the others, but the curl of the wind stress had influence. The currents were more barotropic than predicted by full Ekman layer current models. Tidal-like currents driven by the longitudinal seiches of Lake Erie dominate the island-filled passages between the western and central Lake Erie basins, with currents across the whole island chain very closely in phase. Processes of hypolimnion volume entrainment are suggested from the central basin temperature recordings. Large volume water exchanges between the central and eastern basins occurred after the water mass in the vicinity of the shallow ridge that separates them had become unstratified. These and other topics are discussed as the large data set generated from the experiment is explored.     

Spatiotemporal trends of polychlorinated biphenyls (PCBs) in surface and suspended sediments from the Lake Ontario Canadian nearshore 1994–2018: A fish consumption advisory perspective

Spatiotemporal trends for polychlorinated biphenyls (PCBs) were examined in surface and suspended sediments collected between 1994 and 2018 from over twenty nearshore stations on the Canadian side of Lake Ontario and the St. Lawrence River. In 2018, PCB concentrations ranged over an order-of-magnitude in surface (<10 ng/g?357 ng/g) and suspended sediments (<10 ng/g–330 ng/g), illustrating the presence of legacy hotspots as well as diffuse urban inputs. PCB concentrations in both surface and suspended sediments were consistently highest in Hamilton Harbour, but were also elevated at other stations around the perimeter of the Niagara basin as well as near Trenton and Kingston, Ontario. Generally, higher PCB concentrations were found in surface sediment relative to paired suspended sediment samples suggesting temporal improvements in nearshore sediment quality. However, many stations demonstrated temporal variability in PCB concentrations, a likely factor in the lack of an overall nearshore declining trend. Given that PCBs are listed as a consumption-limiting contaminant for all fish sampling blocks in the Canadian waters of Lake Ontario and are responsible for 78% of restricted advisories, sediment quality benchmarks that account for bioaccumulation potential should be considered over toxicologically-derived guidelines to gauge severity of PCB contamination of nearshore sediments. Relatively higher TOC-normalized PCB concentrations in the western end of Lake Ontario suggests additional research on PCB bioavailability from nearshore sediment would be beneficial in understanding basinwide trends of PCBs in fish, and that an adaptive approach to sediment management may be needed in the context of consumption advisories.     

Reproductive Potential and Fecundity of Lake Trout Strains in Southern and Eastern Waters of Lake Ontario, 1977–1994

We assessed the reproductive potential of various genetic strains of hatchery lake trout (Salvelinus namaycush) in southern and eastern Lake Ontario from indices of fecundity and indices of male abundance. Indices were constructed from catches of mature lake trout in gill nets during September 1980 to 1994 after correcting for mortality from sea lampreys (Petromyzon marinus) which occurred between September sampling and late fall spawning. Strain and age were assigned to individual lake trout based on clipped fins and maxillary bones or coded wire tags. Fecundity-length relationships for fish of the same age, determined from mature females collected in 1977 to 1981 and 1994, were not different (P > 0.05) among genetic strains. For all strains combined, fecundity-length relationships in 1977 to 1981 were not different among fish of various ages but in 1994, age-5 and -6 fish had fewer eggs (P < 0.003) than age-7 fish, and age-7 fish had fewer eggs (P < 0.003) than fish of age 8, 9, or 10. Annual indices of fecundity varied 19 fold and indices of mature males varied 11 fold; both indices were low in the early 1980s, increased sharply in the mid 1980s, and peaked in 1993. The strain which dominated fecundity and mature male indices shifted during the study from Seneca Lake strain to Lake Superior strain and then back to Seneca Lake strain. However, changes in either reproductive potential or genotypes do not appear responsible for the abrupt appearance of naturally-produced yearling lake trout throughout southern and eastern Lake Ontario in 1994–1995, the first widespread occurrence of juveniles produced by hatchery lake trout in Lake Ontario.