Pendant les transitions glaciaire-interglaciaires,on observe une augmentation en partie abrupte de près de 100 ppm du CO2atmosphérique, indiquant une redistribution majeure entre les réservoirs de carbone des continents, de l'océan et de l'atmosphère.Expliquer les flux de carbone associés à ces transitions est un défi scientifique, qui nécessite une meilleure compréhension du stock de carbone ‘initial’ dans la biosphère terrestre au cours de la période glaciaire. L’objectif de cette thèse est d’améliorer la compréhension du fonctionnement des écosystèmes terrestres et des stocks de carbone au cours du dernier maximum glaciaire (LGM, il y a environ21.000 ans), à travers plusieurs nouveaux développements dans le modèle global de végétation ORCHIDEE-MICT, pour améliorer la représentation de la dynamique de la végétation, la dynamique du carbone dans le sol du pergélisol et les interactions entre les grands herbivores et la végétation dans le modèle de la surface terrestre.Pour la première partie, la représentation de la dynamique de la végétation dans ORCHIDEEMICT pour les régions des moyennes et hautes latitudes, a été calibrée et évaluée avec un ensemble de données spatiales de classes de végétation, production primaire brute, et de biomasse forestière pour la période actuelle.Des améliorations sont obtenues avec la nouvelle version du modèle dans la distribution des groupes fonctionnels de végétation. Ce modèle a ensuite été appliqué pour simuler la distribution de la végétation au cours de laLGM, montrant un accord général avec les reconstructions ponctuelles basées sur des données de pollen et de macro-fossiles de plantes.Une partie du pergélisol (sols gelés en permanence) contient des sédiments épais,riches en glace et en matières organiques appelés Yedoma, qui contiennent de grandes quantités de carbone organique, et sont des reliques des stocks de carbone du Pléistocène.Ces sédiments ont été accumulés sous des climats glaciaires. Afin de simuler l'accumulation du carbone dans les dépôts de Yedoma, j’ai proposé une nouvelle paramétrisation de la sédimentation verticale dans le module de carbone dans le sol de ORCHIDEE-MICT. L'inclusion de ce processus a permis de reproduire la distribution verticale de carbone observée sur des sites de Yedoma. Une première estimation du stock de carbone dans le pergélisol au cours du LGM est obtenue, de l’ordre de ~ 1550 PgC, dont 390 ~446 PgC sous forme de Yedoma encore intacts aujourd’hui (1,3 millions de km2).Potentiellement, une plus grande surface de Yedoma pourrait être présente pendant leLGM, qui a disparue lors de la déglaciation.Pour la troisième partie, à la lumière des impacts écologiques des grands animaux, et le rôle potentiel des méga-herbivores comme une force qui a maintenu les écosystèmes steppiques pendant les périodes glaciaires, j'ai incorporé un modèle de d’herbivores dans ORCHIDEE-MICT, basé sur des équations physiologiques pour l'apport énergétique et les dépenses, le taux de natalité, et le taux de mortalité pour les grands herbivores sauvages.Le modèle a montré des résultats raisonnables de biomasse des grands herbivores en comparaison avec des observations disponibles aujourd’hui sur des réserves naturelles. Nous avons simulé un biome de prairies très étendu pendant le LGM avec une densité importante de grands herbivores. Les effets des grands herbivores sur la végétation et le cycle du carbone du LGM ont été discutés, y compris la réduction de la couverture forestière, et la plus grande productivité des prairies.Enfin, j’ai réalisé une estimation préliminaire du stock total de carbone dans le permafrost pendant le LGM, après avoir tenu compte des effets des grands herbivores et en faisant une extrapolation de l'étendue spatiale des sédiments de type Yedoma basée sur des analogues climatiques et topographiques qui sont similaires à la région de Yedoma actuelle. / During the repeated glacialinterglacialtransitions, there has been aconsistent and partly abrupt increase of nearly100 ppm in atmospheric CO2, indicating majorredistributions among the carbon reservoirs ofland, ocean and atmosphere. A comprehensiveexplanation of the carbon fluxes associatedwith the transitions is still missing, requiring abetter understanding of the potential carbonstock in terrestrial biosphere during the glacialperiod. In this thesis, I aimed to improve theunderstanding of terrestrial carbon stocks andcarbon cycle during the Last Glacial Maximum(LGM, about 21,000 years ago), through aseries of model developments to improve therepresentation of vegetation dynamics,permafrost soil carbon dynamics, andinteractions between large herbivores andvegetation in the ORCHIDEE-MICT landsurface model.For the first part, I improved theparameterization of vegetation dynamics inORCHIDEE-MICT for the northern mid- tohigh-latitude regions, which was evaluatedagainst present-day observation-based datasetsof land cover, gross primary production, andforest biomass. Significant improvements wereshown for the new model version in thedistribution of plant functional types (PFTs),including a more realistic simulation of thenorthern tree limit and of the distribution ofevergreen and deciduous conifers in the borealzone. The revised model was then applied tosimulate vegetation distribution during theLGM, showing a general agreement with thepoint-scale reconstructions based on pollen andplant macrofossil data.Among permafrost (perennially frozen) soils,the thick, ice-rich and organic-rich siltysediments called yedoma deposits hold largequantities of organic carbon, which areremnants of late-Pleistocene carbonaccumulated under glacial climates. In order tosimulate the buildup of the thick frozen carbonin yedoma deposits, I implemented asedimentation parameterization in the soilcarbon module of ORCHIDEE-MICT. Theinclusion of sedimentation allowed the modelto reproduce the vertical distribution of carbonobserved at the yedoma sites, leading toseveral-fold increase in total carbon. Simulatedpermafrost soil carbon stock during the LGMwas ~1550 PgC, among which 390~446 PgCwithin today’s known yedoma region (1.3million km2). This result was still anunderestimation since the potentially largerarea of yedoma during the LGM than todaywas not yet taken into account.For the third part, in light of the growingevidence on the ecological impacts of largeanimals, and the potential role of megaherbivoresas a driving force that maintainedthe steppe ecosystems during the glacialperiods, I incorporated a dynamic grazingmodel in ORCHIDEE-MICT, based onphysiological equations for energy intake andexpenditure, reproduction rate, and mortalityrate for wild large grazers. The model showedreasonable results of today’s grazer biomasscompared to empirical data in protected areas,and was able to produce an extensive biomewith a dominant vegetation of grass and asubstantial distribution of large grazers duringthe LGM. The effects of large grazers onvegetation and carbon cycle were discussed,including reducing tree cover, enhancinggrassland productivity, and increasing theturnover rate of vegetation living biomass.Lastly, I presented a preliminary estimation ofpotential LGM permafrost carbon stock, afteraccounting for the effects of large grazers, aswell as extrapolations for the spatial extent ofyedoma-like thick sediments based on climaticand topographic features that are similar to theknown yedoma region. Since these results werederived under LGM climate and constantsedimentation rate, a more realistic simulationwould need to consider transient climate duringthe last glacial period and sedimentation ratevariations in the next step.

No description available.

No description available.

Au cours des dernières décennies, les observations de glace de glacier enfouie exposée dans les falaises côtières et les glissements de terrain causés par le dégel du pergélisol arctique démontrent que des quantités importantes de glace de glacier ont survécu à la déglaciation et sont toujours préservées dans le pergélisol. Le premier volet de cette étude visait à caractériser des expositions de la glace massive observée à l’île Bylot, au Nunavut, afin de connaitre l’origine de la glace. Puisque la glace de glacier enfouie peut jouer un rôle important dans l’évolution des paysages périglaciaires, cette thèse s’intéresse également au rôle joué par la glace de glacier enfouie dans l'initiation et l’évolution de lacs de thermokarst. Nos résultats démontrent que le pergélisol de l'île Bylot contient des restes de glace de glacier du Pléistocène qui ont survécu aux dernières déglaciations. Dans la vallée Qarlikturvik, des masses de glace intraglaciaire (dérivée du névé) sont associées à un courant de glace de l’inlandsis laurentidien qui recouvrait une partie de la plaine sud de l’île vers la fin du Pléistocène. Cette masse de glace formait une zone de convergence avec les glaciers alpins locaux s'écoulant de la calotte glaciaire centrée sur les monts Byam Martin. Sur l’un des plateaux bordant la vallée Qarlikturvik, une masse de glace de glacier enfouie est associée à la partie basale d'un glacier dont l’âge minimal est estimé à environ 0.77 Ma, mais pourrait être aussi vieux que 2.6 Ma. En raison de sa localisation, à environ 500 m d’altitude, la glace proviendrait vraisemblablement d’une avancée glaciaire régionale et pourrait être associée à la Glaciation de Baffin, soit la plus vieille avancée glaciaire régionale reconnue sur l'île Bylot. De plus, cette glace représente le plus vieux reste de glacier connu en Amérique du Nord et l’une des premières indications de glaciations dans l’est de l’Arctique canadien. La persistance et la fonte tardive de ces épaisses couches de glace de glacier datant du Pléistocène ont eu des effets importants sur le paysage de l'île Bylot, notamment sur les lacs. En effet, nos résultats démontrent que l'initiation des lacs profonds (> 5 m) est liée à la fonte de la glace de glacier enfouie. Ces lacs de thermokarst glaciaire continueront d’évoluer dans un contexte périglaciaire par la fonte de la glace intrasédimentaire (p. ex., glace de ségrégation) et des coins de glace formés ultérieurement dans les sédiments encaissants lors de l’aggradation du pergélisol suivant le retrait glaciaire. Alors qu’une grande partie des paysages arctiques est encore fortement déterminée par leur héritage glaciaire, la fonte de ces masses de glace aura un impact important sur la dynamique des géosystèmes et écosystèmes arctiques. / Over the past decades, observations of buried glacier ice exposed in coastal bluffs and headwalls of retrogressive thaw slumps of the Arctic have indicated that considerable amounts of late Pleistocene glacier ice survived the deglaciation and are still preserved in permafrost. The first phase of this project aimed to characterize two exposures of massive ice observed on Bylot Island (Nunavut) to infer their origins. Since buried glacier ice can play a significant role in reshaping periglacial landscapes, this study also investigates the initiation and development of thermokarst lakes in a tundra valley in response to the melting of buried glacier ice. Our results show that the permafrost of Bylot Island contains remnants of Pleistocene glacier ice that survived the past deglaciations. In the Qarlikturvik valley, bodies of englacial ice (firn-derived) originated from an ice stream flowing from the Laurentide Ice Sheet, which covered part of the southern plain of the island towards the end of the Pleistocene. These glacier ice bodies formed a convergence zone with local alpine glaciers flowing from the ice cap centred over the Byam Martin Mountains. On the edge of a flat plateau bordering the Qarlikturvik Valley, a buried glacier ice body is associated with the basal part of a glacier whose minimum age is estimated at 0.77 Ma, but could be as old as 2.6 Ma. Due to its location on a 500-m a.s.l. plateau, the ice likely originates from a regional glacial advance and could be associated with the Baffin Glaciation, which is the oldest known glaciation on Bylot Island. In addition, this buried glacier ice represents the oldest glacier ice preserved in ice-free Arctic landscapes, and the earliest evidence of a Pleistocene glaciation in the eastern Canadian Arctic Archipelago. The persistence and delayed melting of these thick beds of buried Pleistocene glacier ice had wide-ranging effects on the landscape of Bylot Island. Our results suggest that the initiation of deeper thermokarst lakes (> 5 m) was triggered by the melting of buried glacier ice in our study area, while shallow thermokarst lakes were triggered from the melting of intrasedimental ice and ice wedges. These glacial thermokarst lakes will continue to evolve in a periglacial context through the melting of intrasedimental ice (e.g. segregation ice) and ice wedges subsequently formed in the surrounding sediments during permafrost aggradation following the glacial retreat. As most of the glaciated Arctic landscapes are still strongly determined by their glacial legacy, the melting of these large ice bodies will have significant impacts on Arctic ecosystems and geosystems.

Ce mémoire a pour objectif de mieux comprendre la tomodensitométrie par rayons-x appliquée à la géocryologie. La majorité des études ayant pour but d’analyser le pergélisol en utilisant la tomodensitométrie sont qualitative. De ce fait, il est principalement question de décrire les cryostructures et la disposition géométrique de la glace dans le pergélisol. Il existe peu d’études quantitatives où le volume de glace est estimé avec la tomodensitométrie. De ce fait les méthodologies quantitatives sont peu développées et leur degré de confiance est mal connu. Des échantillons de pergélisol ont été scannés puis analysés avec le logiciel de visualisation Dragonfly©. Une méthode d’analyse quantitative a été développée et combinée à des mesures en laboratoire afin de déterminer l’erreur entre la proportion en glace estimée sur les Ct-Scans et celle mesurée sur l’échantillon physique. Ensuite, une équation pour essayer de corriger cette erreur a été appliquée. Celle-ci a été déterminée à partir de comparaisons des estimations des volumes de glace sur les Ct-Scans et leurs valeurs mesurées en laboratoire. Ces comparaisons ont également permis d’estimer la quantité de matière organique présente dans un échantillon. Les cryostructures du pergélisol ont été modélisées en 3D. Cela a permis de bonifier la classification Nord-Américaine établie par Murton et French en 1994 et de quantifier les volumes respectifs de glace et de sédiments. / This thesis aims to better understand x-ray computed tomography in geosciences and more particularly its use with frozen ground science. Most studies in permafrost science using x-ray computed tomography were qualitative. Therefore, they mainly describe the permafrost cryostructures and the geometric arrangement of the ice in permafrost. There were few quantitative studies where the volumetric ice content of permafrost was estimated with Ct-Scans. As a result, quantitative methodologies are poorly developed, and their degree of confidence is also poorly assessed. Samples from different locations in the Canadian Arctic were scanned and analyzed using Dragonfly©, a visualization software. A quantitative analysis method has been developed and combined with laboratory measurements to calculate the error between the proportion of ice estimated on the Ct-Scans and the one measured on the physical sample. A correction method was developed and applied to the Ct-scans dataset to reduce the error in volumetric ice content estimations. This equation also allows to give an estimation of the organic matter proportion inside a sample. The permafrost cryostructures were modeled in 3D. This improved the visualization of the cryostructures of the North American classification established by Murton and French in 1994, and it allowed to quantify their volumetric ice and sediment contents.

<p>Examinator Professor Peter Jansson, Institutionen för naturgeografi</p> / Research school for teachers on climate and environment grant 2017-06037

Le Haut-Arctique est caractérisé par de nombreux phénomènes climatiques et géomorphologiques extrêmes. Les températures très froides, les faibles précipitations et la couverture végétale disparate permettent une pénétration en profondeur du froid dans le sol. Ceci, combiné à un grand nombre de cycles de gel et de dégel, donne naissance à différentes formes de surface qui affectent les dynamiques d’évolution du pergélisol et induisent une variabilité verticale et latérale de la distribution des teneurs en glace dans le sol. Ce mémoire porte sur l’influence d’un lobe de solifluxion de type rocheux sur les dynamiques d’évolution du pergélisol sur l’île Ward Hunt, Nunavut (Canada). En utilisant une approche cryostratigraphique, les objectifs sont de caractériser la cryostratigraphie d’un lobe de solifluxion, ce qui permettra de proposer un modèle d’évolution du pergélisol en lien avec le déplacement de matériel par solifluxion. L’analyse cryostratigraphique a révélé que la formation du lobe de solifluxion a mené au développement syngénétique d’une couche de pergélisol avec une teneur en glace variant en fonction de la morphologie du lobe, et à l’enfouissement et la préservation d’un corps de glace massive préexistant en bas de pente. La succession verticale et latérale des cryofaciès fait état du déplacement du lobe et de l’impact que ce dernier a sur la variabilité spatiale et temporelle du pergélisol, et ce, tant sur les aspects d’aggradation liés à l’accumulation de matériel en surface que sur les aspects de dégradation associés aux réchauffements climatiques actuels. / The high Arctic is characterized by many extreme climatic and geomorphologic phenomena. Very cold temperatures, low precipitation and sparse vegetation cover, permit a deep penetration of cold in the soil. The latter, combined with a great number of freeze-thaw cycles, give rise to different surface features which in turn affect permafrost evolution dynamics and induce a vertical and lateral variability of the ice content distributions in the soil. This thesis focusses on the influence of a stone-banked solifluction lobe on permafrost evolution dynamics on Ward Hunt Island, Nunvut (Canada). By using a cryostratigraphic approach, the objectives are first, to characterize the cryostratigraphy of a solifluction lobe and second, to propose an evolution model of permafrost in reaction to the downslope displacement of material by solifluction. The cryostratigraphic analysis revealed that the solifluction lobes formation lead to the development of a syngenetic layer of permafrost with an ice content that varied according to the morphology of the lobe, and to the burial and preservation of a pre-existing body of massive ice at the base of the slope. The vertical and lateral sequence of the cryofacies presents the displacement of the lobe and its impact on spatial and temporal variability of the permafrost, that being, the aspects of aggradation related to the accumulation of material at the surface as well as degradation due to the actual climatic warming.

Les régions nordiques à pergélisol seront largement affectées par l'augmentation prévue des températures. Un nombre croissant d’infrastructures qui étaient autrefois construites avec confiance sur des sols gelés en permanence commencent déjà à montrer des signes de détérioration. Les processus engendrés par la dégradation du pergélisol peuvent causer des dommages importants aux infrastructures et entrainer des coûts élevés de réparation. En conséquence, le contexte climatique actuel commande que la planification des projets dans les régions nordiques s’effectue en tenant compte des impacts potentiels de la dégradation du pergélisol. Ce mémoire porte sur l’utilisation de systèmes d’information géographique (SIG) appliqués à l’évaluation du potentiel d’aménagement des territoires situés en milieu de pergélisol. En utilisant une approche SIG, l’objectif est d’élaborer une méthodologie permettant de produire des cartes d'évaluation des risques afin d’aider les collectivités nordiques à mieux planifier leur environnement bâti. Une analyse multi-échelle du paysage est nécessaire et doit inclure l'étude des dépôts de surface, la topographie, ainsi que les conditions du pergélisol, la végétation et les conditions de drainage. La complexité de l'ensemble des interactions qui façonnent le paysage est telle qu'il est pratiquement impossible de rendre compte de chacun d'eux ou de prévoir avec certitude la réponse du système suite à des perturbations. Ce mémoire présente aussi certaines limites liées à l’utilisation des SIG dans ce contexte spécifique et explore une méthode innovatrice permettant de quantifier l'incertitude dans les cartes d'évaluation des risques. / Northern regions underlain by permafrost will largely be affected by the projected increase in air temperature. A growing number of structures that were once built with great confidence on perennially frozen soils are already starting to show signs of deterioration. Processes caused by permafrost degradation can cause significant damages to infrastructure and require high costs of repair. The current climatic context therefore commands that the implementation of projects in permafrost regions follows a well-thought planning in order to account for the potential impacts of permafrost degradation. This thesis focuses on the use of geographic information systems (GIS) applied to the identification of the development potential of communities located in permafrost regions. Using a GIS approach, the goal is to develop a methodology to produce risk-assessment maps to help northern communities better plan their built environment. A multi-scale analysis of the landscape is necessary and should include the investigation of surficial deposits, topography, as well as permafrost, vegetation and drainage conditions. The complexity of all the interactions that shape the landscape is such that it is virtually impossible to account for all of them or to predict with certainty the response of the system following disturbances. This research also presents some of the limitations to the use of GIS in this specific context and explores an innovative method for quantifying uncertainty in risk-assessment maps.

Au Nunavik, l’augmentation des températures atmosphériques en lien avec les changements climatiques influence la capacité portante du sol en modifiant le régime thermique de la couche active et du pergélisol. Selon la nature du sol et des cryostructures, les processus de fluage et tassement différentiel du sol peuvent être considérables. La forte croissance démographique au Nunavik nécessite la construction de plusieurs nouvelles unités d’habitation. La présente étude a produit des outils pour permettre un aménagement durable du territoire par l’interprétation géomorphologique du paysage et l’étude des propriétés géotechniques des cryofaciès qui compose la cryostratigraphie des unités géomorphologiques de Puvirnituq afin d’estimer le comportement du pergélisol au dégel. Trois unités géomorphologiques principales ont été distinguées sur le territoire. Le socle rocheux est stable au dégel ainsi que les dépôts minces qui le recouvrent puisqu’ils sont contenus dans la couche active (unité P1). Le till remanié est dominé par du sable pauvre en glace, lui conférant une stabilité, mais les moraines de De Geer pourraient incorporer des dépôts marins riches en glace à l’intérieur d’un till remanié (unité P2). Les dépressions comblées de dépôts marins et littoraux sont probablement situées sur du till remanié mais contiennent une zone de transition riche en glace entre ~0,9 et 1,7 m de profondeur, excluant les possibilités de construction (unité P3). La principale contribution à l’aménagement du territoire découlant de cette étude a été la production de deux cartes. La première présente la distribution des unités géomorphologiques sur l’ensemble du territoire, permettant d’identifier rapidement le potentiel d’aménagement de l’environnement de Puvirnituq; la seconde carte montre l’emplacement des moraines de De Geer au-delà du village. Ces données trouveront des applications dans l’évaluation de substrats pour l’expansion de l’environnement bâti de Purvinituq et de granulats non-gélifs pour la construction des infrastructures construites sur remblais. Cette carte souligne l’importance de l’histoire Quaternaire sur les conditions actuelles de pergélisol. / In Nunavik, a sustained increase of surface temperatures linked to climate change, affects the bearing capacity of the permafrost by modifying its thermal regime. Depending on the nature of the soil as well as its cryostructures, creep and differential thaw-settlement processes can be significant. The rapid demographic growth in Nunavik requires the construction of several new housing units. In order to estimate the behaviour of permafrost to thaw, this study has contributed to a sustainable land-use planning by the geomorphological interpretation of the landscape and the study of geotechnical properties of the different cryofacies that compose the geomorphological units of Puvirnituq. Three main geomorphological units were distinguished within the territory. The bedrock is thaw-stable as well as the veneer of surficial deposits since they are contained within the active layer (P1 unit). The reworked till is dominated by ice-poor sand, providing a certain stability, but the De Geer moraines could incorporate ice-rich till and ice-rich thaw sensitive marine deposits (P2 unit). The depressions, filled with marine and littoral deposits, probably contain stable reworked till at the bottom but the presence of an ice-rich transition zone between ~0,9 and 1,7 m depth, makes this geomorphological unit improper for construction (P3 unit). The principal contribution to landuse planning of this study is the production of two maps. The first presents the distribution of the geomorphological units over the entire territory, allowing for a quick identification of potential areas to support the development of Puvirnituq; the second map shows the location of De Geer moraines beyond the study area to extrapolate the result of this study to the landscape scale.

Les régions nordiques à pergélisol seront largement affectées par l'augmentation prévue des températures. Un nombre croissant d’infrastructures qui étaient autrefois construites avec confiance sur des sols gelés en permanence commencent déjà à montrer des signes de détérioration. Les processus engendrés par la dégradation du pergélisol peuvent causer des dommages importants aux infrastructures et entrainer des coûts élevés de réparation. En conséquence, le contexte climatique actuel commande que la planification des projets dans les régions nordiques s’effectue en tenant compte des impacts potentiels de la dégradation du pergélisol. Ce mémoire porte sur l’utilisation de systèmes d’information géographique (SIG) appliqués à l’évaluation du potentiel d’aménagement des territoires situés en milieu de pergélisol. En utilisant une approche SIG, l’objectif est d’élaborer une méthodologie permettant de produire des cartes d'évaluation des risques afin d’aider les collectivités nordiques à mieux planifier leur environnement bâti. Une analyse multi-échelle du paysage est nécessaire et doit inclure l'étude des dépôts de surface, la topographie, ainsi que les conditions du pergélisol, la végétation et les conditions de drainage. La complexité de l'ensemble des interactions qui façonnent le paysage est telle qu'il est pratiquement impossible de rendre compte de chacun d'eux ou de prévoir avec certitude la réponse du système suite à des perturbations. Ce mémoire présente aussi certaines limites liées à l’utilisation des SIG dans ce contexte spécifique et explore une méthode innovatrice permettant de quantifier l'incertitude dans les cartes d'évaluation des risques. / Northern regions underlain by permafrost will largely be affected by the projected increase in air temperature. A growing number of structures that were once built with great confidence on perennially frozen soils are already starting to show signs of deterioration. Processes caused by permafrost degradation can cause significant damages to infrastructure and require high costs of repair. The current climatic context therefore commands that the implementation of projects in permafrost regions follows a well-thought planning in order to account for the potential impacts of permafrost degradation. This thesis focuses on the use of geographic information systems (GIS) applied to the identification of the development potential of communities located in permafrost regions. Using a GIS approach, the goal is to develop a methodology to produce risk-assessment maps to help northern communities better plan their built environment. A multi-scale analysis of the landscape is necessary and should include the investigation of surficial deposits, topography, as well as permafrost, vegetation and drainage conditions. The complexity of all the interactions that shape the landscape is such that it is virtually impossible to account for all of them or to predict with certainty the response of the system following disturbances. This research also presents some of the limitations to the use of GIS in this specific context and explores an innovative method for quantifying uncertainty in risk-assessment maps.