Geology topics

Research about Lake Superior

Source-linked reports with geographic coverage including Lake Superior.

At least 19 recordsLinked to original sources

USGS2026

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2026

Status and trends in the Lake Superior fish community, 2025

The U.S. Geological Survey has conducted annual trawl surveys across Lake Superior since 1978 that describe trends in fish species occurrence and relative abundance to support fisheries science and management. In 2025, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 72 nearshore stations in May and June and 36 offshore locations in July. Nearshore bottom trawls collected 63,005 fish represented by 30 species or morphotypes. The number of species collected at each location ranged from 1 to 13, with a median of 7.0 species. Estimated fish biomass density at individual stations ranged from <0.1 to 125.7 kg per ha with a lakewide mean of 7.6 kg per ha. Offshore bottom trawls collected 30,342 fish represented by 13 species or morphotypes. Estimated fish biomass density at individual stations ranged from 1.6 to 34.7 kg per ha with a lakewide mean of 9.1 kg per ha, which was the highest since the offshore survey began in 2011. Lakewide average numerical densities (fish per ha) of age-1 fish were 0.67 per ha for Bloater, 0.01 per ha for Cisco, 1.15 per ha for Kiyi, 0.95 per ha for Lake Whitefish, and 336.33 per ha for Rainbow Smelt. Surface trawling collected 1,562 larval Coregonus individuals which was the fewest Coregonus larvae collected in a whole lake survey since the larval fish survey began in 2014. Nearshore mean larval Coregonus numerical densities were 156 fish per ha in May and June 2025 and 21 fish per ha in July 2025. May mean surface water temperatures (6.1°C) were near the warmest for the period-of-record, while June (6.8°C) and July (10.0°C) were near average or below.

Lake Superior↗

USGS2025

Surface ice effects on water levels, wave dynamics and wave runup

Surface ice-induced changes in hydrodynamics may have significant implications on coastal hazards such as erosion and flooding. The objective of this study was to investigate how Lake Superior's water levels (WLs) and wave dynamics, that is, hydrodynamics, as well as the wave power and wave runup are influenced by the formation of surface ice. This study utilized a coupled circulation and spectral wave model, ADvanced CIRCulation Model (ADCIRC) + Simulating WAves Nearshore (SWAN). In this study, we considered six different ice-on, incorporating either ice thickness or concentration, and one hypothetical ice-free, ignoring any form of ice forcing, scenarios. The impact of surface ice on lake's WLs was found to be fairly insignificant. In contrast, surface ice-induced modifications in wave spectra resulted in smaller significant wave heights but longer peak wave periods. These observations were attributed to less energetic and narrower wave spectra of ice-on models compared to that of the ice-free model. The power loss caused by surface ice in a month was found to have enough capacity to meet the electrical demands of 55–100 households for an entire year in the United States. According to wave runup comparisons of ice-on and ice-free models at different regions, the effect of surface ice on wave runup could potentially exacerbate the adverse effect of coastal flooding by increasing the wave runup more than 20 cm. Such a geographic variation in the wave runup was observed to be influenced by the distance traveled by the waves within the ice field and duration of ice-wave interaction.

Lake Superior↗

USGS2025

Lake Superior fish community and fisheries, 2001–2022: An era of stability

Lake Superior is the least anthropogenically impacted of the Laurentian Great Lakes ecosystems, yet dramatic changes to the fish community are evident. Previous published works chronicled those changes and the efforts to rehabilitate the fish community through the year 2000. Here, we review through the year 2022, where post-rehabilitation stability was driven by lean lake trout ( Salvelinus namaycush namaycush) as the most abundant piscivore in nearshore waters, siscowet lake trout ( Salvelinus namaycush siscowet) as the most abundant piscivore in offshore waters, and a healthy, intact assemblage of native prey species, which created ecological redundancies and helped stabilize the food web. Stocking of non-native salmonines was reduced 74%, and populations of Chinook salmon ( Oncorhynchus tshawytscha) and coho salmon ( Oncorhynchus kisutch) were maintained through natural reproduction. Despite reduced stocking, yield from recreational fisheries was stable. Likewise, developments in population modeling led to evaluations and refinement of management strategies that helped create stability for lake trout, lake whitefish ( Coregonus clupeaformis ), and cisco ( Coregonus artedi ) fisheries. With lake trout rehabilitation achieved, focus shifted toward rehabilitation of native brook trout ( Salvelinus fontinalis ), lake sturgeon ( Acipenser fulvescens ), and walleye ( Sander vitreus ). Despite continued control efforts, sea lamprey ( Petromyzon marinus ) abundance increased considerably, and estimates of fish killed by lampreys averaged 2.65 million kg annually. Environmental changes have benefited sea lampreys and fostered thermal habitats more suitable to non-native organisms, posing new challenges for managers and researchers. Nevertheless, the post-rehabilitation stability in the contemporary fish community will help provide resilience to future perturbations in the ecosystem.

Lake Superior↗

USGS2025

Connecting tributary mercury loads to nearshore and offshore sediments in Lake Superior

Lake Superior has a vast and largely undeveloped watershed in comparison to the other Great Lakes, which makes it challenging to study mercury (Hg) sources and cycling. To examine Hg inputs to Lake Superior, we conducted an expansive binational assessment in 40 watersheds from a diverse range of landcover types. We further paired tributary Hg data to sediment source portfolios in the nearshore and offshore zones of Lake Superior through partnership with the Great Lakes Sediment Surveillance Program. We observed that total Hg loads were highest in the spring driven by the combination of elevated Hg concentrations and increased water discharge from snowmelt. In addition, total Hg concentrations in tributaries from remote, heavily forested regions, such as Pukaskwa National Park and the Minnesota Northshore, were higher than the Southshore and Thunder Bay regions. Methylmercury concentrations and loads were more spatially dependent, often corresponding to regions with more wetlands (e.g., Michigan Upper Peninsula). We estimated that the total Hg tributary load to Lake Superior in 2021 was 126 kg per year. To further examine the fate of watershed Hg sources, we examined sediments from 28 sites in Lake Superior using Hg stable isotopes. At open water sites, precipitation was the primary Hg source to sediments, but within nearshore sites Hg originated predominantly from watershed runoff. This work further defines the sources and fate of Hg within Lake Superior and highlights how Hg delivery is intrinsically tied to varying hydrologic regimes.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗

USGS2025

Food web structure of the Lake Superior fish community in 2021–2022

The trophic linkages and ecological requirements of the Lake Superior fish community have not been assessed on a whole lake scale in over a decade. Here, we investigated the trophic dynamics across multiple species and habitat zones of Lake Superior. From April to October of 2021 and 2022, a total of five piscivore, four planktivore, and six benthivore species were collected by region and length class during bottom-trawl, standard gillnet, commercial cisco gillnet, and recreational angler surveys. To assess trophic linkages, stomach contents were measured to estimate biomass consumed and a multivariate analysis was used to assess diet composition by species, length, and region. We found a high degree of interconnectedness in the Lake Superior food web, with Mysis as a critical diet item for most fishes. Native piscivore diets varied by region and with ontogeny. Lake charr were important habitat couplers in the lake, exhibiting a diverse diet and opportunistic foraging strategy. Conversely, Pacific salmon were more restricted in their foraging. Planktivores and benthivores primarily consumed Mysis , with less reliance on Diporeia compared to previous studies. Lake whitefish and cisco were the exception to this pattern, with broader bathymetric depth distributions represented in their diets. We found the food web to be supported by a predominantly native species assemblage, with redundancies at all trophic levels.

Lake Superior↗

USGS2025

Partial diel vertical migration and niche partitioning in Mysis revealed by stable isotopes

Diel vertical migration (DVM) is critical for moving energy and nutrients between surface and deep waters. Mysis sp. (Crustacea: Mysidae) facilitates this process by serving as predator and prey in both benthic and pelagic habitats. Mysis can also exhibit partial DVM (pDVM), where some individuals do not migrate into the pelagia at night or to the benthos during the day. However, whether Mysis pDVM is a fixed (i.e., same individuals migrate) or random (i.e., random individuals migrate) behavior remains unclear. To evaluate that, we tested whether Mysis exhibit pDVM and niche partitioning in Lake Superior by collecting Mysis from benthic and pelagic habitats day and night across depths ranging from 50 to 250 m and estimating their isotopic niche size (δ 13 C, δ 15 N; Corrected Standard Ellipse Area − SEA c ) and overlap among three life stages. At the population level, Mysis exhibited fixed pDVM structured by life stage. Benthic Mysis (12.9 ± 3.2 mm, mean ± SD) were larger than pelagic Mysis during night (9.6 ± 3.6 mm) and day (8.5 ± 3.6 mm). Adult Mysis (> 15 mm) had larger SEA c (1.8 ± 0.4 ‰ 2 ) compared to juveniles (< 10 mm; 0.3 ± 0.1 ‰ 2 ) and sub-adults (10–15 mm; 0.7 ± 0.2 ‰ 2 ), and their isotopic niche did not overlap with smaller life stages. Adults exhibited random pDVM (i.e., high isotopic niche overlap), whereas juveniles exhibited fixed pDVM (i.e., low overlap). Our observations indicate complex behaviors across and within Mysis life stages, likely due to varying pressures including size-selective predation, light and temperature thresholds, and nutritional requirements. Consequently, the benthic behavior of adult Mysis needs to be considered in monitoring programs.

Minnesota, Wisconsin↗

USGS2025

Status and trends in the Lake Superior fish community, 2024

The U.S. Geological Survey has conducted annual fishery surveys across Lake Superior since 1978 that describe trends in fish species occurrence and relative abundance to inform fisheries management and ecosystem health. In 2024, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 72 nearshore locations in June and 36 offshore locations in July. Nearshore bottom trawls collected 22,190 fish represented by 27 species or morphotypes. The number of species collected at each location ranged from 1 to 12, with a median of 5.5 species. Estimated fish biomass at individual locations ranged from <0.1 to 62.9 kg per ha with a lakewide mean of 3.7 kg per ha. Offshore bottom trawls collected 33,634 fish represented by 12 species or morphotypes. Estimated fish biomass at individual locations ranged from 0.6 to 25.8 kg per ha with a lakewide mean of 8.3 kg per ha, which was the second highest for the period-of-record. Lakewide average densities (fish per ha) of age-1 fish were 1 per ha for Bloater, 5 per ha for Cisco, 1 per ha for Lake Whitefish, 60 per ha for Rainbow Smelt, and 19 per ha for Kiyi. Surface trawling collected 5,177 larval Coregonus individuals which was the third fewest Coregonus larvae collected in a whole lake survey since the larval fish survey began in 2014. Nearshore mean larval Coregonus densities were 176 fish per ha in June 2024 and offshore densities were 7 fish per ha in July 2024. June and July surface water temperatures were near the warmest for the period-of-record.

Lake Superior↗

USGS2024

Trophic ecology of juvenile lean and siscowet lake charr (Salvelinus namaycush) in Lake Superior: Assessing for potential competition

We investigated the spatial overlap, diet, isotopic niche, and growth of juvenile lean and siscowet lake charr ( Salvelinus namaycush ) in Lake Superior to address concerns of potential competition with implications to the study of resource polymorphism. Catch data revealed the greatest levels of sympatry in waters from 40 to 60 m. Juvenile lean and siscowet diet changed ontogenetically with Mysis being the dominant prey item for the smallest lake charr but differentiating with onset of piscivory. As ecotypes increased in size, lean diets became dominated by pelagic prey, whereas siscowets had equal proportions of benthic and pelagic prey. Isotopic niche overlap declined between ecotypes coincident with siscowet lake charr shifting to deeper habitats around 400 mm. Lean and siscowet exhibited different growth trajectories. However, length at age-4 declined in parallel for both ecotypes with no trend in condition suggesting that lake charr growth is sensitive to prey biomass and unlikely related to competition. Our findings indicate minimal evidence of competition and support the concept that multiple sympatric ecotypes of lake charr in Lake Superior are maintained by resource polymorphism.

Lake Superior↗

USGS2024

Paired comparisons with quiet surface drones show evidence of fish behavioral response to motorized vessels during acoustic surveys in Lake Superior

Acoustic surveys are important for fish stock assessments, but fish responses to survey vessels can bias acoustic estimates. We leveraged quiet uncrewed surface vessels (USVs) to characterize potential bias in acoustic surveys. Five conventional motorized ships overtook USVs from astern over 2 km transects at night in Lake Superior in 2022. We examined the difference in acoustic backscatter, average target depth, and average target strength (TS) between USV and motorized vessels. Although sound level measurements from the motorized vessels sometimes exceeded recommendations for scientific vessels, we did not detect differences in acoustic measures among survey vessels. However, the USVs recorded 2 dB higher acoustic backscatter and TS than motorized vessels, leading to ~15% higher fish densities with drones when using in situ TS and echo integration. Differences in fish density would increase to 30-60% if a standard TS value was applied. Target depth did not differ between USVs and motorized ships. These results are consistent with a change in orientation but not depth of insonified fish and limited horizontal avoidance of motorized survey vessels.

Lake Superior↗

USGS2024

Revisiting geophysical interpretations of the Midcontinent Rift below Lake Superior— Insights from GLIMPCE seismic-reflection line C

The 1.1 Ga Midcontinent Rift System (MRS) has been investigated in the Lake Superior region for more than a century. The most influential geophysical data for modern paradigms has come from seismic-reflection profiles collected by the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) in the late 1980s. We are revisiting seismic-reflection GLIMPCE line C by constructing a detailed velocity model for conversion of the seismic data measured in two-way travel time to a section plotted versus depth. Preliminary results from the velocity modeling, depth conversion, and ties to onshore geology have led to a significantly different view of Line C as primarily a sag basin rather than a half-graben. Narrow intervals of high-velocity material that emanate upwards along both sides of the sag basin are interpreted as mafic intrusions connected to the mantle. This new interpretation suggests that basin subsidence was the dominant process in the development of rift stage troughs rather than major half-graben structures.

Conference Paper↗

USGS2023

Compressional-wave seismic velocity, bulk density, and their empirical relations for geophysical modeling of the Midcontinent Rift System in the Lake Superior region

Compressional-wave seismic velocity (velocity) and bulk density (density) data were compiled from published sources for rock suites and earth materials that are significant for geophysical modeling of the Mesoproterozoic Midcontinent Rift System in the Lake Superior region. The data include laboratory measurements of outcrop and drill core samples, seismic refraction studies, and a sonic log from a 1.5-kilometer-deep exploration well. Rock suites of the Midcontinent Rift System include basalts of the Mesoproterozoic Keweenawan Supergroup, Oronto Group sedimentary rocks (divided into arenaceous versus argillaceous units), and several sedimentary formations overlying the Oronto Group that have been correlated across the area. Intrusive units include diabase, gabbro, and felsic igneous rocks. Other geologic units important for geophysical modeling in the Lake Superior region include Archean crystalline crust, Paleoproterozoic metasedimentary and crystalline rocks, lower Mesoproterozoic sedimentary rocks, and Holocene to Pleistocene surficial deposits. Empirical velocity-density relations for each rock suite were determined by comparing the compiled data to published relations, such as the Nafe-Drake curve, Gardner’s relation, and best-fit equations developed for different rock types from laboratory studies. Graphical representations of these velocity-density relations provide a way to easily understand how velocity and density differ between tectonic settings and by rock type. Overlaps in velocity and density ranges for different geologic units are significant and have especially important implications for geologic interpretation of seismic data. Important examples include similar velocities but differing densities for argillaceous Oronto Group versus units overlying the Oronto Group and arenaceous Oronto Group versus basalt of the Keweenawan Supergroup. Similar densities but differing velocities were found for diabase versus gabbro. In addition, expected velocity ranges by rock type show that high-velocity intervals (6.9–7.1 kilometers per second) interpreted as basalt in previous seismic-reflection studies more likely indicate diabase or gabbro, suggesting that these interpretations may warrant additional consideration.

Lake Superior↗

USGS2023

Status and trends in the Lake Superior fish community, 2022

In 2022, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 71 nearshore locations in May-June and 35 offshore locations in July, and at 51 Coordinated Science and Monitoring Initiative (CSMI) locations in July-October with bottom trawls, surface trawls, mid-water trawls and acoustics that were previously sampled in 2011 and 2016. Nearshore bottom trawls collected 11,603 fish from 25 species or morphotypes. Nearshore mean biomass was 1.6 kg per ha which was one of the lowest biomass estimates over survey’s 45-year history. Offshore bottom trawls collected 13,876 fish from 11 species or morphotypes. Offshore mean biomass was 5.1 kg per ha, which was less than the annual average since 2011 of 6.5 kg per ha. Recruitment, as measured by age-1 densities, was near zero for Bloater ( Coregonus hoyi ), Cisco ( C. artedi ), and Kiyi ( C. kiyi ), 2 age-1 fish per ha for Lake Whitefish ( C. clupeaformis ) and 77 age-1 fish for Rainbow Smelt ( Osmerus mordax ). All were less than the long-term averages. Sampling at the CSMI locations collected 26 species and morphotypes. The most abundant species’ lakewide were Deepwater Sculpin (all years), young-of-year ciscoe (Bloater, Cisco, and Kiyi, 2022), and Rainbow Smelt (2011 and 2016). Cisco had the highest estimated lakewide biomass in 2011 and 2022 and siscowet Lake Trout had the highest estimated lakewide biomass in 2016. Native species were more abundant than invasive species by numbers (80, 65, and 92%) and biomass (94, 93, 96%) in 2011, 2016, and 2022, respectively. Total lakewide benthic fish biomass declined from 47 thousand metric tons in 2011 to 29 thousand metric tons in 2016 and increased to 33 thousand metric tons in 2022. Total lakewide pelagic fish biomass declined from 61 thousand metric tons in 2011 to 25 thousand metric tons in 2016 and increased to 54 thousand metric tons in 2022. The most unexpected result from our sampling in 2022 was the 2 billion age-0 ciscoe estimate from the mid-water trawl and acoustic sampling in August-October. These fish were broadly distributed across the lake, being collected at 53 of the 54 locations, and their population estimates were highest in the depths >100 m. The factors underlying the survival of these ciscoes into late summer in 2022 as compared to previous years have not been identified, but our annual population surveys of larval ciscoes suggests that lake conditions in June and July may have differed from previous years and enhanced survival. In 2022, ciscoe larval densities in May were lower than average (likely due to a cold winter and spring that delayed hatching), June densities were similar to previous years, and July density estimates were more than double that of any previous year’s estimate.

Lake Superior↗

USGS2023

Exploring the geology of the Midcontinent Rift under western Lake Superior using a preliminary velocity model of seismic line GLIMPCE C

Seismic-reflection data were collected in the 1980s as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) to investigate the 1.1 Ga Midcontinent Rift System (MRS). GLIMPCE Line C crosses western Lake Superior from north to south shores (Fig. 1 inset). Many previous workers have interpreted the MRS in Line C as an asymmetric central graben filled with 10–20 km of subaerial basalt flows, overlain by 7-10 km of sedimentary section, and underlain by magmatic underplating. The central graben was interpreted to have formed from extensional normal faults, later reactivated as high-angle reverse faults. The northern part of Line C crosses over a prominent gravity low called the Grand Marais Ridge (GMR; Fig. 1 inset), previously interpreted as an Archean granitic basement high.

Lake Superior↗

USGS2023

Status and trends in the Lake Superior fish community, 2023

The U.S. Geological Survey annually conducts fishery surveys across Lake Superior that describe trends in fish species occurrence and relative abundance to inform fisheries management and large lake ecology. In 2023, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 51 nearshore locations in June and 31 offshore locations in July. Nearshore bottom trawls collected 157,804 fish from 25 species or morphotypes. Nearshore mean biomass was 18.3 kg per ha which was the second highest biomass estimate over the survey’s 46-year history. Offshore bottom trawls collected 15,458 fish from 10 species or morphotypes. Offshore mean biomass was 5.2 kg per ha, which was less than the annual average of 6.3 kg per ha. Recruitment, as measured by age-1 densities, was the highest recorded for Bloater, Cisco, and Rainbow Smelt in the nearshore and for Kiyi in the offshore survey’s period-of-records. Lakewide average densities (fish per ha) of age-1 fish were 140 for Bloater, 1,019 for Cisco, 616 for Rainbow Smelt, and 54 for Kiyi, which were the highest estimates for the survey’s period-of-record. Period-of-record averages for these species were 10, 67, and 9 age-1 fish per ha, respectively. Age-1 Lake Whitefish averaged 9 fish per ha which was similar to the long-term average of 8 age-1 Lake Whitefish per ha. If the future can be predicted by past large Bloater, Cisco, and Kiyi (collectively, ciscoe) year-class events, the unprecedented survival of the 2022 ciscoe yearclass will influence the Lake Superior ecosystem for the next 10 to 20-years.

Lake Superior↗

USGS2022

Lake Superior Kiyi reproductive biology

Objective The Lake Superior Kiyi Coregonus kiyi is an understudied species being considered for reintroduction into Laurentian Great Lakes where it no longer occurs. Herein, we provide descriptions of Kiyi reproductive biology with the intention of guiding potential gamete collections for propagation. Methods Data were collected on Kiyi spawning timing, spawning locations, spawning season catch rates, length at sexual maturity, sex ratios, fecundity, egg size, and larval occurrences in Lake Superior from 1996–2021. These data were compared to observations made a century prior in Lakes Michigan, Ontario, and Superior. Result Contemporary Kiyi spawning occurred between late December and late January when surface water temperatures cooled to <4°C. Spawning Kiyi were caught almost exclusively in 38.1-mm stretch mesh, as compared to larger meshes (50.8–76.2 mm). Capture depths for developing, ripe, running, and spent female Kiyi were similar and ranged from 82 to 221 m. Fifty percent of female and male Kiyi were classified as sexually mature at ~150 mm total length. Fecundity estimates ranged from 1,578 to 6,720 eggs/female. Mean diameter of unfertilized eggs was 1.7 mm. Recently hatched larval Kiyi were collected at the surface during May–July at 62 of the 113 locations sampled throughout the lake in 2019. Conclusion Our work suggests that Kiyi gamete collection efforts from mid-December through January using 38.1-mm gill-net panels set at bathymetric depths of at least 100 m would maximize the collection of spawning Kiyi and reduce the bycatch of other Coregonus species. Future research questions include the following: (1) “Do Kiyi form spawning aggregations at specific spawning areas, or do they spawn indiscriminately across the lake?”; (2) “Do Kiyi spawn near the bottom or up in the water column?”; (3) “What is the relationship between fall lake overturn and Kiyi spawn timing?”; and (4) “Could summer larval and age-0 Kiyi collections provide an opportunity for establishing a captive broodstock?”

Lake Superior↗

USGS2022

Quantifying status and trends from monitoring surveys: Application to Pygmy Whitefish (Prosopium coulterii) in Lake Superior

Population assessments of fish species often rely on data from surveys with different objectives, such as measuring biodiversity or community dynamics. These surveys often contain spatial–temporal dependencies that can greatly influence conclusions drawn from analyses. Pygmy whitefish (PWF, Prosopium coulterii ) populations in Lake Superior were recently assessed as Threatened by the Committee on the Status of Endangered Species in Canada, which motivated a thorough analysis of available data to improve our understanding of its population status. The US Geological Survey conducts annual bottom trawl surveys in Lake Superior that commonly capture PWF. We used these data (1989–2018) to model temporal trends in PWF biomass density and make lake-wide population projections. We used a Bayesian approach, integrated nested Laplace approximation (INLA), and compared the impact of including different random structures on model fit. Inclusion of spatial structure improved model fit and conclusions differed from models omitting random effects. PWF populations have experienced periodic fluctuations in biomass density since 1989, though 2018 may represent the lowest density in the 30-year time series. Lake-wide biomass was estimated to be 71.5 t.

Lake Superior↗