BioRisk e 243-25 | (20 | 0) Mb achest Es emcees doi: 10.3897/biorisk.5.853 RESEARCH ARTICLE & B lO R IS kK http://biorisk-journal.com/ The local species richness of Dragonflies in mountain waterbodies: an indicator of climate warming? Beat Oertli University of Applied Sciences Western Switzerland, Ecole d’Ingénieurs HES de Lullier, 150 route de Presinge, CH — 1254 Jussy / Geneva, Switzerland Corresponding author: Beat Oertli (beat.oertli@hesge.ch) Academic editor: /iirgen Ott | Received 29 July 2010 | Accepted 22 September 2010 | Published 30 December 2010 Citation: Oertli B (2010) The local species richness of Dragonflies in mountain waterbodies: an indicator of climate warming? In: Ott J (Ed) (2010) Monitoring Climatic Change With Dragonflies. BioRisk 5: 243-251. doi: 10.3897/ biorisk.5.853 Abstract With climate warming, many Odonata species are extending their geographical area. In Switzerland, as in many parts of the world, this phenomenon may lead to a regional increase in species richness. The local richness (the richness of individual waterbodies) is also expected to increase, particularly in the alpine or subalpine areas where the waterbodies are particularly species—poor. Based on the species richness recorded in 109 waterbodies scattered all across Switzerland, a model is presented here relating the local species richness (adult dragonflies) to environmental variables, including the mean annual air temperature. This model predicts a sharp increase in species richness for alpine or subalpine waterbodies, which is expected to double or even treble before the end of this century. This increase would mainly be the consequence of the immigration of eurythermal species extending their geographical range, together with potential local extinctions of the cold stenothermal species. Keywords biodiversity, Odonata, alpha richness, boreo-alpine species, alpine ponds, colonisation, extinction Introduction Global changes, and particularly climate warming, are likely to have various impacts on Odonata (see other chapters of this book), such as a shift in species’ geographical distribution, an earlier timing of emergence, or a shorter duration of the aquatic live Copyright Beat Oertli. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 244 Beat Oertli/ BioRisk 5: 243-251 (2010) stage. In the northern hemisphere, the northerly shift in the geographical distribu- tion of several species will lead in many areas to an enrichment of the regional species pool, mainly resulting from an increase in the number of eurythermal species. Such an increase will largely encompass the associated decrease in the number of cold stenother- mal species. A regional increase in species richness has already been reported in many European states such as Switzerland, Germany, Netherlands, Belgium, or UK. Similar changes are also expected to occur with altitude. Indeed, the elevation gradient is often claimed to mirror the latitudinal gradient (Rahbek 1995), hence conditions at higher altitudes resemble conditions at higher latitudes (Ricklefs 1990). If the regional increase in species richness is now an evidence, changes in the local richness (richness of a given ecosystem, e.g. pond, stream, wetland) were at the moment poorly documented (but see the example of permanent vegetation plots, in Pauli et al. 2007). Nevertheless, such changes would undoubtedly occur, especially in areas where the local richness is low. This is particularly the case in altitude (or high latitude), where local species richness is clearly lower than in lowlands (or lower latitude) (e.g. Oertli et al. 2008). Many alpine ponds in Switzerland do not presently host any Odonata species (or only a few), whereas it is common to observe ten to twenty species in lowland ponds. The present case study aimed to predict the expected changes in the local species richness as a consequence of various scenarios of temperature elevation. For this pur- pose, the species richness of adults Odonata was assessed in 109 ponds distributed in all altitudinal regions (and therefore thermal regions) of Switzerland. This data set is used to relate the species richness of the ponds with the mean annual air temperature. Material and methods Measure of local species richness and environmental variables A set of 109 ponds scattered in the four altitudinal belts of Switzerland (collinean, montane, subalpine, and alpine; Figure 1) was investigated between 1996 and 2002. This set represents a broad range of thermal conditions with mean annual air tempera- ture ranging from —2.2°C to 12.1°C (Figure 2). The sampling of adult Odonata and the assessment of species richness were standardised and followed the PLOCH method (Oertli et al. 2005), for which a representative species list is gathered during two sam- pling days (at the end of spring and summer). To account for sample-size inadequacies (i.e. it is unlikely to perform the same sampling effort on each pond), the measured species richness is further corrected by means of Chao estimator of true richness (see Oertli 2008; Oertli et al. 2002; 2005). About a hundred of environmental variables were also measured, at the local and regional scales, including water physico-chemistry (e.g. pH, nutrients, conductibility, transparency), morphometry (e.g. pond area, depth), watershed characteristics (geol- ogy, land use), connectivity (isolation from other waterbodies), and climatic data (air Local species richness as indicator of climate warming 245 A alpine A subalpine O montane © collinean Figure |. Distribution of the 109 sampled ponds in Switzerland, with indication of their altitudinal belt. 25 20 15 10 number of ponds <0 Oto1 1to2 2to3 3to4 4to5 5106 6to7 7to8 8to9 Yto >10 é 10 annual air temperature (°C) Figure 2. Range of mean annual air temperature covered by the 109 sampled ponds. Annual mean air temperature was calculated on the basis of monthly values from 115 climate stations and of a digital terrain model with a 25-m grid (data from the Swiss Federal Institute for Forest, Snow and Landscape Research). temperature, quantity of precipitations, cloud cover, evapotranspiration, solar radia- tion). For additional information on sampling procedure and pond characteristics, see Oertli et al. (2002, 2005). 246 Beat Oertli/ BioRisk 5: 243-251 (2010) Statistical analyses > «99 In a first descriptive step, correlation coefhcients (Spearman's “o”) served to assess the relationships between species richness and mean annual air temperature, and each en- vironmental variable. A stepwise linear regression procedure (LR) was then used to model the relation- ship between estimated Odonata species richness and the annual mean air temperature (also including the other significant environmental variables), by setting the prob- ability to 0.05. The model was used to predict the changes in species richness for two “virtual” alpine and subalpine ponds (each one with a surface area of 3300 m7”) in response to global change. ‘The typical alpine pond had a mean annual air temperature of 0.6°C and the typical subalpine pond of 3.4°C. The temperature increases ranged between +1.8°C to +4.0°C, and were based on seven emission scenarios: the six IPCC scenarios for 2090-2100 (IPCC 2007) and the CH2070 scenario for Switzerland in 2070 (Hohmann et al. 2007). Results Relation between pond species richness and mean annual air temperature From the correlation analysis, it resulted that species richness best correlated with, first, mean annual air temperature (r= 0.72), and then with altitude (r= - 0.67). These two variables were highly correlated (r= 0.99). Other important variables were identified by calculating their correlation with the residuals of the linear regression between species richness and annual air temperature. This analysis highlighted a significant relation with pond area (r= 0.47), proportion of the environment covered by forest (r= - 0.32), fish presence (r= 0.30), and mean pond depth (r= 0.23). Based on these preliminary analyses, a selection of 15 variables was used to build the model (stepwise LR) relating the Odonata species richness to the air tempera- ture: mean annual air temperature, pond area, pond age, mean depth, shore devel- opment, water conductivity, water transparency, pond eutrophication state, propor- tion of pond shaded, fish presence, proportion of pond area covered by floating vegetation or by submerged vegetation, proportion of agricultural landcover in the catchment, and proportion of the environment covered by forest and the connectiv- ity in a radius of 1000 m. The final model selected two variables, mean annual temperature and pond area, and showed that species richness (S) increased with an increase of temperature (T, in °C) or/and pond area (A, in m’): S = -4.65 + 1.34*T° + 1.39 * log(A) r= 0.51 p<0.0001 Local species richness as indicator of climate warming 247 Prediction of the local species richness of mountain ponds as a function of differ- ent scenarios of temperature increase The resulting model was used for the simulation of the impact of temperature increases on local species richness. For a given alpine or subalpine pond (fixed area of 3300 m’), we tested seven scenarios (Figure 3) representative of the trends predicted for the next 65 to 95 years. All seven scenarios evidenced a clear increase in Odonata species richness. The magnitude of this increase is particularly high, as the richness is likely to double, or even to treble in the case of alpine ponds. Whilst the species rich- ness currently observed in an alpine pond is one species, this is expected to increase to three to six species during the next decades. In the case of a subalpine pond, species richness may rise from six (mean value currently observed) to 8-11 species. Discussion The predictions of climate warming presented here point out future drastic changes in the species richness of mountain waterbodies. For alpine ponds, the local richness of Odonata is likely to double or treble. This increase in local species richness of moun- tain waterbodies is the consequence of immigration events that will largely exceed ex- tinction events. Indeed, in Switzerland, the present Odonata species pool is composed of 72 indigenous species (Gonseth and Monnerat 2002), the majority being euryther- mal species associated with lowland waterbodies. Only seven species are restricted to altitude areas and can be considered as “cold stenothermal species”. With the warming of the climate, many eurythermal species will extend their geographic distribution to higher altitudes, and will therefore be able to colonise mountain ponds from which they are currently absent. The extension of the geographical area to higher latitude has already been reported for European Odonata by Ott (2007). This phenomenon is likely to have also occurred along the altitudinal gradient through the highest eleva- tions, though this has not been yet reported. Such altitudinal shift in geographical species distribution is also observed for plant species in the Alps (e.g. Pauli et al. 2007). Beside the immigration of eurythermal species in waterbodies, it is noteworthy that extinction of currently established species will probably also occur. As a conse- quence of warming, the seven cold stenothermal species could disappear from their living waterbodies. The lower altitudinal limit of their geographical distribution could increase and could consequently lead to an upward shift of their habitats; in the long term, this could lead to species loss at a regional scale. In Switzerland, the geographical area (and therefore the altitudinal distribution) of Odonata species is particularly well documented by the databases of the Swiss Centre for Fauna Cartography (CSCF) (see also the Swiss Atlas; Wildermuth et al. 2005). Odonata are absent of almost all high-elevation waterbodies (higher than 2500 m). Sixteen species are frequent above 1500 m (Figure 4), from which a first set of seven species can be considered as cold stenothermal due to their predominance in the upper 248 Beat Oertli/ BioRisk 5: 243-251 (2010) IPCC A‘FI S$ = -4.65 + 1.34 * T° + 1.39 * log(area) r IPCC A2 CH2070.. @ IPCC A1B_~@ IPCC B2/A1T. uw o IPcc Bg re) o jon actual e subalpine i) A 5 & IPCC AIF 2 IPCC A2 E A*® cH2070 > A IPCC A1B IPCC B2/A1T IPCC Bi actual alpine A 0 1 2 3 4 5 6 7 8 Annual air temperature (°C) Figure 3. Potential change in the Odonata species richness (S) for an alpine and a subalpine pond as predicted by the stepwise LR equation according to seven climate change scenarios (six IPCC emission scenarios and the CH2070 scenario). “T°” is the annual mean air temperature. The two virtual ponds are typical of subalpine and alpine altitudinal belts, with a fixed area of 3300 m”. For the alpine pond, present “T°” is 0.6°C, while it is 3.4°C for the subalpine pond elevation areas: Aeshna caerulea, A. juncea, A. subarctica, Coenagrion hastulatum, Leu- corrhinia dubia, Somatochlora alpestris, S. arctica. The most elevated altitudinal ranges are those of the boreo-alpine species Aeshna caerulea and Somatochlora alpestris, and are particularly narrow. Species belonging to this set are those presenting the highest risk of extinction at the local scale (waterbody), but also on the long range at the regional scale. A second set is composed by nine eurythermal species, characterised by a broad altitudinal distribution, mainly located below 1000 m. ‘This set includes the candidate species to the colonisation of mountain waterbodies, and will be responsible for the increase of the local species richness. The predictions of changes in local richness presented in this paper are potential values. They present a global trend that provides a baseline for describing biotic re- sponses to warming. The magnitude of the changes, however, may be slightly different than that described here, being either higher or lower. Indeed, the predictive models are based on the change of only one parameter affected by climate change, i.e. mean annual air temperature. Other variables could interact with temperature and either diminish or increase the magnitude of the predicted changes. For example, variables such as the seasonal timing of warming (i.e. winter or summer), the quantity and frequency of precipitations, the number of days of ice cover, or the radiation (e.g. UVs), are factors potentially able to interfere directly with temperature, or to lead to secondary changes (e.g. hydrology, productivity, water chemistry). Land use, which could have many secondary consequences, is a relevant factor that is likely to greatly vary in mountain area (Maurer et al. 2006). Besides environmental factors, biotic vari- ables should also be accounted for. The dynamics of the colonisation processes in both Local species richness as indicator of climate warming 249 Aeschna caerulea Aeschna juncea Aeschna subarctica Coenagrion hastulatum Leucorrhinia dubia 2000-2500m 51% 23% 10% 32% Ez ' ee —_—— —— —__— — —_ — ——_—__§ @ 100-2000 30% 4% TE a% Ed . = 1000-1500m 9% 18% 2% 51% 22% a c & 2 500-1000m 1% 11% 4% 8% = = 250-500m n= 123 n= 785 2% n=ig nadQ 2% n= 148 ‘Observation frequency [%] Somatochlora alpestris Somatochlora arctica 2000-2500m 44 5% E e@ 1500-2000m 48% 25% 2 (4 se 1000-1500m 5% 50% a a 2 500-1000m 17% a = 250-500m ns2ni n= 125 Observation frequency [%) Aeschna cyanea Aeschna grandis Cordulegaster bidentata Cordulia aenea Enallagma cyathigerum 2000-2500m Pac E @ 1800-2000m B% 3% 5% a% 1% a = £ — 1000-1500m 18% 13% 5g, z 3 a E 500-1000m 45% 41% 74% 25% 28% = . 250-500m n=7io 31% n= 107 44% n= 7965 | 97% n= 453 65% n= 349 65% Observation frequency [%] Lestes sponsa Libellula quadrimaculata Somatochiora metallica Sympeotrum danae 2000-2500m @ =«1500-2000m 3% 1% 20% 4% a € - = 1000-1500m 18% 5% a 15% z — oe a 2 500-1000 36% 40% 14% 40% b = 250-500m m= 125 45% n= 785 54% n= 299 58% n= 180 42% Observation frequency [%] Figure 4. Altitudinal distribution in Switzerland of a set of 16 species frequently observed at high alti- tude (above 1500 m). 4a The seven cold stenothermal species, expected to exhibit a decrease in their geo- graphical area (at risk of extinction on the long range). 4b The nine eurythermal species, likely to become more frequent at higher altitude. The data report observations of exuviae and subaldults only, two life «cD stages attesting the reproduction at the observed altitude. “n” indicates the number of observation. Data were gathered before 2003 (mainly between 1990 and 2002) by the Swiss Centre for Fauna Cartography. space and time could be very complex, as they are species-specific. The time at which the upward dispersal begins, or the speed at which a species extend its distributional area, depend on its own biological or ecological characteristics. At the assemblage scale, this could result in a space-time overlapping and changing pattern of species distribu- tions. As an example, Anisoptera are much more efficient fliers and active colonisers than Zygoptera. The latter group can nevertheless compensate this miss by an usually elevated number of individuals per population and by the passive dispersal by wind. Other biotic interactions, such as predation, may occur when an immigrant colonises a new pond. Local processes are therefore likely to present some resistance to immigra- tion and to delay the increase in species richness. With their small size and their relative simple community structure, ponds con- stitute ideal sentinel and early warning systems (De Meester et al. 2005). ‘This is par- ticularly true for alpine or subalpine ponds, characterised by species-poor communities (e.g. Oertli et al. 2007). Such systems should therefore be used for monitoring the 250 Beat Oertli/ BioRisk 5: 243-251 (2010) biotic impacts of climate changes. Monitoring should be conducted with a set of ponds situated at various altitudes and covering a large range of latitudinal scales, both in Switzerland and in other countries. Among other invertebrates, Odonata is certainly one of the most suitable groups for conducting monitoring, either alone or, better, in conjunction with another indica- tor group (i.e. Oertli 2008). Indeed, the regional species pool of Odonata in altitude areas is small (e.g. Figure 4), but nevertheless large enough for conducting long range monitoring. Here, the local species richness has been shown to be a particular sensitive metric, hence being a good candidate to join the set of metrics used for long range monitoring. Local extinctions of cold stenothermal species, or colonisation of lowland species, are early warning events for mountain waterbodies that should be monitored. Acknowledgement The author is grateful to the Swiss Centre for Fauna Cartography (CSCF) for provid- ing data on species distribution and to the Swiss Federal Institute for Forest, Snow and Landscape for Research for the climatic data. Thanks to Michael de la Harpe and to Véronique Rosset for producing the figures 1, 3 and 4. Jane O’Rourke and Franck Cat- taneo provided helpful linguistic corrections. References De Meester L, Declerck S, Stoks R, Louette G, van de Meutter F, de Bie T, Michels E, Bren- donck L (2005) Ponds and pools as model systems in conservation biology, ecology and evolutionary biology. Aquatic Conservation 15: 715-725. Gonseth Y, Monnerat C (2002) Liste rouge des espéces menacées en Suisse. Libellules. Centre Suisse de Cartographie de la Faune. Neuchatel. Hohmann R, Volken E, Miller-Ferch G, Neu U, Ritz Ch, Kull Ch (Eds) (2007) Climate change and Switzerland in 2050. Expected Impacts on Environment, Society and Eco- nomics. OcCC & ProClim, Bern, 168pp. [in German] Intergovernmental Panel on Climate Change (2007) Climate change 2007: the physical science basis Summary for policymakers. Working group I fourth assessment report. Cambridge University Press, Cambridge. Maurer K, Weyand A, Fischer M, Stocklin J (2006) Old cultural traditions, in addition to land use and topography, are shaping plant diversity of grasslands in the Alps. Biological Conservation 130: 438-446. Oertli B (2008) Dragonflies in the Assessment and the Monitoring of Aquatic Habitats. In: Cérdoba-Aguilar A (Ed) Dragonflies: Model Organisms for Ecological and Evolutionary Research: in press. Oxford University Press, Oxford, New York. Local species richness as indicator of climate warming 251 Oertli B, Auderset Joye D, Castella E, Juge R, Cambin D, Lachavanne J-B (2002) Does size matter? The relationship between pond area and biodiversity. Biological Conservation 104: 59-70. Oertli B, Auderset Joye D, Castella E, Juge R, Lehmann A, Lachavanne J-B (2005) PLOCH: a standardised method for sampling and assessing the biodiversity in ponds. Aquatic Con- servation 15: 665-679. Oertli B, Indermuehle N, Angélibert $, Hinden H, Stoll A (2008) Macroinvertebrate assem- blages in 25 high alpine ponds of the Swiss National Park (Cirque of Macun) and relation to environmental variables. Hydrobiologia, 597: 29-41. Ott J (2007) The expansion of Crocothemis erythraea (Brullé, 1832) in Germany — an indicator of climatic changes. In: Tyagi BK (Ed) Biology of Dragonflies: 201-222. Scientific Pub- lishers, Jodhpur, India, 210-222. Pauli H, Gottfried M, Reiter K, Klettner C, Grabherr G (2007) Signals of range expansions and contractions of vascular plants in the high Alps: observations (1994—2004) at the GLORIA* master site Schrankogel, Tyrol, Austria. Global Change Biology 13: 147-156. Rahbek C (1995) The elevational gradient of species richness: a uniform pattern? Ecography 18: 200-205. Ricklefs R (1990) Ecology, 3 edn. W.H. Freeman and Company, New York. Samways M (2009) Dragonflies as Focal Organisms in Contemporary Conservation Biology. In: Cérdoba-Aguilar A (Ed) Dragonflies: Model Organisms for Ecological and Evolution- ary Research. Oxford University Press, Oxford, New York. Wildermuth H, Gonseth Y, Maibach A (2005) Odonata. Les Libellules de Suisse. Centre Suisse de Cartographie de la Faune, Neuchatel.