Niche, distribution and global change : modeling insights into biogeography and conservation biology Olivier Broennimann

RESUME Les évidences montrant que les changements globaux affectent la biodiversité s'accumulent. Les facteurs les plus influant dans ce processus sont les changements et destructions d'habitat, l'expansion des espèces envahissantes et l'impact des changements climatiques. Une évaluation pertinente de la réponse des espèces face à ces changements est essentielle pour proposer des mesures permettant de réduire le déclin actuel de la biodiversité. La modélisation de la répartition d'espèces basée sur la niche (NBM) est l'un des rares outils permettant cette évaluation. Néanmoins, leur application dans le contexte des changements globaux repose sur des hypothèses restrictives et demande une interprétation critique. Ce travail présente une série d'études de cas investiguant les possibilités et limitations de cette approche pour prédire l'impact des changements globaux. Deux études traitant des menaces sur les espèces rares et en danger d'extinction sont présentées. Les caractéristiques éco-géographiques de 118 plantes avec un haut degré de priorité de conservation sont revues. La prévalence des types de rareté sont analysées en relation avec leur risque d'extinction UICN. La revue souligne l'importance de la conservation à l'échelle régionale. Une évaluation de la rareté à échelle globale peut être trompeuse pour certaine espèces car elle ne tient pas en compte des différents degrés de rareté que présente une espèce à différentes échelles spatiales. La deuxième étude test une approche pour améliorer l'échantillonnage d'espèces rares en incluant des phases itératives de modélisation et d'échantillonnage sur le terrain. L'application de l'approche en biologie de la conservation (illustrée ici par le cas du chardon bleu, Eryngium alpinum), permettrait de réduire le temps et les coûts d'échantillonnage. Deux études sur l'impact des changements climatiques sur la faune et la flore africaine sont présentées. La première étude évalue la sensibilité de 227 mammifères africains face aux climatiques d'ici 2050. Elle montre qu'un nombre important d'espèces pourrait être bientôt en danger d'extinction et que les parcs nationaux africains (principalement ceux situé en milieux xériques) pourraient ne pas remplir leur mandat de protection de la biodiversité dans le futur. La seconde étude modélise l'aire de répartition en 2050 de 975 espèces de plantes endémiques du sud de l'Afrique. L'étude propose l'inclusion de méthodes améliorant la prédiction des risques liés aux changements climatiques. Elle propose également une méthode pour estimer a priori la sensibilité d'une espèce aux changements climatiques à partir de ses propriétés écologiques et des caractéristiques de son aire de répartition. Trois études illustrent l'utilisation des modèles dans l'étude des invasions biologiques. Une première étude relate l'expansion de la laitue sáuvage (Lactuca serriola) vers le nord de l'Europe en lien avec les changements du climat depuis 250 ans. La deuxième étude analyse le potentiel d'invasion de la centaurée tachetée (Centaures maculosa), une mauvaise herbe importée en Amérique du nord vers 1890. L'étude apporte la preuve qu'une espèce envahissante peut occuper une niche climatique différente après introduction sur un autre continent. Les modèles basés sur l'aire native prédisent de manière incorrecte l'entier de l'aire envahie mais permettent de prévoir les aires d'introductions potentielles. Une méthode alternative, incluant la calibration du modèle à partir des deux aires où l'espèce est présente, est proposée pour améliorer les prédictions de l'invasion en Amérique du nord. Je présente finalement une revue de la littérature sur la dynamique de la niche écologique dans le temps et l'espace. Elle synthétise les récents développements théoriques concernant le conservatisme de la niche et propose des solutions pour améliorer la pertinence des prédictions d'impact des changements climatiques et des invasions biologiques. SUMMARY Evidences are accumulating that biodiversity is facing the effects of global change. The most influential drivers of change in ecosystems are land-use change, alien species invasions and climate change impacts. Accurate projections of species' responses to these changes are needed to propose mitigation measures to slow down the on-going erosion of biodiversity. Niche-based models (NBM) currently represent one of the only tools for such projections. However, their application in the context of global changes relies on restrictive assumptions, calling for cautious interpretations. In this thesis I aim to assess the effectiveness and shortcomings of niche-based models for the study of global change impacts on biodiversity through the investigation of specific, unsolved limitations and suggestion of new approaches. Two studies investigating threats to rare and endangered plants are presented. I review the ecogeographic characteristic of 118 endangered plants with high conservation priority in Switzerland. The prevalence of rarity types among plant species is analyzed in relation to IUCN extinction risks. The review underlines the importance of regional vs. global conservation and shows that a global assessment of rarity might be misleading for some species because it can fail to account for different degrees of rarity at a variety of spatial scales. The second study tests a modeling framework including iterative steps of modeling and field surveys to improve the sampling of rare species. The approach is illustrated with a rare alpine plant, Eryngium alpinum and shows promise for complementing conservation practices and reducing sampling costs. Two studies illustrate the impacts of climate change on African taxa. The first one assesses the sensitivity of 277 mammals at African scale to climate change by 2050 in terms of species richness and turnover. It shows that a substantial number of species could be critically endangered in the future. National parks situated in xeric ecosystems are not expected to meet their mandate of protecting current species diversity in the future. The second study model the distribution in 2050 of 975 endemic plant species in southern Africa. The study proposes the inclusion of new methodological insights improving the accuracy and ecological realism of predictions of global changes studies. It also investigates the possibility to estimate a priori the sensitivity of a species to climate change from the geographical distribution and ecological proprieties of the species. Three studies illustrate the application of NBM in the study of biological invasions. The first one investigates the Northwards expansion of Lactuca serriola L. in Europe during the last 250 years in relation with climate changes. In the last two decades, the species could not track climate change due to non climatic influences. A second study analyses the potential invasion extent of spotted knapweed, a European weed first introduced into North America in the 1890s. The study provides one of the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. Models fail to predict the current full extent of the invasion, but correctly predict areas of introduction. An alternative approach, involving the calibration of models with pooled data from both ranges, is proposed to improve predictions of the extent of invasion on models based solely on the native range. I finally present a review on the dynamic nature of ecological niches in space and time. It synthesizes the recent theoretical developments to the niche conservatism issues and proposes solutions to improve confidence in NBM predictions of the impacts of climate change and species invasions on species distributions.

Spatial distribution and biomass of aquatic rooted macrophytes and their relevance in the metabolism of a Mediterranean coastal lagoon

This work aims to characterise the current autotrophic compartment of the Albufera des Grau coastal lagoon (Menorca, Balearic Islands) and to assess the relationship between the submerged macrophytes and the limnological parameters of the lagoon. During the study period the submerged vegetation was dominated by the macrophyte Ruppia cirrhosa, which formed dense extensive meadows covering 79% of the surface. Another macrophyte species, Potamogeton pectinatus, was also observed but only forming small stands near the rushing streams. Macroalgae were only occasionally observed. Macrophyte biomass showed a clear seasonal trend, with maximum values in July. The biomass of R. cirrhosa achieved 1760 g DW m-2, the highest biomass ever reported for this species in the literature. The seasonal production-decomposition cycle of the macrophyte meadows appears to drive the nutrient dynamics and carbon fluxes in the lagoon. Despite the significant biomass accumulation and the absence of a washout of nutrients and organic matter to the sea, the lagoon did not experience a dystrophic collapse. These results indicate that internal metabolism is more important than exchange processes in the lagoon.

Spatial pattern and process at the plant neighbourhood scale: insights from communities dominated by the clonal grass Elymus repens (L.) Gould.

Question: How do clonal traits of a locally dominant grass (Elymus repens (L.) Gould.) respond to soil heterogeneity and shape spatial patterns of its tillers? How do tiller spatial patterns constrain seedling recruitment within the community?Locations: Artificial banks of the River Rhone, France.Material and Methods: We examined 45 vegetation patches dominated by Elymus repens. During a first phase we tested relationships between soil variables and three clonal traits (spacer length, number of clumping tillers and branching rate), and between the same clonal traits and spatial patterns (i.e. density and degree of spatial aggregation) of tillers at a very fine scale. During a second phase, we performed a sowing experiment to investigate effects of density and spatial patterns of E. repens on recruitment of eight species selected from the regional species pool.Results: Clonal traits had clear effects - especially spacer length - on densification and aggregation of E. repens tillers and, at the same time, a clear response of these same clonal traits as soil granulometry changed. The density and degree of aggregation of E. repens tillers was positively correlated to total seedling cover and diversity at the finest spatial scales.Conclusions: Spatial patterning of a dominant perennial grass responds to soil heterogeneity through modifications of its clonal morphology as a trade-off between phalanx and guerrilla forms. In turn, spatial patterns have strong effects on abundance and diversity of seedlings. Spatial patterns of tillers most probably led to formation of endogenous gaps in which the recruitment of new plant individuals was enhanced. Interestingly, we also observed more idiosyncratic effects of tiller spatial patterns on seedling cover and diversity when focusing on different growth forms of the sown species.

Spatial Exploration of Age Distribution in Catalan Municipalities

This paper takes the shelf and digs into the complex population’s age structure of Catalan municipalities for the year 2009. Catalonia is a very heterogeneous territory, and age pyramids vary considerably across different areas of the territory, existing geographical factors shaping municipalities’ age distributions. By means of spatial statistics methodologies, this piece of research tries to assess which spatial factors determine the location, scale and shape of local distributions. The results show that there exist different distributional patterns across the geography according to specific local determinants. Keywords: Spatial Models. JEL Classification: C21.

Predicting current and future spatial community patterns of plant functional traits

Community-level patterns of functional traits relate to community assembly and ecosystem functioning. By modelling the changes of different indices describing such patterns - trait means, extremes and diversity in communities - as a function of abiotic gradients, we could understand their drivers and build projections of the impact of global change on the functional components of biodiversity. We used five plant functional traits (vegetative height, specific leaf area, leaf dry matter content, leaf nitrogen content and seed mass) and non-woody vegetation plots to model several indices depicting community-level patterns of functional traits from a set of abiotic environmental variables (topographic, climatic and edaphic) over contrasting environmental conditions in a mountainous landscape. We performed a variation partitioning analysis to assess the relative importance of these variables for predicting patterns of functional traits in communities, and projected the best models under several climate change scenarios to examine future potential changes in vegetation functional properties. Not all indices of trait patterns within communities could be modelled with the same level of accuracy: the models for mean and extreme values of functional traits provided substantially better predictive accuracy than the models calibrated for diversity indices. Topographic and climatic factors were more important predictors of functional trait patterns within communities than edaphic predictors. Overall, model projections forecast an increase in mean vegetation height and in mean specific leaf area following climate warming. This trend was important at mid elevation particularly between 1000 and 2000 m asl. With this study we showed that topographic, climatic and edaphic variables can successfully model descriptors of community-level patterns of plant functional traits such as mean and extreme trait values. However, which factors determine the diversity of functional traits in plant communities remains unclear and requires more investigations.

Effects of alternative sets of climatic predictors on species distribution models and associated estimates of extinction risk: a test with plants in an arid environment

Aim To evaluate the effects of using distinct alternative sets of climatic predictor variables on the performance, spatial predictions and future projections of species distribution models (SDMs) for rare plants in an arid environment. . Location Atacama and Peruvian Deserts, South America (18º30'S - 31º30'S, 0 - 3 000 m) Methods We modelled the present and future potential distributions of 13 species of Heliotropium sect. Cochranea, a plant group with a centre of diversity in the Atacama Desert. We developed and applied a sequential procedure, starting from climate monthly variables, to derive six alternative sets of climatic predictor variables. We used them to fit models with eight modelling techniques within an ensemble forecasting framework, and derived climate change projections for each of them. We evaluated the effects of using these alternative sets of predictor variables on performance, spatial predictions and projections of SDMs using Generalised Linear Mixed Models (GLMM). Results The use of distinct sets of climatic predictor variables did not have a significant effect on overall metrics of model performance, but had significant effects on present and future spatial predictions. Main conclusion Using different sets of climatic predictors can yield the same model fits but different spatial predictions of current and future species distributions. This represents a new form of uncertainty in model-based estimates of extinction risk that may need to be better acknowledged and quantified in future SDM studies.

Spatial statistical analysis and selection of genotypes in plant breeding

The objective of this study was to evaluate the efficiency of spatial statistical analysis in the selection of genotypes in a plant breeding program and, particularly, to demonstrate the benefits of the approach when experimental observations are not spatially independent. The basic material of this study was a yield trial of soybean lines, with five check varieties (of fixed effect) and 110 test lines (of random effects), in an augmented block design. The spatial analysis used a random field linear model (RFML), with a covariance function estimated from the residuals of the analysis considering independent errors. Results showed a residual autocorrelation of significant magnitude and extension (range), which allowed a better discrimination among genotypes (increase of the power of statistical tests, reduction in the standard errors of estimates and predictors, and a greater amplitude of predictor values) when the spatial analysis was applied. Furthermore, the spatial analysis led to a different ranking of the genetic materials, in comparison with the non-spatial analysis, and a selection less influenced by local variation effects was obtained.

Are niche-based species distribution models transferable in space?

Aim To assess the geographical transferability of niche-based species distribution models fitted with two modelling techniques. Location Two distinct geographical study areas in Switzerland and Austria, in the subalpine and alpine belts. Methods Generalized linear and generalized additive models (GLM and GAM) with a binomial probability distribution and a logit link were fitted for 54 plant species, based on topoclimatic predictor variables. These models were then evaluated quantitatively and used for spatially explicit predictions within (internal evaluation and prediction) and between (external evaluation and prediction) the two regions. Comparisons of evaluations and spatial predictions between regions and models were conducted in order to test if species and methods meet the criteria of full transferability. By full transferability, we mean that: (1) the internal evaluation of models fitted in region A and B must be similar; (2) a model fitted in region A must at least retain a comparable external evaluation when projected into region B, and vice-versa; and (3) internal and external spatial predictions have to match within both regions. Results The measures of model fit are, on average, 24% higher for GAMs than for GLMs in both regions. However, the differences between internal and external evaluations (AUC coefficient) are also higher for GAMs than for GLMs (a difference of 30% for models fitted in Switzerland and 54% for models fitted in Austria). Transferability, as measured with the AUC evaluation, fails for 68% of the species in Switzerland and 55% in Austria for GLMs (respectively for 67% and 53% of the species for GAMs). For both GAMs and GLMs, the agreement between internal and external predictions is rather weak on average (Kulczynski's coefficient in the range 0.3-0.4), but varies widely among individual species. The dominant pattern is an asymmetrical transferability between the two study regions (a mean decrease of 20% for the AUC coefficient when the models are transferred from Switzerland and 13% when they are transferred from Austria). Main conclusions The large inter-specific variability observed among the 54 study species underlines the need to consider more than a few species to test properly the transferability of species distribution models. The pronounced asymmetry in transferability between the two study regions may be due to peculiarities of these regions, such as differences in the ranges of environmental predictors or the varied impact of land-use history, or to species-specific reasons like differential phenotypic plasticity, existence of ecotypes or varied dependence on biotic interactions that are not properly incorporated into niche-based models. The lower variation between internal and external evaluation of GLMs compared to GAMs further suggests that overfitting may reduce transferability. Overall, a limited geographical transferability calls for caution when projecting niche-based models for assessing the fate of species in future environments.

Spatial Polarization of the Ecological Footprint distribution

The international allocation of natural resources is determined, not by any ethical or ecological criteria, but by the dominance of market mechanisms. From a core-periphery perspective, this allocation may even be driven by historically determined structural patterns, with a core group of countries whose consumption appropriates most available natural resources, and another group, having low natural resource consumption, which plays a peripheral role. This article consists of an empirical distributional analysis of natural resource consumption (as measured by Ecological Footprints) whose purpose is to assess the extent to which the distribution of consumption responds to polarization (as opposed to mere inequality). To assess this, we estimate and decompose different polarization indices for a balanced sample of 119 countries over the period 1961 to 2007. Our results points toward a polarized distribution which is consistent with a core-periphery framework.

Do stacked species distribution models reflect altitudinal diversity patterns?

The objective of this study was to evaluate the performance of stacked species distribution models in predicting the alpha and gamma species diversity patterns of two important plant clades along elevation in the Andes. We modelled the distribution of the species in the Anthurium genus (53 species) and the Bromeliaceae family (89 species) using six modelling techniques. We combined all of the predictions for the same species in ensemble models based on two different criteria: the average of the rescaled predictions by all techniques and the average of the best techniques. The rescaled predictions were then reclassified into binary predictions (presence/absence). By stacking either the original predictions or binary predictions for both ensemble procedures, we obtained four different species richness models per taxa. The gamma and alpha diversity per elevation band (500 m) was also computed. To evaluate the prediction abilities for the four predictions of species richness and gamma diversity, the models were compared with the real data along an elevation gradient that was independently compiled by specialists. Finally, we also tested whether our richness models performed better than a null model of altitudinal changes of diversity based on the literature. Stacking of the ensemble prediction of the individual species models generated richness models that proved to be well correlated with the observed alpha diversity richness patterns along elevation and with the gamma diversity derived from the literature. Overall, these models tend to overpredict species richness. The use of the ensemble predictions from the species models built with different techniques seems very promising for modelling of species assemblages. Stacking of the binary models reduced the over-prediction, although more research is needed. The randomisation test proved to be a promising method for testing the performance of the stacked models, but other implementations may still be developed.

Revisiting spatial distribution and biochemical composition of calcium-containing crystals in human osteoarthritic articular cartilage.

INTRODUCTION: Calcium-containing (CaC) crystals, including basic calcium phosphate (BCP) and calcium pyrophosphate dihydrate (CPP), are associated with destructive forms of osteoarthritis (OA). We assessed their distribution and biochemical and morphologic features in human knee OA cartilage. METHODS: We prospectively included 20 patients who underwent total knee replacement (TKR) for primary OA. CaC crystal characterization and identification involved Fourier-transform infra-red spectrometry and scanning electron microscopy of 8 to 10 cartilage zones of each knee, including medial and lateral femoral condyles and tibial plateaux and the intercondyle zone. Differential expression of genes involved in the mineralization process between cartilage with and without calcification was assessed in samples from 8 different patients by RT-PCR. Immunohistochemistry and histology studies were performed in 6 different patients. RESULTS: Mean (SEM) age and body mass index of patients at the time of TKR was 74.6 (1.7) years and 28.1 (1.6) kg/m², respectively. Preoperative X-rays showed joint calcifications (chondrocalcinosis) in 4 cases only. The medial femoro-tibial compartment was the most severely affected in all cases, and mean (SEM) Kellgren-Lawrence score was 3.8 (0.1). All 20 OA cartilages showed CaC crystals. The mineral content represented 7.7% (8.1%) of the cartilage weight. All patients showed BCP crystals, which were associated with CPP crystals for 8 joints. CaC crystals were present in all knee joint compartments and in a mean of 4.6 (1.7) of the 8 studied areas. Crystal content was similar between superficial and deep layers and between medial and femoral compartments. BCP samples showed spherical structures, typical of biological apatite, and CPP samples showed rod-shaped or cubic structures. The expression of several genes involved in mineralization, including human homolog of progressive ankylosis, plasma-cell-membrane glycoprotein 1 and tissue-nonspecific alkaline phosphatase, was upregulated in OA chondrocytes isolated from CaC crystal-containing cartilages. CONCLUSIONS: CaC crystal deposition is a widespread phenomenon in human OA articular cartilage involving the entire knee cartilage including macroscopically normal and less weight-bearing zones. Cartilage calcification is associated with altered expression of genes involved in the mineralisation process.

Spatial polarization of the ecological footprint distribution

The international allocation of natural resources is determined, not by any ethical or ecological criteria, but by the dominance of market mechanisms. From a core-periphery perspective, this allocation may even be driven by historically determined structural patterns, with a core group of countries whose consumption appropriates most available natural resources, and another group, having low natural resource consumption, which plays a peripheral role. This article consists of an empirical distributional analysis of natural resource consumption (as measured by Ecological Footprints) whose purpose is to assess the extent to which the distribution of consumption responds to polarization (as opposed to mere inequality). To assess this, we estimate and decompose different polarization indices for a balanced sample of 119 countries over the period 1961 to 2007. Our results points toward a polarized distribution which is consistent with a core-periphery framework. Keywords: Polarization, Core-Periphery, Ecological Footprint

Comparative population genetic structure in a plant-pollinator/seed predator system.

Comparative analyses of spatial genetic structure of populations of plants and the insects they interact with provide an indication of how gene flow, natural selection and genetic drift may jointly influence the distribution of genetic variation and potential for local co-adaptation for interacting species. Here, we analysed the spatial scale of genetic structure within and among nine populations of an interacting species pair, the white campion Silene latifolia and the moth Hadena bicruris, along a latitudinal gradient across Northern/Central Europe. This dioecious, short-lived perennial plant inhabits patchy, often disturbed environments. The moth H. bicruris acts both as its pollinator and specialist seed predator that reproduces by laying eggs in S. latifolia flowers. We used nine microsatellite markers for S. latifolia and eight newly developed markers for H. bicruris. We found high levels of inbreeding in most populations of both plant and pollinator/seed predator. Among populations, significant genetic structure was observed for S. latifolia but not for its pollinator/seed predator, suggesting that despite migration among populations of H. bicruris, pollen is not, or only rarely, carried over between populations, thus maintaining genetic structure among plant populations. There was a weak positive correlation between genetic distances of S. latifolia and H. bicruris. These results indicate that while significant structure of S. latifolia populations creates the potential for differentiation at traits relevant for the interaction with the pollinator/seed predator, substantial gene flow in H. bicruris may counteract this process in at least some populations.

Evaluation of spatial distribution of flatfish by laser scannig

Peer reviewed

Spatial predictions of land-use transitions and associated threats to biodiversity: the case of forest regrowth in mountain grasslands

Question Can we predict where forest regrowth caused by abandonment of agricultural activities is likely to occur? Can we assess how it may conflict with grassland diversity hotspots? Location Western Swiss Alps (4003210m a.s.l.). Methods We used statistical models to predict the location of land abandonment by farmers that is followed by forest regrowth in semi-natural grasslands of the Western Swiss Alps. Six modelling methods (GAM, GBM, GLM, RF, MDA, MARS) allowing binomial distribution were tested on two successive transitions occurring between three time periods. Models were calibrated using data on land-use change occurring between 1979 and 1992 as response, and environmental, accessibility and socio-economic variables as predictors, and these were validated for their capacity to predict the changes observed from 1992 to 2004. Projected probabilities of land-use change from an ensemble forecast of the six models were combined with a model of plant species richness based on a field inventory, allowing identification of critical grassland areas for the preservation of biodiversity. Results Models calibrated over the first land-use transition period predicted the second transition with reasonable accuracy. Forest regrowth occurs where cultivation costs are high and yield potential is low, i.e. on steeper slopes and at higher elevations. Overlaying species richness with land-use change predictions, we identified priority areas for the management and conservation of biodiversity at intermediate elevations. Conclusions Combining land-use change and biodiversity projections, we propose applied management measures for targeted/identified locations to limit the loss of biodiversity that could otherwise occur through loss of open habitats. The same approach could be applied to other types of land-use changes occurring in other ecosystems.