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
Prediction of lake depth across a 17-state region in the United States
Lake depth is an important characteristic for understanding many lake processes, yet it is unknown for the vast majority of lakes globally. Our objective was to develop a model that predicts lake depth using map-derived metrics of lake and terrestrial geomorphic features. Building on previous models that use local topography to predict lake depth, we hypothesized that regional differences in topography, lake shape, or sedimentation processes could lead to region-specific relationships between lake depth and the mapped features. We therefore used a mixed modeling approach that included region-specific model parameters. We built models using lake and map data from LAGOS, which includes 8164 lakes with maximum depth (Zmax) observations. The model was used to predict depth for all lakes ≥4 ha (n = 42 443) in the study extent. Lake surface area and maximum slope in a 100 m buffer were the best predictors of Zmax. Interactions between surface area and topography occurred at both the local and regional scale; surface area had a larger effect in steep terrain, so large lakes embedded in steep terrain were much deeper than those in flat terrain. Despite a large sample size and inclusion of regional variability, model performance (R2 = 0.29, RMSE = 7.1 m) was similar to other published models. The relative error varied by region, however, highlighting the importance of taking a regional approach to lake depth modeling. Additionally, we provide the largest known collection of observed and predicted lake depth values in the United States.
more »
« less
- Award ID(s):
- 1065786
- NSF-PAR ID:
- 10025551
- Date Published:
- Journal Name:
- Inland waters
- Volume:
- 6
- ISSN:
- 2044-2041
- Page Range / eLocation ID:
- 314-324
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Depth regulates many attributes of aquatic ecosystems, but relatively few lakes are measured, and existing datasets are biased toward large lakes. To address this, we used a large dataset of maximum (
Z max;n = 16,831) and mean (Z mean;n = 5,881) depth observations to create new depth models, focusing on lakes < 1,000 ha. We then used the models to characterize patterns in lake basin shape and volume. We included terrain metrics, water temperature and reflectance, polygon attributes, and other predictors in a random forest model. Our final models generally outperformed existing models (Z max; root mean square error [RMSE] = 8.0 m andZ mean; RMSE = 3.0 m). Our models show that lake depth followed a Pareto distribution, with 2.8 orders of magnitude fewer lakes for an order of magnitude increase in depth. In addition, despite orders of magnitude variation in surface area, most size classes had a modal maximum depth of ~ 5 m. Concave (bowl‐shaped) lake basins represented 79% of all lakes, but lakes were more convex (funnel‐shaped) as surface area increased. Across the conterminous United States, 9.8% of all lake water was within the top meter of the water column, and 48% in the top 10 m. Excluding the Laurentian Great Lakes, we estimate the total volume in the conterminous United States is 1,057–1,294 km3, depending on whetherZ maxorZ meanwas modeled. Lake volume also exhibited substantial geographic variation, with high volumes in the upper Midwest, Northeast, and Florida and low volumes in the southwestern United States. -
null (Ed.)Abstract. The Pleistocene sand sea on the Arctic Coastal Plain (ACP) ofnorthern Alaska is underlain by an ancient sand dune field, a geologicalfeature that affects regional lake characteristics. Many of these lakes,which cover approximately 20 % of the Pleistocene sand sea, are relativelydeep (up to 25 m). In addition to the natural importance of ACP sand sealakes for water storage, energy balance, and ecological habitat, the needfor winter water for industrial development and exploration activities makeslakes in this region a valuable resource. However, ACP sand sea lakes havereceived little prior study. Here, we collect in situ bathymetric data totest 12 model variants for predicting sand sea lake depth based on analysisof Landsat-8 Operational Land Imager (OLI) images. Lake depth gradients weremeasured at 17 lakes in midsummer 2017 using a Humminbird 798ci HD SI Comboautomatic sonar system. The field-measured data points were compared tored–green–blue (RGB) bands of a Landsat-8 OLI image acquired on 8 August2016 to select and calibrate the most accurate spectral-depth model for eachstudy lake and map bathymetry. Exponential functions using a simple bandratio (with bands selected based on lake turbidity and bed substrate)yielded the most successful model variants. For each lake, the most accuratemodel explained 81.8 % of the variation in depth, on average. Modeled lakebathymetries were integrated with remotely sensed lake surface area toquantify lake water storage volumes, which ranged from 1.056×10-3 to 57.416×10-3 km3. Due to variations in depthmaxima, substrate, and turbidity between lakes, a regional model iscurrently infeasible, rendering necessary the acquisition of additional insitu data with which to develop a regional model solution. Estimating lakewater volumes using remote sensing will facilitate better management ofexpanding development activities and serve as a baseline by which toevaluate future responses to ongoing and rapid climate change in the Arctic.All sonar depth data and modeled lake bathymetry rasters can be freelyaccessed at https://doi.org/10.18739/A2SN01440 (Simpson and Arp, 2018) andhttps://doi.org/10.18739/A2HT2GC6G (Simpson, 2019), respectively.more » « less
-
Abstract The physical processes of heat exchange between lakes and the surrounding atmosphere are important in simulating and predicting terrestrial surface energy balance. Latent and sensible heat fluxes are the dominant physical process controlling ice growth and decay on the lake surface, as well as having influence on regional climate. While one-dimensional lake models have been used in simulating environmental changes in ice dynamics and water temperature, understanding the seasonal to daily cycles of lake surface energy balance and its relationship to lake thermal properties, atmospheric conditions, and how those are represented in models is still an open area of research. We evaluated a pair of one-dimensional lake models, Freshwater Lake (FLake) and the General Lake Model (GLM), to compare modeled latent and sensible heat fluxes against observed data collected by an eddy covariance tower during a 1-yr period in 2017, using Lake Mendota in Madison, Wisconsin, as our study site. We hypothesized transitional periods of ice cover as a leading source of model uncertainty, and we instead found that the models failed to simulate accurate values for large positive heat fluxes that occurred from late August into late December. Our results ultimately showed that one-dimensional models are effective in simulating sensible heat fluxes but are considerably less sensitive to latent heat fluxes than the observed relationships of latent heat flux to environmental drivers. These results can be used to focus future improvement of these lake models especially if they are to be used for surface boundary conditions in regional numerical weather models. Significance Statement While lakes consist of a small amount of Earth’s surface, they have a large impact on local climate and weather. A large amount of energy is stored in lakes during the spring and summer, and then removed from lakes before winter. The effect is particularly noticeable in high latitudes, when the seasonal temperature difference is larger. Modeling this lake energy exchange is important for weather models and measuring this energy exchange is challenging. Here we compare modeled and observed energy exchange, and we show there are large amounts of energy exchange happening in the fall, which models struggle to capture well. During periods of partial ice coverage in early winter, lake behavior can change rapidly.more » « less
-
Abstract Global change may contribute to ecological changes in high-elevation lakes and reservoirs, but a lack of data makes it difficult to evaluate spatiotemporal patterns. Remote sensing imagery can provide more complete records to evaluate whether consistent changes across a broad geographic region are occurring. We used Landsat surface reflectance data to evaluate spatial patterns of contemporary lake color (2010–2020) in 940 lakes in the U.S. Rocky Mountains, a historically understudied area for lake water quality. Intuitively, we found that most of the lakes in the region are blue (66%) and were found in steep-sided watersheds (>22.5°) or alternatively were relatively deep (>4.5 m) with mean annual air temperature (MAAT) <4.5°C. Most green/brown lakes were found in relatively shallow sloped watersheds with MAAT ⩾4.5°C. We extended the analysis of contemporary lake color to evaluate changes in color from 1984 to 2020 for a subset of lakes with the most complete time series ( n = 527). We found limited evidence of lakes shifting from blue to green states, but rather, 55% of the lakes had no trend in lake color. Surprisingly, where lake color was changing, 32% of lakes were trending toward bluer wavelengths, and only 13% shifted toward greener wavelengths. Lakes and reservoirs with the most substantial shifts toward blue wavelengths tended to be in urbanized, human population centers at relatively lower elevations. In contrast, lakes that shifted to greener wavelengths did not relate clearly to any lake or landscape features that we evaluated, though declining winter precipitation and warming summer and fall temperatures may play a role in some systems. Collectively, these results suggest that the interactions between local landscape factors and broader climatic changes can result in heterogeneous, context-dependent changes in lake color.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5268/IW-6.3.957
