Page 1
RESPONSES OF FISH AND ZOOPLANKTON TO CLIMATE VARIATION ON
THE PRAIRIES, AND THEIR SENSITIVITY TO CLIMATE CHANGE
A Thesis
Submitted to the Faculty of Graduate Studies and Research
In Partial Fulfillment of the Requirements
For the Degree of
Master of Science
in
Biology
University of Regina
By
Elizabeth Raye Starks
Regina, Saskatchewan
December, 2012
©2012, Elizabeth Starks
Page 2
UNIVERSITY OF REGINA
FACULTY OF GRADUATE STUDIES AND RESEARCH
SUPERVISORY AND EXAMINING COMMITTEE
Elizabeth Raye Starks, candidate for the degree of Master of Science in Biology, has presented a thesis titled, Responses of Fish and Zooplankton to Climate Variation on the Prairies, and Their Sensitivity to Climate Change, in an oral examination held on December 14, 2012. The following committee members have found the thesis acceptable in form and content, and that the candidate demonstrated satisfactory knowledge of the subject material. External Examiner: Dr. Jeff Hudson, University of Saskatchewan
Supervisor: Dr. Björn Wissel, Department of Biology
Committee Member: Dr. Peter R. Leavitt, Department of Biology
Committee Member: Dr. Christopher Somers, Department of Biology
Chair of Defense: Dr. Kyle Hodder, Department of Geography
Page 3
i
ABSTRACT
Climate change is anticipated to significantly increase temperatures and alter
current rainfall patterns, which will have important ramifications for aquatic habitats and
their biological communities. Current observations indicate that climate effects will vary
depending on region and lake type, and some lacustrine areas, such as the Great Plains,
are particularly sensitive to these effects. Variations in local climate and lake
morphometry create different habitats, which each have unique environmental controls.
The anticipated impacts of climate change on aquatic biota can be difficult to evaluate
because of potentially contrasting effects of temperature and hydrology on lake
ecosystems, particularly in closed-basin lakes within semi-arid regions.
To address these challenges, I quantified decade-scale changes in chemical and
biological properties of 20 endorheic lakes in central North America in response to a
pronounced transition from a drought to a pluvial period during the early 21st century.
Lakes exhibited marked changes in chemical characteristics and formed two discrete
clusters corresponding to periods of substantially differing effective moisture (as per
Palmer Drought Severity Index, PDSI). Discriminant function analysis (DFA) explained
90% of variability in fish assemblage composition and showed that fish communities
were predicted best by environmental conditions during the arid interval (PDSI < -2).
DFA also predicted that lakes could support more fish species during pluvial periods, but
their realized occurrences may be limited by periodic stress due to recurrent droughts and
physical barriers to colonization. Zooplankton taxonomic compositions in fishless lakes
were resilient to short-term changes in meteorological conditions, and did not vary
between drought and deluge periods. Conversely, zooplankton taxa that were exposed to
Page 4
ii
fish decreased substantially in biomass during the wet interval, likely due to increased
zooplanktivory by fish.
Based on my results, climate change is expected to alter fish species distributions,
but it is less clear to what extent non-lethal environmental effects will influence physical
health of populations in fish-habitable lakes. To address this question, I investigated the
environmental controls of body condition and parasite load in walleye (Sander vitreus),
northern pike (Esox lucius) and yellow perch (Perca flavescens) in seven lakes from the
prior study. Over a two-year observation period (2009 vs. 2010), I observed large
differences in the number of days within the favorable temperature range for ambient fish
species. Surprisingly, environmental variables such as lake morphometry and nutrient
levels had little relevance, despite their importance in previous studies conducted in
boreal lakes. Instead, temperature and salinity were important correlates of fish health.
In regard to species-specific effects, walleye was most sensitive to interannual
temperature differences, as well as salinity, while yellow perch and northern pike
exhibited temperature sensitivity to a lesser degree. Apparently, temperature increases
are of particular concern in prairie lakes, as their polymictic nature deprives fishes of a
hypolimnetic thermal refuge.
Together these findings suggest that semi-arid lakes provide a useful model
system for anticipating the effects of global climate change on aquatic communities in
closed-basin lakes of semi-arid regions. The particular importance of temperature and
salinity indicates that the interaction of global climate change and local hydrology may
have particularly detrimental effects not only on the health but also the survival of
established fish populations of the Great Plains.
Page 5
iii
ACKNOWLEDGMENTS
This work was made possible through the contributions and support of several
individuals. In particular I thank my supervisor Dr. Björn Wissel for welcoming me to
Saskatchewan and providing his mentorship, funding support, guidance and expertise. I
also thank the members of my academic committee, Drs. Peter Leavitt and Chris Somers,
for their constructive feedback and supervision throughout my graduate program. Dr.
Chris Wilson contributed knowledgeable proofreading and training. Several individuals
contributed to this research effort in the laboratory or the field, including Samantha Pham,
Sydney Chow, David Scott, Matt Bogard, Dr. Zoraida Quiñones-Rivera, Peter Dowdy,
David Braun and John Kalyn. This research was funded by NSERC, the Saskatchewan
Ministry of Environment, the Saskatchewan Fish and Wildlife Development Fund, the
Prairie Adaptation Research Collaborative, and the University of Regina Faculty of
Graduate Studies and Research. Finally, I’d like to thank my spouse Peter Dowdy,
Denise Brooks, my parents and siblings for their love and support.
Page 6
iv
TABLE OF CONTENTS
Chapter Page
ABSTRACT…………………………………………………………………………….…i
ACKNOWLEDGMENTS……………………………………………………………..…iii
TABLE OF CONTENTS………………………………………………………..….…….iv
LIST OF TABLES………………………………………………………………….……vii
LIST OF FIGURES………………………………………………………………………ix
CHAPTER 1: INTRODUCTION………………………………………………………...1
1.1.Climate Effects on Lakes………………………………………………...……1
1.2.Prairie Lake Climate Response…………………………………………..……2
1.3.Prior Research on Prairie Biota……………………………………………….3
1.4.Research Needs………………………………………………………………..4
1.5.Objectives and Relevance………………………………………………..……6
CHAPTER 2: EFFECTS OF DROUGHT AND PLUVIAL PERIODS ON FISH AND
ZOOPLANKTON COMMUNITIES IN PRAIRIE LAKES: DETERMINISTIC VS.
STOCHASTIC RESPONSES……………………………………………………….……8
2.1.INTRODUCTION…………………………………………………………....8
2.2.METHODS……………………………………………………………..……11
2.2.1. Study Area…………………………………………………….…11
2.2.2. Lake Sampling……………………………………………….…..14
Page 7
v
2.2.3. Laboratory Analysis……………………………………………...15
2.2.4. Statistical Analyses………………………………………………16
2.3.RESULTS……………………………………………………………………18
2.3.1. Climate………………………………………………...…………18
2.3.2. Zooplankton…………………………………………………...…23
2.3.3. Fishes……………………………………………………….……23
2.4.DISCUSSION……………………………………………………………….28
2.4.1. Climate………………………………………………………..…28
2.4.2. Water Chemistry…………………………………………………28
2.4.3. Zooplankton…………………………………………………..…29
2.4.4. Fishes………………………………………………………….…31
2.4.5. System Responses……………………………………………..…32
2.4.6. Conclusions………………………………………………………33
CHAPTER 3. CONDITION AND PARASITE LOAD OF FISHES IN RESPONSE TO
CLIMATE-INDUCED CHANGES IN TEMPERATURE AND SALINITY……….…34
3.1.INTRODUCTION………………………………………………………..…34
3.2.METHODS………………………………………………………………..…38
3.2.1. Study Sites and Field Sampling……………………………….…38
3.2.2. Statistical Analysis…………………………………………….…41
3.3.RESULTS……………………………………………………………………43
3.3.1. Climate………………………………………………………...…43
3.3.2. Body Condition…………………………………………..………43
3.3.3. Stomach Contents……………………………………………..…49
Page 8
vi
3.3.4. Parasite Load…………………………………………………..…49
3.4.DISCUSSION……………………………………………………………..…53
3.4.1. Body Condition………………………………………………..…53
3.4.2. Parasite Load…………………………………………………..…55
3.4.3. Species-specific Differences…………………………………..…56
3.4.4. Anticipated Future Effects of Climate Change………………..…56
CHAPTER 4. CONCLUSION………………………………………………………..…58
4.1.SYNTHESIS……………………………………………………………..…..58
4.2.FUTURE RESEARCH DIRECTIONS…………………………………...…59
4.3.GENERAL CONCLUSIONS……………………………………………..…61
REFERENCES………………………………………………………………………..…63
Page 9
vii
LIST OF TABLES
Table Page
Table 2.1. Characteristics of the 20 study lakes utilized for the survey from 2002 to 2010.
Ecoregions are abbreviated as MG (mixed grassland), AP (aspen parkland) and MMG
(moist mixed grassland). Fish communities are abbreviated as Fishl. (fishless), Plankt.
(planktivorous) and Pisc. (piscivorous). Zooplankton species clusters are identified as 1
through 4 (see Figure 4 for details). Maximum lake depth (m), salinity (g L-1
), DOC (mg
L-1
) and TKN (µg L-1
) are calculated averages across the 9-year study period by lake…13
Table 2.2. Comparison of average water chemistry parameters for fishless, planktivorous
and piscivorous lake communities. Water chemistry records shown here include
maximum lake depth (Zmax, m), salinity (g L-1
), TKN (µg L-1
), TDP (µg L-1
) and DOC
(mg L-1
). For each community type, means are reported for the drought and fluvial
periods and mid-range values are reported for the drought-based DFA model……….…20
Table 3.1. Environmental characteristics of the seven lakes included in the study from
2009 to 2010. Gamefish recorded include walleye (W), northern pike (N) and yellow
perch (Y). Recorded values of salinity (as total dissolved solids, TDS), Total Kjeldahl
Nitrogen (TKN), Total Phosphorous (TDP) and Dissolved Organic Carbon (DOC) are
means from June and August of 2009 and 2010.Table 3.2. Climate conditions from 2008
to 2010…………………………………………………………………………………...39
Table 3.2. Climate conditions from 2008 to 2010. Temperatures are recorded from the
Saskatoon weather station and ice cover data were taken from the Buffalo Pound water
treatment facility. Ice duration data are recorded for the entire winter leading up to the
Page 10
viii
warm season of the indicated year. Mean highs and lows are recorded for the full ice off
period………………………………………………………………………………….…44
Table 3.3. Wr for walleye, pike and perch over the full study period from 2009 to 2010.
NA = sample size not adequate for calculating an average…………………………...…45
Page 11
ix
LIST OF FIGURES
Figure Page
Figure 2.1. Locations of the 20 study lakes in southern Saskatchewan, Canada in the
northern Great Plains. Sampling sites are coded by fish community type
(piscivorous=black, planktivorous=gray, fishless=white). Dashed lines represent annual
precipitation deficits. The solid line represents the transition from southern grassland to
northern forest……............................................................................................................12
Figure 2.2. Palmer Drought Severity Index (PDSI) for Saskatoon in August of 2000-2010.
Negative values indicate effective precipitation deficit, while positive values indicate
excess………………………………………………………………………………….…19
Figure 2.3. (a) PCA of water chemistry in the 20 study lakes. Labeled gray arrows
represent water chemistry variables used to generate the ordination space (TDS, TDP,
TKN, DOC, Zmax, Secchi depth). Each black arrow represents an individual lake,
originating at environmental conditions of this lake during the drought and extending to
conditions during the deluge. Generally, directional changes trend toward increased depth
and water clarity and reduced nutrient load. (b) Circular plot identifying directional
changes in water chemistry from drought to deluge. Labeled arrows represent the
superimposed PCA water chemistry variables. Each unlabeled arrow represents the
direction and magnitude of change in water chemistry for one of the 20 study lakes. With
some exceptions, most lakes increased in depth and water clarity and decreased in
nutrient load and temperature. The directional change was statistically significant
(p<0.05)…….......22
Page 12
x
Figure 2.4. Principal Components Analysis of zooplankton taxonomic composition
during drought- and deluge-influenced periods. Labeled lines represent taxa used to
generate the ordination space. Vectors of individual lakes within taxonomic space are
plotted between drought and deluge periods, with arrows terminating in the deluge period.
Of the four zooplankton assemblages, the most change occurred in the piscivores group
with small cladocerans in the lower right segment (Ceriodaphnia, Simocephalus,
Diaphanosoma)…………………………………………………………………………..24
Figure 2.5. Discriminant Function Analysis of the 20 study lakes using drought-
influenced water chemistry data (a) and deluge-influenced chemistry data (b). The
discriminant function classified lakes as fishless (white), planktivorous (gray) or
piscivorous (black) using conductivity, dissolved organic carbon (DOC), depth (Zmax) and
nutrient load (TDP and TKN) as variables. The components of the discriminant functions
are listed on each axis in order of contribution. During the drought and deluge, nitrogen
and conductivity were the most important factors on axis 1. Axis 2 was not significant
for either scenario. Classifications based on the drought model were 90% accurate, with
a high level of distinction between all community types. Deluge-based classifications
were only 60% accurate, with most misclassifications arising due to overlap in the
chemistry of planktivorous and piscivorous lakes……………………………………….26
Figure 2.6. Principal Components Analysis of the 20 study lakes during the drought- and
deluge-influenced periods (a and b, respectively), coded by level of winterkill risk
(Barica and Mathias 1979). White indicates low risk, gray indicates moderate risk and
Page 13
xi
dark gray indicates high risk. Moderate to high winterkill risk affected 9 lakes during the
drought period, which was reduced to 4 lakes during the deluge. Risk is negatively
associated with depth and positively associated with nutrient load……………………...27
Figure 3.1. Location of the seven study lakes (filled circles) in the grassland ecoregion of
southern Saskatchewan, Canada. Ice dates were recorded from Buffalo Pound Lake
(filled square), and air temperature data were recorded from Saskatoon………………..40
Figure 3.2. Regressions of length by weight revealed a significant decrease in walleye
body weight by length from June to August 2009 (above), and an increase in 2010
(below)…………………………………………………………………………………...47
Figure 3.3. Regressions of length by weight revealed no significant change in perch body
weight for the full lake set from early to late 2009, and an increase in weight from early
to late 2010……………………………………………………………………………….48
Figure 3.4. Female walleye (gray bar) exhibited a significantly higher mean tapeworm
count than males (hatched bar) in both years. Overall tapeworm counts were
significantly higher in 2009 than 2010. No males were infected in 2010………………51
Figure 3.5. Proportion of walleye populations infected with tapeworms was significantly
associated with salinity in 2009 (r2 = 0.934), but not in 2010. Shannon Lake, excluded
from the analysis, is shown in white……………………………………………………..52
Page 14
1
1. INTRODUCTION
1.1. Climate Effects On Lakes
The IPCC has stated that climate change effects on freshwater resources are a
critical concern (Bates et al. 2008). Changes in climate are expected to include regional
shifts in temperature, precipitation and frequency of extreme weather events, such as
drought or flooding, all of which have ramifications for sustainability of freshwater
resources (Gleick 1989, Millet et al. 2009). Climate change research on lake ecosystems
is broad in scope, including assessments of hydrology (Gleick 1989), water chemistry
(Hammer 1990, Jeppesen et al. 2009), and biota (Chu et al. 2005, O’Reilly et al. 2003,
Winder and Schindler 2004).
Climate impacts on lake ecosystems occur primarily by changing inputs of solar
energy and mass (water and dissolved/particulate substances) (Leavitt et al. 2009). Solar
energy enhances regional evapotranspiration rates, not only reducing water volume but
also concentrating solutes. Mass inputs are usually more localized and can come from
tributaries, surface runoff, seasonal meltwater and groundwater, at rates which are
determined by precipitation, seasonal thermal variation and local hydrology.
While most lacustrine regions are expected to experience warming and reduced
precipitation (Adrian et al. 2009, Bates et al. 2008), research indicates that climate
responses will differ depending on ecoregion and lake type (Williamson et al. 2009). As
a consequence, small boreal lakes are expected to stratify earlier in the season and with
greater intensity (De Stasio et al. 1996), while lakes in areas with substantial groundwater
contributions might show asynchronous responses (Webster et al. 2000), and saline lakes
could be at risk of complete desiccation (Williams 2002). Of all lake types, endorheic
Page 15
2
prairie lakes are considered among the most vulnerable to climate change due to their
closed hydrology and high dependency on winter precipitation (Hammer 1990).
1.2. Prairie Lake Climate Response
The regional climate of the northern Great Plains is extremely seasonal,
experiencing a thermal range of up to 85˚C throughout the year, ranging from -45˚C to
40˚C (Hammer et al. 1978). Lakes are frozen for approximately 155 days per year due to
long and extreme winters (Buffalo Pound Water Treatment Plant report). The grassland
region of southern Saskatchewan is classified as sub-humid to semi-arid, based on
precipitation deficits. Yet, precipitation rates are highly variable in time and space and
can result in a range from dry to humid conditions (Sauchyn et al. 2002).
For endorheic lakes, seasonal meltwater and groundwater are key contributors to
lake water budgets. Nevertheless, seasonal meltwater contributions vary substantially
from year to year depending on winter snowfall, and a year with low winter precipitation
can result in drought-like conditions for lakes during the subsequent summer (McGowan
et al. 2005). Groundwater structure on the prairies is complex and spatially variable, and
depending on climate and the state of the aquifer, direction of ground water flow shift
between sink and source in respect to lake water budgets (Winter and Rosenberry 1998).
Concentrations of dissolved substances in lake waters are closely linked to
climate either through evaporative concentration in arid periods, or dilution in periods of
high precipitation. As a consequence of evaporative deficits and abundance of closed-
basin lakes in this region, many lakes have elevated salinities and nutrient concentrations,
which could be intensified by saline groundwater contributions (Last and Ginn 2005).
Additionally, the naturally high concentrations of phosphorus and nitrogen in the region
Page 16
3
are augmented by fertilizer run-off, such that most lakes are classified as eutrophic or
hypereutrophic (Biehuizen and Prepas 1985). However, saline lake nutrients are often
less bioavailable (Waiser and Robarts 1995).
1.3. Prior Research On Prairie Biota
Aquatic communities in prairie lakes are diverse in species richness and trophic
complexity, supporting between two and four trophic levels and ranging from simple
invertebrate communities to diverse communities of fishes, zooplankton and benthic
invertebrates (Cooper and Wissel 2012b, Hammer 1978). Although trophic state and
salinity levels are higher in southern Saskatchewan lakes than many other lake ecoregions,
Saskatchewan has historically supported a variety of fish populations, with 26 species
identified in a 1944 survey (Rawson and Moore 1944).
In respect to environmental controls of community structure, salinity has been
identified as a master variable in prairie lakes, reducing species richness and
preferentially selecting halophiles (Hammer 1978, Rawson and Moore 1944). Among
major taxonomic groups, fishes were identified as most sensitive to salinity, followed by
zooplankton and finally benthic invertebrates (Hammer 1978). Surprisingly, recent
studies showed that most fishes are now eliminated from lakes at salinities well below
their physiological thresholds, indicating that controls other than salinity may be
important for structuring these lake communities (Cooper and Wissel 2012b). Seasonal
or spatial oxygen depletion is commonly responsible for loss of fishes and other sensitive
organisms, and is likely an important community structuring mechanism in these lakes
(Robarts et al. 2005). Due to the polymictic nature of many prairie lakes, summer
hypoxia is usually not a major concern except in meromictic lakes. In contrast, oxygen
Page 17
4
depletion in winter could pose a serious threat to established fish populations on the
prairies as low water depth and high eutrophic state can significantly deplete oxygen
reserves under ice during the prolonged winter (Robarts et al. 2005). Climate change
scenarios predict a reduction in winter hypoxia as winter temperatures rise, although
summer temperatures may rise beyond comfortable limits for fishes in some regions,
leaving the net effect of a temperature increase somewhat ambiguous (Fang et al. 2004).
Such lethal effects clearly determine the presence and absence of species in
particular ecosystems. Nevertheless, environmental parameters can also impact aquatic
species at a sub-lethal level. For example, lake depth is often positively associated with
body condition for larger fishes (Cena et al. 2006), and a high rate of primary
productivity reduces foraging efficiency of visual predators (Craig and Babaluk 1989),
while optimal temperature ranges are associated with better growth rates (Quist et al.
2002). Additionally, environmental stressors can reduce body condition and increase
parasite load by weakening the immune system (Bly et al. 1997, Lewis et al. 2003).
1.4. Research Needs
Climate and other environmental conditions in prairie lakes are distinct from
better studied boreal lakes. So far, very few studies have tested the impacts of climate and
water quality on fish populations in northern polymictic lakes. Effects of temperature
and water chemistry on these populations are as yet not fully characterized. Shallow
temperate reservoirs are the nearest analogue for temperature effects, with studies
showing a large effect of temperature changes on fishes in these systems due to the lack
of thermal stratification (Quist et al. 2002). It is currently assumed that high temperatures
do not exceed critical thresholds for northern fish populations, but based on projected
Page 18
5
temperature increases, this may no longer be the case in the future (Sharma et al. 2007).
Fish populations in these lakes currently experience water chemistry that is unusual (high
nutrient content and high levels of uncommon ions such as Mg+ and SO4-) relative to
other regions (Rawson and Moore 1944), and yet the sub-lethal effects of the unique
chemical properties have not been fully characterized (Pollock et al. 2010).
Climate change scenarios for moderate climates often predict increases in
seasonal variability, as well as increases in ice-free time (Bates et al. 2008). In contrast,
the climate in central Canada is already considered extreme, with very pronounced
seasonal shifts in temperature, a brief growing season and high interannual variation in
precipitation (Hammer et al. 1978). This area also experiences cyclical drought and
deluge events that may further construct species assemblages (Michels et al. 2007). Due
to the extreme variability in short-term (seasonal, annual) and long-term (decadal,
centennial) weather this region experiences, it is highly suited for climate research, and
the periodic droughts and deluges experienced by Saskatchewan could serve well as a
proxy for climate change conditions.
Few studies have had the ability to evaluate environmental and community effects
simultaneously in order to get a more comprehensive view of the biological response to
climate change. Studies of stratified lakes imply that zooplankton sensitivity to
environmental change, particularly shifts in thermal stratification and timing of seasonal
events, could reduce their overlap with predators (Hampton et al. 2006, Winder and
Schindler 2004). However, in the polymictic lakes of Saskatchewan such habitat overlap
effects would not occur, so it is less clear how predator/prey interactions would be altered
in a climate-change scenario. Due to polymixis, it is possible that community outcomes
Page 19
6
and any resulting shifts in relative abundance will be dependent on the abilities of each
species to thrive under future conditions.
1.5. Objectives and Relevance
The primary goal of the second chapter of this study was to determine how
climate variability, particularly during a transition between deluge and drought, structures
aquatic communities in prairie lakes. Because climate may also have indirect effects
through altering species interactions, environmental and community effects were
evaluated simultaneously using zooplankton species abundances and fish habitat
suitability models. These interactions were evaluated within the context of
environmental parameters of known relevance to aquatic species composition, such as
lake morphometry and nutrient loads.
The third chapter addressed non-lethal environmental effects on game fish. Non-
lethal effects of environment and climate on fish are of importance due to their relevance
to total fisheries yields and long-term population sustainability, and are therefore an
important consideration beyond habitability. Because Saskatchewan lakes have such
variable water chemistry, I compared the importance of inter-lake chemical variation
versus a two-year difference in seasonal temperature that was similar across the study
lakes. All parameters identified as important to species composition were evaluated as
potential sublethal stressors in order to determine whether these parameters were relevant
to physical health, and could have detrimental effects increasing with the degree of
environmental change.
Finally, I evaluated potential global climate change effects on local biota using
what was learned from chapters 2 and 3. Projected reductions in winter precipitation and
Page 20
7
increasing temperatures could induce conditions similar to regional droughts, with
extreme drought conditions acting as a habitat filter (Chase 2007), while anticipated lake
surface-temperature increases would surpass thermal thresholds for many cool- and cold-
water fishes (Sharma et al. 2007). This study therefore demonstrates that the inherent
climate variability of this region makes it a useful model system for climate change
projections.
Page 21
8
2. EFFECTS OF DROUGHT AND PLUVIAL PERIODS ON FISH AND
ZOOPLANKTON COMMUNITIES IN PRAIRIE LAKES: DETERMINISTIC VS.
STOCHASTIC RESPONSES
2.1. INTRODUCTION
The Intergovernmental Panel on Climate Change (IPCC) has identified arid
regions and endorheic basins as particularly vulnerable to climate change effects due to
low water availability and high temporal variation in precipitation and evaporation (Bates
et al. 2008). Lakes in these areas do not receive continuous inputs through streams or
rivers and frequently lack surface outlets. Across the northern Great Plains, endorheic
lakes are sustained by spring snow melt, while evaporation during summer represents the
predominant loss of water from lakes (McGowan et al. 2005, Pham et al. 2008).
Consequently, these lakes are susceptible to partial or total desiccation during droughts,
particularly those lasting multiple years. Conversely, the low relief of the landscape can
result in significant flooding during spring snow melt or intense summer deluges (Winter
and Rosenberry 1998), resulting in dilution of dissolved substances (e.g., major ions,
nutrients, dissolved carbon). These strong effects of hydrologic variability on
limnological conditions is further intensified as many lakes in central North America are
shallow, exhibit high surface to depth ratios, and are often polymictic. In particular, the
lack of continuous thermal stratification during summer can increase vulnerability of the
biota to rising water temperatures as cool hypolimnetic refuges are absent. Ultimately,
changes in hydrology and the associated changes in nutrient concentrations and salinity
have the potential to strongly alter biological diversity, taxonomic composition and
Page 22
9
trophic interactions within aquatic food webs (Cooper and Wissel 2012a, Hammer 1990,
Wissel et al. 2011).
Endorheic lakes on the semi-arid Great Plains represent an important model to
study the long-term effects of climate on lake ecosystems (Pham et al 2008, 2009, Wissel
et al. 2011). At present, global circulation models predict a uniform increase of about
4 °C in all seasons of 2050 across the Canadian prairies (Barrow 2009), further raising
the probability of drought conditions due to intensified evaporation. In conjunction,
extreme precipitation events are expected to become more common as elevated
atmospheric energy in summer favors increased evapotranspiration rates and atmospheric
vapor content (Bates et al. 2008), while warmer winters may reduce snowpack and
subsequent runoff during spring melt periods. In both cases, the sustainability of prairie
lakes may be compromised by rapid changes in lake level during severe summer droughts
(Van der Kamp et al. 2008) or pluvial intervals (Winter and Rosenberry 1998).
Furthermore future climate change may exhibit high geographic variability, with
increased precipitation in the southeastern prairie of Minnesota and the Dakotas, and
increased temperature combined with static or declining precipitation in the northwestern
prairies including southern Alberta and Saskatchewan (Millet et al. 2009). If such trends
continue, the prairie region could diverge into two distinct ecozones, with central regions
in Canada experiencing a more arid climate. Such changes may have profound socio-
economic implications for central North America, a region which encompassed a high
proportion of continental cropland (e.g., 40% of arable farmland in Canada).
During drought periods when energy inputs (radiation) dominate, a higher
synchrony of water chemistry is expected across lakes, while wet periods with high mass
Page 23
10
inputs (hydrological flux) result in asynchrony of lake parameters (Leavitt et al. 2009,
Pham et al. 2009). A synchronous response of water chemistry during drought could
shift environmental parameters in a common direction that restricts the diversity of
species assemblages. In instances where species composition is controlled primarily
through deterministic mechanisms (environmental filters), community similarity should
increase during droughts, as inferred from analysis of zooplankton communities in
experimental ponds (Chase 2007), and plankton in hydrologically connected prairie lakes
(Vogt et al. 2011). Conversely, if water chemistry responds asynchronously to a
prolonged deluge, increased differences in water chemistry reduce the effect of
environmental filters and favor elevated biological diversity. Unfortunately, the lack of
long-term monitoring of endorheic lakes combined with limited ranges of observed
meteorological conditions has often prevented explicit tests of these hypotheses.
In this paper, I describe a recent decade-scale shift from drought conditions to an
above-average pluvial interval that provided a unique opportunity to study the effects of
extreme climate on water chemistry and biota of this vulnerable prairie region. During
2000-2003, central Canada experienced a severe multi-year drought (Sauchyn et al. 2005)
which was followed immediately by several years of unusually high rain fall in 2005-
2010 (Van der Kamp et al. 2008). Effects of these hydrological extremes on lake
ecosystems were evaluated by comparing lake properties and biological communities
monitored during 2002-2005 with those recorded during 2007-2010.
Previously, it has been shown that landscape differences in taxonomic
compositions in endorheic lakes are regulated by pronounced spatial variation in water
chemistry and lake morphology, whereas short-term meteorological variability (2-3 yrs.)
Page 24
11
had little significant effects on food-web composition (Cooper and Wissel 2012a, Wissel
et al. 2011). Here I take advantage of the wide range of chemical and morphological
conditions associated with decadal scale shifts in effective moisture to quantify the
effects of climate change on fish and zooplankton communities within endorheic lakes
and to evaluate if the observed community changes are associated with deterministic
mechanisms or are stochastic. The results of this study are important not only to better
understand the responses of prairie lake communities to climate change, but should also
be applicable to lakes in other arid regions world-wide (Hughes 2003).
2.2. METHODS
2.2.1. Study Area
Meteorological conditions in southern Saskatchewan are extremely variable.
Annual temperatures vary between extremes of -40 ºC and +30 ºC, and precipitation
averages 350 mm year-1
(Saskatoon) but ranges from < 250 to > 600 mm year-1
. Relative
humidity rarely exceeds 50% during the growing season, which together with the high
number of sunshine hours (> 2000) and low precipitation, results in evaporative deficits
in most years (Last and Ginn 2005). Sensitivity to evaporative deficits is commonly
expressed as Palmer Drought Severity Index (PDSI), accounting for precipitation and
temperature to determine effective precipitation. PDSI values are centered on 0 under
conditions of normal effective precipitation, while negative or positive values represent
dry or wet periods, respectively. PDSI values below -2 are considered droughts.
To evaluate the impacts of climate on lake ecosystems, 20 study lakes were sampled over
a nine-year observation period (Table 2.1). The study lakes were located in the prairie
ecozone of southern Saskatchewan in central Canada (Fig. 2.1), and were sampled
Page 25
12
Table 2.1. Characteristics of the 20 study lakes utilized for the survey from 2002 to 2010.
Ecoregions are abbreviated as MG (mixed grassland), AP (aspen parkland) and MMG
(moist mixed grassland). Fish communities are abbreviated as Fishl. (fishless), Plankt.
(planktivorous) and Pisc. (piscivorous). Zooplankton species clusters are identified as 1
through 4 (see Figure 2.4 for details). Maximum lake depth (m), salinity (g L-1
), DOC
(mg L-1
) and TKN (µg L-1
) are calculated averages across the 9-year study period by lake.
Lake
Ecoregion
Salinity
TKN
Depth
DOC
Fish
Zooplankton
Charron AP 6 2350 7 92 Fishl. 2
Success MG 22 2837 16 65 Fishl. 4
Arthur AP 12 3414 5 68 Fishl. 3
L. Manitou MMG 43 7106 5 142 Fishl. 4
Antelope MG 13 4498 5 92 Fishl. 3
Rabbit AP 5 2174 4 32 Fishl. 3
Snakehole MG 64 11357 3 372 Fishl. 4
Middle AP 23 6182 5 79 Fishl./Plankt 4
Edouard AP 0.3 1511 6 21 Plankt. 2
Deadmoose AP 2.0 1199 13 31 Plankt. 3
Redberry AP 11 2081 13 39 Plankt. 3
Waldsea AP 16 1688 12 61 Plankt. 3
Clair AP 1.4 1481 4 29 Plankt. 2
Pelletier MG 0.4 713 9 16 Pisc. 1
Wakaw AP 1.9 1023 9 18 Pisc. 1
Shannon AP 2.6 1582 8 21 Pisc. 1
Humboldt AP 1.3 1875 7 29 Pisc. 1
Kipabiskau AP 0.4 1167 8 31 Pisc. 1
Fishing AP 2.0 1199 13 31 Pisc. 1
Lenore AP 1.3 870 9 33 Pisc. 1
Page 26
13
Fig. 2.1: Locations of the 20 study lakes in southern Saskatchewan, Canada in the
northern Great Plains. Sampling sites are coded by fish community type
(piscivorous=black, planktivorous=gray, fishless=white). Dashed lines represent annual
precipitation deficits. The solid line represents the transition from southern grassland to
northern forest.
Page 27
14
over a nine-year observation period (Table 2.1). Four lakes were situated in the semi-arid
southwestern mixed grasslands ecoregion and the remaining sixteen existed in the sub-
humid mixed grasslands/aspen parkland ecoregions of central Saskatchewan. The lakes
have long retention times, exceeding ten years when known, (Pham et al. 2009), and most
lakes are classified as endorheic. Salinity (as total dissolved solids, TDS) ranged from
fresh to hypersaline, maximum lake depth was between 3 m and 30 m, and trophic status
covered mesotrophic to hypereutrophic conditions (Wissel et al. 2011). Of the 20 lakes,
seven had piscivorous, planktivorous and benthivorous fishes (hereafter referred to as
piscivorous), six had only planktivorous and benthivorous fishes (planktivorous), and
seven were fishless (fishless). Further analyses of the lakes were based on these three
fish community classes.
2.2.2. Lake Sampling
The 20 lakes were sampled at central locations over their maximum depth in May,
July and September of 2002 to 2005, and in June and August of 2007 to 2010. Lakes
were not sampled in 2006. I relied on annual averages for all subsequent analyses, as
monthly sampling of a sub-set of eight lakes (between May and September from 2008 to
2010) identified no significant differences between the two sampling periods (analysis
not shown).
Vertical profiles of temperature, oxygen and salinity were recorded with a YSI
model 556 probe at 1 meter depth intervals. Changes in lake depth were recorded at
established GPS coordinates with an electronic depth finder (model), and water clarity
was measured using a 20 cm diameter Secchi disk. Small and large pelagic zooplankton
were collected using 30 cm and 50 cm diameter plankton nets with mesh sizes of 80 µm
Page 28
15
and 500 µm, respectively, towed vertically to the surface from 0.5 m above the sediment
or above the anoxic monimolimnion in meromictic lakes. Zooplankton samples were
stored in 50% ethanol sucrose, identified to species (Hanley et al. 2010, Hudson and
Lesko 2003) and enumerated using a dissecting microscope. Pelagic fish were collected
using overnight sets of multi-panel 1.5x30 meter gillnets with 55 mm mesh size. Littoral
fishes were caught during daylight using a beach seine (2 m x 30 m, 10 mm mesh) from
shore to a depth of two meters (Cooper and Wissel 2012a). All fish were identified to
species (Scott and Crossman 1973) and preserved frozen for further analysis.
2.2.3. Laboratory Analysis
Nutrient analyses were conducted on integrated water samples taken with a 5 cm
diameter tygon tubing sampler deployed from the surface to the bottom in shallow
polymictic lakes, from 0 to 6 meters in deeper polymictic lakes, or to the depth of the
monimolimnion in meromictic lakes. Water samples were prescreened in 80-μm mesh
then filtered through sterile 0.45 μm glass fiber filters and analyzed for total dissolved
phosphorus (TDP), total Kjeldahl nitrogen (TKN), nitrate (NO3), nitrite (NO2),
ammonium (NH4), dissolved organic carbon (DOC) dissolved inorganic carbon (DIC)
and chlorophyll α (Chl α) using standard methods described by Cooper and Wissel
(2012a). Nutrient analyses were conducted at the University of Alberta water chemistry
laboratory following the procedures of Stainton et al. (1977). DOC and DIC were
quantified using a Shimadzu TOC5000 analyzer at the Environmental Quality Analysis
Laboratory (EQAL) at the University of Regina. To estimate algal biomass (Chl α), an
aliquot of 1 liter was filtered onto GF/C filters and stored at -10 °C until extraction with
Page 29
16
acetone:methanol:water (80:15:5), and analyzed with a spectrophotometer (Hewlett-
Packard model 84452A) (Wetzel and Likens 1991).
2.2.4. Statistical Analyses
To test if the chemical characterization of lakes by year would be associated with
drought and pluvial-influenced periods I used a Multi-Response Permutation Procedure
(MRPP, PC-ORD version 5.0) (McCune and Mefford 1999). MRPP tests for clustering of
group members by comparing intragroup average distances with the average distances of
all possible data permutations into groups. Included parameters (NH4, TKN, TDP,
salinity, and maximum lake depth (Zmax)) were selected based on prior analysis of the
lakes (Cooper and Wissel 2012a, Wissel et al. 2011). All data were z-transformed around
a common origin of zero and a standard deviation of one in order to standardize water
chemistry among lakes and isolate temporal differences. Because the transformed dataset
contained negative values, relative Euclidean distance was used as distance measure, and
the distance matrix was weight-transformed due to the high heterogeneity of the data.
Subsequently, I conducted a Principal Components Analysis (PCA) to evaluate
the differences in water chemistry between drought and pluvial periods using averaged
log10-transformed water chemistry data from the two periods. The polar coordinates of
drought and deluge water chemistry data for each lake were exported to Oriana (version 2,
2003) and converted to vectors for a quantitative characterization of the strength and
direction of change in water-chemistry between the two time periods. Oriana calculates
directional statistics by quantifying the distance (r) and direction (μ) that data points shift
from a common origin, testing for a significant directional trend.
Page 30
17
Previous zooplankton cluster analysis results (Wissel et al. 2011) showed that
zooplankton taxa fall into four different groups that were largely defined by the impacts
of salinity, fish assemblage and lake depth. Subsequently, PCA was used to generate
coordinates for each taxon during drought- and deluge-influenced periods.
Environmental parameters used for the PCA were identical to those used for the MRPP,
and log10-transformed abundance data were used to ordinate the zooplankton taxa in
environmental space. The coordinates of the zooplankton taxa during drought and deluge
periods were used as inputs for Oriana to analyze the strength and directionality of
changes in the zooplankton communities between these two periods (see above). The
four identified zooplankton clusters were tested independently for significant directional
change in Oriana.
Fish assemblage data were only collected during the deluge-influenced time
period (2007-2010). Therefore, climate effects on fish assemblages were evaluated
indirectly using discriminant functions based on water chemistry variables. Discriminant
functions were developed for both drought and deluge water chemistry data in order to
determine which climate period best predicted the complexity of fish community as
indicated by the three classes: fishless, planktivorous or piscivorous. Predictive abilities
of averaged water chemistry during drought and deluge periods were quantified using
Discriminant Function Analysis (DFA, PASW Version 18.0). DFA develops a predictive
function for a categorical variable (fish trophic class) using one or more continuous
independent variables (water chemistry). The input variables for DFA (TDP, salinity,
lake depth and DOC) were selected based on ecological relevance and statistical
performance, with variables selected to minimize covariance and maximize effect size.
Page 31
18
Due to eutrophic conditions in combination with shallow water depth, prairie
lakes are often prone to winterkill (Barica and Mathias 1979), and the susceptibility of
individual lakes to winterkill can be altered by climatic conditions (Robarts et al. 2005).
To investigate the potential impacts of drought vs. deluge conditions on the occurrence of
winterkill I applied the model of Barica and Mathias (1979) to the 20 study lakes. This
model was developed particularly for eutrophic prairie lakes, and incorporates initial
oxygen storage, stratification, lake depth and eutrophic status to estimate lake oxygen
depletion for up to 100 days of ice cover. Model outcomes are categorical and lakes are
classified as low, moderate or high risk. Subsequently, I conducted two PCAs using
average water chemistry data for either drought or wet periods. The results of the
winterkill models were then superimposed onto the PCAs of water chemistry during
drought and deluge, respectively, to evaluate the potential role of additional
environmental parameters not included in the model.
2.3. RESULTS
2.3.1. Climate
The most intense drought conditions across southern Saskatchewan were reported
to occur during 2001, a year prior to the study, however, extremely arid conditions
extended until 2003 (Bonsal and Wheaton 2005, Marchildon et al. 2007). The Palmer
Drought Severity Index (PDSI) exceeded -5 in Saskatoon during 2001, and surpassed the
limit for drought in both 2002 and 2003 before rising to neutral values in 2004 (Bonsal et
al. 2011) (Fig. 2.2). Pluvial conditions (PDSI ≥ 3) were recorded in 4 of 6 years
thereafter, and flooding occurred in many lakes starting in 2007 (Saskatchewan
Watershed Authority 2007 runoff report).
Page 32
19
Fig. 2.2: Palmer Drought Severity Index (PDSI) for Saskatoon in August of 2000-2010.
Negative values indicate effective precipitation deficit, while positive values indicate
excess.
Page 33
20
Table 2.2. Comparison of average water chemistry parameters for fishless, planktivorous
and piscivorous lake communities. Lake characteristic records shown here include
maximum lake depth (Zmax, m), salinity (g L-1
), TKN (µg L-1
), TDP (µg L-1
) and DOC
(mg L-1
). For each community type, means are reported for the drought and fluvial
periods and mid-range values are reported for the drought-based DFA model.
Community Fishless Planktivorous Piscivorous
Parameter Drought Deluge Model Drought Deluge Model Drought Deluge Model
Zmax 5.5 7.0 3.0 8.4 10.9 7.5 8.7 9.4 >7
Salinity 24.5 22.0 26.0 11.1 6.6 10.0 1.7 1.2 <8
TKN 4752 4549 4927 2078 1731 2505 1329 1332 <2400
TDP 240 216 612 51 50 88 52 74 102
DOC 80 111 165 34 41 42 22 27 41
Page 34
21
The effects of drought and deluge periods were clearly reflected in the changes in
water chemistry (Table 2.2). Hierarchical cluster analysis revealed two primary temporal
groups based on lake water chemistry, one ranging from 2002-2005 (drought-influenced)
and the other from 2007-2010 (deluge-influenced). A Principal Components Analysis
explained 62.5% of total variation on axis one and two, and showed good separation
between climatic groups along PCA axis 2, mainly reflecting differences of lake depth
and salinity. PCA axis 1 largely separated lakes on the basis of nutrient concentrations,
which differed to a lesser degree between time periods. The MRPP test of these two
groups across all lakes generated a within-group agreement of A = 0.14, with values
above 0.1 in ecological data suggesting a strong separation between groups (McCune et
al. 2002). All further analyses were conducted on the mean values for drought- and
pluvial-influenced temporal groups (2002-2005 and 2007-2010, respectively).
Analysis of survey data with circular statistics revealed that lake depth increased
during the deluge period, while salinity and nutrient concentrations decreased (Fig. 2.3).
This pattern was further supported by the significant differences in mean water chemistry
and morphometric values for the 20 lakes between the two time periods (Table 2.2).
Average lake deepening was at least 1 m for lakes of all fish community types. Total
dissolved solids decreased by 1 g L-1
on average in piscivorous lakes, and by 6 g L-1
in
planktivorous lakes. Reductions in TKN and TDP also occurred in fishless and
planktivorous lakes but were not substantial enough to change the trophic classifications
of individual basins. Fishless lake water chemistry changed less over the observation
period than in the other lake types.
Page 35
22
Fig. 2.3: (a) PCA of water chemistry in the 20 study lakes. Labeled gray arrows
represent water chemistry variables used to generate the ordination space (TDS, TDP,
TKN, DOC, Zmax, Secchi depth). Each black arrow represents an individual lake,
originating at environmental conditions of this lake during the drought and extending to
conditions during the deluge. Generally, directional changes trend toward increased depth
and water clarity and reduced nutrient load.
(b) Circular plot identifying directional changes in water chemistry from drought to
deluge. Labeled arrows represent the superimposed PCA water chemistry variables.
Each unlabeled arrow represents the direction and magnitude of change in water
chemistry for one of the 20 study lakes. With some exceptions, most lakes increased in
depth and water clarity and decreased in nutrient load and temperature. The directional
change was statistically significant (p<0.05).
Page 36
23
Synchronous regional changes in water chemistry resulted in an increase in lake
similarity, forming a tighter cluster when ordinated on environmental axes during periods
of high inputs of water. In contrast similarity of lake chemistry was reduced during
periods of higher solar energy input.
2.3.2. Zooplankton
Four zooplankton taxonomic groups were identified using PCA, explaining 42%
of the total variability on axis 1 and 2 (Fig. 2.4). The four separate groups coincided with
piscivorous communities, planktivorous communities, mesosaline fishless communities
and hypersaline fishless communities, respectively (Table 2.1, see Wissel et al. 2011 for
more detailed taxonomic resolution). Overall, the freshwater piscivorous component was
composed mainly of small daphnids, non-daphnid cladocerans and cyclopoid copepods,
whereas the planktivorous group contained amphipods, calanoid copepods and daphnids.
Large daphnids and calanoid copepods characterized mesosaline fishless zooplankton,
while hypersaline fishless communities contained Artemia franciscana and harpacticoid
copepods. Based on analysis with circular statistics, of the four zooplankton clusters,
only the freshwater piscivore-influenced cluster showed a significant shift in community
composition and species abundances during the transition from drought to deluge (p <
0.05). The remaining three communities showed a stochastic response as they shifted
unpredictably in composition and abundance without significant directional trends.
2.3.3. Fishes
The fish communities in prairie lakes were relatively simple, not exceeding ten
species per lake (Cooper and Wissel 2012b). The most common piscivorous fishes were
walleye (Sander vitreus), northern pike (Esox lucius) and yellow perch (Perca flavescens).
Page 37
24
Fig. 2.4 Principal Components Analysis of zooplankton taxonomic composition during
drought- and deluge-influenced periods. Labeled lines represent taxa used to generate the
ordination space. Vectors of individual lakes within taxonomic space are plotted between
drought and deluge periods, with arrows terminating in the deluge period. Of the four
zooplankton assemblages, the most change occurred in the piscivores group with small
cladocerans in the lower right segment (Ceriodaphnia, Simocephalus, Diaphanosoma).
Page 38
25
Common benthivorous species included white suckers and lake whitefish. Species with
mixed planktivorous / benthivorous preference included spottail shiner and brook
stickleback, as well as juvenile piscivores and a few less prevalent species. Fish
assemblages were most effectively predicted by a discriminant function model that was
based on water chemistry data during drought conditions (Fig. 2.5). The drought-based
model showed 90% classification success relative to the 60% success of the model that
was based on water chemistry during the deluge. Of the included variables, TKN was
most important, and differed notably between the three lake types during this time period
(Table 2.2). In addition to TKN, salinity, DOC and lake depth also distinguished
between fish-habitable and fishless lakes on axis 1.
During the pluvial interval, total nitrogen levels were lower and more similar
across lakes (Table 2.2), and salinity replaced total nitrogen as the most influential
variable separating communities along axis 1. When the drought-based model was
applied to water chemistry values during the deluge period, many formerly fishless,
borderline fresh and mesosaline lakes were identified as fish-habitable. In the case of
Middle Lake, water chemistry values shifted from fishless to planktivorous classification,
consistent with the appearance of planktivore/benthivore communities during the 2007 to
2010 observation period. While the DFA classified two other historically fishless lakes
as fish-supporting during the wet interval, no other fishless lakes were colonized by fish
over the study period.
The risk of winterkill in prairie lakes may have been alleviated during the deluge
due to changes in water chemistry associated with increased precipitation. Based on the
model of Barica and Mathias (1979), nine lakes were at risk of winterkill during the
Page 39
26
Fig. 2.5: Discriminant Function Analysis of the 20 study lakes using drought-influenced
water chemistry data (a) and deluge-influenced chemistry data (b). The discriminant
function classified lakes as fishless (white), planktivorous (gray) or piscivorous (black)
using conductivity, dissolved organic carbon (DOC), depth (Zmax) and nutrient load (TDP
and TKN) as variables. The components of the discriminant functions are listed on each
axis in order of contribution. During the drought and deluge, nitrogen and conductivity
were the most important factors on axis 1. Axis 2 was not significant for either scenario.
Classifications based on the drought model were 90% accurate, with a high level of
distinction between all community types. Deluge-based classifications were only 60%
accurate, with most misclassifications arising due to overlap in the chemistry of
planktivorous and piscivorous lakes.
Page 40
27
Fig. 2.6: Principal Components Analysis of the 20 study lakes during the drought- and
deluge-influenced periods (a and b, respectively), coded by level of winterkill risk
(Barica and Mathias 1979). White indicates low risk, gray indicates moderate risk and
dark gray indicates high risk. Moderate to high winterkill risk was predicted for 9 lakes
during the drought period, which was reduced to 4 lakes during the deluge. Risk is
negatively associated with depth and positively associated with nutrient load.
Page 41
28
drought, seven of which were fishless. Only four of these nine at-risk lakes retained
moderate or high risk during the deluge (Fig. 2.6). No lakes exhibited an increase in risk
of winterkill during the pluvial period, but six decreased in risk ranking. At-risk lakes
shared morphological and chemical properties that formed a cluster in the PCAs. In
particular, shallow lake depth and high nitrogen and phosphorus concentrations were
strongly associated with elevated risk of winterkill. Similarly parameters included in the
Barica and Mathias (1979) model (lake depth and primary productivity (Chl α)), were
prominent factors in the PCA. However, nutrient concentration was also a critical and
non-redundant covariate in the PCAs.
2.4 DISCUSSION
2.4.1. Climate
In 2001, Saskatchewan experienced one of the driest years recorded in a century,
followed by a persistent period of above-average effective precipitation from 2005 to
2010. This climate variation occurred over a relatively brief period of time, yet
comparable drought and fluvial periods have been encountered repeatedly over the last
several centuries (Michels et al. 2007). Longer and sometimes more extreme droughts
and pluvial periods have been occurring on millennial scales in the past, but projections
imply that climate extremes could occur over shorter periods in the future (Michels et al.
2007). Consequently, the observation period may provide a glimpse of processes
representing the future "climate-normal" as a consequence of climate change.
2.4.2. Water Chemistry
The extreme weather shift demonstrably impacted the water chemistry of lakes,
but with a short temporal delay. This delay may reflect the relatively slow flow of
Page 42
29
groundwater into and out of lakes. The groundwater acting as a capacitor could have
mitigated the initial effects of drought by supplying groundwater to the lakes, and storing
excess water during the deluge that was then slowly released into the lakes (Bredehoeft
and Durbin 2009, Webster et al. 2000).
During the observation period, lake water chemistries were more disparate during
the drought interval, but became more similar during the pluvial period, possibly because
the increasing addition of precipitation resulted in an asymptotic approach toward a
freshwater, oligotrophic endpoint. In particular, inputs of water to parkland lakes during
the wet interval had a homogenizing effect compared to the drought period, in contrast to
prior observations that mass inputs tend to reduce intra-annual synchrony of lake
ecosystems (Pham et al. 2008, Vogt et al. 2011). However, both studies evaluated
synchrony rather than homogeneity or directionality of change as performed here, making
the discrepancy less surprising. These apparently contrasting observations may also be
reflected by the sampling frequency, as Pham et al. (2008) and considered intra-annual
changes in synchrony whereas my study focused on a pronounced step change in
effective moisture resulting from the end of a multi-year drought and influx of above
average precipitation (Fig. 2.2). In particular buffering effects of groundwater may mute
the effects of short-term changes in precipitation, yet may increase ecosystem coherence
over decadal time scales. Webster et al. (2000) determined that lakes within geographic
regions with simple surface runoff hydrology responded synchronously to precipitation
changes, whereas groundwater and lake-specific factors resulted in an unstructured
response to short-term changes in water flux. Consistent with this hypothesis, once the
Page 43
30
systems became overwhelmed with mass influx toward later years of the deluge, lakes
increased in similarity, possibly because groundwater ceased to regulate mass flux.
2.4.3. Zooplankton
Zooplankton species compositions are known to vary in response to strong
differences in environmental conditions across the northern Great Plains. Previous
studies by Wissel et al. (2011) and Cooper and Wissel (2012b) show that zooplankton
communities in prairie lakes can be classified into four main functional groups based on
water chemistry and predatory effects of fishes. Governing variables identified here were
similar to prior studies, as zooplankton composition was controlled primarily by salinity,
total nitrogen, lake depth and fish assemblage. While some differences in taxonomic
composition occurred across the studies, the four groups were functionally similar with
substantial taxonomic overlap.
Zooplankton assemblages were more differentiated among lakes in drought years
than in pluvial periods, as indicated by a higher taxonomic diversity among lakes (ß-
diversity). These findings contrast with those of Chase (2007), who found that ß-
diversity of zooplankton decreased during drought years in fishless ephemeral ponds.
This comparison suggests that size and permanence of a water body, as well as trophic
complexity, are important factors to consider before generalizations can be made across
ecosystems. Nevertheless, increased ß-diversity likely occurs in lakes of the northern
Great Plains due to the highly variable hydrology of endorheic basins, resulting in diverse
responses to a uniform climate influence (Last and Ginn 2005). In contrast, prolonged
high precipitation would partially negate differences in local hydrology and create more
uniform (dilute) water chemistry across ecosystems.
Page 44
31
Zooplankton communities in piscivorous freshwater lakes increased in similarity
and declined in population density during the pluvial interval. Potentially, these changes
may reflect a greater sensitivity of freshwater species to changes in environmental
conditions; however, I discount this explanation because absolute changes in water
chemistry were small relative to documented tolerances of most species (Table 2.2).
Similarly, dilution did not reduce Chl α, as nutrients likely did not dilute sufficiently to
limit algal growth, and in fact, mesotrophic lakes exhibited elevated Chl α content during
the wet interval, while Chl α decreased in historically eutrophic and hypereutrophic
systems. Instead I hypothesize that changes in abundances of freshwater zooplankton
taxa were due to increased young-of-the-year recruitment and planktivory during the
deluge when fishes encountered more suitable environmental conditions. Fish
communities are key regulators of overall lake communities, and their loss or
reintroduction has drastic impacts on zooplankton composition (Srinivasan et al. 2007).
In contrast, zooplankton communities in fishless mesohaline and hypersaline lake seemed
robust to changes, and exhibited little response to climate-related water chemistry
changes over the study period.
2.4.4. Fishes
Discriminant function analysis suggested that the level of habitability for a
complex fish community depended on a functional interaction of nitrogen, phosphorus,
lake depth, salinity and DOC. Although abundances of fish populations could not be
quantified, past studies have observed strong associations between fish recruitment and
climate-related environmental variables. Fry recruitment can vary considerably between
years due to changes in lake depth and water chemistry (Zalewski 1990), with late-
Page 45
32
summer water temperature and disease prevalence being particularly good predictors of
recruitment of some species (Paxton et al. 2004).
The long-term habitability of lakes for gamefish appeared to be defined best by
conditions associated with drought intervals, in particular osmotic stress (salinity) and
winterkill risk. The factors influencing winterkill risk also appeared to be interactive
because multiple variables contribute to the rate of under-ice oxygen depletion. Risk of
winterkill coincided with community classifications based on DFA, as moderate to high
risk lakes only occurred in the fishless and planktivorous lakes, while all piscivorous
lakes were classified as low risk. Piscivorous communities only existed in lakes with
salinities below 2 g L-1
, despite being documented in lakes with salinity up to 7 g L-1
in
recent history (Bayly 1972, Burnham and Peterka 1975, Rawson and Moore 1944).
Consistent with these observations, the DFA identified a conductivity threshold of < 8
mS m-1
for piscivorous communities (approximately equivalent to 5 g L-1
). The
interactive effects of depth, nutrient levels and primary productivity (leading to winterkill)
seemed, however, to limit fish populations in several of the 20 study lakes before
conductivity thresholds were reached. The importance of DOC as a predictive variable is
likely associated with in situ DOC production due to high primary productivity (Baines
and Pace 1991), rather than its effects on visual and thermal lake properties (Williamson
et al. 1999), as in these systems DOC is colorless and has no associations with light,
temperature or oxygen profiles (Cooper and Wissel 2012a). It is also correlated with
Chl-α, yet more coherent between measurements and therefore a more reliable indicator
for this sampling regime.
2.4.5. Ecosystem responses
Page 46
33
The decadal transition from drought to deluge intervals was accompanied by
pronounced changes in water chemistry in prairie lakes, with lake-depth increases and
reduced concentrations of nutrients and major ions (Fig. 2.3). Concomitantly, fish
communities were largely determined by lake habitability due to changes in water
chemistry during droughts. In contrast zooplankton showed a predominantly asystematic
response to meteorological variability during the study period, as most taxa were resilient
to short-term drought impacts, the notable exception being piscivore lakes where declines
in abundances during the deluge appeared to reflect more planktivory by fish. It is
probable that zooplanktivore lakes did not experience a similar increase in recruitment
and planktivory due to their differing habitat requirements from piscivore larvae.
In the future the northern Great Plains could experience substantial fragmentation
and differentiation to ultimately have a wider diversity of community types, but may
support fewer gamefish species and have fewer lakes containing piscivorous fish.
Endorheic basins are particularly sensitive to climate-change impacts, as they respond
rapidly to drought events due to their shallow morphology and long retention times.
Direct and indirect effects of drought and warming, such as reduced lake depth and
increased productivity, can increase the risk of winterkill, potentially extirpating fish
populations. Yet, if precipitation events were to drastically increase, as envisioned in
some climate scenarios (Barrow 2009), saline communities could be lost, resulting in a
reduction of landscape diversity and an increase in gamefish-supporting lakes. In such a
scenario, the thermal stress of climate change could replace water chemistry as the
primary determinant of game fish prevalence.
2.4.6. Conclusions
Page 47
34
This and other studies have shown that endorheic lake systems across the northern
Great Plains and other semi-arid regions are sensitive to the impacts of increasing
temperatures and extreme weather due to future climate change (Adrian et al. 2009, Bates
et al. 2008). In general, the impacts of climate change on endorheic systems will appear
to be more severe than anticipated for boreal temperate lakes, which are expected to be
largely impacted by deepening of epilimnion, stronger stratification and changes in light
penetration (De Stasio et al. 1996, Fee et al. 1996, King et al. 1999, Schindler 1998,
Williamson et al. 2009). The greater resilience of boreal lakes to climate change is in
part due to stronger hydrological connectivity, lower retention times and stability of
thermal stratification (Fee et al. 1996, Heitman 1973, Kling 1988, Schindler et al. 1996,
Webster et al. 2000). Consequently, endorheic lakes should be considered separately in
respect to the impacts of future climate change on lake ecosystems.
Page 48
35
3. CONDITION AND PARASITE LOAD OF FISHES IN RESPONSE TO
CLIMATE-INDUCED CHANGES IN TEMPERATURE AND SALINITY
3.1 INTRODUCTION
Assemblages of fish communities can be impacted by climate-related changes in
environmental conditions, which ultimately determine the presence or absence of
individual species (Bond et al. 2008, Magoulick and Kobza 2003). In particular,
temperature, salinity, lake depth and nutrient levels define the limits of fish habitability
either directly by exceeding physiological limits of fishes (Mehner et al. 2005, Ostrand
and Wilde 2001), or indirectly, for example by causing seasonal hypoxia (e.g., winterkill
or summer hypoxia) (Robarts et al. 2005).
As environmental parameters vary tremendously among lakes, sub-lethal impacts of
these factors should occur within the fish-habitable range, and may manifest as lower
relative body mass or impaired immune function. Laboratory studies revealed sensitivity
of fishes to varying concentrations of dissolved compounds such as salts and metals
(Anderson 1997, Victoria et al. 1992) or varying environmental temperatures (Cherry et
al. 1977). Such approaches provide invaluable information on physiological limits of
fishes under controlled conditions, but they cannot capture longer-term and interactive
impacts under natural conditions. Studies on natural systems have shown negative
impacts of extreme climate events, such as droughts on both lotic and lentic systems
(Bond et al. 2008), which are usually a consequence of changes in hydrology, water
chemistry and temperature.
Lakes in endorheic (closed) drainage basins are particularly sensitive to the impacts
of climate variability, as salinity, nutrients, temperature and water depth in endorheic
Page 49
36
lakes all respond quickly and significantly to altered weather conditions (Adrian et al.
2009, Chapter 2 of manuscript). Hence, lake systems such as prairie lakes in the
Canadian Great Plains likely represent very good model systems to study the potential
impacts of climate variability on fish communities (Hammer 1990, Williams 2002).
Understanding the relationships between climate variability and fish health and
performance is particularly important because the potential impacts of future climate
change raise new questions about the sustainability of freshwater fisheries, particularly in
drought-prone areas (Ficke et al. 2007). In semi-arid systems, such as the northern Great
Plains, lake-specific water chemistry is governed by the relative importance of
evaporation, precipitation and groundwater inputs (Last and Ginn 2005). Under a high-
precipitation climate-change scenario (Barrow et al. 2009), excess precipitation may
compensate for increased evaporation due to elevated temperatures, and may maintain
sufficient water quality to sustain suitable lake habitats for fishes. Yet, a warm but
relatively dry future climate-change scenario (Barrow et al. 2009) would enhance
evaporation and result in increased salinity and dissolved nutrients in lake basins, in
combination with lower water depths (Chapter 2 of manuscript). In the latter scenario,
water chemistry changes could be a significant threat to fish populations. Negative
relationships between body condition of fishes and eutrophic state, stocking density, and
salinity are often reported, as well as positive relationships with water level and access to
prey (Dicenzo et al. 1996, Fisher et al. 1996, Marwitz and Hubert, 1997).
It is foreseeable that the predicted thermal increases of up to 5 ºC by 2050
(Barrow 2009) may have both positive and negative effects. For species with fairly
narrow thermal requirements, preliminary findings suggest a warming climate could
Page 50
37
improve body condition in colder areas of the species range, but reduce condition in
warmer areas (Drinkwater 2005). A water temperature around 22 ºC appears to be
optimal for walleye (Sander vitreus) and yellow perch (Perca flavescens) in captivity as
well as in natural systems, while temperatures significantly above or below may result in
reduced growth rates (Huh et al. 1976, Quist et al. 2002). Hence, temperature increases
may offer a direct benefit of enhanced metabolism and faster weight gain if adequate
food is available, and the associated reduction in winter duration could lower the risk of
winterkill and extend the growing season (Quist et al. 2002). Yet, excessive temperature
increases could be detrimental or lethal once bioenergetic limits are surpassed (Ficke et al.
2007, Quist et al. 2002), especially in isothermal prairie lakes that often lack a cool
hypolimnion. The metabolic changes from increased temperatures can reduce maximum
consumption rates while increasing metabolic rates, and either increase risk of starvation
or increase risk of predation through compensatory feeding activity (Biro et al. 2007).
Additionally, long hot summers are often associated with toxic algal blooms and hypoxia,
threatening fish condition and survival (Heisler et al. 2008).
For this study I chose fish condition and intestinal parasite load as measures of
fish health. Parasite load is a good indicator of immune function and host-pathogen
interactions, is associated with water quality (Lewis et al. 2003, Poulin 1992), and can
significantly respond to environmental change (Marcogliese et al. 2005). At elevated
concentrations, naturally-occurring salts and other inorganic and organic compounds can
function as environmental stressors and increase parasite prevalence (Bly et al. 1997).
By suppressing immune functions or creating a more favorable environment for the
pathogens, water chemistry changes can increase host infection rates (Hauton et al. 2000).
Page 51
38
Relative weight (Wr) of fishes (or plumpness) is another important metric to
express fish performance in a changing environment (Blackwell et al. 2000). This metric
utilizes species-specific standard weight equations derived from a regression-line
percentile technique using representative population data (Murphy et al. 1990). Low
plumpness not only reduces the commercial quality of a fish, but is also a precursor to
starvation (Dutil and Lambert 2000). Fish body weight is a commonly-used indicator of
the health of a fishery, and plumpness metrics such as relative weight have been used as
health indicators, to assess prey availability, to inform adjustment of stocking levels, to
make comparisons between fish populations, and as an indicator of habitat quality
(Blackwell et al. 2000).
For this study, I selected walleye as the primary focal species. As a common top-
predator, it not only regulates many lake functions (Vander Zanden and Vadeboncoeur
2002), but is also a desired species for commercial and recreational fisheries (Fenton et al.
1996). Additionally, many lakes in southern Saskatchewan are supplemented with
walleye coming from a single hatchery (Lake Diefenbaker). Hence, walleye populations
in this region are genetically uniform across lakes (Starks, unpublished data), allowing
for isolation of environmental from genetic effects. Unmanaged co-existing northern
pike (Esox lucius) and yellow perch populations were also evaluated to determine if any
of the observed environmental effects were species-specific.
To be able to analyze the importance of lake-, climate- and species-specific
effects on fish health, I included seven different lakes over a two-year period that was
characterized by contrasting climatic conditions. I expected that lake-specific effects
Page 52
39
would dominate as the study lakes differed strongly in water chemistry and lake
morphometry, followed by species and climate effects.
3.2 METHODS
3.2.1 Study sites and field sampling
All seven study lakes were sampled in June and August of 2009 and 2010 (Table
3.1). Lakes varied in water chemistry and physical characteristics, ranging in salinity
from 0.5 to 3 g L-1
of total dissolved solids (TDS). Trophic state ranged from
mesotrophic to hypereutrophic, with a large variability in both nitrogen and phosphorus
concentrations. All lakes were located in the prairie ecoregion of southern Saskatchewan
in central Canada (Fig. 3.1). Of these seven lakes, all had walleye populations, and six
were inhabited by yellow perch and northern pike.
Pelagic fish were netted overnight using a single multi-panel 1.5 x 30 meter
gillnet with 55 mm mesh size, set perpendicular to shore at a starting depth of 3 meters.
Littoral fishes were captured with a single 2 x 20 meter beach seine during daytime with
5 mm mesh size from shore to a depth of 2 meters. Fish were identified to species (Scott
and Crossman 1973) and frozen prior to further analysis. Standard length and weight
were recorded for all fish samples, as well as gonad and liver weights for fishes
exceeding a length of 10 cm. Stomach contents were identified to family, and recorded
in order of prevalence. Internal parasites were recorded for all three piscivorous fish
species (walleye, northern pike and yellow perch) in adults exceeding 10 cm.
Sampling for water chemistry was conducted at the location of maximum depth,
which was generally in the central area of each lake. Dissolved oxygen (mg L-1
), total
dissolved solids (g L-1
), conductivity (mS cm-1
), water temperature (°C) and pH were
Page 53
40
Table 3.1. Environmental characteristics of the seven lakes included in the study from
2009 to 2010. Gamefish recorded include walleye (W), northern pike (N) and yellow
perch (Y). Recorded values of max depth in meters, surface area in km2, salinity in g L
-1
(as total dissolved solids, TDS), Total Kjeldahl Nitrogen in ug L-1
(TKN), Total
Phosphorous in ug L-1
(TDP) and Dissolved Organic Carbon in mg L-1
(DOC) are
means from June and August of 2009 and 2010.
Lake Max Depth Surface Area TDS TKN TDP DOC Gamefish
Humboldt 8.0 19 0.80 2159 249 29.4 WNY
Fishing 19.5 32 1.58 1194 9.8 32.5 WNY
Pelletier 10 2 0.41 757 11.7 16.4 WNY
Kipabiskau 8.5 5 0.44 1413 162 32.5 WNY
St. Brieux 11 0.5 1.40 1110 22.0 37.5 WNY
Wakaw 11 11 1.61 990 6.3 18.1 WNY
Shannon 10 1 1.70 1700 56.0 20.4 W
Page 54
41
Fig. 3.1. Location of the seven study lakes (filled circles) in the grassland ecoregion of
southern Saskatchewan, Canada. Ice dates were recorded from Buffalo Pound Lake
(filled square), and air temperature data were recorded from Saskatoon.
Page 55
42
measured throughout the water column during daytime at 1 m intervals using a YSI multi
probe (model 556). Water transparency was measured with a 20 cm black and white
Secchi disk. Using a tube sampler, I collected integrated, prescreened (80-μm mesh)
water samples for water chemistry (total Kjeldahl nitrogen (TKN), nitrate (NO3-), total
phosphorus (TDP), soluble reactive phosphorus (SRP), ammonium (NH4+), dissolved
inorganic carbon (DIC), dissolved organic carbon (DOC), calcium (Ca) and chlorophyll
a concentration (Chl a)). The tube sampler was suspended in the water column either
down to 6 m for deeper polymictic lakes, or down to 0.5 m above bottom sediments for
shallower polymictic lakes. Water samples were filtered through sterile 0.45 μm glass
fiber filters and analyses of TDP, TKN, NO3, NH4, were conducted at the University of
Alberta water chemistry laboratory after Stainton et al. (1977). DOC and DIC were
quantified using a Shimadzu TOC5000 analyzer at the Environmental Quality Analysis
Laboratory (EQAL) at the University of Regina. To estimate algal biomass, (Chl a), an
aliquot of 1 liter was filtered onto GF/C filters and stored in -10 ºC until extraction with
acetone:methanol:water (80:15:5) and analyzed using trichromatic methods (Wetzel and
Likens, 1991).
3.2.2 Statistical analysis
Daily air temperatures were obtained from the Environment Canada weather
station in Saskatoon for 2008 through 2010 (Fig. 3.1). The dates of ice-off to ice-on were
recorded for nearby Buffalo Pound Lake, which serves as a drinking water reservoir for
the city of Regina (Fig. 3.1), and used as a proxy for the length of the growing season
across the study area. This was deemed appropriate due to the high coherence of ice
coverage dates between lakes in the region (Vogt et al. 2011).
Page 56
43
To evaluate the response of fishes to environmental conditions and community
factors, relative weight (Wr) indices were calculated for each gamefish species at each
lake, and averaged for individual sampling periods (Murphy et al. 1991). Wr scores were
calculated with species-specific length-weight functions based on a representative sample
of North American fish populations. Species-specific functions were used for walleye,
northern pike and yellow perch (Blackwell et al. 2000). Scores between 95 and 105 are
considered optimal for the species. Scores near or below 80 imply depleted fat reserves
due to thermally-induced weight loss and/or lack of food, while extremely high scores
can be indicative of hypertrophy or obesity (Murphy et al. 1991).
To compare growth patterns of fishes over time, linear regressions of weight by
cubed length were calculated for each species and sampling period. Because individual
lake populations were not significantly different within sampling periods (1-way
ANOVA P>0.05), length and weight data from all lakes were pooled for each period and
species. Regression lines were tested for significant differences in slopes and intercepts
between sampling periods. Samples were restricted to similar size ranges (> 100 mm) for
each comparison to avoid confounding effects of sub-adult growth patterns.
Due to low parasite densities in yellow perch and northern pike, tapeworm
infection levels could only be sufficiently evaluated for walleye. Gut tapeworm parasite
loads for walleye were compared between lakes and over time, and used as an indicator
of immune performance. Number of tapeworms per host and proportion of lake
population infected were recorded for each lake and time interval. These values were
tested against other environmental factors such as water chemistry using linear
regressions; and the effects of gender, lake, year and season were quantified using t-tests
Page 57
44
or 1-way ANOVAs. Other parasites, such as nematodes and Ergasilus sp. were also
recorded but did not occur in sufficient abundances to be further analyzed.
3.3. RESULTS
3.3.1. Climate
The ice-free season of 2008 was slightly longer than long-term average conditions,
and had many days within the optimal thermal range for cool water species, such as
walleye (10-22 ˚C) (Cherry et al. 1977, Huh et al. 1976). 2009 had a short ice-free period
and was characterized by many days outside the optimal metabolic range for walleye
(Table 3.2). In contrast, the 2010 ice-free season was particularly good for walleye, with
the vast majority of days in the optimal thermal range (197 days in 2010, relative to 146
in 2009).
3.3.2. Body condition
Throughout the study period, I analyzed 129, 121, and 36 individuals for walleye,
perch and pike, respectively. The range of specimens per sampling date varied from 21
to 36 for walleye, while for pike and perch the range was 0 to 25 and 11 to 48,
respectively. Walleye, pike and perch did not differ significantly in weight and size class
between lakes due to high intrapopulation variability, and were therefore pooled by
species for all subsequent analyses. When averaged across all seven lakes and both years,
walleye exhibited an intermediate mean relative weight (107) in comparison to perch (97)
and pike (119). Walleye and perch achieved the Wr target proposed by Murphy et al.
(1990) (Table 3.3), while northern pike were exceeding the target range, despite
performing poorly in some lakes (data not shown).
Page 58
45
Table 3.2. Climate conditions from 2008 to 2010. Temperatures are recorded from the
Saskatoon weather station and ice cover data were taken from the Buffalo Pound water
treatment facility. Ice duration data are recorded for the entire winter leading up to the
warm season of the indicated year. Mean highs and lows are recorded for the full ice off
period.
Year Mean Low (ºC)
Mean High (ºC)
Ice Duration (Days)
Ice Free (Days)
Ice Off Date
Days > 22 ºC
Days > 27 ºC
Days < 10 ºC
2008 4.2 18.5 150 214 19 April 87 26 13
2009 3.6 17.2 163 206 1 May 83 24 36
2010 4.9 16.5 137 224 9 April 61 9 18
Norm 3.9 16.8 155 210 18 April NA NA 44.7
Page 59
46
Table 3.3. Wr for walleye, pike and perch over the full study period from 2009 to 2010.
NA = sample size not adequate for calculating an average.
Sampling Period Walleye Perch Pike
June 2009 156 (n=21) 91 (n=48) 157 (n=8)
Aug 2009 95 (n=49) 92 (n=35) NA
June 2010 81 (n=23) 89 (n=27) 84 (n=25)
Aug 2010 105 (n=36) 123 (n=11) 117 (n=3)
Grand Mean 107 ± 1.7 97 ± 4.7 119± 3.6
Page 60
47
Walleye length-weight regressions performed on pooled monthly data (all lakes
combined) were significantly different between June and August of both years,
decreasing in intercept in 2009, but increasing in intercept in 2010 (P < 0.01, no
difference in slopes) (Fig. 3.2). While Wr of walleye was significantly different between
June of 2009 and June of 2010, there was no difference between Wr in August of both
years. Yellow perch length-weight regressions did not change significantly from early to
late summer of 2009, but intercepts differed significantly in 2010 due to an increase in Wr
from June to August (Fig. 3.3). Unfortunately, data for pike were insufficient for
regression analysis of length-weight relationships.
Walleye and northern pike were in good condition during early 2009, and
substantial body fat reserves were observed upon dissection (Sydney Chow, personal
observation). Despite their high fitness in early summer, these two species were slightly
underweight in late summer of 2009. Conditions of yellow perch were somewhat
different, as Wr in June of 2009 was only 91 (slightly below the target range and well
below the high values recorded for the other species), and values did not further increase
until August of 2009 (Fig. 3.3). Over the warmer and longer ice-free season of 2010,
body mass generally increased from June to August for all three species. In June of 2010,
walleye condition was poor, with a low relative weight of 81, but by August mean
walleye condition across lakes reached a relative weight of 105. The seasonal increases in
Wr values for pike and perch were even higher with increases from 85 to 123 and 84 to
117, respectively (Table 3.3).
The majority of walleye populations achieved the target weight within 2010, but
in Shannon, Kipabiskau and Fishing Lakes body conditions remained poor. Overall,
Page 61
48
Fig. 3.2. Regressions of length by weight revealed a significant decrease in walleye body
weight by length from June to August 2009 (above), and an increase in 2010 (below).
Page 62
49
Fig. 3.3. Regressions of length by weight revealed no significant change in perch body
weight for the full lake set from early to late 2009, and an increase in weight from early
to late 2010.
Page 63
50
Fishing and Kipabiskau Lakes were below target for all three species (not significant),
while stocked walleye represented the only gamefish species in Shannon Lake.
Potentially, the lower performance of walleye in these lakes may have been associated
with poor water quality, as these lakes did not represent an optimal habitat for walleye.
While Shannon and Fishing Lakes had elevated salinities, Kipabiskau and Shannon Lakes
were eutrophic to hypereutrophic, based on integrated water column samples of total
phosphorus and chlorophyll a. However, no significant differences in fitness were
detected among lakes, and linear regressions revealed no predictive relationship between
Wr and water chemistry or lake morphometry parameters.
Liver somatic index (LSI) averaged approximately 1 for walleye, 1.15 for
northern pike, and 1.3 for yellow perch, indicating that liver weight was relatively normal
in these populations (Adams and McLean 1984). LSI did not differ between 2009 and
2010 for any of the species, or between June and August. Furthermore, LSI did not
covary with Wr or any of the measured water chemistry or climate parameters.
3.3.3. Stomach contents
Stomach fullness was not significantly different between 2009 and 2010 for any
of the three species. Common stomach contents of walleye and pike included small
fishes (brook stickleback and fathead minnows) and benthic prey (amphipods and larval
insects). Yellow perch consumed zooplankton, benthic prey and small fishes.
Approximately 30% of sampled walleye had empty stomachs, and 35% of pike and <1%
of perch stomachs were empty. Some differences in stomach fullness between lakes were
observed over the study period, but did not correlate with inter-lake differences in fitness
or any environmental parameters. Pike and perch in Kipabiskau Lake exhibited
Page 64
51
significantly higher benthivory than in other lakes, while walleye in Shannon Lake
exhibited a higher degree of benthivory (P < 0.05).
3.3.4. Parasite load
Tapeworm levels were highest in walleye, and significantly more prevalent in
2009 than in 2010 (Fig. 3.4). Population-level parasite prevalence was significantly
correlated with salinity in 2009 (r2=0.93), but not in 2010 (Fig. 3.5). Additionally,
salinities were significantly higher in 2009 (P<0.05). No other environmental parameters
showed any significant relationships with tapeworm prevalence. Larger walleye (total
length > 10 cm) were significantly more likely to host tapeworms (P<0.05).
Subsequently, tapeworm loads were also tested as weight-based loads against
environmental variables, but this adjustment did not change the general findings. In 2009,
female and male walleye had similar rates of tapeworm infection, but when infected,
females hosted significantly more tapeworms (p<0.05). In 2010, only female walleye
were infected.
Page 65
52
Fig. 3.4. Female walleye (gray bar) exhibited a significantly higher mean tapeworm
count than males (hatched bar) in both years. Overall tapeworm counts were
significantly higher in 2009 than 2010. No males were infected in 2010.
Page 66
53
Fig. 3.5. Proportion of walleye populations infected with tapeworms was significantly
associated with salinity in 2009 (r2 = 0.934), but not in 2010. Shannon Lake, excluded
from the analysis, is shown in white.
Page 67
54
3.4. DISCUSSION
Game fish fitness and parasite load, among walleye in particular, correlated
significantly with changes in salinity and temperature over the two year study.
Accordingly, temperature and precipitation were related to the fitness of ambient fish
populations. Furthermore, the changes in fish condition and parasite load were even
detectable on short-term (annual) time scales. Surprisingly, no other factors known to
predict fish habitability were influential in this study. Potentially, environmental
parameters such as nutrient concentrations, lake depth, Secchi depth and algal biomass
established the habitability thresholds of the lakes by increasing winterkill risk, but had
no observable effects on ambient fish populations during summer. In this case, no effects
would be observable below the threshold for winterkill. This is in contrast to other
gamefish populations, which have shown seasonally reduced condition along with lower
Secchi depth (Craig and Babaluk 1989), poor growth with small lake surface area (Cena
et al. 2006), and changes in recruitment relative to lake depth, volume and Chl-α (Mehner
et al. 2005). Nevertheless, the impacts of lake morphometry and trophic state on a
landscape scale are critical, but are more likely to control the inhabitability of lakes for
fishes, rather than controlling their condition in those lakes that have established fish
communities. Therefore, established gamefish of this region are likely to be more
strongly impacted by natural variability of weather conditions and climate change, as
both temperature and precipitation are directly linked to climate and hydrology.
3.4.1. Body condition
The poor fitness values in June of 2010 were likely a result of the shorter growing
season and less optimal temperatures during 2009. In 2009, many more days were above
Page 68
55
and below the optimal thermal range for walleye growth, which likely resulted in the
decline of body condition from June to August. Temperature has a direct effect on fish
bioenergetics, increasing the metabolic rate of fish (Kitchell et al. 1991). Assuming that
adequate prey is available, high temperatures will result in weight gain up to a critical
threshold, beyond which it is physically impossible to sustain body mass. High
temperatures have not only been shown to induce thermal stress but also reduce overlap
of walleye with warm water prey species at temperatures above 22˚C, even in polymictic
lakes (Kocovsky and Carline 2001, Quist et al. 2002). However, a warmer spring and fall
could compensate for summer growth losses if summer temperatures remain sublethal
(Quist et al. 2002).
In 2010, the number of days in the optimal thermal range were well above-
average, resulting in a substantial improvement in condition over the 2010 growing
season. These trends were uniform across the target size range, despite findings in other
studies that temperature effects can be size-dependent (Carginelli and Gross 1997,
Gabelhouse 1991). It is noteworthy that mostly young of the year (YOY) fish, which
were not included in this study, may experience high over-winter mortality and condition
decline relative to older conspecifics (Pangle et al. 2004). Hence, the lack of size-
dependency was not a concern.
While dietary proportions of fish and benthic prey varied among lakes, this did
not have an apparent effect on body condition. Paradis et al. (2006) found that walleye
maintained similar body condition regardless of dietary composition, but grew fastest on
a piscivorous diet. Unfortunately, this effect could not be tested here, as not all fishes
could be aged. Furthermore, food availability and bottom-up effects were not likely to
Page 69
56
have impacted the condition of fishes. Almost no perch were found with empty stomachs,
indicating that food was not limiting (Knight et al. 1984). For walleye and pike, only 30
to 35% of specimens had no detectable food in their stomach, which is fairly common for
species that have a higher degree of piscivory and therefore rely on less frequent meals of
larger prey items (Chapman et al. 1988).
3.4.2. Parasite load
While salinity did not correlate with body condition, there was a significant
increase in parasite load and prevalence in walleye corresponding with increasing salinity
in 2009. The insignificant relationship between salinity and parasite load in 2010 may
have been due to the concurrent freshening of all lakes between years. In 2009 salinities
exceeded 2.0 g L-1
, whereas in 2010 they were all below 1.5 g L-1
in all lakes. Tapeworm
levels were correspondingly low across all lakes in 2010. The strong correlation between
parasite load and an environmental stressor (salinity) was not surprising. Environmental
stress and pathogenicity are believed to act synergistically (Marcogliese et al. 2005);
however, fitness effects of this parasite and environment interaction were beyond the
scope of this study. If the here identified threshold of 1.5 g L-1
salinity is universal or
region-specific will have to be evaluated in future studies. Koel and Peterka (1995)
identified a threshold between 1 and 2 g L-1
sodium sulfate for walleye hatching success;
however, thresholds for effects of salinity on immune function in walleye are not
documented.
Several alternative explanations for the discrepancy in parasite load between years
could not be directly evaluated here. The interannual difference in temperature may have
had some effect on cestode life cycles; however, this seems unlikely as eggs are often
Page 70
57
stimulated to hatch by warmer temperatures (Wright and Curtis 2000), and 2010
experienced warmer temperatures overall. It is also plausible that selective mortality of
infected hosts occurred between years, although this would not fully explain the apparent
lack of new infections in 2010. I therefore hypothesize that salinity-induced
immunocompromisation of the host is the most likely explanation for the salinity and
parasite load correlation in 2009 and the significant interannual difference in parasite load.
3.4.3. Species-specific differences
Despite the similar general trends across species, I detected several differences
among walleye, pike and perch. Walleye were a significantly greater definitive host for
tapeworms in this system as in others (Muzzall and Haas 1998), and was the only species
that exhibited a change in tapeworm load and prevalence in response to environmental
salinity. Pike and perch exhibited such low parasite load in both years that this data
could not be analyzed. Most likely, the cause for different parasite prevalences is due to
dietary preferences, which include a larger benthic composition (Beaudoin et al. 1999,
Hayward and Margraf 1987). Pike and perch are therefore less exposed to vertical
transfer of parasites from pelagic prey (Marcogliese 2002). Yellow perch were
apparently least selective in regard to prey, as their stomachs were full significantly more
often than those of pike and walleye, which is like associated with their lower trophic
position (Cooper and Wissel 2012a).
Overall, walleye seemed to be the most vulnerable species to observed
environmental changes based on its greater intra-annual and inter-annual changes in
fitness and parasite load, followed by pike and perch. Hence, while stocking efforts
maintain walleye populations in these lakes, their overall resilience is low relative to
Page 71
58
populations of naturally occurring fish species. In fact, the single-source stocking
approach for walleye may have detrimental effects on these populations by inviting larger
parasite loads (Poole and Dick 1985) and reducing mean body weight (Li et al. 1996).
3.4.4. Anticipated future effects of climate change
Overall growth rates of fishes will likely improve if climate change would only
prolong the duration of summer without increasing maximum temperatures. Yet, the
predicted increase of 5˚C by 2050 relative to current maximum temperatures (Barrow
2009) will likely result in severe metabolic deficits in walleye and other fish species.
Observations of walleye in more southern lakes of the Great Plains indicate that increases
in summer temperatures result in growth declines; yet, simultaneous warming in spring
and fall may actually increase total growth over the whole year (Quist et al. 2002).
Simulations predict that maximum lake temperatures in southern Saskatchewan
will increase to about 25 to 30 °C by the year 2100, relative to current maxima of 20 to
25 °C (Sharma et al. 2007). Current temperature maxima bracket the thermal optimum
for walleye (22 °C), while future conditions would not only go beyond the thermal
optimum, but would also reach or even exceed the bioenergetic (lethal) limit of adults
(27 °C). While walleye and yellow perch have similar optimal thermal ranges (Kitchell
et al. 1977), perch aren’t as detrimentally affected by suboptimal temperatures (Huh et al.
1976). Based on this information, perch should not be as negatively impacted by climate
warming. Unfortunately, pike samples were not adequate in 2009 to measure the impact
of temperature on growth, but since their thermal optimum is also similar to walleye
(Christie and Regier 1988), this species is also like to be natively impacted in the future,
though it is unclear to what extent.
Page 72
59
Among the three piscivorous species, walleye appeared most vulnerable to
climate change in these lakes. In addition to their fitness correlations to temperature and
salinity changes, only walleye exhibited a high parasite load. Walleye body condition
was overall less variable among lakes and over time than for the other two game fish
species, potentially due to the low inter-lake genetic variability of this species as a result
of single-source stocking (Starks, unpublished data). Because of their popularity as game
and commercial fish species, identifying and prioritizing more optimal habitats for their
maintenance should be a priority. Small lakes such as those studied here may not be
suitable in the long term, and more potentially resilient species such as northern pike and
yellow perch are likely to replace walleye in small lakes (which would largely also be the
current scenario in the absence of stocking). Larger and deeper freshwater lakes that
would be more suited to support walleye are rare in the prairie region (Last and Ginn
2005), but might be a prime habitat in northern, boreal parts of Saskatchewan, especially
as temperature increase.
In systems other than endorheic prairie lakes, lake size, genetics and diet, vary
much more than in this dataset (Cena et al. 2006), while hydrology is less variable (Pham
et al. 2008, Winter and Rosenberry 1998). Temperature and salinity were relevant to
fitness of game fishes in Saskatchewan prairie lakes, but lakes that are thermally stratified
and less prone to salinization are unlikely to be negatively impacted by climate change to
the same degree as prairie lakes (Adrian et al. 2009). Nevertheless, this study clearly
indicates that endorheic systems, such as prairie lakes have a high capacity as model
system to study early impacts of climate change on lake ecosystems. This is particularly
true, since the two study years fell into a relatively wet period and the observed changes
Page 73
60
in climate and water chemistry did not reflect the typical gradients that can occur in this
system. Hence, longer-term studies will likely identify even stronger effects on fish
fitness.
Page 74
61
4. CONCLUSION
4.1 SYNTHESIS
Species assemblages in endorheic lakes of southern Saskatchewan are highly
diverse due to the extreme inter-lake variability in environmental parameters and trophic
interactions (Cooper and Wissel 2012b). This study revealed that regular drought events
such as the recent subdecadal drought functioned as an environmental filter to define fish
community composition in prairie lakes, and divided lakes into fairly discrete clusters of
fishless, planktivorous (plus benthivorous) and piscivorous (plus planktivorous) fish
communities. Environmental variables most relevant to fish community composition
were nutrient concentrations and lake depth, which are strongly associated with winterkill
risk. Zooplankton responded to fish predation pressure with a reduction in overall
quantity and biomass during wet periods, rather than a change in species composition.
These observations offer an interesting contrast to other lakes which may stratify to
reduce habitat thermal habitat overlap of zooplanktivorous fishes with their zooplankton
prey, thereby reducing predation pressure (De Stasio et al. 1996).
While environmental variables, such as lake morphometry, nutrient levels and
salinity have been described to be important for defining fish community composition
(Cooper and Wissel 2012b, Mehner et al. 2005), surprisingly few of these variables
produced visible sub-lethal effects for walleye, perch and pike in my study system. For
prairie lakes, seasonal and interannual temperature variation significantly affected body
condition of these species, while the large environmental differences among lakes had no
influences. Furthermore, salinity had a significant positive correlation with parasite load,
but this effect was only detected for walleye populations. In climate change scenarios
Page 75
62
(Barrow 2009), variables that induce sub-lethal effects (temperature and salinity) are
predicted to change more uniformly across lakes than most other variables, indicating
that climate change effects on fish health would be fairly consistent between lakes in this
ecoregion.
4.2 FUTURE RESEARCH DIRECTIONS
This study was designed to assess the effects of a subdecadal climate shift on
adult fish population and zooplankton assemblages in small endorheic lakes, and evaluate
the potential of these lakes as a model system for ecosystem responses to climate change.
While the environmental conditions in my study lakes generally encompasses a large
variability for the most important parameters, a larger sample size with more evenly
distributed values for depth, surface area and water chemistry parameters could establish
that factors other than temperature and salinity are also relevant for determining fish
health in these prairie lakes. For example, lakes with sufficient depth could eliminate
both winterkill risk and thermal variation that dominate smaller lakes, potentially
revealing more subtle underlying effects of water chemistry or seasonal thermal
stratification.
Further, the sub-lethal effects on fish identified here were only examined for
mature adults. Sub-lethal and partially lethal effects are known to vary across life stages
and are often more pronounced at earlier life stages (Mélard et al. 1996). In particular,
salinity, suspended sediments or altered temperatures can substantially reduce egg
hatching rates (Koel and Peterka 1995, Paxton et al. 2004). Tapeworm parasite load is
likely an exception, as parasite loads tend to increase with age (Zelmer and Arai 1998).
Page 76
63
Also, a study of all age-classes over multiple years could reveal generational effects that
may cause significant demographic shifts over time (Ficke et al. 2007).
This research is helpful for understanding the nature of climate response in
climate-stressed lakes. Yet, for better lake management, it would be important to develop
quantitative models the future responses of prairie lake communities to climate change.
This will require a longer data set as well as an assessment of fish biomass to
complement the existing more detailed zooplankton data. Although it was apparent that
drought conditions more accurately predicted species composition, the absence of reliable
biomass and abundances of fished prevented a more detailed analysis of population
response to climate. Shifts in relative abundances of species, as well as shifts in biomass
and increases or declines in recruitment numbers are well-documented in other systems
(Christie and Regier 1988, Pangle et al. 2004, Paxton et al. 2004), but little information is
currently available on this topic for these ecosystems. Once a better data structure exists,
a temporal predator-prey modeling approach such as MAR could be an effective option
to predict trophic interactions between piscivores, planktivores and zooplankton for
climate-change scenarios (Hampton et al. 2006).
4.3 GENERAL CONCLUSIONS
Endorheic prairie lakes served as an excellent model system due to their high
hydrologic sensitivity to climate and tendency to regularly experience cyclical climate
extremes (Last and Ginn 2005). This allowed for the detection of significant climate
responses over a relatively brief (intradecadal) period of time. Similarly, interannual
variation in temperature and winter length was substantial enough to reveal its sub-lethal
Page 77
64
effects on fish populations. This high degree of natural variability provides valuable
information for predicting effects of future environmental change.
Furthermore, this study successfully encapsulated environmental and trophic
effects at the community to intraspecific levels. Fish communities were assessed based
on trophic position, while zooplankton were evaluated in terms of species abundances,
and sub-lethal effects were assessed for three piscivorous fish species. Such a
comprehensive analysis facilitates evaluation of the relative importance of these factors
and how they may act in concert.
This study revealed that observations about ecosystem response to disturbances
may not be universal. In this system, only fishes were directly affected by water
chemistry changes, while zooplankton were primarily regulated by predation. Further,
lake chemistry parameters differed more during drought, while many other systems
responded synchronously in times of environmental stress (Vogt et al. 2011, Webster et
al. 2000). Also, lake communities differed more during drought, suggesting that
landscape-level diversity is heightened by disturbance in this system, despite typically
being lowered by systematic disturbance in others (Balata et al. 2007, Chase 2007).
Aside from these differences, it is clear that this system, like many others, would
likely be detrimentally affected by changes in climate, resulting in reduced trophic
complexity and lower species diversity within fish and zooplankton communities.
Similarly, fish communities would experience physical stress in more extreme
environments. Conservation of these valued and threatened species is therefore essential
across all systems; however, differences in the underlying mechanisms of change require
that model systems be identified.
Page 78
65
REFERENCES
Adams, S. M., and R. B. McLean. 1985. Estimation of largemouth bass, Micropterus
salmoides Lacepede, growth using the liver somatic index and physiological
variables. Journal of Fish Biology 26:111–126.
Adrian, R., C. M. O’Reilly, H. Zagarese, S. B. Baines, D. O. Hessen, W. Keller, D. M.
Livingstone, R. Sommaruga, D. Straile, E. Van Donk, G. A. Weyhenmeyer, and M.
Winder. 2009. Lakes as sentinels of climate change. Limnology and Oceanography
54:2283–2297.
Anderson, D. P. 1997. 7 Environmental factors in fish health: Immunological aspects.
Fish Physiology 15:289–310.
Baines, S., and M. Pace. 1991. The production of dissolved organic matter by
phytoplankton and its importance to bacteria: Patterns across marine and freshwater
systems. Limnology and Oceanography 36:1078–1090.
Balata, D., L. Piazzi and L. Benedetti-Cecchi. 2007. Sediment disturbance and loss of
beta diversity on subtidal rocky reefs. Ecology 88: 2455-2461.
Barica, J., and J. A. Mathias. 1979. Oxygen depletion and winterkill risk in small prairie
lakes under extended ice cover. Journal of the Fisheries Research Board of Canada
36:980–986.
Barrow, E. 2009. Climate scenarios for Saskatchewan. 131 pp. Prairie Adaptation
Research Collaborative (PARC), Regina.
Bates, B., Z. W. Kundzewicz, S. Wu, and J. Palutikof. 2008. Climate change and water.
200 pp. Intergovernmental Panel on Climate Change (IPCC), Geneva.
Battarbee, R., J.-A. Grytnes, and R. Thompson et al. 2002. Comparing paleolimnological
and instrumental evidence of climate change for remote mountain lakes over the last
200 years. Journal of Paleolimnology 28:161–179.
Bayly, I. A. E. 1972. Salinity tolerance and osmotic behavior of animals in athalassic
saline and marine hypersaline waters. Annual Review of Ecology and Systematics
3:233–268.
Beaudoin, C., W. Tonn, E. Prepas and L. Wassenaar. 1999. Individual specialization
and trophic adaptability of northern pike (Esox lucius): an isotope and dietary
analysis. Oecologia 120: 386-396.
Blackwell, B. G., M. L. Brown, and D. W. Willis. 2000. Relative weight (Wr) status and
current use in fisheries assessment and management. Reviews in Fisheries Science
8:1–44.
Page 79
66
Blenckner, T., A. Omstedt, and M. Rummukainen. 2002. A Swedish case of
contemporary and possible future consequences of climate change on lake function.
Aquatic Sciences 64:171–184.
Bly, J. E., S. M. Quiniou, and L. W. Clem. 1997. Environmental effects on fish immune
mechanisms. Developments in biological standardization 90:33–43.
Bond, N. R., P. S. Lake, and A. H. Arthington. 2008. The impacts of drought on
freshwater ecosystems: an Australian perspective. Hydrobiologia 600:3–16.
Bonsal, B. R., and E. E. Wheaton. 2005. Atmospheric circulation comparisons between
the 2001 and 2002 and the 1961 and 1988 Canadian prairie droughts. Atmosphere-
Ocean 43:163–172.
Bonsal, B. R., E. E. Wheaton, A. C. Chipanshi, C. Lin, D. J. Sauchyn, and L. Wen. 2011.
Drought research in Canada: A review. Atmosphere-Ocean 49:303–319.
Bredehoeft, J., and T. Durbin. 2009. Ground water development--the time to full capture
problem. Ground Water 47:506–514.
Burnham, B. L., and J. J. Peterka. 1975. Effects of saline water from North Dakota lakes
on survival of fathead minnow (Pimephales promelas) embryos and sac fry. Journal
of the Fisheries Research Board of Canada 32:809–812.
Cargnelli, L. M., and M. R. Gross. 1997. Notes: Fish energetics: Larger individuals
emerge from winter in better condition. Transactions of the American Fisheries
Society 126:153–156.
Cena, C. J., G. E. Morgan, M. D. Malette, and D. D. Heath. 2006. Inbreeding,
outbreeding and environmental effects on genetic diversity in 46 walleye (Sander
vitreus) populations. Molecular Ecology 15:303–320.
Chapman, L. J., W. C. Mackay, and C. W. Wilkinson. 1989. Feeding flexibility in
northern pike (Esox lucius): Fish versus Invertebrate Prey. Canadian Journal of
Fisheries and Aquatic Sciences 46:666–669.
Chase, J. M. 2007. Drought mediates the importance of stochastic community assembly.
Proceedings of the National Academy of Sciences of the United States of America
104:17430–4.
Cherry, D. S., K. L. Dickson, J. Cairns Jr., and J. R. Stauffer. 1977. Preferred, avoided,
and lethal temperatures of fish during rising temperature conditions. Journal of the
Fisheries Research Board of Canada 34:239–246.
Page 80
67
Christie, G. C., and H. A. Regier. 1988. Measures of optimal thermal habitat and their
relationship to yields for four commercial fish species. Canadian Journal of Fisheries
and Aquatic Sciences 45:301–314.
Chu, C., N. E. Mandrak, and C. K. Minns. 2005. Potential impacts of climate change on
the distributions of several common and rare freshwater fishes in Canada. Diversity
and Distributions 11:299–310.
Cooper, R., and B. Wissel. 2012a. Loss of trophic diversity in saline prairie lakes as
indicated by stable-isotope based community metrics. Aquatic Biosystems 8:6.
Cooper, R. and B. Wissel. 2012b. Interactive effects of chemical and biological controls
on food-web composition in saline prairie lakes. Aquatic Biosystems 8: 29.
Craig, J. F., and J. A. Babaluk. 1989. Relationship of condition of walleye (Stizostedion
vitreum) and northern pike (Esox lucius) to water clarity, with special reference to
Dauphin Lake, Manitoba. Canadian Journal of Fisheries and Aquatic Sciences
46:1581–1586.
Dicenzo, V. J., M. J. Maceina, and W. C. Reeves. 1995. Factors related to growth and
condition of the Alabama subspecies of spotted bass in reservoirs. North American
Journal of Fisheries Management 15:794–798.
Drinkwater, K. 2005. The response of Atlantic cod (Gadus morhua) to future climate
change. ICES Journal of Marine Science 62:1327–1337.
Dutil, J.-D., and Y. Lambert. 2011. Natural mortality from poor condition in Atlantic cod
(Gadus morhua). Canadian Journal of Fisheries and Aquatic Sciences 57:826–836.
Eaton, J. G., and R. M. Scheller. 1996. Effects of climate warming on fish thermal habitat
in streams of the United States. Limnology and Oceanography 41:1109–1115.
Fang, X., H. G. Stefan, J. G. Eaton, J. H. McCormick, and S. R. Alam. 2004. Simulation
of thermal/dissolved oxygen habitat for fishes in lakes under different climate
scenarios. Ecological Modelling 172:13–37.
Fee, E. J., R. E. Hecky, S. E. M. Kasian, and D. R. Cruikshank.1996. Effects of lake size,
water clarity, and climatic variability on mixing depths in Canadian Shield lakes.
Limnology and oceanography 41:912–920.
Fenton, R., J. A. Mathias, and G. E. E. Moodie. 1996. Recent and future demand for
walleye in North America. Fisheries 21:6–12.
Ficke, A. D., C. A. Myrick, and L. J. Hansen. 2007. Potential impacts of global climate
change on freshwater fisheries. Reviews in Fish Biology and Fisheries 17:581–613.
Page 81
68
Fisher, S. J., D. W. Willis, and K. L. Pope. 1996. An assessment of burbot ( Lota lota )
weight–length data from North American populations. Canadian Journal of Zoology
74:570–575.
Fritz, S. C. 1996. Paleolimnological records of climatic change in North America.
DigitalCommons@University of Nebraska, Lincoln.
Gabelhouse, D. W. 1991. Seasonal changes in body condition of white crappies and
relations to length and growth in Melvern Reservoir, Kansas. North American
Journal of Fisheries Management 11:50–56.
Gleick, P. H. 1989. Climate change, hydrology, and water resources. Reviews of
Geophysics 27:329.
Hall, R. I., P. R. Leavitt, R. Quinlan, A. S. Dixit, and J. P. Smol. 1999. Effects of
agriculture, urbanization, and climate on water quality in the northern Great Plains.
Limnology and Oceanography 44:739–756.
Hammer, U. 1990. The effects of climate change on the salinity, water levels and biota of
Canadian prairie saline lakes. Internationale Vereinigung fuer Theoretische und
Angewandte Limnologie 24:321–326.
Hammer, U. T. 1978. The Saline Lakes of Saskatchewan III. Chemical Characterization.
Internationale Revue der gesamten Hydrobiologie und Hydrographie 63:311–335.
Hammer, U. T. 1986. Saline Lake Ecosystems of the World, Volume 59 (Google eBook).
Page 632. . Dr. W. Junk Publishers, Dordrecht.
Hampton, S. E., M. D. Scheuerell, and D. E. Schindler. 2006. Coalescence in the Lake
Washington story: Interaction strengths in a planktonic food web. Limnology and
Oceanography 51:2042–2051.
Hanley, J. et al. 2010. An-image-based Key to the Zooplankton of the Northeast.
Hauton, C., L. E. Hawkins, and S. Hutchinson. 2000. The effects of salinity on the
interaction between a pathogen (Listonella anguillarum) and components of a host
(Ostrea edulis) immune system. Comparative Biochemistry and Physiology Part B:
Biochemistry and Molecular Biology 127:203–212.
Hayward, R. and F. Margraf. 1987. Eutrophication effects on prey size and food
available to yellow perch in Lake Erie. Transactions of the American Fisheries
Society 116: 210-223.
Heisler, J., P. M. Glibert, J. M. Burkholder, D. M. Anderson, W. Cochlan, W. C.
Dennison, Q. Dortch, C. J. Gobler, C. a. Heil, E. Humphries, a. Lewitus, R. Magnien,
H. G. Marshall, K. Sellner, D. a. Stockwell, D. K. Stoecker, and M. Suddleson. 2008.
Page 82
69
Eutrophication and harmful algal blooms: A scientific consensus. Harmful Algae
8:3–13.
Heitmann, M. L. 1973. Stability of thermal stratification in a lake. Hydrology of Lakes:
321–326.
Hudson, P. L., and L. T. Lesko. 2003. Free-living and Parasitic Copepods of the
Laurentian Great Lakes: Keys and Details on Individual Species. Ann Arbor, MI:
Great Lakes Science Centre Home Page. URL:
http//:www.glsc.usgs.gov/greatlakescopepods/
Hughes, L. 2003. Climate change and Australia: Trends, projections and impacts. Austral
Ecology 28:423–443.
Huh, T. H., H. E. Calbert, and D. A. Stuiber. 1976. Effects of temperature and light on
growth of yellow perch and walleye using formulated feed. Transactions of the
American Fisheries Society 105:254–258.
Jeppesen, E., B. Kronvang, M. Meerhoff, M. Søndergaard, K. M. Hansen, H. E.
Andersen, T. L. Lauridsen, L. Liboriussen, M. Beklioglu, A. Ozen, and J. E. Olesen.
2007. Climate change effects on runoff, catchment phosphorus loading and lake
ecological state, and potential adaptations. Journal of environmental quality
38:1930–41.
van der Kamp, G., D. Keir, and M. S. Evans. 2008. Long-term water level changes in
closed-basin lakes of the Canadian prairies. Canadian Water Resources Journal
33:23–38.
King, J., B. Shuter, and A. Zimmerman. 1999. Empirical links between thermal habitat,
fish growth, and climate change. Transactions of the American Fisheries Society
128:656–665.
Kitchell, J., and D. J. Stewart. 1977. Applications of a bioenergetics model to yellow
perch (Perca flavescens) and walleye (Stizostedion vitreum vitreum). Journal of the
Fisheries Research Board of Canada 34:1922–1935.
Kling, G. 1988. Comparative transparency, depth of mixing, and stability of stratification
in lakes of Cameroon, West Africa. Limnology and Oceanography 33:27–40.
Knight, R. L., F. J. Margraf, and R. F. Carline. 1984. Piscivory by walleyes and yellow
perch in western Lake Erie. Transactions of the American Fisheries Society
113:677–693.
Kocovsky, P. and R. Carline. Influence of extreme temperatures on consumption and
condition of walleyes in Pymatuning Sanctuary, Pennsylvania. North American
Journal of Fisheries Management 21: 198-207.
Page 83
70
Koel, T. M., and J. J. Peterka. 1995. Survival to hatching of fishes in sulfate-saline waters,
Devils Lake, North Dakota. Canadian Journal of Fisheries and Aquatic Sciences
52:464–469.
Larocque, I., R. I. Hall, and E. Grahn. 2001. Chironomids as indicators of climate change :
a 100-lake training set from a subarctic region of northern Sweden ( Lapland ).
Journal of Paleolimnology 26:307–322.
Last, W. M., and F. M. Ginn. 2005. Saline systems of the Great Plains of western Canada:
an overview of the limnogeology and paleolimnology. Saline systems 1:1–10.
Leavitt, P. R., S. C. Fritz, N. Anderson, P. Baker, T. Blenckner, L. Bunting, J. Catalan, D.
Conley, W. Hobbs, E. Jeppesen, A. Korhola, S. McGowan, K. Ruhland, J. Rusak, G.
Simpson, N. Solovieva, and J. Werne. 2009. Paleolimnological evidence of the
effects on lakes of energy and mass transfer from climate and humans. Limnology
and Oceanography 54:2330–2348.
Lewis, J., D. Hoole, and L. H. Chappell. 2003. Parasitism and environmental pollution:
parasites and hosts as indicators of water quality. Parasitology 126:S1–S3.
Li, J., Y. Cohen, D. Schupp, and I. Adelman. 1996. Effects of walleye stocking on
population abundance and fish size. North American Journal of Fisheries
Management 16: 830-839.
Magnuson, J. J., J. D. Meisner, and D. K. Hill. 1990. Potential changes in the thermal
habitat of Great Lakes fish after global climate warming. Transactions of the
American Fisheries Society 119:254–264.
Magoulick, D. D., and R. M. Kobza. 2003. The role of refugia for fishes during drought:
a review and synthesis. Freshwater Biology 48:1186–1198.
Marchildon, G. P., S. Kulshreshtha, E. Wheaton, and D. Sauchyn. 2007. Drought and
institutional adaptation in the Great Plains of Alberta and Saskatchewan, 1914–1939.
Natural Hazards 45:391–411.
Marcogliese, D. J. 2002. Food webs and the transmission of parasites to marine fish.
Parasitology 124:83–99.
Marcogliese, D. J., L. G. Brambilla, F. Gagné, and A. D. Gendron. 2005. Joint effects of
parasitism and pollution on oxidative stress biomarkers in yellow perch Perca
flavescens. Diseases of aquatic organisms 63:77–84.
Marwitz, T. D., and W. A. Hubert. 1997. Trends in relative weight of walleye stocks in
wyoming reservoirs. North American Journal of Fisheries Management 17:44–53.
Page 84
71
McCune, B., and M. Mefford. 1999. PC-ORD: Multivariate analysis of ecological data. .
MJM Software Design, Gleneden Beach, OR.
McCune, B., J. B. Grace, and D. L. Urban. 2002. Analysis of ecological communities. .
MJM Software Design, Gleneden Beach, OR.
McGowan, S., P. R. Leavitt, and R. I. Hall. 2005. A whole-lake experiment to determine
the effects of winter droughts on shallow lakes. Ecosystems 8:694–708.
Mehner, T., M. Diekmann, U. Bramick, and R. Lemcke. 2005. Composition of fish
communities in German lakes as related to lake morphology, trophic state, shore
structure and human-use intensity. Freshwater Biology 50:70–85.
Mélard, C., P. Kestemont, and J. C. Grignard. 1996. Intensive culture of juvenile and
adult Eurasian perch (P. fluviatilis): effect of major biotic and abiotic factors on
growth. Journal of Applied Ichthyology 12:175–180.
Michels, A., K. R. Laird, S. E. Wilson, D. Thomson, P. R. Leavitt, R. J. Oglesby, and B.
F. Cumming. 2007. Multidecadal to millennial-scale shifts in drought conditions on
the Canadian prairies over the past six millennia: implications for future drought
assessment. Global Change Biology 13:1295–1307.
Millett, B., W. C. Johnson, and G. Guntenspergen. 2009. Climate trends of the North
American prairie pothole region 1906–2000. Climatic Change 93:243–267.
Murphy, B., M. Brown and T. Springer. 1990. Evaluation of relative weight (Wr) index,
with new applications to walleye. North American Journal of Fisheries
Management 10: 85-97.
Murphy, B., D. Willis, and T. Springer. 1991. The relative weight index in fisheries
management: status and needs. Fisheries 16:30–38.
Muzzall, P. M., and R. C. Haas. 1998. Parasites of Walleyes, Stizostedion vitreum, from
Saginaw Bay, Lake Huron, and the Other Great Lakes. Journal of Great Lakes
Research 24:152–158.
Nielsen, D. L., M. a. Brock, G. N. Rees, and D. S. Baldwin. 2003. Effects of increasing
salinity on freshwater ecosystems in Australia. Australian Journal of Botany 51:655.
O’Reilly, C. M., S. R. Alin, P.-D. Plisnier, A. S. Cohen, and B. a McKee. 2003. Climate
change decreases aquatic ecosystem productivity of Lake Tanganyika, Africa.
Nature 424:766–8.
Ostrand, K. G., and G. R. Wilde. 2001. Temperature, dissolved oxygen, and salinity
tolerances of five prairie stream fishes and their role in explaining fish assemblage
patterns. Transactions of the American Fisheries Society 130:742–749.
Page 85
72
Pangle, K. L., T. M. Sutton, R. E. Kinnunen, and M. H. Hoff. 2004. Overwinter survival
of juvenile lake herring in relation to body size , physiological condition , energy
stores , and food ration. Transactions of the American Fisheries Society 133:1235–
1246.
Parmesan, C., G. Yohe, and others. 2003. A globally coherent fingerprint of climate
change impacts across natural systems. Nature 421:37–42.
Patz, J. a, T. K. Graczyk, N. Geller, and a Y. Vittor. 2000. Effects of environmental
change on emerging parasitic diseases. International journal for parasitology
30:1395–405.
Paxton, C. G. M., I. J. Winfield, J. M. Fletcher, D. G. George, and D. P. Hewitt. 2004.
Biotic and abiotic influences on the recruitment of male perch in Windermere, U.K.
Journal of Fish Biology 65:1622–1642.
Pham, S. V., P. R. Leavitt, S. McGowan, and P. Peres-Neto. 2008. Spatial variability of
climate and land-use effects on lakes of the northern Great Plains. Limnology and
Oceanography 53:728–742.
Pham, S. V., P. R. Leavitt, S. McGowan, B. Wissel, and L. I. Wassenaar. 2009. Spatial
and temporal variability of prairie lake hydrology as revealed using stable isotopes
of hydrogen and oxygen. Limnology and Oceanography 54:101–118.
Pinel-Alloul, B., T. Niyonsenga, and P. Legendre. 1995. Spatial and environmental
components of freshwater zooplankton structure. Ecoscience 2:1–19.
Pollock, M., Z. Hoover, M. Ferrari, D. Chivers, and G. McMaster. 2010. Tolerance of
juvenile northern pike (Esox lucius), juvenile walleye (Sander vitreus), adult fathead
minnows (Pimephales promelas) and zooplankton to representative environment
levels of total dissolved solids. Page 51.
Poole, B., and T. Dick. 1985. Parasite recruitment by stocked walleye, Stizostedion
vitreum vitreum (Mitchill), fry in a small boreal lake in central Canada. Journal of
Wildlife Diseases 21:371–376.
Poulin, R. 1992. Toxic pollution and parasitism in freshwater fish. Parasitology Today
8:58–61.
Quinlan, R., P. R. Leavitt, A. S. Dixit, R. I. Hall, and J. P. Smol. 2002. Landscape effects
of climate, agriculture, and urbanization on benthic invertebrate communities of
Canadian prairie lakes. Limnology and Oceanography 47:378–391.
Quist, M. 2002. Seasonal variation in condition, growth and food habits of walleye in a
Great Plains reservoir and simulated effects of an altered thermal regime. Journal of
Fish Biology 61:1329–1344.
Page 86
73
Rawson, D. S., and J. E. Moore. 1944. The saline lakes of Saskatchewan. Canadian
Journal of Research 22d:141–201.
Robarts, R. D., M. J. Waiser, M. T. Arts, and M. S. Evans. 2005. Seasonal and diel
changes of dissolved oxygen in a hypertrophic prairie lake. Lakes and Reservoirs:
Research and Management 10:167–177.
Sharma, S., D. A. Jackson, C. K. Minns, and B. J. Shuter. 2007. Will northern fish
populations be in hot water because of climate change? Global Change Biology
13:2052–2064.
De Stasio, B. T., D. K. Hill, J. M. Kleinhans, N. P. Nibbelink, and J. J. Magnuson. 1996.
Potential effects of global climate change on small north-temperate lakes : Physics,
fish, and plankton. Limnology and oceanography 41:1136–1149.
Sauchyn, D., M. Khandekar, E. R. Garnett, P. M. August, and A. Pennylegion. 2005. The
Science , Impacts and Monitoring of Drought in Western Canada. The Prairie
Adaptation Research Collaborative ( PARC ) and University of Regina, Regina, SK.
Saxon, E., B. Baker, W. Hargrove, F. Hoffman, and C. Zganjar. 2005. Mapping
environments at risk under different global climate change scenarios. Ecology
Letters 8:53–60.
Schindler, D., S. Bayley, and B. Parker et al. 1996. The effects of climatic warming on
the properties of boreal lakes and streams at the Experimental Lakes Area,
Northwestern Ontario. Limnology and Oceanography 41:1004–1017.
Schindler, D. W. 1998. A dim future for boreal waters and landscapes. BioScience
48:157–164.
Schindler, D. W. 2001. The cumulative effects of climate warming and other human
stresses on Canadian freshwaters in the new millennium. Canadian Journal of
Fisheries and Aquatic Sciences 58:18–29.
Scott, W., and E. Crossman. 1973. Freshwater fishes of Canada. Fisheries Research
Board of Canada Bulletin 184.
Srinivasan, U. T., J. A. Dunne, J. Harte, and N. D. Martinez. 2007. Response of complex
food webs to realistic extinction sequences. Ecology 88:671–682.
Stainton, M., C. MJ, and F. Armstrong. 1974. The chemical analysis of fresh water. Page
120. . Research and Development Directorate, Freshwater Institute, Winnipeg.
Vander Zanden, M., and Y. Vadeboncoeur. 2002. Fishes as integrators of benthic and
pelagic food webs in lakes. Ecology 83:2152–2161.
Page 87
74
Victoria, C. J., B. S. Wilkerson, R. J. Klauda, and E. S. Perry. 1992. Salinity tolerance of
yellow perch eggs and larvae from coastal plain stream populations in maryland,
with comparison to a pennsylvania lake population. Copeia 1992:859–865.
Vogt, R. J., J. A. Rusak, A. Patoine, and P. R. Leavitt. 2011. Differential effects of
energy and mass influx on the landscape synchrony of lake ecosystems. Ecology
92:1104–14.
Waiser, M. and R. Robarts. 1995. Microbial nutrient limitation in prairie saline lakes
with high sulfate concentration. Limnology and Oceanography 40: 566-574.
Webster, K. E., P. A. Soranno, S. B. Baines, T. K. Kratz, C. J. Bowser, P. J. Dillon, P.
Campbell, E. J. Fee, and R. E. Hecky. 2000. Structuring features of lake districts:
landscape controls on lake chemical responses to drought. Freshwater Biology
43:499–515.
Wetzel, R. G., and G. E. Likens. 1996. Limnological Analyses. 456 pp., 2nd edition.
Springer-Verlag, New York.
Williams, W. D. 2002. Environmental threats to salt lakes and the likely status of inland
saline ecosystems in 2025. Environmental Conservation 29:154–167.
Williamson, C. E., D. P. Morris, M. L. Pace, and O. G. Olson. 1999. Dissolved organic
carbon and nutrients as regulators of lake ecosystems: resurrection of a more
integrated paradigm. Limnology and Oceanography 44:759–803.
Williamson, C. E., J. E. Saros, W. F. Vincent, and J. P. Smol. 2009. Lakes and reservoirs
as sentinels, integrators, and regulators of climate change. Limnology and
Oceanography 54:2273–2282.
Winder, M., and D. E. Schindler. 2004. Climate change uncouples trophic interactions in
an aquatic ecosystem. Ecology 85:2100–2106.
Winter, T., and D. Rosenberry. 1998. Hydrology of prairie pothole wetlands during
drought and deluge: a 17-year study of the cottonwood lake wetland complex in
north dakota in the perspective of longer term measured and proxy hydrological
records. Climatic Change 40:189–209.
Wissel, B., R. N. Cooper, P. R. Leavitt, and S. V. Pham. 2011. Hierarchical regulation of
pelagic invertebrates in lakes of the northern Great Plains: a novel model for
interdecadal effects of future climate change on lakes. Global Change Biology
17:172–185.
Wright, M. and M. Curtis. 2000. Temperature effects on embryonic development and
the life cycle of Diphyllobothrium dendriticum. International Journal of
Parasitology 30: 849-852.
Page 88
75
Zalewski, M., B. Brewinska-Zaraś, P. Frankiewicz, and S. Kalinowski. 1990. The
potential for biomanipulation using fry communities in a lowland reservoir:
concordance between water quality and optimal recruitment. Hydrobiologia
200:549–556.
Zelmer, D. A., and H. P. Arai. 1998. The contributions of host age and size to the
aggregated distribution of parasites in yellow perch, Perca flavescens, from Garner
Lake, Alberta, Canada. The Journal of parasitology 84:24–8.