Geology topics

Charles P. Madenjian

Publications and source records attributed to Charles P. Madenjian.

At least 19 recordsLinked to original sources

USGS2026

The Great Lakes aquatic tissue analysis repository (GLATAR)

Empiricists and modellers use information on energy density, proximate composition, stable isotopes, fatty acids, thiamine, and bioaccumulative tracers (e.g., PCBs and mercury) to understand the state and inter-relationships of aquatic food webs. Data exist in many published and unpublished sources, but are not consolidated in an easily accessible database that would serve as a vital resource to i) provide basic estimates of these diet-derived measures of body composition, ii) understand sources of variation in the underlying data, iii) facilitate exploration and development of data proxies, and iv) assist in study design. We designed GLATAR (Great Lakes Aquatic Tissue Analysis Repository, glatar.org ) to address this need. GLATAR is an open-access, searchable database linked to a web-based toolbox to visualise and generate user-defined summaries on diet-derived ecological metrics, with a focus on taxa of importance to the Great Lakes. GLATAR currently contains over 50,000 records on energy density, chemical tracers, and proximate body composition from 67 species of fish and 72 invertebrate taxa. We hope this user-friendly interface will entice others to upload their published and unpublished data to the repository, enriching the breadth of data accessible to researchers. In this way, GLATAR will become a ‘living’ and interactive resource for empiricists and modellers working in freshwater ecosystems as they make critical decisions related to growth, production, and consumption across a diverse group of economically and ecologically important aquatic species.

Journal of Great Lakes Research↗

USGS2026

Status and trends of pelagic and benthic prey fish populations in Lake Michigan, 2025

Lake Michigan↗

USGS2025

Assessing fish growth changes under an ecosystem regime shift: The approach of linear mixed-effects modeling with application to Lake Huron lake trout

Linear mixed-effects modeling has been used to characterize fish growth changes over time, but early studies did not fully consider the implicit relationship among the three major factors of fish age, year-class, and sampling year in model development. Our objective is to resolve this issue and develop a linear mixed-effects model to clearly assess how lake trout ( Salvelinus namaycush ) growth has responded to the ecosystem regime shift in Lake Huron. We found that, from the late 1970s through the early 2000s, lake trout weight-at-age in northern Lake Huron had large variation with no apparent temporal trend, whereas the weight-at-age in southern Lake Huron steadily declined but remained larger than that in northern Lake Huron. The difference between northern and southern Lake Huron almost disappeared prior to the 2003 collapse of alewives, which accelerated the weight-at-age declines during the 2000s in both regions. By 2010, lake trout growth had a modest recovery and then stabilized. Our model fitted and explained variation in lake trout weight-at-age between years, regions, and recruitment origins, with the fixed age effect representing the general growth trajectory, the fixed year effect as the year-specific growth index, and the random year-class effect as the cohort-specific growth index.

Lake Huron↗

USGS2025

Comparing year-class strength indices from longitudinal analysis of catch-at-age data with those from catch-curve regression: Application to Lake Huron lake trout

Fish year-class strength (YCS) has been estimated via longitudinal analysis of catch-at-age data and via catch-curve regression, but no study has compared the two approaches. The objective of this study was to compare YCS estimates derived from both approaches applied to catch-at-age data for the lake trout (Salvelinus namaycush) population in the main basin of Lake Huron, one of the Laurentian Great Lakes of North America. YCS was reconstructed for both hatchery-stocked and wild lake trout. Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used to compare 14 linear mixed-effects models for longitudinal analysis of catch-at-age data, and three linear mixed-effects models for catch-curve regression. From the best models based on AIC or BIC comparisons, YCS estimates with year-class as a fixed effect were consistent with those estimated with year-class as a random effect. Patterns and trends in the YCS estimates were also the same or similar between the longitudinal analysis of catch-at-age data approach and the catch-curve regression approach, suggesting that both modeling approaches are applicable to a variety of fish populations. indicating that both approaches provide robust measures of YCS. Potential bias in using the approach of catch-curve regression could be caused by abrupt changes in adult mortality. It is also critical to recognize multiple recruitment origins for using the approach of longitudinal analysis of catch-at-age data.

Lake Huron↗

USGS2025

Burbot population dynamics in northern Lake Michigan, 2001–2023, and predation on newly stocked lake trout by burbot

Although the worldwide status of burbot ( Lota lota ) across its Holarctic distribution has recently been assessed, adverse effects of dreissenid mussels on burbot abundance in waters invaded by mussels have not been investigated. We analyzed data from an annual spring gillnet survey to characterize burbot population dynamics in northern Lake Michigan during 2001–2023. Adult burbot abundance in northern Lake Michigan trended neither significantly upward nor downward during 2001–2023. Thus, expansion of the quagga mussel ( Dreissena rostriformis bugensis ) population into deeper waters of Lake Michigan during 2003–2008 had no detectable negative effects on burbot abundance. Burbot continued to feed on lake trout ( Salvelinus namaycush ) stocked in the Northern Refuge of Lake Michigan through 2023, which supported a hypothesis that the burbot population potentially impeded lake trout restoration in the refuge.

Michigan↗

USGS2025

Status and trends of pelagic and benthic prey fish populations in Lake Michigan, 2024

Fall bottom trawl (fall BT) and lakewide acoustic (AC) surveys are conducted annually to generate indices of pelagic and benthic prey fish densities in Lake Michigan. The fall BT survey has been conducted each fall since 1973 using 12-m trawls at depths ranging from 9 to 110 m at fixed locations distributed across seven transects; this survey estimates densities of seven prey fish species [i.e., Alewife ( Alosa pseudoharengus ), Bloater ( Coregonus hoyi ), Rainbow Smelt ( Osmerus mordax ), Deepwater Sculpin ( Myoxocephalus thompsonii ), Slimy Sculpin ( Cottus cognatus ), Round Goby ( Neogobius melanostomus ), Ninespine Stickleback ( Pungitius pungitius )] as well as age-0 Yellow Perch ( Perca flavescens ) and large (> 350 mm) Burbot ( Lota lota ). In recent years, wild juvenile (<400 mm) Lake Trout ( Salvelinus namaycush ) have also become more common in the fall bottom trawl. The AC survey has been conducted each late summer/early fall since 2004 (except 2020). The 2024 AC survey consisted of 24 transects [468 km total (291 miles)] covering bottom depths ranging from 16 to 173 m and 38 midwater trawl tows at 4 to 72 m; this survey estimates densities of three prey fish species (i.e., Alewife, Bloater, and Rainbow Smelt). The data generated from these surveys are used to estimate various population parameters that are, in turn, used by state and tribal agencies in managing Lake Michigan fish stocks. In spring of 2024, an additional spring bottom trawl survey (spring BT) was implemented across six of the transects sampled in the fall and sites ranged in depth from 9 to 237 m. The goal of the spring BT, conducted annually since 2021 with differing levels of effort, was to explore seasonal differences in biomass density and distributions of key prey species, most notably Alewife. Total prey fish biomass density from the spring BT was 5.7 kg/ha. For the AC survey, total biomass density of prey fish equaled 10.8 kg/ha, more than double the long-term average (20042023) of 5.1 kg/ha but 4.0 kg/ha lower than the 2023 estimate. For the fall BT, total biomass density of prey fish equaled 2.1 kg/ha, the lowest value since 2020 and 69% lower than the average from 2004-2023 (6.8 kg/ha). The 2024 fall BT biomass density was only 6.3% of the average over the entirety of the time series (1973-2023; 33.1 kg/ha). Over the period both surveys have been conducted (2004-2024), total biomass density has trended downward in the fall BT (despite a high 2022 estimate) and remained relatively stable in the AC survey. Deepwater Sculpin and Bloater were the most common species (by biomass) among prey fishes in the spring BT while the AC survey and fall BT reported co-dominance of Bloater and Alewife. Mean biomass of yearling and older (YAO) Alewife was 1.30 kg/ha in the spring BT, 4.7 kg/ha in the AC survey, and 0.68 kg/ha in the fall BT. Since 2014, annual survey results suggest that the catchability of YAO Alewives for the fall BT is substantially lower than the AC survey. Like previous spring surveys, Alewives were aggregated in deeper habitats, with 93% of biomass collected between 110 and 201 m. Results of the 2024 spring BT align with past spring surveys and do not suggest that spring bottom trawling provides a better index of age-2 and older Alewives than fall bottom trawling, even with adjustments for differences in habitat use. However, the spring BT does appear to index age-1 Alewives more effectively than the fall BT. The 2024 AC survey YAO Alewife biomass density estimate was 77% higher than the average from 2004-2023. The Alewife population of Lake Michigan appears to be composed mostly of young fish and the proportion of age-4 and older Alewives was <1.8% in each of the three surveys. Age-0 Alewife numeric density from the AC survey was 510 fish/ha in 2024, slightly higher than the long-term mean (486 fish/ha). Biomass density of large (≥120 mm) Bloater was 5.2 kg/ha in the AC survey and 0.76 kg/ha in the fall BT, while total Bloater biomass in the spring BT was 1.8 kg/ha - all three estimates were much lower than what was estimated by the fall BT between 1981 and 1998. The density of small (<120 mm) Bloater was 456 fish/ha in the AC survey, the second highest value in the time series and potentially reflective of an above average 2024 year-class. Meanwhile, small Bloater density estimated in the fall BT was only 16 fish/ha. Biomass density of large Rainbow Smelt (≥90 mm) was 0.21 kg/ha in the AC survey and 0.03 kg/ha in the fall BT survey, continuing the trend of low large Rainbow Smelt biomass observed since 2001. Numeric density of small (<90 mm) Rainbow Smelt was 31 fish/ha in the AC survey and 143 fish/ha in the fall BT. All four prey fish species indexed only by the fall BT had below-average biomass densities regardless of trawling season. Deepwater Sculpin biomass density was 0.26 kg/ha, which makes 14 of the past 15 years with biomass <1 kg/ha. Spring BT Deepwater Sculpin biomass density (2.0 kg/ha) was higher than any fall BT estimate since 2006, likely reflective of including bottom trawls at greater depths in the spring than the fall. Slimy Sculpin was estimated to be < 0.04 kg/ha in the spring and fall BT, an order of magnitude lower than the long-term average from the fall BT. Round Goby biomass density estimates were low and similar across seasons (0.43 kg/ha in the spring and 0.10 kg/ha in the fall). Ninespine Stickleback density was 3.9 fish/ha in the fall BT and no fish were collected in the spring BT.

Lake Michigan↗

USGS2025

Evaluating how changing climate and water clarity can affect restoration of native coregonine fishes in midwestern lakes

In the upper midwestern United States, coregonine fishes (e.g., ciscoes and whitefishes) have provided a key food source for Indigenous peoples since time immemorial. In the last century, however, several anthropogenic stressors including overfishing, declining quality of key habitat (e.g., dams, algal blooms), and negative interactions with invasive species have led to declines in the fisheries that coregonines support. Hence, fishery managers have undertaken steps to conserve existing diversity or re-introduce species that have been locally extirpated. This study assessed the extent to which changes in water clarity and climate pose challenges for conservation of coregonines and sought to provide information to inform management decisions about which populations would be best for reintroduction. Through experiments, we found that a Cisco ( Coregonus artedi ) population at the southern edge of its range had no marked advantages in respiration or thermal tolerance, relative to a northern Great Lake population. This result indicates that the latter population would be a reasonable candidate for reintroduction in Lake Erie, where warming water could be a potential impediment. Other laboratory experiments revealed that Cisco could be vulnerable to increasing ultraviolet radiation (UV; 280 – 320 nm) penetration in lakes due to reduced ice cover and increased water clarity owing to water quality improvements and/or invasive mussel filtration. Our experiments found that UV exposure causes earlier hatching, more irregular pigmentation near the head and yolk-sac region, and reduced heart rate. Therefore, UV could be a heretofore unrealized threat to Cisco populations. Through the completion of additional final products from this project and additional research connecting laboratory and field data together, the extent to which the eggs and larvae of coregonine species could be vulnerable to increased UV associated with increasing water clarity or reduced ice cover will help inform restoration or conservation planning for these fish.

Cooperator Report↗

USGS2024

Fish bioenergetics modeling

A bioenergetics model for a fish is defined as a quantitative description of the fish's energy budget. Bioenergetics modeling can be applied to free-ranging fish in rivers, lakes, and oceans to estimate annual food consumption by fish populations; such applications have proven useful in managing fisheries. In addition, bioenergetics models have been used to better understand fish growth and consumption in ecosystems, to guide conservation management of endangered or threatened fish populations, to determine the importance of the role of fish in cycling nutrients within ecosystems, to predict fish growth and consumption under climate change , and to identify the important factors regulating contaminant accumulation in free-ranging fish. Over the past 60 years, fish bioenergetics modeling has led to significant advancements in both aquatic ecology, especially fish ecology, and fisheries management.

Book chapter↗

USGS2024

Spatial distribution and diet of Lake Michigan juvenile lake trout

Most studies of Lake Michigan lake trout ( Salvelinus namaycush ) have focused on adults, with scant attention to juveniles (<400 mm). We explored the spatial distribution and diet of juvenile lake trout using U.S. Geological Survey September bottom trawl data (2015–2022) and stomach content information opportunistically collected since 2012 by various agencies using multiple gear types. Most juvenile lake trout in the September bottom trawl survey were caught at 37–64 m bottom depths. Length frequency data from the bottom trawl survey identified three size classes likely associated with wild juvenile lake trout age: < 85 mm (∼age-0), 85–170 mm (∼age-1) and > 170 mm (∼age-2+). Largest catches of wild lake trout < 170 mm occurred along a northeastern transect (near Frankfort, Michigan), whereas most > 170 mm were collected along southern transects. Mysis diluviana was the dominant prey for juvenile lake trout < 170 mm, and > 250 mm were primarily piscivorous, while 170–250 mm appeared to be a transitional period of switching from Mysis to fish. Species composition of prey fishes consumed by lake trout varied spatially and we found evidence of seasonal and annual diet variation within Grand Traverse Bay. Diporeia , once an important component of juvenile lake trout diet, appears to no longer be consumed by juvenile lake trout in Lake Michigan to any measurable degree. Continued research on the ecology of juvenile lake trout may provide insight into the effects of a changing ecosystem on juvenile lake trout diet and growth, thereby contributing to the effort to rehabilitate the Lake Michigan lake trout population.

Lake Michigan↗

USGS2024

Comparison of cisco (Coregonus artedi) aerobic scope and thermal tolerance between two latitudinally-separated populations

The cisco Coregonus artedi is a coldwater fish that is distributed throughout much of Canada and the northern United States, including the Laurentian Great Lakes. Cisco historically supported large commercial fisheries in the Great Lakes during the late 1800s and early 1900s, but many populations declined and never recovered. Restoration efforts focusing on re-establishing cisco in the Great Lakes are underway, but increasing water temperatures may hinder these efforts. Therefore, we examined aerobic scope and thermal tolerance of allopatric cisco populations from different latitudes and habitats to determine if a southern latitude population (Crooked Lake, Indiana, USA) near the southern edge of cisco distribution was better adapted to withstand warmer water temperatures than a northern latitude population (Les Cheneaux Islands, Michigan, USA; Lake Huron). As expected, both stocks demonstrated increases in metabolic rates and absolute aerobic scope with increased temperature. Northern cisco had significantly lower aerobic scope compared to southern cisco at both treatment temperatures of 10 and 13 °C. Both cisco stocks had high thermal tolerances when challenged by temperatures increased to 20 and 23 °C but low tolerances at 26 °C. Cisco thermal tolerances increased with acclimation temperature, but we did not detect a difference in thermal tolerances between northern and southern cisco. Although southern cisco had higher capacity for aerobic metabolism, both stock sources had high thermal tolerances at the upper end of their thermal limits. Therefore, either population would be likely suitable for reintroduction into Great Lakes habitats, even with expected warming in the future.

Indiana, Michigan↗

USGS2024

Synchrony of alewife, Alosa pseudoharengus, year-class strength in the Great Lakes region

Fish recruitment is interannually variable and challenging to predict. Annual recruitment is often regionally synchronized among populations and identifying drivers of such synchrony may help shed light on recruitment dynamics. We investigated interannual variation of alewife Alosa pseudoharengus recruitment by estimating year-class strength for populations from three of the Laurentian Great Lakes (lakes Huron, Michigan, and Ontario) using annual assessment survey data from 1968 to 2022. We first determined the age when year-class strength was set for each dataset. We then used a mixed-modeling approach incorporating age, year-class, and sampling year to estimate year-class strength for each study lake. Thereafter, we evaluated regional synchrony in year-class strength across the three lakes and evaluated potential climatic drivers of synchrony. Our results suggest that alewife year-class strength is set by age-1. Our model-derived indices of alewife year-class strength were synchronized among the three lakes, and we identified spring-summer degree-days as a potential regional driver of synchrony. This analysis highlights the potential for strong influence of climatic conditions on fish recruitment in large freshwater systems.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2024

Effects of temperature, body size, and sex on the standard metabolic rates of a sexually dimorphic freshwater piscivore, walleye (Sander vitreus)

Walleye ( Sander vitreus ) are a sexually dimorphic freshwater piscivore that have long been studied using bioenergetics modeling, yet robust estimates of metabolic rates for use in those models have been lacking. Therefore, we quantified the effects of body mass, water temperature, and sex, on standard metabolic rate (SMR). We estimated SMR across temperatures ranging 0.5–24 °C using intermittent-flow respirometry for male ( n = 54) and female ( n = 53) walleye raised in hatchery conditions. We found a significant interaction between sex and temperature, whereby males had lower SMR below 16 °C, but increased such that SMR was ∼16% higher than females at 22 °C. The mass-scaling coefficient, b , was similar for both sexes at a value around 1. These estimates are the first to be generated using respirometry for adult walleye and differ slightly from those being used by previous bioenergetics models, which may have underestimated the energetic costs of SMR, and did not include sex-specific estimates for metabolism.

Canadian Journal of Fisheries and Aquatic Sciences↗

USGS2024

Status and trends of pelagic and benthic prey fish populations in Lake Michigan, 2023

Fall bottom trawl (fall BT) and lakewide acoustic (AC) surveys are conducted annually to generate indices of pelagic and benthic prey fish densities in Lake Michigan. The fall BT survey has been conducted each fall since 1973 using 12-m trawls at depths ranging from 9 to 110 m at fixed locations distributed across seven transects; this survey estimates densities of seven prey fish species [i.e., Alewife ( Alosa pseudoharengus ), Bloater ( Coregonus hoyi ), Rainbow Smelt ( Osmerus mordax ), Deepwater Sculpin ( Myoxocephalus thompsonii ), Slimy Sculpin ( Cottus cognatus ), Round Goby ( Neogobius melanostomus ), Ninespine Stickleback ( Pungitius pungitius )] as well as age-0 Yellow Perch ( Perca flavescens ) and large (> 350 mm) Burbot ( Lota lota ). The AC survey has been conducted each late summer/early fall since 2004, and the 2023 survey consisted of 20 transects [450 km total (336 miles)] covering bottom depths ranging from 12 to 248 m and 29 midwater trawl tows above bottom depths ranging 16 to 246 m; this survey estimates densities of three prey fish species (i.e., Alewife, Bloater, and Rainbow Smelt). The data generated from these surveys are used to estimate various population parameters that are, in turn, used by state and tribal agencies in managing Lake Michigan fish stocks. A spring bottom trawl survey (spring BT) was implemented across 3 of the transects sampled in the fall and sites ranged in depth from 18 to 164 m. The goal of the spring BT is to explore seasonal differences in biomass density and distributions of key prey species, most notably Alewife. Additionally, we conducted acoustic sampling while bottom trawling to evaluate the vertical distribution of fish relative to the height of the trawl. The abbreviated spring BT survey results indicated that Alewives were primarily offshore with peak biomass density at the 91 m bottom depth. There was no evidence of higher acoustic density above the trawl at depths of 18, 73, 128, and 146 m. At 91 and 164 m, acoustic density above the trawl was >2x that in the trawl path, but sample size at these two depths was low. For the AC survey, total biomass density of prey fish equaled 14.8 kg/ha, 223% higher than the longterm average (2004-2022) of 4.6 kg/ha and 8.8 kg/ha higher than the 2022 estimate. For the fall BT, total biomass density of prey fish equaled 3.6 kg/ha, about 50% lower than the average value from 2004-2022 (6.9 kg/ha). The 2023 fall BT biomass was an order of magnitude lower than the average over the entirety of the time series (1973-2022; 33.7 kg/ha). Bloater was the dominant species (by biomass) among prey fishes in the fall BT, while the AC survey reported dominance of Alewife. Mean biomass of yearling and older (YAO) Alewife was 10.3 kg/ha in the AC survey, and 0.7 kg/ha in the fall BT. Since 2014, catchability of YAO Alewives for the fall BT has been substantially lower than the AC survey. While limited in scope, the results of the 2023 spring BT do not suggest that catchability is substantially higher in the spring than the fall, which aligns with the 2022 survey results. Comparing the AC estimate to previous years, YAO Alewife biomass was 359% higher than the average from 2004-2022. Numeric density of age-0 Alewife from the AC survey was 1,205 fish/ha in 2023, which is the third highest in the time series and well above the long-term mean of 428 fish/ha. Biomass density of large (≥120 mm) Bloater in 2023 was 3.5 kg/ha in the AC survey and 2.1 kg/ha in the fall BT - each at least an order of magnitude lower than what was estimated by the fall BT between 1985 and 1997. Following a record high year in 2021 (1,034 fish/ha), the numeric density of small (<120 mm) Bloater was 142 fish/ha in the AC survey, similar to the long-term mean of 120 fish/ha. Meanwhile, small Bloater density estimated in the fall BT was 2 fish/ha. Biomass density of large Rainbow Smelt (≥90 mm) was <0.05 kg/ha in the AC and fall BT surveys, continuing the trend of low Rainbow Smelt biomass that has been observed since 2001. Numeric density of small (<90 mm) Rainbow Smelt was 119 fish/ha in the AC survey and 7 fish/ha in the fall BT, indicating a weak year-class. All four prey fish species sampled only by the fall BT indicated below average biomass densities. Deepwater Sculpin biomass density was estimated at 0.4 kg/ha, which makes 13 of the past 14 years when biomass was <1 kg/ha. Slimy Sculpin was estimated at 0.02 kg/ha, only 5% of the long-term average. Round Goby was estimated at 0.3 kg/ha, below the average biomass of 0.85 kg/ha since 2008 but similar to intermittent low values observed throughout the dataset. Ninespine Stickleback density was 1 fish/ha. Only 35 small (<100 mm) Yellow Perch were caught, indicating a weak Yellow Perch year-class in 2023.

Lake Michigan↗

USGS2023

Estimates of adult Lake Trout mortality from coded wire tags in a population with developing natural reproduction in southern Lake Michigan

Objective Overfishing and Sea Lamprey Petromyzon marinus predation led to extirpation of Lake Trout Salvelinus namaycush from Lake Michigan in the 1950s. Large populations of hatchery‐reared fish were developed by the 1970s, but natural reproduction was limited until the early 2000s when it began to increase in the southern main basin. Hypothesizing that the relatively low mortality of spawning‐aged fish contributed to this reproductive success, we estimated the total annual mortality rate for this population. Methods We used catch curves to estimate the total instantaneous mortality rate Z using coded wire tags, which provided definitive ages. We made separate estimates from fish collected in three on‐going surveys: a spring gill‐net survey, a fall gill‐net spawning survey, and a sport fishery survey. Result Our estimates of Z ± SE were 0.297 ± 0.019, 0.239 ± 0.009, and 0.205 ± 0.007 for the spring, spawning, and sport fishery surveys, respectively. We suggest that the mean Z ± SE of all survey estimates of 0.247 ± 0.027 would be a reasonable estimate for this population, which equates to a total annual mortality of 22 ± 3%. This estimate is in the low range of rates reported for the species and is in the same range as other populations in the Great Lakes with well‐established natural reproduction. Conclusion We concluded that these low total mortality rates contributed to the reproductive success in southern Lake Michigan through increasing spawning stock density and age structure and that previous estimates of another important population parameter, the instantaneous natural mortality rate M , were too high. Estimates of M ranged from 0.210 to 0.240 and were based on the Pauly equation, a growth‐ and temperature‐based estimator. We suggest maximum‐age‐based estimators of M are more appropriate for Lake Trout. Several alternative maximum‐age‐based estimators produced estimates for M of 0.132–0.058, all of which are more compatible with our estimate of Z .

Lake Michigan↗

USGS2023

Lake Trout population dynamics in the Northern Refuge of Lake Michigan: Importance of stocking rate

Objective The Northern Refuge of Lake Michigan was established under the auspices of the Great Lakes Fishery Commission in 1985 as an additional step toward restoring a self‐sustaining Lake Trout Salvelinus namaycush population to the lake. The overall goal of our study was to assess progress toward Lake Trout rehabilitation in the Northern Refuge through 2021. Methods We conducted annual (1992–2021) gill‐net surveys in the fall to assess the adult population and beam trawl surveys in the spring to assess naturally reproduced age‐0 Lake Trout. Result Spawner abundance averaged 45 fish • 305 m of gill net −1 • day −1 during 1992–1999, just 4 fish • 305 m of gill net −1 • day −1 during 2000–2014, and then 67 fish • 305 m of gill net −1 • day −1 during 2018–2021. Sea Lamprey Petromyzon marinus wounding rate decreased by nearly 90% between the 2000–2014 and 2015–2021 periods. The average annual mortality rate of adult Lake Trout decreased from 43% during the 1990s and 2000s to 22% during 2015–2021. All of these population dynamics were partly driven by changes in Lake Trout stocking rates, which declined by roughly 50% during 1994–1996 but then increased roughly twofold during 2006–2010 and remained at this higher level during 2010–2018. In 2021, early signs of natural recruitment of Lake Trout became evident. Beam trawling during June 2021 yielded a wild age‐0 Lake Trout, the first ever caught since the surveys began in 1992. In addition, 15 (4% of total catch) of the Lake Trout caught in the fall 2021 gill‐net survey were unclipped (wild) fish. Estimated ages of these wild fish indicated that natural recruitment had begun sometime around 2015. Conclusion With Alewife Alosa pseudoharengus abundance in Lake Michigan at record low levels, Sea Lamprey wounding rate greatly reduced, and Lake Trout spawner abundance at relatively high levels, the outlook for Lake Trout rehabilitation in the Northern Refuge holds promise.

Northern Refuge of Lake Michigan↗

USGS2023

A generalized application of the catch-curve regression with comparisons of adult mortality and year-class strength between hatchery-stocked and wild-reared lake trout in US waters of Lake Huron

The recently developed approach to estimating the instantaneous total mortality of coded-wire-tagged lake trout ( Salvelinus namaycush ) is generally applicable to catch-at-age data. We further formalized the technique to objectively incorporate the year-class and year effects into the model structure of catch-curve regression. We used this new method to compare adult mortality and year-class strength between the hatchery-stocked and wild-reared lake trout in US waters of Lake Huron, one of the Laurentian Great Lakes. Model comparisons showed no difference in adult mortality between the hatchery-stocked and wild-reared lake trout. Based on 95% confidence intervals, the estimate of adult mortality using the simple catch-curve regression with average number-at-age was not statistically different from the estimate using the linear mixed model with individual number-at-age of multiple year-classes and sampling years. The linear mixed model, however, also quantified lake trout year-class strength and indicated that since 2003, the increases in recruitment of wild-reared lake trout did not fully compensate for the rapid declines in recruitment of hatchery-stocked lake trout in Lake Huron.

Lake Huron↗

USGS2023

Status and trends of pelagic and benthic prey fish populations in Lake Michigan, 2022

Fall bottom trawl (fall BT) and lakewide acoustic (AC) surveys are conducted annually to generate indices of pelagic and benthic prey fish densities in Lake Michigan. The fall BT survey has been conducted each fall since 1973 using 12-m trawls at depths ranging from 9 to 110 m at fixed locations distributed across seven transects; this survey estimates densities of seven prey fish species [i.e., Alewife (Alosa pseudoharengus), Bloater ( Coregonus hoyi ), Rainbow Smelt ( Osmerus mordax ), Deepwater Sculpin ( Myoxocephalus thompsonii ), Slimy Sculpin ( Cottus cognatus ), Round Goby ( Neogobius melanostomus ), Ninespine Stickleback ( Pungitius pungitius )] as well as age-0 Yellow Perch ( Perca flavescens ) and large (> 350 mm) Burbot ( Lota lota ). The AC survey has been conducted each late summer/early fall since 2004, and the 2022 survey consisted of 26 transects [570 km total (354 miles)] covering bottom depths ranging from 5 to 255 m and 37 midwater trawl tows above bottom depths ranging 5 to 232 m; this survey estimates densities of three prey fish species (i.e., Alewife, Bloater, and Rainbow Smelt). The data generated from these surveys are used to estimate various population parameters that are, in turn, used by state and tribal agencies in managing Lake Michigan fish stocks. In spring of 2022, an additional spring bottom trawl survey (spring BT) was implemented across six of the transects sampled in the fall and sites ranged in depth from 9 to 236 m. The goal of the spring BT was to explore seasonal differences in biomass density and distributions of key prey species, mostly notably Alewife. Total prey fish biomass density from the spring BT was 2.1 kg/ha. For the AC survey, total biomass density of prey fish equaled 6.2 kg/ha, 37% higher than the long-term average (2004-2021) of 4.5 kg/ha and 0.43 kg/ha higher than the 2021 estimate. For the fall BT, total biomass density of prey fish equaled 8.7 kg/ha, the highest value since 2013 and 21% higher than average value from 20042021 (6.8 kg/ha). The 2022 fall BT biomass density was still well below the average over the entirety of the time series (1973-2021; 34.3 kg/ha). Over the period both surveys have been conducted (2004-2021), total biomass density has trended downward in the fall BT (despite a high 2022 estimate) and remained relatively stable in the AC survey. Bloater was the dominant species (by biomass) among prey fishes in both the spring and fall BT, while the AC survey reported co-dominance of Bloater and Alewife. Mean biomass of yearling and older (YAO) Alewife was 0.38 kg/ha in the spring BT, 3.0 kg/ha in the AC survey, and 0.10 in the fall BT. Alewife were aggregated in deepwater habitats in the spring of 2022 (> 110 m). Since 2014, catchability of YAO Alewives for the fall BT has been substantially lower than the AC survey. Results of the 2022 spring BT do not suggest that catchability is substantially higher in the spring than the fall. Comparing the acoustic estimate to previous years, YAO Alewife biomass was 40% higher than the average from 2004-2021. An age-7 fish was recorded for the first time since 2009. Despite the rare catches of older fish, the Alewife age distribution still appears truncated, with age-1 fish as the most represented age class in all three surveys. Numeric density of age-0 Alewife from the AC survey was 7 fish/ha in 2022, which is the third lowest in the time series and well below the longterm mean of 452 fish/ha. Biomass density of large (≥120 mm) Bloater was 2.7 kg/ha in the AC survey and 4.4 kg/ha in the fall BT - each at least an order of magnitude lower than what was estimated by the fall BT between 1981 and 1998. Following a record high year in 2021 (1,037 fish/ha), the numeric density of small (<120 mm) Bloater was only 15 fish/ha in the AC survey. Meanwhile, small Bloater density estimated in the fall BT was 261 fish/ha, the highest value since 1990 and likely partially reflective of a large 2021 year-class. Biomass density of large Rainbow Smelt (≥90 mm) was 0.29 kg/ha in the AC survey and 0.12 kg/ha in the fall BT survey, continuing the trend of low Rainbow Smelt biomass that has been observed since 2001. Numeric density of small (<90 mm) Rainbow Smelt was 21 fish/ha in the AC survey and 2.7 fish/ha in the fall BT, indicating a weak year-class. All four prey fish species sampled only by the fall BT indicated below average biomass densities. Deepwater Sculpin biomass density was estimated at 0.41 kg/ha, which makes 12 of the past 13 years when biomass was <1 kg/ha. Slimy Sculpin was estimated at 0.10 kg/ha, the highest estimate since 2016 but still only 25% of the long-term average. Round Goby was estimated at 1.3 kg/ha, above the average biomass of 0.82 kg/ha since 2008 but similar to intermittent high values observed throughout the dataset. Ninespine Stickleback density was 1.5 fish/ha. Burbot biomass remained near record low levels, and no age-0 Yellow Perch were caught, indicating a weak Yellow Perch year-class in 2022.

Lake Michigan↗

USGS2022

Egg hatchability and egg thiamine concentrations in Lake Trout from the Northern Refuge of Lake Michigan

The Northern Refuge of Lake Michigan was established in 1985 under the auspices of the Great Lakes Fishery Commission as part of an intensely managed effort to restore a self‐sustaining Lake Trout Salvelinus namaycush population to the lake. Despite the stocking of millions of Lake Trout into the refuge and survival of these stocked fish to adulthood, signs of natural recruitment by Lake Trout in the refuge have been very sparse. To ascertain whether low egg viability or poor fry survival could be limiting natural recruitment by Lake Trout in the refuge, hatching rates of eggs taken from Lake Trout caught in the Northern Refuge during October 2019 were determined in the laboratory. Moreover, survival rates of the hatched‐out fry were determined along with the thiamine concentrations in the fry. In addition, thiamine concentrations in eggs of ovulating Lake Trout caught in the Northern Refuge during 2019–2021 were determined. Egg hatching rate in the laboratory averaged 93%. Survival rate of Lake Trout fry up to 5 weeks after hatching averaged 97%. Mean (SE) thiamine concentrations in eggs of Lake Trout caught in 2019, 2020, and 2021 were 19.1 (3.6), 5.9 (1.3), and 6.2 (0.4) nmol/g, respectively; these levels are well above the threshold concentration of 2.3 nmol/g, above which the likelihood of experiencing mortality due to thiamine deficiency complex (TDC) is low. Thiamine concentrations in the Lake Trout fry at the end of the laboratory experiment averaged 11.2 nmol/g, and abnormal behavior in the fry was not observed. We conclude that TDC does not appear to represent a major impediment to Lake Trout rehabilitation in the Northern Refuge at present, although TDC may have contributed to the lack of detectable natural recruitment in the past.

Lake Michigan, Northern Refuge↗