An official website of the United States government
The phenology of critical biological events in aquatic ecosystems are rapidly shifting due to climate change. Growing variability in phenological cues can increase the likelihood of trophic mismatches, causing recruitment failures in commercially, culturally, and recreationally important fisheries. We tested for changes in spawning phenology of regionally important walleye (Sander vitreus) populations in 194 Midwest US lakes in Minnesota, Michigan, and Wisconsin spanning 1939-2019 to investigate factors influencing walleye phenological responses to climate change and associated climate variability, including ice-off timing, lake physical characteristics, and population stocking history. Data from Wisconsin and Michigan lakes (185 and 5 out of 194 total lakes, respectively) were collected by the Wisconsin Department of Natural Resources (WDNR) and the Great Lakes Indian Fish and Wildlife Commission (GLIFWC) through standardized spring walleye mark-recapture surveys and spring tribal harvest season records. Standardized spring mark-recapture population estimates are performed shortly after ice-off, where following a marking event, a subsequent recapture sampling event is conducted using nighttime electrofishing (typically AC – WDNR, pulsed-DC – GLIFWC) of the entire shoreline including islands for small lakes and index stations for large lakes (Hansen et al. 2015) that is timed to coincide with peak walleye spawning activity (G. Hatzenbeler, WDNR, personal communication; M. Luehring, GLIFWC, personal communication; Beard et al. 1997). Data for four additional Minnesota lakes were collected by the Minnesota Department of Natural Resources (MNDNR) beginning in 1939 during annual collections of walleye eggs and broodstock (Schneider et al. 2010), where date of peak egg take was used to index peak spawning activity. For lakes where spawning location did not match the lake for which the ice-off data was collected, the spawning location either flowed into (Pike River) or was within 50 km of a lake where ice-off data were available (Pine River) and these ice-off data were used. Following the affirmation of off-reservation Ojibwe tribal fishing rights in the Ceded Territories of Wisconsin and the Upper Peninsula of Michigan in 1987, tribal spearfishers have targeted walleye during spring spawning (Mrnak et al. 2018). Nightly harvests are recorded as part of a compulsory creel survey (US Department of the Interior 1991). Using these records, we calculated the date of peak spawning activity in a given lake-year as the day of maximum tribal harvest. Although we were unable to account for varying effort in these data, a preliminary analysis comparing spawning dates estimated using tribal harvest to those determined from standardized agency surveys in the same lake and year showed that they were highly correlated (Pearson’s correlation: r = 0.91, P < 0.001). For lakes that had walleye spawning data from both agency surveys and tribal harvest, we used the data source with the greatest number of observation years. Ice-off phenology data was collected from two sources – either observed from the Global Lake and River Ice Phenology database (Benson et al. 2000)t, or modeled from a USGS region-wide machine-learning model which used North American Land Data Assimilation System (NLDAS) meteorological inputs combined with lake characteristics (lake position, clarity, size, depth, hypsography, etc.) to predict daily water column temperatures from 1979 - 2022, from which ice-off dates could be derived (https://www.sciencebase.gov/catalog/item/6206d3c2d34ec05caca53071; see Corson-Dosch et al. 2023 for details). Modeled data for our study lakes (see (Read et al. 2021) for modeling details), which performed well in reflecting ice phenology when compared to observed data (i.e., highly significant correlation between observed and modeled ice-off dates when both were available; r = 0.71, p < 0.001). Lake surface area (ha), latitude, and maximum depth (m) were acquired from agency databases and lake reports. Lake class was based on a WDNR lakes classification system (Rypel et al. 2019) that categorized lakes based on temperature, water clarity, depth, and fish community. Walleye stocking history was defined using the walleye stocking classification system developed by the Wisconsin Technical Working Group (see also Sass et al. 2021), which categorized lakes based on relative contributions of naturally-produced and stocked fish to adult recruitment by relying heavily on historic records of age-0 and age-1 catch rates and stocking histories. Wisconsin lakes were divided into three groups: natural recruitment (NR), a combination of stocking and natural recruitment (C-ST), and stocked only (ST). Walleye natural recruitment was indexed as age-0 walleye CPE (number of age-0 walleye captured per km of shoreline electrofished) from WDNR and GLIFWC fall electrofishing surveys (see Hansen et al. 2015 for details). We excluded lake-years where stocking of age-0 fish occurred before age-0 surveys to only include measurements of naturally-reproduced fish.
more »
« less
This content will become publicly available on December 1, 2025
Diminishing productivity and hyperstable harvest in northern Wisconsin walleye fisheries
Managing fisheries in a changing socio-ecological environment may require holistic approaches for identifying and adapting to novel ecosystem dynamics. Using 32 years of Ceded Territory of Wisconsin (CTWI) walleye ( Sander vitreus) data, we estimated production ( P), biomass ( B), biomass turnover ( P/B), yield ( Y), and yield over production ( Y/P) and tested for hyperstability in walleye yield. Most CTWI walleye populations showed low P and B, and Y/P < 1 . Yet, production overharvest ( Y/P > 1) was prevalent among Wisconsin walleye recruitment-based management approaches (natural recruitment (NR), sustained only by stocking, combination). Production, B, and P/B have declined in NR populations, while Y and Y/P have remained constant. Walleye Y was hyperstable along a production gradient among all management approaches and fishery types (i.e., angling only, angling/tribal harvest combined). Diminishing productivity and hyperstable yield may be jointly contributing to observed walleye declines. We classified lakes into management groups of low, moderate, or high vulnerability to harvest based on Y/P and P/B dynamics and identify that harvest may benefit from declines to maintain or increase the adaptive capacity of CTWI walleye.
more »
« less
- Award ID(s):
- 2025982
- PAR ID:
- 10571748
- Publisher / Repository:
- Canadian Science Publishing
- Date Published:
- Journal Name:
- Canadian Journal of Fisheries and Aquatic Sciences
- Volume:
- 81
- Issue:
- 12
- ISSN:
- 0706-652X
- Page Range / eLocation ID:
- 1650 to 1665
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Biological reference points for fishery management depend on estimates of current stock status relative to unfished biomass (depletion). The ratio of fish density outside to inside a marine reserve, the density ratio, could serve as a proxy for depletion for data-poor management. However, transient dynamics associated with time lags in returning to the unfished state following reserve implementation make that proxy inaccurate on short time scales. We assessed density ratio management rules using an age-structured, spatially explicit model of four US west coast nearshore fishes following reserve implementation, with scenarios encompassing sampling error, recruitment variability, and uncertainty in natural mortality. In deterministic simulations, management incorporating time lags generally resulted in a higher mean and lower variability in biomass over 20 years, but lower mean yield compared to management that did not. However, when stochastic recruitment was included, differences among simulations due to stochasticity were much greater than any difference in performance between management strategies. Nonetheless, in certain cases, accounting for time lags could help avoid unwarranted increases in harvest effort after reserve implementation.more » « less
-
Abstract Yield improvement in cotton could be accelerated through selection for functional yield drivers such as interception of cumulative photosynthetically active radiation (∑IPAR), radiation use efficiency (RUE), and harvest index (HI). However, information on the extent to which these traits vary in cotton in the southeastern United States is limited. It was hypothesized that functional yield drivers would vary significantly within a diverse cotton collection. This study was conducted in Tifton and Athens, GA, and included a total of 4 site‐years. Lint yield, total biomass production, ∑IPAR, RUE, and HI were all affected by genotype. Biomass was more strongly correlated with RUE than ∑IPAR. Even among the highest yielding genotypes, values for functional yield drivers (biomass and harvest index) differed significantly, indicating that high yields could be achieved by differentially manipulating these underlying traits. However, when considered for all genotypes, only HI exhibited a significant positive correlation with yield. Boll production and intra‐boll yield components were also affected by genotype. When considered across upland genotypes, lint per boll, lint per seed, and lint percent were strongly associated with HI and lint yield, whereas boll mass and seed number per boll were not. We conclude that the genotypes evaluated in the current study achieve high lint production per boll and lint yields by manipulating different yield drivers. However, lint yield was primarily maximized through an increase in HI due to increases in boll production and within‐boll distribution of biomass to fiber, not due to increases in total biomass production or boll size.more » « less
-
Abstract Size‐based harvest limits or gear regulations are often used to manage fishing mortality and ensure the spawning biomass of females is sufficiently protected. Yet, management interactions with species’ mating systems that affect fishery sustainability and yield are rarely considered. For species with obligate male care, it is possible that size‐specific harvest of males will decrease larval production. In order to examine how size‐based management practices interact with mating systems, we modelled fisheries of two species with obligate care of nests, corkwing wrasse (Symphodus melops, Labridae) and lingcod (Ophiodon elongatus,Hexigrammidae) under two management scenarios, a minimum size limit and a harvest slot limit. We simulated the population dynamics, larval production and yield to the fishery under a range of fishing mortalities. We also modelled size‐dependent male care to determine its interaction with management. In both species, the slot limit decreased yield by <12% (relative to minimum size limits) at low fishing mortalities; at higher mortalities, individuals rarely survived to outgrow the slot and spawning potential decreased substantially relative to unfished levels, similar to minimum size limits. Spawning potential decreased less when managed with a slot limit if we included a positive feedback between male size, care and hatching success, but the benefit of implementing the slot depended both on the relative proportions of each sex selected by the fishery and on our assumptions regarding male size and care. This work highlights that the effects of size‐ and sex‐selective fisheries management can be nuanced and produce counter‐intuitive results.more » « less
-
Abstract Marine Protected Areas (MPAs) are being implemented worldwide, yet there are few cases where managers make specific predictions of the response of previously harvested populations to MPA implementation.Such predictions are needed to evaluate whether MPAs are working as expected, and if not, why. This evaluation is necessary to perform adaptive management, identifying whether and when adjustments to management might be necessary to achieve MPA goals.Using monitoring data and population models, we quantified expected responses of targeted species to MPA implementation and compared them to monitoring data.The model required two factors to explain observed responses in MPAs: (a) pre‐MPA harvest rates, which can vary at local spatial scales, and (b) recruitment variability before and after MPA establishment. Low recruitment years before MPA establishment in our study system drove deviations from expected equilibrium population size distributions and introduced an additional time lag to response detectability.Synthesis and applications. We combined monitoring data and population models to show how (a) harvest rates prior to Marine Protected Area (MPA) implementation, (b) variability in recruitment, and (c) initial population size structure determine whether a response to MPA establishment is detectable. Pre‐MPA harvest rates across MPAs plays a large role in MPA response detectability, demonstrating the importance of measuring this poorly known parameter. While an intuitive expectation is for response detectability to depend on recruitment variability and stochasticity in population trajectories after MPA establishment, we address the overlooked role of recruitment variability before MPA establishment, which alters the size structure at the time of MPA establishment. These factors provide MPA practitioners with reasons whether or not MPAs may lead to responses of targeted species. Our overall approach provides a framework for a critical step of adaptive management.more » « less
