Unfolding nature "in silico" : overcoming practical and metodological limitations of species distribution models

Abstract : The existence of a causal relationship between the spatial distribution of living organisms and their environment, in particular climate, has been long recognized and is the central principle of biogeography. In turn, this recognition has led scientists to the idea of using the climatic, topographic, edaphic and biotic characteristics of the environment to predict its potential suitability for a given species or biological community. In this thesis, my objective is to contribute to the development of methodological improvements in the field of species distribution modeling. More precisely, the objectives are to propose solutions to overcome limitations of species distribution models when applied to conservation biology issues, or when .used as an assessment tool of the potential impacts of global change. The first objective of my thesis is to contribute to evidence the potential of species distribution models for conservation-related applications. I present a methodology to generate pseudo-absences in order to overcome the frequent lack of reliable absence data. I also demonstrate, both theoretically (simulation-based) and practically (field-based), how species distribution models can be successfully used to model and sample rare species. Overall, the results of this first part of the thesis demonstrate the strong potential of species distribution models as a tool for practical applications in conservation biology. The second objective this thesis is to contribute to improve .projections of potential climate change impacts on species distributions, and in particular for mountain flora. I develop and a dynamic model, MIGCLIM, that allows the implementation of dispersal limitations into classic species distribution models and present an application of this model to two virtual species. Given that accounting for dispersal limitations requires information on seed dispersal, distances, a general methodology to classify species into broad dispersal types is also developed. Finally, the M~GCLIM model is applied to a large number of species in a study area of the western Swiss Alps. Overall, the results indicate that while dispersal limitations can have an important impact on the outcome of future projections of species distributions under climate change scenarios, estimating species threat levels (e.g. species extinction rates) for a mountainous areas of limited size (i.e. regional scale) can also be successfully achieved when considering dispersal as unlimited (i.e. ignoring dispersal limitations, which is easier from a practical point of view). Finally, I present the largest fine scale assessment of potential climate change impacts on mountain vegetation that has been carried-out to date. This assessment involves vegetation from 12 study areas distributed across all major western and central European mountain ranges. The results highlight that some mountain ranges (the Pyrenees and the Austrian Alps) are expected to be more affected by climate change than others (Norway and the Scottish Highlands). The results I obtain in this study also indicate that the threat levels projected by fine scale models are less severe than those derived from coarse scale models. This result suggests that some species could persist in small refugias that are not detected by coarse scale models. Résumé : L'existence d'une relation causale entre la répartition des espèces animales et végétales et leur environnement, en particulier le climat, a été mis en évidence depuis longtemps et est un des principes centraux en biogéographie. Ce lien a naturellement conduit à l'idée d'utiliser les caractéristiques climatiques, topographiques, édaphiques et biotiques de l'environnement afin d'en prédire la qualité pour une espèce ou une communauté. Dans ce travail de thèse, mon objectif est de contribuer au développement d'améliorations méthodologiques dans le domaine de la modélisation de la distribution d'espèces dans le paysage. Plus précisément, les objectifs sont de proposer des solutions afin de surmonter certaines limitations des modèles de distribution d'espèces dans des applications pratiques de biologie de la conservation ou dans leur utilisation pour évaluer l'impact potentiel des changements climatiques sur l'environnement. Le premier objectif majeur de mon travail est de contribuer à démontrer le potentiel des modèles de distribution d'espèces pour des applications pratiques en biologie de la conservation. Je propose une méthode pour générer des pseudo-absences qui permet de surmonter le problème récurent du manque de données d'absences fiables. Je démontre aussi, de manière théorique (par simulation) et pratique (par échantillonnage de terrain), comment les modèles de distribution d'espèces peuvent être utilisés pour modéliser et améliorer l'échantillonnage des espèces rares. Ces résultats démontrent le potentiel des modèles de distribution d'espèces comme outils pour des applications de biologie de la conservation. Le deuxième objectif majeur de ce travail est de contribuer à améliorer les projections d'impacts potentiels des changements climatiques sur la flore, en particulier dans les zones de montagnes. Je développe un modèle dynamique de distribution appelé MigClim qui permet de tenir compte des limitations de dispersion dans les projections futures de distribution potentielle d'espèces, et teste son application sur deux espèces virtuelles. Vu que le fait de prendre en compte les limitations dues à la dispersion demande des données supplémentaires importantes (p.ex. la distance de dispersion des graines), ce travail propose aussi une méthode de classification simplifiée des espèces végétales dans de grands "types de disperseurs", ce qui permet ainsi de d'obtenir de bonnes approximations de distances de dispersions pour un grand nombre d'espèces. Finalement, j'applique aussi le modèle MIGCLIM à un grand nombre d'espèces de plantes dans une zone d'études des pré-Alpes vaudoises. Les résultats montrent que les limitations de dispersion peuvent avoir un impact considérable sur la distribution potentielle d'espèces prédites sous des scénarios de changements climatiques. Cependant, quand les modèles sont utilisés pour évaluer les taux d'extinction d'espèces dans des zones de montages de taille limitée (évaluation régionale), il est aussi possible d'obtenir de bonnes approximations en considérant la dispersion des espèces comme illimitée, ce qui est nettement plus simple d'un point dé vue pratique. Pour terminer je présente la plus grande évaluation à fine échelle d'impact potentiel des changements climatiques sur la flore des montagnes conduite à ce jour. Cette évaluation englobe 12 zones d'études réparties sur toutes les chaines de montages principales d'Europe occidentale et centrale. Les résultats montrent que certaines chaines de montagnes (les Pyrénées et les Alpes Autrichiennes) sont projetées comme plus sensibles aux changements climatiques que d'autres (les Alpes Scandinaves et les Highlands d'Ecosse). Les résultats obtenus montrent aussi que les modèles à échelle fine projettent des impacts de changement climatiques (p. ex. taux d'extinction d'espèces) moins sévères que les modèles à échelle large. Cela laisse supposer que les modèles a échelle fine sont capables de modéliser des micro-niches climatiques non-détectées par les modèles à échelle large.

Understanding the ecological niche to elucidate spatial strategies of the southernmost Tupinambis lizards

Understanding factors that shape ranges of species is central in evolutionary biology. Species distribution models have become important tools to test biogeographical, ecological and evolutionary hypotheses. Moreover, from an ecological and evolutionary perspective, these models help to elucidate the spatial strategies of species at a regional scale. We modelled species distributions of two phylogenetically, geographically and ecologically close Tupinambis species (Teiidae) that occupy the southernmost area of the genus distribution in South America. We hypothesized that similarities between these species might have induced spatial strategies at the species level, such as niche differentiation and divergence of distribution patterns at a regional scale. Using logistic regression and MaxEnt we obtained species distribution models that revealed interspecific differences in habitat requirements, such as environmental temperature, precipitation and altitude. Moreover, the models obtained suggest that although the ecological niches of Tupinambis merianae and T. rufescens are different, these species might co-occur in a large contact zone. We propose that niche plasticity could be the mechanism enabling their co-occurrence. Therefore, the approach used here allowed us to understand the spatial strategies of two Tupinambis lizards at a regional scale.

Plant species richness and environmental heterogeneity in a mountain landscape: effects of variability and spatial configuration

The loss of biodiversity has become a matter of urgent concern and a better understanding of local drivers is crucial for conservation. Although environmental heterogeneity is recognized as an important determinant of biodiversity, this has rarely been tested using field data at management scale. We propose and provide evidence for the simple hypothesis that local species diversity is related to spatial environmental heterogeneity. Species partition the environment into habitats. Biodiversity is therefore expected to be influenced by two aspects of spatial heterogeneity: 1) the variability of environmental conditions, which will affect the number of types of habitat, and 2) the spatial configuration of habitats, which will affect the rates of ecological processes, such as dispersal or competition. Earlier, simulation experiments predicted that both aspects of heterogeneity will influence plant species richness at a particular site. For the first time, these predictions were tested for plant communities using field data, which we collected in a wooded pasture in the Swiss Jura mountains using a four-level hierarchical sampling design. Richness generally increased with increasing environmental variability and "roughness" (i.e. decreasing spatial aggregation). Effects occurred at all scales, but the nature of the effect changed with scale, suggesting a change in the underlying mechanisms, which will need to be taken into account if scaling up to larger landscapes. Although we found significant effects of environmental heterogeneity, other factors such as history could also be important determinants. If a relationship between environmental heterogeneity and species richness can be shown to be general, recently available high-resolution environmental data can be used to complement the assessment of patterns of local richness and improve the prediction of the effects of land use change based on mean site conditions or land use history.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Rarity types among plant species with high conservation priority in Switzerland

We investigated the ecogeographic characteristics of 118 Swiss plant species listed as those deserving highest conservation priority in a national conservation guide and classified them into the seven Rabinowitz' rarity types, taking geographic distribution, habitat rarity and local population size into account. Our analysis revealed that species with high conservation priority in Switzerland mostly have a very restricted geographic distribution in Switzerland and generally occur in rare habitats, but do not necessarily constitute small populations and are generally not endemics on a global scale. Moreover, species that are geographically very restricted on a regional scale are not generally restricted on a global scale. By analysing relationships between rarity and IUCN extinction risks for Switzerland, we demonstrated that species with the highest risk of extinction are those with the most restricted geographic distribution; whereas species with lower risk of extinction (but still high conservation priority) include many regional endemics. Habitat rarity and local population size appeared to be of minor importance for the assessment of extinction risk in Switzerland, but the total number of fulfilled rarity criteria still correlated positively with the severity of extinction risk. Our classification is the first preliminary assessment of the relative importance of each rarity type among endangered plant species of the Swiss flora and our results underline the need to distinguish between a regional and a global responsibility for the conservation of rare and endangered species.

Evaluation des pratiques d'irrigation, de fertilisation et d'application de pesticides dans l'agriculture périurbaine de Dakar, Sénégal

RESUMEL'agriculture urbaine et périurbalne - nommée ci-après AU - est un thème fort de recherche transversale, au vu des nombreux enjeux économiques, sociaux et environnementaux. L'objectif de cette recherche était de contribuer à une meilleure connaissance des processus de transfert de polluants et du cycle des nutriments à l'échelle locale, afin de déterminer sous quelles conditions l'AU de Dakar peut être pratiquée sans porter atteinte à la santé et à l'environnement.Une approche basée sur l'étude des processus géochimiques dans ie sol jusqu'à la nappe a été choisie, à l'échelle de la parcelle cultivée et à une échelle un peu plus large de la zone périurbaine de Dakar pour déterminer les influences du type d'occupation du sol.L'évaluation des impacts de l'irrigation avec des eaux usées brutes et des eaux de nappe saumâtres sur la qualité des sols (chapitre 2) a montré que l'alcalinité et les teneurs en calcium élevées des eaux saumâtres induisent la précipitation de CaC03 dans l'horizon superficiel du sol. Na remplace consécutivement Ca sur le complexe argilo-humique du sol et les bicarbonates diminuent dans la solution du sol. Le carbone organique dissout (COD) augmente significativement dans la solution du sol et dans la nappe sous-jacente. Malgré l'alcalinité et les teneurs très élevées en calcium des eaux usées, il y a peu de précipitation de CaC03 dans l'horizon superficiel du sol et une faible augmentation du sodium échangeable ESP. La nitrification de l'ammonium des eaux usées (moy 190mg/L à Pikine) produit des protons, qui ne sont plus tamponnés par les bicarbonates exportés hors du profil. Il y a alors une nette baisse de pH des sols irrigués par des eaux usées non traitées. Les sols irrigués par des eaux usées et saumâtres stockent moins de C et Ν que les sols de référence.L'évaluation de l'influence de l'occupation des sols en zone périurbaine sur à la nappe phréatique peu profonde (chapitre 3) a permis de déterminer les traceurs représentatifs de l'occupation du sol, à savoir Br/CI, NO3/CI et δ180-Ν03 pour l'irrigation par des eaux usées, pH et δ15Ν-Ν03 pour l'irrigation par des eaux de nappe, et Rb+Cr et Κ pour les lixiviats de fosses septiques. Ce chapitre a mis en évidence des points importants de la dynamique de l'azote en zone périurbaine sous deux occupations du sol : (1) la dénitrification est un processus important dans l'agrosystème périurbain de Dakar en bas de dune, dans les gleysols où l'on trouve des conditions temporairement réduites et un substrat organique favorables aux microorganismes de la dénitrification. Les teneurs en nitrates sont presque nulles avec irrigation d'eau de nappe. (2) en bas de pente, mais avec irrigation quotidienne par les eaux usées, l'apport continu d'ammonium inhibe probablement la dénitrification, mais favorise la volatilisation. (3) la nitrification de l'ammonium dans la nappe lors de la lixiviation des fosses septiques se distingue de la nitrification de l'ammonium dans la zone non saturée dans la zone d'agriculture périurbaine par la composition isotopique de l'oxygène de l'eau. Une comparaison des flux d'azote entre l'agrosystème et les quartiers périurbains de Dakar (chapitre 4) ont révélé que ces derniers étaient du même ordre de grandeur par unité de surface, à savoir 2-4 tonnes Ν /ha/an.L'évaluation des flux de pesticides dans l'agrosystème et des risques induits pour les eaux souterraines (chapitre 5) a révélé un fiux total de pesticides de 60kg/ha/an, totalisant 15 matières actives. Seules deux de ces matières actives sont autorisées par le comité des pesticides sahélien. Les pesticides les plus utilisés par les producteurs sont l'organochloré dicofol, les organophosphorés methamidophos, dimethoate et fenithrotion ainsi que le cabamate methomyl. Les flux les plus importants sont de 9 à 7 kg/ha/an (methomyl, methamidophos, ethoprophos et dicofol). Les pesticides qui présentent un risque élevé de contamination des eaux souterraines et qui devraient être prioritaires pour un suivi analytique sont : le carbofuran, le dimethoate, l'ethoprophos et le methomyl.En conclusion, une meilleure gestion de la fertilisation est nécessaire dans la zone d'AU de Dakar, afin de (1) réduire les pertes gazeuses qui contribuent à l'effet de serre, (2) de ralentir la minéralisation du carbone et de l'azote organiques pour créer un stock de C et Ν dans ces sols, (3) de limiter le lessivage dans la nappe et enfin, 4) d'augmenter l'efficacité d'utilisation de Ν par les plantes. Une optimisation de l'irrigation devrait limiter l'alcalinisation secondaire. Enfin, la mise en place d'une lutte intégrée ou biologique contre les ravageurs est indispensable afin de minimiser les risques pour les eaux souterraines et les mares permanentes.ABSTRACTUrban and periurban agriculture (UA) is an important issue in southern countries, because of its key role in their social and economical development and its environmental concern. The goal of this study was to contribute to a better understanding of pollutant transfer and nutrient cycling at the local scale, in order to implement the necessary improvements to guarantee the sustainability of this practice.An approach based on geochemical processes occurring in the vadose zone from the surface down to the groundwater level was chosen, at the scale of cultivated plots and at the regional scale of Dakar periurban areas, to determine the influence of land use.The assessment of irrigation with untreated domestic wastewater and brackish water on soil quality (chapter 2) showed: (1) that the high alkalinity and calcium contents of brackish water induce CaC03 precipitation in the top layer of the soil and therefore a replacement of Ca by Na on the clay- humic complexes, strongly marked during the dry season. Dissolved organic carbon (DOC) increased significantly in the soil solution and in the underlying groundwater. (2) in spite of the similarly high alkalinity and Ca contents of waste water, there is only little CaC03 precipitation and a low increase of the percentage of exchangeable sodium (ESP) in the soil top layer. The nitrification of the ammonium of wastewater (mean 190 mg/L in Pikine) produces protons, which are not any more buffered by bicarbonates exported out of the soil profile, which leads to a net decline of soil pH. Both soils irrigated with untreated wastewater and brackish water store less of C and Ν than soils irrigated with non saline groundwater.The assessment of the impact of land use on the shallow groundwater (chapter 3) allowed determining representative tracers of the land use. Low Br/CI ratio, high NO3/CI ratio and low δ1βΟ- nitrate indicated the influence of wastewater; high pH and high 515N-nitrates indicated the influence of brackish water together with high amendments of organic fertilizers; high Rb+Cr and Κ indicated the influence of poor sanitation facilities in periurban districts (septic tank leakage). This chapter also pointed out the following facts about the nitrogen dynamics : (1) denitrification is a key-process in the Dakar UA agrosystem in the gleysols irrigated with groundwater. The underlying groundwater is almost nitrate free. (2) in the Gleysols irrigated with waste water, ammonium inhibits denitrification but facilitate ammoniac volatilization. A comparison of nitrogen balance between the UA agrosystem and the periurban districts of Dakar (chapter 4) revealed similar flows per surface unit, namely 2-4 tons Ν / ha / year.The evaluation of pesticides use in the UA agrosystem and the risk assessment for the groundwater (chapter 5) revealed a total flow of pesticides of 60kg / ha / year, totalizing 15 active substances. Only two of these are authorized by the Sahelian Pesticides Committee. The most used pesticides are dicofol (organochlorinated), methamidophos, dimethoate and fenithrotion (organophosphate) as well as methomyl. (carbamate). The most important flows vary between 9 to 7 kg / ha / year. Pesticides with a high risk of groundwater contamination - according to SIRIS and EPRIP 2 indicators - are: carbofuran, dimethoate, ethoprophos and methomyl. These substances should be established as a priority for an analytical follow-up in the different environmental compartments.In conclusion, a better management of the fertilization is necessary in the Dakar UA, (1) to reduce the gaseous losses which contribute to greenhouse emissions (2) to slow down the mineralization of the organic carbon and the nitrogen, in order to enhance the C and Ν stock in these soils, (3) to limit the nitrate leaching in the groundwater and finally, 4) to increase the N-use efficiency of plants. An optimization of the irrigation scheme should limit the secondary sodisation if coupled with an increase the stable organic matter of the soil. An integrated or biologic crop pest strategy is urgently needed to minimize risks with respect to ground and surface water (ponds used for fishing).RESUME LARGE PUBLICL'agriculture mondiale connaît actuellement une crise majeure, affectée par les changements climatiques, la sécurité alimentaire et les dégradations de l'environnement. Elle n'a plus le rôle unique de produire, mais devient un élément essentiel de la protection des ressources naturelles et du paysage. Les politiques agricoles basées sur les marchés mondiaux devront se réorienter vers une agriculture locale basée sur le développement durable.La production alimentaire située dans l'enceinte des villes, nommée agriculture urbaine ou périurbaine (AU ci-après) joue un rôle important dans le contexte actuel d'accroissement de la population et de la pauvreté urbaines. L'AU concerne en effet la majorité des mégapoies du monde, fait vivre plus de 200 millions de personnes dans les pays du Sud, fournit jusqu'à 80% de la demande urbaine en certains produits frais, fait barrière à l'extension urbaine et permet un recyclage de certains déchets urbains. L'AU a pour particularité d'être à cheval entre des politiques rurales et urbaines, d'où un délaissement ce cette activité au secteur informel. Ce qui a développé de nombreuses stratégies à risques, comme à Dakar, où les petits producteurs périurbains irriguent quotidiennement avec des eaux usées domestiques par manque d'accès à une eau de bonne qualité et pour raccourcir les cycles de production. L'extrême précarité foncière des acteurs de l'AU de Dakar les empêchent d'investir à long terme et induit des pratiques inadéquates d'irrigation, d'usage de pesticides et de fertilisation de ces sols sableux.L'objectif de cette recherche était de contribuer à une meilleure connaissance des processus de transfert de polluants et du cycle des nutriments à l'échelle des parcelles cultivées par des eaux usées et des eaux saumâtres, afin de déterminer sous quelles conditions l'AU de Dakar peut être pratiquée et surtout maintenue sans porter atteinte à la santé et à l'environnement. Pour cela, une approche basée sur l'étude des processus géochimiques dans le sol jusqu'à la nappe a été choisie, à l'échelle de la parcelle cultivée et à une échelle un peu plus large de la zone périurbaine de Dakar pour déterminer les influences du type d'occupation du sol.Les résultats principaux de cette étude ont montré que (1) il y a un processus de salinisation anthropique des sols (sodisation) lors d'irrigation avec des eaux de nappe saumâtres, un processus accentué en saison sèche et lors d'années à pluviométrie déficitaire. Bien que les eaux usées soient aussi salines que les eaux de nappe, la salinisation des sols irrigués' par des eaux usées est limitée par l'ammonium présent dans les eaux usées (moy 190mg NH4/L à Pikine) qui produit de l'acidité lors de la transformation en nitrates dans le sol (nitrification). (2) les sols irrigués par des eaux usées (EU) stockent moins de C et Ν que les sois de référence, ce qui montrent bien que l'azote des eaux usées n'est pas disponible pour les plantes, mais est lessivé dans la nappe (100 à 450 mg/L N03 sous irrigation par EU, alors que la limite de OMS est de 50mg/L). (3) l'utilisation des isotopes stables des nitrates et des éléments traces, notamment le bore et le brome, ont permis de distinguer l'influence de l'irrigation par des eaux usées, de l'irrigation par des eaux de nappe et des lixiviats de fosses septiques sur les propriétés de la nappe. (4) Le processus de la dénitrification (atténuation naturelle des concentrations en nitrates de la nappe par biotransformation en azote gazeux) est important dans les zones basses de l'agrosystème périurbain de Dakar, sous irrigation par eaux naturelles (ΝΟ3 < 50mg/L). Tandis que sous habitat sans assainissement adéquat, les nitrates atteignent 300 à 700 mg/L. (5) Le flux total de pesticides dans l'AU est énorme (60kg/ha/an) totalisant 15 pesticides, dont deux seulement sont autorisés. Les pesticides les plus utilisés sont des insecticides organophosphorés et organochlorés classés extrêmement dangereux à dangereux par l'OMS, appliqués à des doses de 2 à 9 kg/ha/an. Les pesticides qui ont montré un risque élevé de contamination des eaux souterraines avec les indicateurs SIRIS et EPRIP2 sont : le carbofuran, le dimethoate, l'ethoprophos et le methomyl.En conclusion, nous recommandons la reconstitution d'un horizon superficiel des sols riche en matière organique stable et structuré par production locale de compost. Cette mesure réduira les pertes gazeuses contribuant à l'effet de serre, augmentera le stock de Ν dans ces sols, alors utilisable par les plantes et permettra de diminuer l'irrigation car la capacité de rétention de l'eau dans le sol sera accru, ce qui limitera le lessivage des nitrates dans la nappe et l'alcalinisation secondaire. Enfin, la mise en place d'une lutte intégrée ou biologique contre les ravageurs est indispensable afin de minimiser les risques pour les eaux souterraines et lesmares permanentes.

Rockfall characterisation and structural protection: A review

Rockfall is an extremely rapid process involving long travel distances. Due to these features, when an event occurs, the ability to take evasive action is practically zero and, thus, the risk of injury or loss of life is high. Damage to buildings and infrastructure is quite likely. In many cases, therefore, suitable protection measures are necessary. This contribution provides an overview of previous and current research on the main topics related to rockfall. It covers the onset of rockfall and runout modelling approaches, as well as hazard zoning and protection measures. It is the aim of this article to provide an in-depth knowledge base for researchers and practitioners involved in projects dealing with the rockfall protection of infrastructures, who may work in the fields of civil or environmental engineering, risk and safety, the earth and natural sciences.

What drives invasibility? A multi-model inference test and spatial modelling of alien plant species richness patterns in northern Portugal

Understanding and anticipating biological invasions can focus either on traits that favour species invasiveness or on features of the receiving communities, habitats or landscapes that promote their invasibility. Here, we address invasibility at the regional scale, testing whether some habitats and landscapes are more invasible than others by fitting models that relate alien plant species richness to various environmental predictors. We use a multi-model information-theoretic approach to assess invasibility by modelling spatial and ecological patterns of alien invasion in landscape mosaics and testing competing hypotheses of environmental factors that may control invasibility. Because invasibility may be mediated by particular characteristics of invasiveness, we classified alien species according to their C-S-R plant strategies. We illustrate this approach with a set of 86 alien species in Northern Portugal. We first focus on predictors influencing species richness and expressing invasibility and then evaluate whether distinct plant strategies respond to the same or different groups of environmental predictors. We confirmed climate as a primary determinant of alien invasions and as a primary environmental gradient determining landscape invasibility. The effects of secondary gradients were detected only when the area was sub-sampled according to predictions based on the primary gradient. Then, multiple predictor types influenced patterns of alien species richness, with some types (landscape composition, topography and fire regime) prevailing over others. Alien species richness responded most strongly to extreme land management regimes, suggesting that intermediate disturbance induces biotic resistance by favouring native species richness. Land-use intensification facilitated alien invasion, whereas conservation areas hosted few invaders, highlighting the importance of ecosystem stability in preventing invasions. Plants with different strategies exhibited different responses to environmental gradients, particularly when the variations of the primary gradient were narrowed by sub-sampling. Such differential responses of plant strategies suggest using distinct control and eradication approaches for different areas and alien plant groups.

Réflexions sur les céramiques fines en Gaule interne de La Tène finale au début du Haut-Empire : comparaisons d'ensembles et essai de caractérisation des phénomènes d'acculturation

Le concept de romanisation ne revêt d'utilité que s'il peut être qualifié, et donc analysé en différents « niveaux », afin de comparer les rythmes et les formes du phénomène dans les régions de l'Empire. Avant même de chercher à tenir un discours global sur la romanisation, il faut donc en étudier les multiples aspects de manière spécifique. C'est dans ce cadre que s'inscrit cette thèse de doctorat dont le but est de quantifier et de qualifier l'acculturation des céramiques fines en Gaule interne (200 BC - 50 AD). Fondée sur un corpus de quelques 130 ensembles, provenant de 21 sites urbains localisés entre le centre ouest de la France, l'Auvergne, le Plateau suisse et le Luxembourg, cette étude se caractérise par la mise en place de nouvelles approches méthodologiques, s'appuyant sur le calcul d'un « indice de romanisation » (quantifiant la romanisation d'un ensemble sur une échelle de 1 à 100) et sur l'analyse des modalités du phénomène. L'examen diachronique et cartographique des résultats obtenus montrant de fortes disparités dans les processus d'acculturation des céramiques fines, des niveaux de romanisation ont pu être définis, dont le principal intérêt est d'ouvrir la comparaison avec les autres aspects de l'acculturation des Gaules (instrumentum culinaire, technique de construction, romanisation de l'anthroponymie...). The concept of Romanisation can only be considered useful if it can be qualified and analysed at different levels in order to compare the acculturation rates, frequency and form of the phenomenon among the regions in the Empire. But before an overview of Romanisation can be considered, its multiple aspects must be examined in a specific manner. Among the French speaking areas, it is only recently that research concerning Gaulish and Roman periods has risked expanding the scope to a supra regional scale. Against this background, the purpose of my doctoral thesis is to examine the qualitative and quantitative aspects of acculturation among Gaulish fine wares (200 BCE - 50 CE). Based on a corpus of 130 assemblages from 21 urban sites from the western centre of France, Auvergne, the Swiss plateau and Luxembourg, the study is characterised by new methodological techniques relying on the calculation of a 'Romanisation index' (quantifying the level of Romanisation of an assemblage on a scale of 1 to 100) and the analysis of the methods of the phenomenon. The diachronic and cartographic examination of the results obtained showed large disparities in the acculturation processes concerning fine wares. The level of Romanisation could be defined, thus allowing a comparison with other aspects of the acculturation process among other aspects of Gaulish culture such as tools, cuisine, construction techniques and Romanisation of anthroponyms.

Comparative phylogeography in a specific and obligate pollination antagonism.

In specific and obligate interactions the nature and abundance of a given species can have important effects on the survival and population dynamics of associated organisms. In a phylogeographic framework, we therefore expect that the fates of organisms interacting specifically are also tightly interrelated. Here we investigate such a scenario by analyzing the genetic structures of species interacting in an obligate plant-insect pollination lure-and-trap antagonism, involving Arum maculatum (Araceae) and its specific psychodid (Diptera) visitors Psychoda phalaenoides and Psycha grisescens. Because the interaction is asymmetric (i.e., only the plant depends on the insect), we expect the genetic structure of the plant to be related with the historical pollinator availability, yielding incongruent phylogeographic patterns between the interacting organisms.Using insect mtDNA sequences and plant AFLP genome fingerprinting, we inferred the large-scale phylogeographies of each species and the distribution of genetic diversities throughout the sampled range, and evaluated the congruence in their respective genetic structures using hierarchical analyses of molecular variances (AMOVA). Because the composition of pollinator species varies in Europe, we also examined its association with the spatial genetic structure of the plant.Our findings indicate that while the plant presents a spatially well-defined genetic structure, this is not the case in the insects. Patterns of genetic diversities also show dissimilar distributions among the three interacting species. Phylogeographic histories of the plant and its pollinating insects are thus not congruent, a result that would indicate that plant and insect lineages do not share the same glacial and postglacial histories. However, the genetic structure of the plant can, at least partially, be explained by the type of pollinators available at a regional scale. Differences in life-history traits of available pollinators might therefore have influenced the genetic structure of the plant, the dependent organism, in this antagonistic interaction.

Measuring the relative effect of factors affecting species distribution model predictions

1. Species distribution models are increasingly used to address conservation questions, so their predictive capacity requires careful evaluation. Previous studies have shown how individual factors used in model construction can affect prediction. Although some factors probably have negligible effects compared to others, their relative effects are largely unknown. 2. We introduce a general "virtual ecologist" framework to study the relative importance of factors involved in the construction of species distribution models. 3. We illustrate the framework by examining the relative importance of five key factors-a missing covariate, spatial autocorrelation due to a dispersal process in presences/absences, sample size, sampling design and modeling technique-in a real study framework based on plants in a mountain landscape at regional scale, and show that, for the parameter values considered here, most of the variation in prediction accuracy is due to sample size and modeling technique. Contrary to repeatedly reported concerns, spatial autocorrelation has only comparatively small effects. 4. This study shows the importance of using a nested statistical framework to evaluate the relative effects of factors that may affect species distribution models.

Événements et déformations tardi-métamorphiques dans le segment Ossola-Ticino (Val Vigezzo-Centovalli, Italie-Suisse)

RESUME: Une zone tectonique large et complexe, connue sous le nom de ligne des Centovalli, traverse le secteur des Alpes Centrales compris entre Domodossola et Locarno. Cette région, formée par le Val Vigezzo et la vallée des Centovalli, constitue la terminaison méridionale du dôme Lepontin et représente une portion de la zone des racines des nappes alpines. Elle fait partie d’une grande et complexe zone de cisaillement, en partie associée à des phénomènes hydrothermaux d’âge alpin (<20 Ma), qui comprend le système tectonique Insubrien et celui du Simplon. Le Val Vigezzo et les Centovalli constituent un vrai carrefour entre les principaux accidents tectoniques des Alpes ainsi qu'une zone de juxtaposition du socle Sudalpin avec la zone des racines de l’Austroalpin et du Pennique. Les phases de déformation et les structures géologiques qui peuvent être étudiées s'étalent sur une période comprise entre environ 35 Ma et l'actuel. L’étude détaillée de terrain a mis en évidence la présence de nombreuses roches et structures de déformation de type ductile et cassant tels que des mylonites, des cataclasites, des pseudotachylites, des kakirites, des failles minéralisées, des gouges de faille et des plis. Sur le terrain on a pu distinguer au moins quatre générations de plis liés aux différentes phases de déformation. Le nombre et la complexité de ces structures indiquent une histoire très compliquée, selon plusieurs étapes distinctes, parfois liées, voire même superposées. Une partie de ces structures de déformation affectent aussi les dépôts sédimentaires d’âge quaternaire, notamment des limons et des sables lacustres. Ces sédiments constituent les restes d'un bassin lacustre attribué à l'époque interglaciaire Riss/Würm (éemien, 67.000-120.000 ans) et ils affleurent dans la partie centrale de la zone étudiée, à l'Est de la plaine de Santa Maria Maggiore. Ces sédiments montrent en leur sein toute une série de structures de déformation tels que des plans de faille inverses, des structures conjuguées de raccourcissement et des véritables plis. Ces failles et ces plis représenteraient les évidences de surface d’une déformation probablement active en époque quaternaire. Une autre formation rocheuse a retenu tout notre attention; il s'agit d'un corps de brèches péridotitiques monogéniques qui affleure en discontinuité le long du versant méridional et le long du fond de la vallée Vigezzo sur environ 20 km. Ces brèches se posent indifféremment sur le socle (unités Finero, Orselina) ou sur les sédiments lacustres. Elles sont traversées par des plans de failles qui développent des véritables stries de faille et des gouges de faille; l’orientation de ces plans est la même que celle affectant les failles à gouges du socle. La genèse de cette brèche est liée à l'altération et au modelage glacier (rock-glaciers) d'une brèche tectonique originelle qui borde la partie externe du Corps de Finero. Les structures de déformation de cette brèche, pareillement à celles des sédiments lacustres, ont été considérées comme les évidences de surface d'une tectonique quaternaire active dans la région. La dernière phase de déformation cassante qui affecte cette région peut donc être considérée comme active en époque quaternaire. Une vue d’ensemble de la région étudiée nous permet de reconnaître à l’échelle régionale une zone de cisaillement complexe orientée E-W, parallèlement à l’axe de la vallée Centovalli-Val Vigezzo. Les données de terrain, indiquent que cette zone de cisaillement débute sous conditions ductiles et évolue en plusieurs étapes jusqu’à des conditions de failles cassantes de surface. La reconstruction de l'évolution géodynamique de la région a permis de définir trois étapes distinctes qui marquent le passage, de ce secteur de socle cristallin, de conditions P-T profondes à des conditions de surface. Dans ce contexte, on a reconnu trois phases principales de déformation à l’échelle régionale qui caractérisent ces trois étapes. La phase la plus ancienne est constituée par des mylonites en faciès amphibolite, associées à des mouvements de cisaillement dextre, qui sont ensuite remplacés par des mylonites en faciès schistes verts et des plis rétrovergentes liés au rétrocharriage des nappes alpines. Une deuxième étape est identifiée par le développement d’une phase hydrothermale liée à un système de failles extensives et décrochantes dextres à direction principale E-W, NE-SW et NW-SE. Leur caractérisation minéralogique a permis la mise en évidence des phases cristallines de néoformation liées à cet événement constituées par : K-feldspath (microcline), chlorites (Fe+Mg), épidotes, prehnite, zéolites (laumontite), sphène, calcite. Dans ce contexte, pour obtenir une meilleure caractérisation de cet événement hydrothermal on a utilisé des géothermomètres sur chlorites, sensible aussi à la pression et a la a(H2O), qui ont donné des valeurs descendantes comprises entre 450-200°C. Les derniers mouvements sont mis en évidence par le développement d’une série de plans majeurs de failles à gouge, qui forment une structure en sigmoïdes d’épaisseur kilométrique reconnaissable à l’échelle de la vallée et caractérisée par des mouvements transpressifs avec une composante décrochante dextre toujours importante. Cette phase de déformation forme un système conjugué de failles avec direction moyenne E-W qui coupent la zone des racines des nappes alpines, la zone du Canavese et le corps ultramafique de Finero. Ce système se déroule de manière subparallèle à l'axe de la vallée le long de plusieurs dizaines de kilomètres. Une analyse complète et détaillée des gouges de faille par XRD a montré que la fraction argileuse (<2 µm) de ces gouges contient une partie de néoformation très importante constituée par, des illites, des chlorites et des interstratifiés de type illite/smectite ou chlorite/smectite. Des datations avec méthode K-Ar sur ces illites ont donné des valeurs comprises entre 12 et 4 Ma qui représentent l'âge de cette dernière déformation cassante. L'application de la méthode de la cristallinité de l'illite (C.I.) a permis d'évaluer les conditions thermiques qui caractérisent le déroulement de cette dernière phase tectonique qui se produit sous conditions de température caractéristiques de l'anchizone et de la diagenèse. L'ensemble des structures de déformation qu'on vient de décrire s'insère parfaitement dans le contexte de convergence oblique entre la plaque adriatique et celle européenne qui à produit l'orogène alpin. On peut considérer les structures tectoniques du Val Vigezzo-Centovalli comme l'expression d'une zone majeure de cisaillement "Simplo-Insubrienne". L'empilement structural et les structures tectoniques affleurantes dans la région sont le résultat de l'interaction entre un régime tectonique transpressif et un régime transtensif. Ces deux champs de tension sont antagonistes entre eux mais sont reliés, de toute façon, à une seule phase décrochante dextre principale, due à une convergence oblique entre deux plaques. À l'échelle de l'évolution géodynamique on peut distinguer différentes étapes au cours desquelles les structures de ces deux régimes tectoniques interagissent en manière différente. En accord avec les données géophysiques et les reconstructions paléodynamiques prises dans la littérature on considère que la ligne Rhône-Simplon-Centovalli représente l'évidence de surface de la suture majeure profonde entre la plaque Adriatique et celle Européenne. Les vitesses de soulèvement qui ont été calculées dans cette étude pour cette région des Alpes donnent une valeur moyenne de 0.8 mm/a qui est tout à fait comparable avec les données proposées par la littérature sur cette zone. La zone Val Vigezzo-Centovalli peut être donc considérée comme un carrefour géologique où se croisent différentes phases tectoniques qui représentent les évidences de surface d'une suture profonde majeure entre deux plaques dans un contexte de collision continentale. ABSTRACT: A wide and complex tectonic zone known as Centovalli line, crosses the Central Alps sector between Domodossola and Locarno. This area, formed by the Vigezzo Valley and Centovalli valley, constitutes the southernmost termination of the Lepontin dome and represents a portion of the alpine nappes root zone. It belongs to a large and complex shear-zone, partly associated with hydrothermal phenomena of alpine age (<20 My), which includes the Insubric Line and the Simplon fault zone. Vigezzo Valley and Centovalli constitute a real crossroads between the mains alpines tectonics lines as well as a zone of juxtaposition of the Southalpine basement with the Austroalpin and Pennique root zone. The deformation phases and the geological structures that can be studied between approximately 35 My and the present. The detailed field study showed the presence of many brittle and ductile deformation structures and fault rocks such as mylonites, cataclasites, pseudotachylites, kakirites, mineralized faults, fault gouges and folds. In the field we could distinguish at least four folds generations related to the various deformation phases. The number and the complexity of these structures indicate a very complicated history, comprising several different stages, that sometimes are related and even superimposed. Part of these deformation structures affect also the sedimentary deposits of quaternary age, in particular the silts and sands lake deposit. These sediments constitute the remainders of a lake basin ascribed to the interglacial Riss/Würm (Eemien, 67.000-120.000 years) and outcroping in the central part of the studied area, in the Eastern part of Santa Maria Maggiore plain. These sediments show a whole series of deformation structures such as inverse fault planes, combined shortening structures and true folds. These faults and folds would represent the surface evidence of a probably active tectonic deformation in quaternary time. Another rock formation attracted all our attention. It is a body of monogenic peridotite breccia which outcrops in discontinuity along the southernmost slope and the bottom of the Vigezzo valley on approximately 20 km. This breccia lies indifferently on the basement (Finero and Orselina units) or on the lake sediments. They are crossed by fault planes which developed slikenside and fault gouges whose orientation is the same of the faults gouges in the alpine basement. This breccia results from the weathering and the surface modelling of an original tectonic breccia which borders the external part of Finero peridotite body. This breccia deformation structures, like those of the lake sediments, were regarded as the surface interaction of active quaternary tectonics in the area. So the last brittle deformation phases which affects this area seems to be actives in quaternary time. Theoverall picture of the studied area on a regional scale enables us to point out a complex shear-zone directed E-W, parallel to the axis of the Centovalli and Vigezzo Valley. The field analysis indicates that this shear-zone began under ductile conditions and evolved in several stages to brittle faulting under surface conditions. The analysis of the geodynamic evolution of the area allows to define three different stages which mark the transition of this alpine basement root zone, from deep P-T conditions to P-T surface conditions. In this context on regional scale three principal deformation phases, which characterize these three stages can be distinguished. The oldest phase consisted of the amphibolitie facies mylonites, associated to dextral strikeslip movements. They are then replaced by green-schists facies mylonites and backfolds related to the backthrusting of the alpines nappes. A second episode is caracterized by the development of an hydrothermal phase bound to an extensive fault and dextral strike-slip fault system, with E-W, NW-SE and SE-NW principal directionsThe principal neoformed mineral phases related to this event are: K-feldspar (microcline), chlorites (Fe+Mg), epidotes prehnite, zéolites (laumontite), sphene and calcite. In this context, to obtain a better characterization of this hydrothermal event, we have used an chlorite geothermometer, sensitive also to the pressure and has the a(H2O), which gave downward values ranging between 450-200°C. The last movements are caracterized by the development of important gouge fault plans, which form a sigmoid structure of kilometric thickness which is recognizable at the valley scale, and is characterized by transpressive movements always with a significant dextral strike-slip component. This deformation phase forms a combined faults system with an average E-W direction, which cuts trough the alpine root zone, the Canavese zone and the Finero ultramafic body. This fault system takes place subparallel to the axis of the valley over several tens of kilometers. A complete and detailed XRD analysis of the gouges fault showed that the clay fraction (<2µm) contains a very significant neo-formation of illite, chlorites and mixed layered clays such as illite/smectite or chlorite/smectite. The K-Ar datings of the illite fraction <2µm gave values ranging between 12 and 4 My and the illite fraction <0.2µm gave more recents values until to 2,4-0 My.This values represent the age of this last brittle deformation. The application of the illite crystallinity method (C.I.) allowed evaluating the thermal conditions which characterize this tectonic phase that occured under temperature conditions of the anchizone and diagenesis. The whole set of deformation structures which we just described, perfectly fit the context of oblique convergence between the Adriatic and the European plate that produced the alpine orogen. We can regard the Vigezzo valley and Centovalli tectonic structures as the expression of a major "Simplo-Insubric" shear-zone. Structural stacking and tectonic structures that outcrop in the studied area, are the result of the interaction between a transpressive and a transtensve tectonic phases. These two tension fields are antagonistic but they are also connected, in any event, with only one principal dextral strike-slip movement, caused by an oblique convergence between two plates. On the geodynamic evolution scale we can distinguish various stages during which these two tectonic structures fields interact in various ways. In agreement with the geophysical data and the paleodynamic recostructions taken in the literature we considers that the Rhone-Simplon-Centovalli line are the surface feature of the major collision between the Adriatique and the European plate at depth. The uplift speeds we calculated in this study for this Alpine area give an average value of 0.8 mm/a, which is in good agreement with the data suggested by the literature on this zone. TheVigezzo Valley and Centovalli zone can therefore be regarded as a geological crossroad where various tectonic phases are superimposed. They represent the evidences of a major and deeper suture between two plates in a continental collision context.

Debris flow modeling for susceptibility mapping at regional to national scale in Norway

Debris flows and related landslide processes occur in many regions all over Norway and pose a significant hazard to inhabited areas. Within the framework of the development of a national debris flows susceptibility map, we are working on a modeling approach suitable for Norway with a nationwide coverage. The discrimination of source areas is based on an index approach, which includes topographic parameters and hydrological settings. For the runout modeling, we use the Flow-R model (IGAR, University of Lausanne), which is based on combined probabilistic and energetic algorithms for the assessment of the spreading of the flow and maximum runout distances. First results for different test areas have shown that runout distances can be modeled reliably. For the selection of source areas, however, additional factors have to be considered, such as the lithological and quaternary geological setting, in order to accommodate the strong variation in debris flow activity in the different geological, geomorphological and climate regions of Norway.

Towards Integrated Governance of Landscape Development : The Swiss Model of Regional Nature Parks

Coherent regulation of landscape as a resource is a major challenge. How can the development interests of some actors (eg cable car operators and property developers) be reconciled with those of others (agriculture, forestry) and with conservation of biodiversity and scenic value? To help understand how the newly introduced Regional Nature Parks (RNPs) can improve the coherence of the regulation regime in Switzerland, we highlight current direct mechanisms for regulation of landscape as a resource (bans, inventories, subsidies) as well as indirect mechanisms (taking place through the regulation of the physical basis of landscapes, eg forest, land, and water planning policies). We show that RNPs are fundamentally innovative because they make it possible to manage and coordinate indirect strategies for appropriate regulation of resources at a landscape scale. In other words, RNPs enable organization of governance of landscape as a resource in a perimeter that is not necessarily restricted to administrative boundaries.