The ghost of past species occurrence: improving species distribution models for presence-only data

1. Model-based approaches have been used increasingly in conservation biology over recent years. Species presence data used for predictive species distribution modelling are abundant in natural history collections, whereas reliable absence data are sparse, most notably for vagrant species such as butterflies and snakes. As predictive methods such as generalized linear models (GLM) require absence data, various strategies have been proposed to select pseudo-absence data. However, only a few studies exist that compare different approaches to generating these pseudo-absence data. 2. Natural history collection data are usually available for long periods of time (decades or even centuries), thus allowing historical considerations. However, this historical dimension has rarely been assessed in studies of species distribution, although there is great potential for understanding current patterns, i.e. the past is the key to the present. 3. We used GLM to model the distributions of three 'target' butterfly species, Melitaea didyma, Coenonympha tullia and Maculinea teleius, in Switzerland. We developed and compared four strategies for defining pools of pseudo-absence data and applied them to natural history collection data from the last 10, 30 and 100 years. Pools included: (i) sites without target species records; (ii) sites where butterfly species other than the target species were present; (iii) sites without butterfly species but with habitat characteristics similar to those required by the target species; and (iv) a combination of the second and third strategies. Models were evaluated and compared by the total deviance explained, the maximized Kappa and the area under the curve (AUC). 4. Among the four strategies, model performance was best for strategy 3. Contrary to expectations, strategy 2 resulted in even lower model performance compared with models with pseudo-absence data simulated totally at random (strategy 1). 5. Independent of the strategy model, performance was enhanced when sites with historical species presence data were not considered as pseudo-absence data. Therefore, the combination of strategy 3 with species records from the last 100 years achieved the highest model performance. 6. Synthesis and applications. The protection of suitable habitat for species survival or reintroduction in rapidly changing landscapes is a high priority among conservationists. Model-based approaches offer planning authorities the possibility of delimiting priority areas for species detection or habitat protection. The performance of these models can be enhanced by fitting them with pseudo-absence data relying on large archives of natural history collection species presence data rather than using randomly sampled pseudo-absence data.

Pattern recognition ecological niche models fit to presence-only and presence-absence data

1. Identifying the boundary of a species' niche from observational and environmental data is a common problem in ecology and conservation biology and a variety of techniques have been developed or applied to model niches and predict distributions. Here, we examine the performance of some pattern-recognition methods as ecological niche models (ENMs). Particularly, one-class pattern recognition is a flexible and seldom used methodology for modelling ecological niches and distributions from presence-only data. The development of one-class methods that perform comparably to two-class methods (for presence/absence data) would remove modelling decisions about sampling pseudo-absences or background data points when absence points are unavailable. 2. We studied nine methods for one-class classification and seven methods for two-class classification (five common to both), all primarily used in pattern recognition and therefore not common in species distribution and ecological niche modelling, across a set of 106 mountain plant species for which presence-absence data was available. We assessed accuracy using standard metrics and compared trade-offs in omission and commission errors between classification groups as well as effects of prevalence and spatial autocorrelation on accuracy. 3. One-class models fit to presence-only data were comparable to two-class models fit to presence-absence data when performance was evaluated with a measure weighting omission and commission errors equally. One-class models were superior for reducing omission errors (i.e. yielding higher sensitivity), and two-classes models were superior for reducing commission errors (i.e. yielding higher specificity). For these methods, spatial autocorrelation was only influential when prevalence was low. 4. These results differ from previous efforts to evaluate alternative modelling approaches to build ENM and are particularly noteworthy because data are from exhaustively sampled populations minimizing false absence records. Accurate, transferable models of species' ecological niches and distributions are needed to advance ecological research and are crucial for effective environmental planning and conservation; the pattern-recognition approaches studied here show good potential for future modelling studies. This study also provides an introduction to promising methods for ecological modelling inherited from the pattern-recognition discipline.

Novel methods improve prediction of species' distributions from occurrence data

Prediction of species' distributions is central to diverse applications in ecology, evolution and conservation science. There is increasing electronic access to vast sets of occurrence records in museums and herbaria, yet little effective guidance on how best to use this information in the context of numerous approaches for modelling distributions. To meet this need, we compared 16 modelling methods over 226 species from 6 regions of the world, creating the most comprehensive set of model comparisons to date. We used presence-only data to fit models, and independent presence-absence data to evaluate the predictions. Along with well-established modelling methods such as generalised additive models and GARP and BIOCLIM, we explored methods that either have been developed recently or have rarely been applied to modelling species' distributions. These include machine-learning methods and community models, both of which have features that may make them particularly well suited to noisy or sparse information, as is typical of species' occurrence data. Presence-only data were effective for modelling species' distributions for many species and regions. The novel methods consistently outperformed more established methods. The results of our analysis are promising for the use of data from museums and herbaria, especially as methods suited to the noise inherent in such data improve.

What matters for predicting spatial distributions of trees: techniques, data, or species' characteristics?

Data characteristics and species traits are expected to influence the accuracy with which species' distributions can be modeled and predicted. We compare 10 modeling techniques in terms of predictive power and sensitivity to location error, change in map resolution, and sample size, and assess whether some species traits can explain variation in model performance. We focused on 30 native tree species in Switzerland and used presence-only data to model current distribution, which we evaluated against independent presence-absence data. While there are important differences between the predictive performance of modeling methods, the variance in model performance is greater among species than among techniques. Within the range of data perturbations in this study, some extrinsic parameters of data affect model performance more than others: location error and sample size reduced performance of many techniques, whereas grain had little effect on most techniques. No technique can rescue species that are difficult to predict. The predictive power of species-distribution models can partly be predicted from a series of species characteristics and traits based on growth rate, elevational distribution range, and maximum elevation. Slow-growing species or species with narrow and specialized niches tend to be better modeled. The Swiss presence-only tree data produce models that are reliable enough to be useful in planning and management applications.

Sensitivity of predictive species distribution models to change in grain size

Predictive species distribution modelling (SDM) has become an essential tool in biodiversity conservation and management. The choice of grain size (resolution) of environmental layers used in modelling is one important factor that may affect predictions. We applied 10 distinct modelling techniques to presence-only data for 50 species in five different regions, to test whether: (1) a 10-fold coarsening of resolution affects predictive performance of SDMs, and (2) any observed effects are dependent on the type of region, modelling technique, or species considered. Results show that a 10 times change in grain size does not severely affect predictions from species distribution models. The overall trend is towards degradation of model performance, but improvement can also be observed. Changing grain size does not equally affect models across regions, techniques, and species types. The strongest effect is on regions and species types, with tree species in the data sets (regions) with highest locational accuracy being most affected. Changing grain size had little influence on the ranking of techniques: boosted regression trees remain best at both resolutions. The number of occurrences used for model training had an important effect, with larger sample sizes resulting in better models, which tended to be more sensitive to grain. Effect of grain change was only noticeable for models reaching sufficient performance and/or with initial data that have an intrinsic error smaller than the coarser grain size.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

Modelling ecological niches with support vector machines

1. The ecological niche is a fundamental biological concept. Modelling species' niches is central to numerous ecological applications, including predicting species invasions, identifying reservoirs for disease, nature reserve design and forecasting the effects of anthropogenic and natural climate change on species' ranges. 2. A computational analogue of Hutchinson's ecological niche concept (the multidimensional hyperspace of species' environmental requirements) is the support of the distribution of environments in which the species persist. Recently developed machine-learning algorithms can estimate the support of such high-dimensional distributions. We show how support vector machines can be used to map ecological niches using only observations of species presence to train distribution models for 106 species of woody plants and trees in a montane environment using up to nine environmental covariates. 3. We compared the accuracy of three methods that differ in their approaches to reducing model complexity. We tested models with independent observations of both species presence and species absence. We found that the simplest procedure, which uses all available variables and no pre-processing to reduce correlation, was best overall. Ecological niche models based on support vector machines are theoretically superior to models that rely on simulating pseudo-absence data and are comparable in empirical tests. 4. Synthesis and applications. Accurate species distribution models are crucial for effective environmental planning, management and conservation, and for unravelling the role of the environment in human health and welfare. Models based on distribution estimation rather than classification overcome theoretical and practical obstacles that pervade species distribution modelling. In particular, ecological niche models based on machine-learning algorithms for estimating the support of a statistical distribution provide a promising new approach to identifying species' potential distributions and to project changes in these distributions as a result of climate change, land use and landscape alteration.

Conservation of a threatened European tree frog (Hyla arborea) metapopulation

Summary: Amphibians are among the most vulnerable animals of the world. One third of all species are currently threatened with extinction. Habitat loss is the major menace to pond- and stream-breeding species in the old world. In highly urbanized landscape like the Swiss Plateau, most species suffer from habitat reduction and fragmentation. Among all indigenous species, the European tree frog (Hyla arborea L., 1758) is one of the most endangered. It experienced an alarming decline during the last century and its regional long-term persistence is not guaranteed. We developed a monitoring framework based on calling male counts which included multiple visits to each wetland during the reproduction period in order to precisely determine its distribution on the Lemanic coast. Our results indicate that visiting populations 3 limes under suitable climatic conditions (temperature >20°C) provides reliable presence/absence data. Based on our monitoring data, we analyzed the species requirements regarding its breeding habitat. It appeared that anthropogenic activities had paradoxical effects on the species. On one hand, urbanization, traffic and intensive agriculture had a strong detrimental effect on tree frog distribution. On the other hand, large tree frog populations were frequently associated with gravel pits and military training grounds. Our results allowed us to create a habitat suitability map taking into account detrimental landscape elements around ponds (>1100m away from urban areas and >500m away from first class roads). In parallel, we developed a metapopulation model of the European tree frog in order to identify the critical threats to the long term persistence of the species. Our results indicated that suitable pond density is at the low end of the species requirements. Pond creation must therefore be considered an essential complementary approach to pond conservation and restoration. Our model also provided a mapping solution permitting the location of the must suitable area for pond creation from a metapopulation perspective. As many other amphibians, the European tree frog is not only exposed to an aquatic habitat (breeding and larval period), but also to a terrestrial stage (summer and overwintering habitats). Unfortunately, animals in their terrestrial phase are less conspicuous and, as a consequence, their terrestrial needs are relatively unknown. Using a recent tracking method (the Harmonic Direction Finder), we followed post-breeding frogs and identified favored terrestrial habitats, thus providing another practical conservation tool. We conclude that only the combination of multiple spatially explicit approaches (landscape-scale habitat suitability, metapopulation dynamics and terrestrial needs) is likely to provide wildlife managers with effective tools for the conservation of highly endangered amphibians. Résumé: Les amphibiens font partie des animaux les plus vulnérables du monde. Un tiers des espèces est actuellement menacé d'extinction. Dans l'ancien monde, la disparition des habitats constitue la principale menace pour les grenouilles, crapauds, tritons et salamandres. Dans les paysages fortement urbanisés comme le Plateau Suisse, la plupart des espèces souffrent d'une réduction et d'une fragmentation de leurs habitats. Parmi toutes les espèces indigènes, la rainette verte (Hyla arborea L., 1758) est l'une des plus menacée. Sa distribution a régressé de manière alarmante durant le siècle passé et sa survie régionale à long terme n'est pas assurée. Nous avons développé une méthode de suivi des populations se basant sur le comptage des mâles chanteurs durant la période de reproduction. Cette méthode requiert plusieurs visites à chaque plan d'eau de manière à déterminer précisément la distribution de l'espèce. Nos résultats démontrent que 3 visites par population dans des conditions climatiques favorable (température >20°C) permettent d'obtenir des données de présence/ absence valables. Sur la base de nos comptages sur la Côte lémanique, nous avons analysé les exigences de l'espèce concernant ses sites de reproduction. Il est apparu que les activités humaines avaient un effet paradoxal sur l'espèce. D'une part, l'urbanisation, le trafic routier et l'intensification de l'agriculture ont un effet fortement préjudiciable, tandis que d'autre part les plus grandes populations sont souvent associées à des gravières et autres places d'armes. Nos résultats ont permis de créer une carte de qualité d'habitat prenant en compte les éléments paysagers préjudiciables à la rainette (situé à plus de 1100m de zones urbaines et à plus de 500m de routes de première classe). En parallèle, nous avons développé un modèle métapopulationnel (incluant l'ensemble des populations) de manière à identifier les menaces prépondérantes sur la survie à long terme de l'espèce. Nos résultats ont permis de déterminer que la densité actuelle de plans d'eau adéquats est à la limite inférieure des exigences de l'espèce. La création d'étangs doit donc être considérée comme une approche indispensable et complémentaire à la protection et à la restauration des sites existants. Notre modèle a également fourni des résultats cartographiables permettant l'identification des sites les plus appropriés dans une perspective métapopulationnelle. Comme de nombreux autres amphibiens, la rainette verte est exposée à un habitat aquatique (reproduction et développement larvaire) ainsi qu'à un habitat terrestre (été et hiver). Les animaux étant particulièrement cryptiques dans cette seconde phase, leurs besoins terrestres sont relativement mal connus. Nous avons donc développé une nouvelle méthode de télémétrie basée sur le goniomètre harmonique. Cette méthode nous a permis de suivre des rainettes dans leurs migrations jusqu'à leurs habitats d'été et d'établir ainsi des recommandations pratiques pour la conservation de la rainette. Nous concluons que la combinaison de multiples approches spatialement explicites (qualité d'habitat, dynamique de métapopulation et habitats terrestres) est seule à même de produire des outils efficaces pour la conservation des espèces menacées d'amphibiens.

Changing landscapes of Switzerland: assessing effects of historic and future land uses and land covers on the occurrence of selected butterfly species

Summary Landscapes are continuously changing. Natural forces of change such as heavy rainfall and fires can exert lasting influences on their physical form. However, changes related to human activities have often shaped landscapes more distinctly. In Western Europe, especially modern agricultural practices and the expanse of overbuilt land have left their marks in the landscapes since the middle of the 20th century. In the recent years men realised that mare and more changes that were formerly attributed to natural forces might indirectly be the result of their own action. Perhaps the most striking landscape change indirectly driven by human activity we can witness in these days is the large withdrawal of Alpine glaciers. Together with the landscapes also habitats of animal and plant species have undergone vast and sometimes rapid changes that have been hold responsible for the ongoing loss of biodiversity. Thereby, still little knowledge is available about probable effects of the rate of landscape change on species persistence and disappearance. Therefore, the development and speed of land use/land cover in the Swiss communes between the 1950s and 1990s were reconstructed using 10 parameters from agriculture and housing censuses, and were further correlated with changes in butterfly species occurrences. Cluster analyses were used to detect spatial patterns of change on broad spatial scales. Thereby, clusters of communes showing similar changes or transformation rates were identified for single decades and put into a temporally dynamic sequence. The obtained picture on the changes showed a prevalent replacement of non-intensive agriculture by intensive practices, a strong spreading of urban communes around city centres, and transitions towards larger farm sizes in the mountainous areas. Increasing transformation rates toward more intensive agricultural managements were especially found until the 1970s, whereas afterwards the trends were commonly negative. However, transformation rates representing the development of residential buildings showed positive courses at any time. The analyses concerning the butterfly species showed that grassland species reacted sensitively to the density of livestock in the communes. This might indicate the augmented use of dry grasslands as cattle pastures that show altered plant species compositions. Furthermore, these species also decreased in communes where farms with an agricultural area >5ha have disappeared. The species of the wetland habitats were favoured in communes with smaller fractions of agricultural areas and lower densities of large farms (>10ha) but did not show any correlation to transformation rates. It was concluded from these analyses that transformation rates might influence species disappearance to a certain extent but that states of the environmental predictors might generally outweigh the importance of the corresponding rates. Information on the current distribution of species is evident for nature conservation. Planning authorities that define priority areas for species protection or examine and authorise construction projects need to know about the spatial distribution of species. Hence, models that simulate the potential spatial distribution of species have become important decision tools. The underlying statistical analyses such as the widely used generalised linear models (GLM) often rely on binary species presence-absence data. However, often only species presence data have been colleted, especially for vagrant, rare or cryptic species such as butterflies or reptiles. Modellers have thus introduced randomly selected absence data to design distribution models. Yet, selecting false absence data might bias the model results. Therefore, we investigated several strategies to select more reliable absence data to model the distribution of butterfly species based on historical distribution data. The results showed that better models were obtained when historical data from longer time periods were considered. Furthermore, model performance was additionally increased when long-term data of species that show similar habitat requirements as the modelled species were used. This successful methodological approach was further applied to assess consequences of future landscape changes on the occurrence of butterfly species inhabiting dry grasslands or wetlands. These habitat types have been subjected to strong deterioration in the recent decades, what makes their protection a future mission. Four spatially explicit scenarios that described (i) ongoing land use changes as observed between 1985 and 1997, (ii) liberalised agricultural markets, and (iii) slightly and (iv) strongly lowered agricultural production provided probable directions of landscape change. Current species-environment relationships were derived from a statistical model and used to predict future occurrence probabilities in six major biogeographical regions in Switzerland, comprising the Jura Mountains, the Plateau, the Northern and Southern Alps, as well as the Western and Eastern Central Alps. The main results were that dry grasslands species profited from lowered agricultural production, whereas overgrowth of open areas in the liberalisation scenario might impair species occurrence. The wetland species mostly responded with decreases in their occurrence probabilities in the scenarios, due to a loss of their preferred habitat. Further analyses about factors currently influencing species occurrences confirmed anthropogenic causes such as urbanisation, abandonment of open land, and agricultural intensification. Hence, landscape planning should pay more attention to these forces in areas currently inhabited by these butterfly species to enable sustainable species persistence. In this thesis historical data were intensively used to reconstruct past developments and to make them useful for current investigations. Yet, the availability of historical data and the analyses on broader spatial scales has often limited the explanatory power of the conducted analyses. Meaningful descriptors of former habitat characteristics and abundant species distribution data are generally sparse, especially for fine scale analyses. However, this situation can be ameliorated by broadening the extent of the study site and the used grain size, as was done in this thesis by considering the whole of Switzerland with its communes. Nevertheless, current monitoring projects and data recording techniques are promising data sources that might allow more detailed analyses about effects of long-term species reactions on landscape changes in the near future. This work, however, also showed the value of historical species distribution data as for example their potential to locate still unknown species occurrences. The results might therefore contribute to further research activities that investigate current and future species distributions considering the immense richness of historical distribution data. Résumé Les paysages changent continuellement. Des farces naturelles comme des pluies violentes ou des feux peuvent avoir une influence durable sur la forme du paysage. Cependant, les changements attribués aux activités humaines ont souvent modelé les paysages plus profondément. Depuis les années 1950 surtout, les pratiques agricoles modernes ou l'expansion des surfaces d'habitat et d'infrastructure ont caractérisé le développement du paysage en Europe de l'Ouest. Ces dernières années, l'homme a commencé à réaliser que beaucoup de changements «naturels » pourraient indirectement résulter de ses propres activités. Le changement de paysage le plus apparent dont nous sommes témoins de nos jours est probablement l'immense retraite des glaciers alpins. Avec les paysages, les habitats des animaux et des plantes ont aussi été exposés à des changements vastes et quelquefois rapides, tenus pour coresponsable de la continuelle diminution de la biodiversité. Cependant, nous savons peu des effets probables de la rapidité des changements du paysage sur la persistance et la disparition des espèces. Le développement et la rapidité du changement de l'utilisation et de la couverture du sol dans les communes suisses entre les années 50 et 90 ont donc été reconstruits au moyen de 10 variables issues des recensements agricoles et résidentiels et ont été corrélés avec des changements de présence des papillons diurnes. Des analyses de groupes (Cluster analyses) ont été utilisées pour détecter des arrangements spatiaux de changements à l'échelle de la Suisse. Des communes avec des changements ou rapidités comparables ont été délimitées pour des décennies séparées et ont été placées en séquence temporelle, en rendrent une certaine dynamique du changement. Les résultats ont montré un remplacement répandu d'une agriculture extensive des pratiques intensives, une forte expansion des faubourgs urbains autour des grandes cités et des transitions vers de plus grandes surfaces d'exploitation dans les Alpes. Dans le cas des exploitations agricoles, des taux de changement croissants ont été observés jusqu'aux années 70, alors que la tendance a généralement été inversée dans les années suivantes. Par contre, la vitesse de construction des nouvelles maisons a montré des courbes positives pendant les 50 années. Les analyses sur la réaction des papillons diurnes ont montré que les espèces des prairies sèches supportaient une grande densité de bétail. Il est possible que dans ces communes beaucoup des prairies sèches aient été fertilisées et utilisées comme pâturages, qui ont une autre composition floristique. De plus, les espèces ont diminué dans les communes caractérisées par une rapide perte des fermes avec une surface cultivable supérieure à 5 ha. Les espèces des marais ont été favorisées dans des communes avec peu de surface cultivable et peu de grandes fermes, mais n'ont pas réagi aux taux de changement. Il en a donc été conclu que la rapidité des changements pourrait expliquer les disparitions d'espèces dans certains cas, mais que les variables prédictives qui expriment des états pourraient être des descripteurs plus importants. Des informations sur la distribution récente des espèces sont importantes par rapport aux mesures pour la conservation de la nature. Pour des autorités occupées à définir des zones de protection prioritaires ou à autoriser des projets de construction, ces informations sont indispensables. Les modèles de distribution spatiale d'espèces sont donc devenus des moyens de décision importants. Les méthodes statistiques courantes comme les modèles linéaires généralisés (GLM) demandent des données de présence et d'absence des espèces. Cependant, souvent seules les données de présence sont disponibles, surtout pour les animaux migrants, rares ou cryptiques comme des papillons ou des reptiles. C'est pourquoi certains modélisateurs ont choisi des absences au hasard, avec le risque d'influencer le résultat en choisissant des fausses absences. Nous avons établi plusieurs stratégies, basées sur des données de distribution historique des papillons diurnes, pour sélectionner des absences plus fiables. Les résultats ont démontré que de meilleurs modèles pouvaient être obtenus lorsque les données proviennent des périodes de temps plus longues. En plus, la performance des modèles a pu être augmentée en considérant des données de distribution à long terme d'espèces qui occupent des habitats similaires à ceux de l'espèce cible. Vu le succès de cette stratégie, elle a été utilisée pour évaluer les effets potentiels des changements de paysage futurs sur la distribution des papillons des prairies sèches et marais, deux habitats qui ont souffert de graves détériorations. Quatre scénarios spatialement explicites, décrivant (i) l'extrapolation des changements de l'utilisation de sol tels qu'observés entre 1985 et 1997, (ii) la libéralisation des marchés agricoles, et une production agricole (iii) légèrement amoindrie et (iv) fortement diminuée, ont été utilisés pour générer des directions de changement probables. Les relations actuelles entre la distribution des espèces et l'environnement ont été déterminées par le biais des modèles statistiques et ont été utilisées pour calculer des probabilités de présence selon les scénarios dans six régions biogéographiques majeures de la Suisse, comportant le Jura, le Plateau, les Alpes du Nord, du Sud, centrales orientales et centrales occidentales. Les résultats principaux ont montré que les espèces des prairies sèches pourraient profiter d'une diminution de la production agricole, mais qu'elles pourraient aussi disparaître à cause de l'embroussaillement des terres ouvertes dû à la libéralisation des marchés agricoles. La probabilité de présence des espèces de marais a décrû à cause d'une perte générale des habitats favorables. De plus, les analyses ont confirmé que des causes humaines comme l'urbanisation, l'abandon des terres ouvertes et l'intensification de l'agriculture affectent actuellement ces espèces. Ainsi ces forces devraient être mieux prises en compte lors de planifications paysagères, pour que ces papillons diurnes puissent survivre dans leurs habitats actuels. Dans ce travail de thèse, des données historiques ont été intensivement utilisées pour reconstruire des développements anciens et pour les rendre utiles à des recherches contemporaines. Cependant, la disponibilité des données historiques et les analyses à grande échelle ont souvent limité le pouvoir explicatif des analyses. Des descripteurs pertinents pour caractériser les habitats anciens et des données suffisantes sur la distribution des espèces sont généralement rares, spécialement pour des analyses à des échelles fores. Cette situation peut être améliorée en augmentant l'étendue du site d'étude et la résolution, comme il a été fait dans cette thèse en considérant toute la Suisse avec ses communes. Cependant, les récents projets de surveillance et les techniques de collecte de données sont des sources prometteuses, qui pourraient permettre des analyses plus détaillés sur les réactions à long terme des espèces aux changements de paysage dans le futur. Ce travail a aussi montré la valeur des anciennes données de distribution, par exemple leur potentiel pour aider à localiser des' présences d'espèces encore inconnues. Les résultats peuvent contribuer à des activités de recherche à venir, qui étudieraient les distributions récentes ou futures d'espèces en considérant l'immense richesse des données de distribution historiques.

Model development for capercaillie (Tetrao urogallus) habitat in the Jura Mountains (Western Switzerland)

Capercaillie, Tetrao urogallus, is a threatened species in central Europe, and Swiss populations declined 40 to 50 % between 1970 and 1985. Capercaillie are sensitive to forest structure, and loss of habitat is a major cause of their decline. Knowledge of habitat characteristics is therefore essential for capercaillie conservation. Here, we present models predicting capercaillie probability of occurrence, based on relevant structural habitat variables. Models were built using multiple logistic regression analyses on capercaillie presence/absence data. Vegetation survey was carried out in July 1999 in a 170-km2 forested area (Jura mountains, canton de Vaud, western Switzerland) inhabited by capercaillie and presence/absence of the species was assessed according to dropping presence/absence. The survey was based on 10-m-radius sample plots each in a 1-km2 forest patch (n = 76 with capercaillie droppings, n = 80 without). A first model included seven out of 27 measured habitat variables and a second model only four. The latter model best represents practical needs. It includes three variables which had a negative impact on capercaillie presence: tree and shrub covers and spruce, Picea excelsa, shrub cover, and one which had a positive effect: bilberry, Vaccinium myrtillus, cover, highlighting that capaercaillie selected open forest with high bilberry abundance. The model can be used to map potential capercaillie habitat distribution and to manage the habitat in favour of capercaillie (protection and adapted forestry practices) in the Swiss Jura mountains.

Application of detection probabilities to the design of amphibian monitoring programs in temporary ponds

Failure to detect a species in an area where it is present is a major source of error in biological surveys. We assessed whether it is possible to optimize single-visit biological monitoring surveys of highly dynamic freshwater ecosystems by framing them a priori within a particular period of time. Alternatively, we also searched for the optimal number of visits and when they should be conducted. We developed single-species occupancy models to estimate the monthly probability of detection of pond-breeding amphibians during a four-year monitoring program. Our results revealed that detection probability was species-specific and changed among sampling visits within a breeding season and also among breeding seasons. Thereby, the optimization of biological surveys with minimal survey effort (a single visit) is not feasible as it proves impossible to select a priori an adequate sampling period that remains robust across years. Alternatively, a two-survey combination at the beginning of the sampling season yielded optimal results and constituted an acceptable compromise between sampling efficacy and survey effort. Our study provides evidence of the variability and uncertainty that likely affects the efficacy of monitoring surveys, highlighting the need of repeated sampling in both ecological studies and conservation management.

Evaluating the ability of habitat suitability models to predict species presences.

Models predicting species spatial distribution are increasingly applied to wildlife management issues, emphasising the need for reliable methods to evaluate the accuracy of their predictions. As many available datasets (e.g. museums, herbariums, atlas) do not provide reliable information about species absences, several presence-only based analyses have been developed. However, methods to evaluate the accuracy of their predictions are few and have never been validated. The aim of this paper is to compare existing and new presenceonly evaluators to usual presence/absence measures. We use a reliable, diverse, presence/absence dataset of 114 plant species to test how common presence/absence indices (Kappa, MaxKappa, AUC, adjusted D-2) compare to presenceonly measures (AVI, CVI, Boyce index) for evaluating generalised linear models (GLM). Moreover we propose a new, threshold-independent evaluator, which we call "continuous Boyce index". All indices were implemented in the B10MAPPER software. We show that the presence-only evaluators are fairly correlated (p > 0.7) to the presence/absence ones. The Boyce indices are closer to AUC than to MaxKappa and are fairly insensitive to species prevalence. In addition, the Boyce indices provide predicted-toexpected ratio curves that offer further insights into the model quality: robustness, habitat suitability resolution and deviation from randomness. This information helps reclassifying predicted maps into meaningful habitat suitability classes. The continuous Boyce index is thus both a complement to usual evaluation of presence/absence models and a reliable measure of presence-only based predictions.

New anthropometry-based age- and sex-specific reference values for urinary 24-hour creatinine excretion based on the adult Swiss population (epublish)

BACKGROUND: Urinary creatinine excretion is used as a marker of completeness of timed urine collections, which are a keystone of several metabolic evaluations in clinical investigations and epidemiological surveys. METHODS: We used data from two independent Swiss cross-sectional population-based studies with standardised 24-hour urinary collection and measured anthropometric variables. Only data from adults of European descent, with estimated glomerular filtration rate (eGFR) ≥60 ml/min/1.73 m2 and reported completeness of the urinary collection were retained. A linear regression model was developed to predict centiles of the 24-hour urinary creatinine excretion in 1,137 participants from the Swiss Survey on Salt and validated in 994 participants from the Swiss Kidney Project on Genes in Hypertension. RESULTS: The mean urinary creatinine excretion was 193 ± 41 μmol/kg/24 hours in men and 151 ± 38 μmol/kg/24 hours in women in the Swiss Survey on Salt. The values were inversely correlated with age and body mass index (BMI). CONCLUSIONS: We propose a validated prediction equation for 24-hour urinary creatinine excretion in the general European population, based on readily available variables such as age, sex and BMI, and a few derived normograms to ease its clinical application. This should help healthcare providers to interpret the completeness of a 24-hour urine collection in daily clinical practice and in epidemiological population studies.

Case management for frequent users of the emergency department: study protocol of a randomised controlled trial.

BACKGROUND: We devised a randomised controlled trial to evaluate the effectiveness and efficiency of an intervention based on case management care for frequent emergency department users. The aim of the intervention is to reduce such patients' emergency department use, to improve their quality of life, and to reduce costs consequent on frequent use. The intervention consists of a combination of comprehensive case management care and standard emergency care. It uses a clinical case management model that is patient-identified, patient-directed, and developed to provide high intensity services. It provides a continuum of hospital- and community-based patient services, which include clinical assessment, outreach referral, and coordination and communication with other service providers. METHODS/DESIGN: We aim to recruit, during the first year of the study, 250 patients who visit the emergency department of the University Hospital of Lausanne, Switzerland. Eligible patients will have visited the emergency department 5 or more times during the previous 12 months. Randomisation of the participants to the intervention or control groups will be computer generated and concealed. The statistician and each patient will be blinded to the patient's allocation. Participants in the intervention group (N = 125), additionally to standard emergency care, will receive case management from a team, 1 (ambulatory care) to 3 (hospitalization) times during their stay and after 1, 3, and 5 months, at their residence, in the hospital or in the ambulatory care setting. In between the consultations provided, the patients will have the opportunity to contact, at any moment, the case management team. Participants in the control group (N = 125) will receive standard emergency care only. Data will be collected at baseline and 2, 5.5, 9, and 12 months later, including: number of emergency department visits, quality of life (EuroQOL and WHOQOL), health services use, and relevant costs. Data on feelings of discrimination and patient's satisfaction will also be collected at the baseline and 12 months later. DISCUSSION: Our study will help to clarify knowledge gaps regarding the positive outcomes (emergency department visits, quality of life, efficiency, and cost-utility) of an intervention based on case management care. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01934322.

Influence of age on retinochoroidal healing processes after argon photocoagulation in C57bl/6j mice.

PURPOSE: To analyze the influence of age on retinochoroidal wound healing processes and on glial growth factor and cytokine mRNA expression profiles observed after argon laser photocoagulation. METHODS: A cellular and morphometric study was performed that used 44 C57Bl/6J mice: 4-week-old mice (group I, n=8), 6-week-old mice (group II, n=8), 10-12-week-old mice (group III, n=14), and 1-year-old mice (group IV, n=14). All mice in these groups underwent a standard argon laser photocoagulation (50 microm, 400 mW, 0.05 s). Two separated lesions were created in each retina using a slit lamp delivery system. At 1, 3, 7, 14, 60 days, and 4 months after photocoagulation, mice from each of the four groups were sacrificed by carbon dioxide inhalation. Groups III and IV were also studied at 6, 7, and 8 months after photocoagulation. At each time point the enucleated eyes were either mounted in Tissue Tek (OCT), snap frozen and processed for immunohistochemistry or either flat mounted (left eyes of groups III and IV). To determine, by RT-PCR, the time course of glial fibrillary acidic protein (GFAP), vascular endothelial growth factor (VEGF), and monocyte chemotactic protein-1 (MCP-1) gene expression, we delivered ten laser burns (50 microm, 400 mW, 0.05 s) to each retina in 10-12-week-old mice (group III', n=10) and 1-year-old mice (group IV', n=10). Animals from Groups III' and IV' had the same age than those from Groups III and IV, but they received ten laser impacts in each eye and served for the molecular analysis. Mice from Groups III and IV received only two laser impacts per eye and served for the cellular and morphologic study. Retinal and choroidal tissues from these treated mice were collected at 16 h, and 1, 2, 3, and 7 days after photocoagulation. Two mice of each group did not receive photocoagulation and were used as controls. RESULTS: In the cellular and morphologic study, the resultant retinal pigment epithelium interruption expanse was significantly different between the four groups. It was more concise and smaller in the oldest group IV (112.1 microm+/-11.4 versus 219.1 microm+/-12.2 in group III) p<0.0001 between groups III and IV. By contrast, while choroidal neovascularization (CNV) was mild and not readily identifiable in group I, at all time points studied, CNV was more prominent in the (1-year-old mice) Group IV than in the other groups. For instance, up to 14 days after photocoagulation, CNV reaction was statistically larger in group IV than in group III ((p=0.0049 between groups III and IV on slide sections and p<0.0001 between the same groups on flat mounts). Moreover, four months after photocoagulation, the CNV area (on slide sections) was 1,282 microm(2)+/-90 for group III and 2,999 microm(2)+/-115 for group IV (p<0.0001 between groups III and IV). Accordingly, GFAP, VEGF, and MCP-1 mRNA expression profiles, determined by RT-PCR at 16 h, 1, 2, 3, and 7 days postphotocoagulation, were modified with aging. In 1-year-old mice (group IV), GFAP mRNA expression was already significantly higher than in the younger (10-12 week) group III before photocoagulation. After laser burns, GFAP mRNA expression peaked at 16-24 h and on day 7, decreasing thereafter. VEGF mRNA expression was markedly increased after photocoagulation in old mice eyes, reaching 2.7 times its basal level at day 3, while it was only slightly increased in young mice (1.3 times its level in untreated young mice 3 days postphotocoagulation). At all time points after photocoagulation, MCP-1 mRNA expression was elevated in old mice, reaching high levels of expression at 16 h and day 3 respectively. CONCLUSIONS: Our results were based on the study of four different age groups and included not only data from morphological observations but also from a molecular analysis of the various alterations of cytokine signaling and expression. One-year-old mice demonstrated more extensive CNV formation and a slower pace of regression after laser photocoagulation than younger mice. These were accompanied by differences in growth factors and cytokine expression profiles indicate that aging is a factor that aggravates CNV. The above results may provide some insight into possible therapeutic strategies in the future.