An exploration of Lake Ontario’s offshore epibenthic cladoceran community

Pelagic cladocerans are relatively well studied within the Great Lakes. However, little is known about cladocerans that more closely interact with benthic substrate. In 2018, 26 sites in Lake Ontario were sampled with benthic grabs to collect meiofauna. The stations ranged in depth from 0.1 to 184 m with sand, silt or hard bottom substrate. Epibenthic cladocerans were recovered from all 12 sites shallower than 20 m, five of eight sites between 20 and 40 m and one (63 m deep) of six sites deeper than 40 m. The density of epibenthic cladocerans were at most 12% of the density of harpacticoid copepods (shallow hard bottom sites). The Lake Ontario epibenthic cladoceran community was represented by 16 species from 8 different genera although one of these species, Illyocryptus cf. sordidus, is likely represented by two cryptic species. The most widespread species was identified as Alona sibirica (Sinev, 2020), a newly described species from Russia formerly of the Alona cf. affinis (Leydig, 1860) species complex. The Great Lakes specimens were 98.3% similar genetically to sequences assigned to A. cf. sibirica. Among substrates, sand had the highest densities (1662 m² in 0–20 m) and species richness (10) of epibenthic Cladocera. After cladocerans were identified, the cytochrome oxidase I (COI) gene of 37 specimens were successfully sequenced through the Barcode of Life (BOLD). Six specimens were assigned to three previously existing Barcode Identification Numbers (BINs), the remaining 31 specimens were assigned to eight new BINs.     

Dreissena in Lake Ontario 30 years post-invasion

We examined three decades of changes in dreissenid populations in Lake Ontario and predation by round goby (Neogobius melanostomus). Dreissenids (almost exclusively quagga mussels, Dreissena rostriformis bugensis) peaked in 2003, 13 years after arrival, and then declined at depths <90 m but continued to increase deeper through 2018. Lake-wide density also increased from 2008 to 2018 along with average mussel lengths and lake-wide biomass, which reached an all-time high in 2018 (25.2 ± 3.3 g AFTDW/m²). Round goby densities were estimated at 4.2 fish/m² using videography at 10 to 35 m depth range in 2018. This density should impact mussel populations based on feeding rates, as indicated in the literature. While the abundance of 0–5 mm mussels appears to be high in all three years with measured length distributions (2008, 2013, 2018), the abundance of 5 to 12 mm dreissenids, the size range most commonly consumed by round goby, was low except at >90 m depths. Although the size distributions indicate that round goby is affecting mussel recruitment, we did not find a decline in dreissenid density in the nearshore and mid-depth ranges where goby have been abundant since 2005. The lake-wide densities and biomass of quagga mussels have increased over time, due to both the growth of individual mussels in the shallower depths, and a continuing increase in density at >90 m. Thus, the ecological effects of quagga mussels in Lake Ontario are likely to continue into the foreseeable future.     

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.     

Expansion of Dreissena into offshore waters of Lake Michigan and potential impacts on fish populations

Lake Michigan was invaded by zebra mussels (Dreissena polymorpha) in the late 1980s and then followed by quagga mussels (D. bugensis) around 1997. Through 2000, both species (herein Dreissena) were largely restricted to depths less than 50 m. Herein, we provide results of an annual lake-wide bottom trawl survey in Lake Michigan that reveal the relative biomass and depth distribution of Dreissena between 1999 and 2007 (although biomass estimates from a bottom trawl are biased low). Lake-wide mean biomass density (g/m²) and mean depth of collection revealed no trend between 1999 and 2003 (mean = 0.7 g/m² and 37 m, respectively). Between 2004 and 2007, however, mean lake-wide biomass density increased from 0.8 g/m² to 7.0 g/m², because of increased density at depths between 30 and 110 m, and mean depth of collection increased from 42 to 77 m. This pattern was confirmed by a generalized additive model. Coincident with the Dreissena expansion that occurred beginning in 2004, fish biomass density (generally planktivores) declined 71% between 2003 and 2007. Current understanding of fish population dynamics, however, indicates that Dreissena expansion is not the primary explanation for the decline of fish, and we provide a species-specific account for more likely underlying factors. Nonetheless, future sampling and research may reveal a better understanding of the potential negative interactions between Dreissena and fish in Lake Michigan and elsewhere.     

Lake-wide,annual status of the Mysis diluviana population in Lake Michigan in 2015

Mysis diluviana is one of the most abundant zooplankton by biomass in the Laurentian Great Lakes of North America, a predator of other zooplankton and an important prey for fishes. Studies of long-term trends in Lake Michigan have shown 2005–2016 densities to be 50–80% lower than 1990s densities, but these observations have been based on annual monitoring that is either spatially or seasonally limited. We combined Lake Michigan Mysis data from three annual programs and the 2015 Cooperative Science and Monitoring Initiative to achieve broad spatial coverage during spring, summer, and fall of 2015 and broad depth coverage during spring 2016. Lake-wide, annual density and biomass were 82 (SE: 10) Mysis/m² and 200 (SE: 36) mg dry mass/m². Density and biomass estimates were highest offshore, generally higher in the north basin, and seasonally highest in summer. Annual lake-wide averages for depths >30 m were better captured by seasonally-extensive annual programs than spatially-extensive annual programs, although spring sampling may bias annual values low. Mysis cohorts grew 0.026 mm/d (age-0) and 0.007 to 0.027 mm/d (age-1). Annual mortality was 81–98%. Reproduction was fall-spring and seasonal lake-wide estimates ranged from 0.6 to 19.1% females brooding, 13–20 embryos/brood, and 3–46 embryos/m². Annual production (423 mg dry mass/m²/yr, SE: 31) was lower than all but one previous estimate from lakes Michigan, Huron, and Ontario. While Mysis tend to persist, low Mysis production may be a concern for prey fishes that feed on Mysis.     

Lake Superior zooplankton biomass: Alternate estimates from a probability-based net survey and spatially extensive LOPC surveys

We conducted a probability-based net tow sampling of Lake Superior in 2006 and compared the zooplankton biomass estimate with an estimate from laser optical plankton counter (LOPC) surveys. The net survey consisted of 52 sites stratified across three depth zones (0–30, 30–150, > 150 m). The LOPC tow surveys were extensive and spatially covered much of Lake Superior (> 1300 km of towing). The LOPC was field calibrated to Lake Superior zooplankton samples collected across the years of 2004 to 2006. The volume-weighted lake-wide zooplankton biomass determined by traditional net tows to 100-m sample depth was 20.1 (± 7.8 SD n = 52) mg dry-weight m^− 3. The estimates varied by depth zones within the lake, where nearshore (0–30 m) estimates were highest and highly variable. Net sites for the LOPC field calibration were removed to allow for LOPC validation with independent nets; the resulting net-based estimate 20.0 (± 9.3 SD n = 38) mg dry-weight m^− 3 and LOPC lake-wide estimate 19.1 (± 3.3 SD) mg dry-weight m^− 3 agreed well. Consistency across survey methods for lake-wide estimates suggested that LOPC survey data provides a comparable assessment tool to traditional nets for collecting zooplankton biomass data. We briefly compare our results with some observed historical patterns. Onshore–offshore trends in zooplankton biomass concentrations were similar to the last major lake-wide survey in 1973. The LOPC provided high resolution data on zooplankton biomass distribution. Using simultaneously collected in situ sensor data, the LOPC zooplankton biomass distributions over horizontal and vertical space can be modeled as a function of temperature and fluorescence.     

Incorporation of non-native species in the diets of cisco (Coregonus artedi) from eastern Lake Ontario

Cisco Coregonus artedi was once an important native fish in Lake Ontario; however, after multiple population crashes, the cisco stock has yet to recover to historic abundances. Rehabilitation of cisco in Lake Ontario is a fish community management objective, but the extent to which recent non-native species and pelagic food web changes have influenced cisco is not well understood. We described cisco diets in contemporary Lake Ontario following the addition and spread of non-native zooplankton species. We collected 618 cisco and processed 178 for full diet analysis in eastern Lake Ontario using mid-water trawls and bottom-set gill nets from 2016 to 2020. We found that Lake Ontario cisco were mostly zooplanktivorous, and non-native zooplankton dominated their diet during July and September. Cisco smaller than 300 mm had a more diverse diet including both native and non-native zooplankton, while cisco larger than 300 mm fed almost exclusively on non-native predatory cladocerans Bythotrephes longimanus and Cercopagis pengoi (98.9% consumed prey dry mass). We also found fish eggs, presumed to be of coregonine origin in 75% of non-empty December-collected cisco diets, suggesting eggs subsidize cisco diets when available. Juvenile round goby Neogobius melanostomus, alewife Alosa pseudoharengus and rainbow smelt Osmerus mordax were found in 2% of all analyzed non-empty stomachs. Lake Ontario cisco diet appears to be more similar to zooplanktivorous Lake Superior cisco than Lake Michigan where piscivory is prevalent. Lake Ontario cisco diets reflected zooplankton community changes indicating that non-native predatory cladocerans are now an important energy source supporting this native species.     

In situ-measured primary production in Lake Superior

Water column primary production is a major term in the organic carbon cycle, particularly in large lakes with relatively reduced shoreline and near-shore influence. Presently, there is a large imbalance in the known inputs vs. outputs of organic carbon in Lake Superior. This study examined primary production in offshore Lake Superior using in situ incubations over a range of conditions representing an annual cycle. Primary producers were dominated by small (< 20 μm) cells and included a relatively large abundance of small, spherical flagellates. During conditions with a warm surface layer, chlorophyll concentrations were two- to three-fold higher within the deep chlorophyll maximum (DCM) than at the surface. Volumetric production (mass L^− 1 d^− 1) was maximal at 2-10 m depth, well above the typical DCM depth. On average, 22% of ¹⁴C label appeared in the dissolved pool at the end of the incubation period with the rest appearing in GF/F-strained particles. A statistical model for volumetric production explained 93% of the variance in individual measurements for depths > 2 m, using temperature and light as predictors. This model was applied to annual fields of temperature and light, and a new estimate for whole-lake annual primary production, 9.73 Tg y^− 1, was derived. This combination of new measurements and modeling results brings the organic carbon cycle of Lake Superior closer to being balanced.     

Measurement of Diporeia respiration rate for Lake Superior

The freshwater amphipod Diporeia is a dominant macroinvertebrate species in Lake Superior’s benthic community and an important prey item for many fish. A capacity to predict growth and production rates of Diporeia using a bioenergetics model requires information on physiological processes of the species. The objective of this study is to quantify oxygen consumption of Lake Superior Diporeia and to determine if respiration rate changes with body length. Diporeia were collected from Lake Superior and kept over natural sediment maintained at 4 °C. Dissolved oxygen levels for groups of immature (2 mm), juvenile (4 mm), and adult (6 mm) Diporeia in 20 ml microcosms were measured using a polarographic microelectrode. Mass-specific respiration rates for Lake Superior Diporeia ranged from 32.0 to 44.7 mg O₂ g DW^− 1 day ⁻¹. A significant relationship between body length and mass-specific respiration rate (p > 0.1) was not found. The estimate of Diporeia respiration presented here is significantly higher (p < 0.05) than previous findings from populations in Lakes Michigan and Ontario. This study provides new data on respiration rates of Lake Superior Diporeia and compares findings to studies for other connecting Great Lakes.     

Hydrology influences generalist–specialist bird-based indices of biotic integrity in Great Lakes coastal wetlands

Marsh bird habitats are influenced by water levels which may pose challenges for interpreting bird-based indices of wetland health. We determined how much fluctuating water levels and associated changes in emergent vegetation influence the Index of Marsh Bird Community Integrity (IMBCI) using data collected in Great Lakes coastal wetlands by participants in Bird Studies Canada's Great Lakes Marsh Monitoring Program. IMBCI scores for 90 wetlands in Lake Erie and 131 wetlands in Lake Ontario decreased with decreasing water levels due to decreasing number of marsh-dependent species in Lake Erie and perhaps also in Lake Ontario. The average magnitude of the decrease in scores between extremely high and low water periods for wetlands with sufficient data was 15% in Lake Erie where water dropped 0.9 m on average (n = 11 wetlands) and 18% in Lake Ontario where water dropped 0.5 m (n = 7). Scores in Lake Erie increased with increasing Typha due to increasing numbers of marsh-dependent species and decreased with increasing Phragmites due to increasing numbers of generalist species. The opposite was observed in Lake Ontario, perhaps due to denser Typha and sparser Phragmites. Scores were explained by the naturally fluctuating water levels of Lake Erie, which favored Phragmites expansion and the regulated water levels of Lake Ontario which promoted Typha expansion. Scores were influenced by fluctuating water levels and associated changes in emergent vegetation. Inter-annual water level fluctuations should be considered when interpreting any indicator of wetland health that is based on marsh-dependent bird species.     

Modelling of nearshore microbial water quality at confluence of a local tributary in Lake St. Clair

Microbial water quality, measured as Escherichia coli (E. coli) concentration, at beaches along the southern shore of Lake St. Clair in Canada, often exceeds public safety guidelines. Belle River, located near a public beach and a drinking water intake, is one of the several smaller tributaries of the lake whose contribution to nearshore microbial water quality is currently unknown. A flexible mesh 3D coupled TUFLOW-FV and Aquatic Ecodynamic (AED2+) model was used to simulate the hydrodynamics and microbial water quality in Lake St. Clair. A higher resolution nested model was developed within the lake-wide TUFLOW-FV model for better spatial and temporal resolution in the local region surrounding Belle River. Regular and up to a factor of four difference in predicted E. coli concentrations were observed with the nested and lake-wide models at the public beach next to Belle River, whereas the difference was marginal at the drinking water intake about a kilometre away from the shore. While the E. coli loading to Lake St. Clair from Belle River is considered negligible, >90% of the predicted daily E. coli concentration at the beach and > 50 % at the water intake were attributed to Belle River from amongst all watershed sources to Lake St. Clair considered in the model. The model results also show that the construction of a new 150 m jetty in 2018, replacing the older 25 m jetty separating Belle River from the public beach, is expected to significantly reduce E. coli concentrations observed at the beach.     

Spatial and temporal variability in lake trout diets in Lake Ontario as revealed by stomach contents and stable isotopes

Lake trout (Salvelinus namaycush) are an ecologically and economically important piscivore with reported differences in diet and feeding behaviour throughout its range. Eleven stomach content and stable isotope-based metrics were used to describe diets of 349 lake trout between two years (2013 and 2018) and among geographic zones (west, central, east, Kingston basin) in Lake Ontario. Using individual (e.g., volumetric, %V) and aggregate (e.g., index of relative importance, %IRI) diet metrics, we found an overwhelming dominance of alewife (Alosa pseudoharengus) in lake trout diets among some zones in 2013 (%V = 23.3 – 92.7; %IRI = 12.2 – 99.5) and all zones in 2018 (%V = 83.9 – 96.7; %IRI = 96.5 – 100). Round goby (Neogobius melanostomus) and rainbow smelt (Osmerus mordax) were secondary lake trout prey items with relative diet percentages only marginally reflected by spatial and temporal variation in prey abundance (round goby: %V = 1.0 – 33.3, %IRI = 0.1 – 13.2; rainbow smelt: %V = 2.5 – 54.0, %IRI = 0.1 – 54.0). Carbon (δ¹³C) and nitrogen (δ¹⁵N) isotopic niche areas and orientations were similar across all year-zone combinations reinforcing temporal and spatial consistency in lake trout diet. The findings of this study advance the time series in describing Lake Ontario lake trout diets and can be used to complement stock assessments and management decisions associated with carrying capacity for the diverse salmonid community.     

Patterns in spatial use of land-locked Atlantic salmon (Salmo salar) in a large lake

Understanding the spatial use of reintroduced fish is useful for fisheries management and evaluating restoration success. Atlantic salmon (Salmo salar) were reintroduced into Lake Ontario in the 1990s; however, the movement ecology of these land-locked fish is unknown. Using acoustic telemetry and Floy tag mark-recaptures, we examined seasonal home range and space use of Atlantic salmon in Lake Ontario. Hatchery-raised adult Atlantic salmon were tagged with acoustic transmitters (n = 14; 8 with depth sensors) or Floy tags (n = 1915) and released. Both acoustic telemetry and Floy tag recaptures (n = 90) indicated cross lake movements, and home ranges encompassed nearly the entire lake in summer but was smaller in winter. Movements were nearshore (<2 km from shore) from spring to summer at ～20 m bathymetric depths, with movements closer to shore in the fall, and further offshore (～5.5 km from shore and 45 m bathymetric depths) in winter. Depth use was relatively shallow (<4 m) with occasional deeper dives (max = 28.5 m), and small diel vertical movements (1–5 m), moving deeper during daytime, consistent with ocean movements of Atlantic salmon. There appears to be spatial segregation among Atlantic salmon and other Lake Ontario salmonids, however, overlap likely occurs in nearshore waters during the spring. Wide-ranging movements of Atlantic salmon in binational (Canada/USA) waters reflects the importance of government agencies collaborating to ensure sustainable fisheries and the coordination of species restoration activities. This is the first study to provide detailed spatial use of Lake Ontario Atlantic salmon to assist in the management of this reintroduced species.     

Six decades of Lake Ontario ecological history according to benthos

The Laurentian Great Lakes have experienced multiple anthropogenic changes in the past century, including cultural eutrophication, phosphorus abatement initiatives, and the introduction of invasive species. Lake Ontario, the most downstream lake in the system, is considered to be among the most impaired. The benthos of Lake Ontario has been studied intensively in the last six decades and can provide insights into the impact of environmental changes over time. We used multivariate community analyses to examine temporal changes in community composition over the last 54 years and to assess the major drivers of long-term changes in benthos. The benthic community of Lake Ontario underwent significant transformations that correspond with three major periods. The first period, termed the pre/early Dreissena period (1964–1990), was characterized by high densities of Diporeia, Sphaeriidae, and Tubificidae. During the next period defined by zebra mussel dominance (the 1990s) the same groups were still prevalent, but at altered densities. In the most recent period (2000s to present), which is characterized by the dominance and proliferation of quagga mussels deeper into the lake, the community has changed dramatically: Diporeia almost completely disappeared, Sphaeriidae have greatly declined, and densities of quagga mussels, Oligochaeta and Chironomidae have increased. The introduction of invasive dreissenids has changed the Lake Ontario benthic community, historically dominated by Diporeia, Oligochaeta and Sphaeriidae, to a community dominated by quagga mussels and Oligochaeta. Dreissenids, especially the quagga mussel, were the major drivers of these changes over the last half century.     

European Valve Snail Valvata piscinalis (Müller) in the Laurentian Great Lakes Basin

Previously reported from the lower Great Lakes basin and St. Lawrence and Hudson rivers, the nonindigenous gastropod Valvata piscinalis was found for the first time in Superior Bay (Minnesota) of Lake Superior, Lake Michigan (Wisconsin), and Oneida Lake (New York) of the Lake Ontario basin. This snail was not abundant in Lakes Superior and Michigan, whereas in eutrophic Oneida Lake it reached a maximum density of 1,690 individuals/m² (mean density = 216 individuals/m²). Human-mediated disturbances could facilitate the range extension of this snail by providing dispersal opportunities (e.g., canals, shipping traffic) or increasing nutrients (e.g., eutrophication). A native of the Palaearctic region, V. piscinalis has colonized sites across the Great Lakes basin, suggesting that it will likely become common in disturbed Great Lakes areas.     

Status of Mysis diluviana in Lake Ontario in 2013: Lower abundance but higher fecundity than in the 1990s

Mysis diluviana is a major component of prey fish diets in the Great Lakes, so annual production of M. diluviana is important for understanding and modeling energy flow through Great Lakes food webs. However, only three lake-wide measurements of M. diluviana annual production in Lake Ontario are currently available (1971, 1990, 1995). During 2013, lake-wide coverage of Lake Ontario was achieved during four periods from April to November. Annual mean density and biomass of M. diluviana in 2013 were 99?#/m² (SE: 8) and 318?mg?dw/m² (SE: 28) – approximately half of values observed in 1990s. M. diluviana comprised 13–30% of offshore zooplankton biomass in each period. Reproduction peaked in fall, with mean brood size of 32 embryos (range: 11–49), at least 10% larger than in 1990s. Generation time was two years from embryo to initial reproduction. Growth rates were 0.052?mm/d for the age-0 cohort and 0.027?mm/d for the age-1 cohort. Age-0 growth rate was significantly higher than in 1980s–90s (0.035?mm/d). Annual production in 2013 was 0.85?g?dw/m²/yr (SE: 0.03) which was 30–40% of values observed in 1990 and 1995 (2.23 and 2.53?g/m²/yr). Annual production to biomass ratio (P/B) in 2013 was 2.65?/yr which was 80–85% of values observed in 1990 and 1995 (3.24 and 3.11?/yr), but this difference was not statistically significant. Our results suggest that changes in annual production over time can be estimated using changes in biomass over time and a mean P/B ratio.     

The Invasive New Zealand Mud Snail (Potamopyrgus antipodarum) in Lake Erie

The New Zealand mud snail (Potamopyrgus antipodarum) is an invasive species in Europe, Japan, Australia, and North America. In the western United States it is a species of special concern where population densities in some rivers and streams are very large (∼300,000 per m²) and considerable ecological effects of its presence have been reported. Much less about the effects of this species is known in the Great Lakes, where the snail was found in Lake Ontario and the St. Lawrence River in 1991. Here we report the occurrence of the snail in Lake Erie. Two P. antipodarum were collected in 18 m deep water (sampling range 5–18 m) in Lake Erie off shore of Presque Isle State Park near Erie, Pennsylvania in the summer of 2005 and others were collected off of Sturgeon Point in Lake Erie (sampling range 5–20 m) south of Buffalo, NY and in the central basin of Lake Erie (18 m) in 2006. This finding demonstrates that this species continues to expand its range in the Great Lakes. The range expansion increases the likelihood that it may become established in rivers and streams emptying into the Great Lakes where higher densities and greater ecological damage may result.     

Changes in mid-summer water temperature and clarity across the Great Lakes between 1968 and 2002

In recent decades, three important events have likely played a role in changing the water temperature and clarity of the Laurentian Great Lakes: 1) warmer climate, 2) reduced phosphorus loading, and 3) invasion by European Dreissenid mussels. This paper compiled environmental data from government agencies monitoring the middle and lower portions of the Great Lakes basin (lakes Huron, Erie and Ontario) to document changes in aquatic environments between 1968 and 2002. Over this 34-year period, mean annual air temperature increased at an average rate of 0.037 °C/y, resulting in a 1.3 °C increase. Surface water temperature during August has been rising at annual rates of 0.084 °C (Lake Huron) and 0.048 °C (Lake Ontario) resulting in increases of 2.9 °C and 1.6 °C, respectively. In Lake Erie, the trend was also positive, but it was smaller and not significant. Water clarity, measured here by August Secchi depth, increased in all lakes. Secchi depth increased 1.7 m in Lake Huron, 3.1 m in Lake Ontario and 2.4 m in Lake Erie. Prior to the invasion of Dreissenid mussels, increases in Secchi depth were significant (p < 0.05) in lakes Erie and Ontario, suggesting that phosphorus abatement aided water clarity. After Dreissenid mussel invasion, significant increases in Secchi depth were detected in lakes Ontario and Huron.     

Alternative Ecological Pathways in the Eastern Lake Ontario Food Web—Round Goby in the Diet of Lake Trout

Round goby (Neogobius melanostomus) range expansion and their possible inclusion in the diet of lake trout (Salvelinus namaycush) were investigated. Fish community index bottom trawls in eastern Lake Ontario (Kingston basin) during summer 2003 and 2004 indicated the presence of the round goby at relatively low densities (3.72 × 10^–2 ± 5.24 × 10^–3 fish/m²) in depths up to 30 m. Lake trout (mean fork length = 585 ± 78 mm and mean weight = 2,770 ± 1,134 g) stomach contents showed round goby to be the second most abundant diet item at almost 20% by number (36% by mass). Round goby ingested by lake trout ranged in total length from 50 to 110 mm. The most important prey species in terms of abundance (68%) and mass (56%) was alewife (Alosa pseudoharengus). Alewives were the most important diet item for all sizes of lake trout sampled, except those in the 550 - 650 mm size class, which ingested more round goby by mass than alewife. Round goby range expansion to deep water and prominence in the diet of lake trout signal significant change in the eastern Lake Ontario food web.