Les apports importants d'eau douce en provenance d'un grand bassin versant constituent l'une des caractéristiques principales de la baie d'Hudson, une mer subarctique située au Canada. Dans la baie d'Hudson, le couplage entre l'eau salée et l'eau douce exerce une grande influence sur l'ensemble des paramètres physiques, biogéochimiques et biologiques. Au cours des dernières décennies, la baie d'Hudson a subi d'importantes modifications associées aux changements climatiques mais aussi à l'installation de grandes infrastructures hydroélectriques qui régulent les apports d'eau douce par les rivières. Alors que ces changements sont susceptibles d'avoir un impact sur les écosystèmes dans la baie, peu d'attention a été portée sur la composition et la distribution des communautés planctoniques qui sont pourtant un maillon essentiel des réseaux trophiques marins. Au sein du projet BaySys qui vise à fournir une base scientifique afin de séparer les impacts liés aux changements climatiques de ceux de la régulation des cours d'eau dans la baie d'Hudson, cette thèse de doctorat a pour objectif d'établir un état des lieux de la diversité et de la distribution spatiale des communautés microbiennes présentes dans la baie d'Hudson. Pour répondre à cet objectif, nous avons utilisé une approche de séquençage haut débit des gènes de l'ARNr et de l'ARNr 18S et 16S afin de déterminer l'influence des paramètres environnementaux mesurés sur le terrain sur les communautés de protistes et de procaryotes. Les résultats ont montré que le débit des rivières modifie la composition des communautés de protistes au niveau des estuaires dans la rivière Nelson, Churchill et Grande Rivière de la Baleine. Dans la rivière Nelson, la forte circulation estuarienne permet la remontée d'eau riches en nutriments qui favorise la présence de diatomées du genre Rhizosolenia. Cependant, à l'échelle de la baie, les apports d'eau douce au printemps et a l'été au niveau des cotes entraine une forte stratification qui limite la croissance du phytoplancton et favorise les taxa hétérotrophes, notamment les Dinoflagellés. Dans ces eaux de surface fortement stratifiées, nous avons pu mettre en évidence la présence de 7 espèces phytoplanctoniques potentiellement toxiques. Le changement de régime de la glace de mer au printemps entraine une succession d'organismes eucaryotes et procaryotes associés à la dégradation des efflorescences de phytoplancton. Ces changements de communautés en surface modifient l'origine des apports de matière organique dans les écosystèmes profonds de la baie. En générant un nouvel état des lieux de la composition des communautés de protistes et de procaryotes à différentes échelles spatiales, les résultats de cette étude permettent de mieux comprendre l'influence des apports d'eau douce sur les communautés microbiennes dans la baie d'Hudson. Cette thèse constitue ainsi un nouvel outil pour mieux appréhender l'influence des perturbations d'origine anthropique sur les réseaux microbiens marins dans la baie d'Hudson. / The Hudson Bay is a sub-Arctic sea in Canada that receives massive freshwater inflows from a large watershed. In Hudson Bay, the freshwater-marine coupling has a strong influence on physical, biogeochemical and biological processes. In recent decades, Hudson Bay has undergone significant changes associated with climate change and the installation of large hydroelectric infrastructures that regulate freshwater inflows from rivers. While these changes are likely to have an impact on ecosystems in the bay, little attention has been paid to the composition and distribution of planktonic communities, which are an essential link in marine food webs. As part of the BaySys project, which aims to provide a scientific basis to separate climate change effects from those of regulation of freshwater on physical, biological and biogeochemical conditions in Hudson Bay, this PhD thesis aims to establish a baseline of the diversity and spatial distribution of microbial communities in Hudson Bay. To address this objective, we used a high-throughput sequencing approach for 18S rRNA and 16S rRNA genes to determine the potential influence of environmental parameters measured in the field on protist and prokaryote communities. The results showed that river flow altered the composition of protist communities in the estuaries of the Nelson, Churchill and Great Whale rivers. In the Nelson River estuary, strong estuarine circulation allows the upwelling of nutrient-rich water which favours the presence of diatoms of the genus Rhizosolenia. At the bay scale in spring and summer, the freshwater inflow leads to a strong stratification which limits the growth of phytoplankton species and favours heterotrophic taxa, particularly dinoflagellates. In these highly stratified surface waters, we were able to identify the presence of 7 potentially toxic phytoplankton species. The changing of the sea ice regime in spring resulted in a succession of eukaryotic and prokaryotic organisms associated with the degradation of phytoplankton blooms. These changes modified the diatom sinking particles reaching the bottom of the Hudson Bay. By generating a new picture of the composition of protists and prokaryotes at different spatial scales, the results of this study provide a better understanding of the influence of freshwater inputs on microbial communities in Hudson Bay. This study thus constitutes a new tool to better understand the influence of anthropogenic disturbances on marine microbial networks in Hudson Bay.

Dans le cadre de la présente étude, un modèle géothermique tridimensionnel du pergélisol a été développé pour la toute première fois à l’échelle locale d’un village nordique afin d’évaluer les impacts des changements climatiques sur le pergélisol. Ce village est celui de Salluit au Québec nordique dans la zone de pergélisol continu. La majeure partie des infrastructures de ce village sont fondées sur un pergélisol salin riche en glace vulnérable aux changements climatiques. Le code numérique d’HydroGeoSphere a été modifié pour simuler la transmission de la chaleur par conduction avec changement de phase et tenir compte de la variation des propriétés thermiques du sol avec la température. Différents scénarios de changements climatiques lors du prochain siècle ont été utilisés pour alimenter le modèle géothermique. La variabilité spatiale du régime thermique du pergélisol simulée à l’échelle d’un village nordique est majeure lorsque l’hétérogénéité des dépôts, leur contenu en glace et les différentes conditions de surface observées sur le terrain sont considérés dans une modélisation numérique tridimensionnelle. L’influence des flux de chaleur latéraux dans cette modélisation produit des résultats qui sont dans certains cas significativement différents de ceux simulés au moyen d’un modèle unidimensionnel, particulièrement pour les sols sous d’importantes congères. Selon les scénarios de changements climatiques jugés pessimistes pour le maintien du pergélisol, un talik permanent se formera à partir de la seconde moitié du 21ième siècle sous les congères, en bordure d’un remblai routier où la neige et l’eau de fonte s’accumule à la surface d’un sol mal drainé et dans l’argile saline du secteur nord-est du village de Salluit. L’épaisseur maximale du talik à la fin du 21ième siècle varie entre 1,3 et 19,3 m respectivement dans l’argile et le roc sous les congères, alors qu’elle est de 2,8 m dans l’argile saline du secteur nord-est. Seuls les tills et le roc recouverts par d’importantes congères seront affectés par la présence d’un talik à la toute fin du 21ième siècle advenant un réchauffement climatique optimiste. Une dégradation plus importante du pergélisol sous les remblais routiers que pour un même sol en terrain naturel sous un couvert nival mince est à prévoir selon les scénarios climatiques pessimistes. Les infrastructures contruites sur des remblais seront potentiellement affectées par des tassements différentiels.

La modélisation de l’impact des différentes méthodes de réhabilitation, en fonction du temps d’initiation des dégradations, est le défi majeur dans le domaine de la durabilité des structures. La durabilité d’une structure est liée directement aux conditions d’exposition influencées par le climat et l’environnement où l’ouvrage est localisé. Selon le degré d’exposition, la vitesse de dégradation des structures peut alors s’accélérer ou être lente. Ce travail développe une étude approfondie sur l’exposition des structures en condition hivernales, qui comprend l’utilisation des sels de déglaçage. Une enquête a été réalisée pour connaître la quantité de sel utilisée au Québec et au Canada. En complément, les trois types de sel de déglaçage les plus utilisés ont été évalués afin de déterminer sa vitesse de progression dans un béton sous conditions hivernales sévères. Une station météorologique mobile MexStUL, munie des plusieurs capteurs de mesures temporelles, a été développé pour monitorer des conditions d’exposition dans deux sites expérimentaux soumis au climat hivernal rigoureux des années 2018/19. Il a été alors possible de prédire les facteurs climatiques (température, humidité, vitesse et direction des vents, ensoleillement et précipitations) et environnementaux (trafic, condition de route, concentration de sel sur la chaussée, aux éclaboussures et au brouillard salin) relatifs aux conditions d’exposition des structures. Des modèles ont été proposés afin de prédire l’exposition à travers du biais d’une station météorologique mobile ainsi que pouvoir évaluer l’influence spatiale de la concentration de sel de déglaçage. Un nouveau capteur a été développé, mesurant la concentration de brouillard salin sous condition hivernale sévère. Les résultats démontrent une importante variabilité spatiale de concentration pour l’exposition stagnante. Un modèle d’influence spatiale des éclaboussures a été proposé afin de délimiter ce type d’exposition. La concentration du brouillard salin pendant l’hiver montre une évolution importante et un modèle de concentration a été proposé et le volume de brouillard salin pendant l’hiver a pu être estimé. / Modeling the impact of the various rehabilitation methods, as a function of the initiation time of the damage, is the major challenge in the field of structural durability. The durability of a structure is directly linked to the exposure conditions influenced by the climate and the environment where the structure is located. Depending on the degree of exposure, the rate of degradation of structures can then accelerate or be slow. This work develops an in-depth study on the exposure of structures in winter conditions, which includes the use of de-icing salts. A survey was carried out to find out the amount of salt used in Quebec and Canada. In addition, the three most used types of de-icing salt were evaluated to determine its speed of progression in concrete under severe winter conditions. A mobile weather station MexStUL, designed with several sensors of temporal measures, has been developed to monitor exposure conditions in two different experimental sites on severe winter conditions during 2018/19. It was then possible to predict the climatic (temperature, humidity, speed and wind direction, sunshine and precipitation) and environmental factors (traffic, road condition, salt concentration on road, splashes and salt laden mist) relating to the conditions of exposure of the structures. Models were proposed to predict exposure condition using a meteorological station and evaluation of spatial influence of de-icing salt on different types of exposure as well. A new sensor has been developed, measuring the concentration of salt laden mist under severe winter conditions. Results demonstrate an important spatial variability on stagnant water. A model about spatial influence of brine splashes is proposed to delimit this type of exposure. Airborne concentration during winter climate shows an important evolution and a concentration model is proposed and the volume of salt laden mist is also determined during winter.

L’objectif de cette étude était d’évaluer les effets du réchauffement climatique (augmentation de température et des nutriments dans le sol) combinés au broutement simulé sur la disponibilité et la qualité du bouleau glanduleux (Betula glandulosa), une des ressources alimentaires estivales du caribou migrateur (Rangifer tarandus). Le broutement modéré et élevé a diminué la biomasse des feuilles de bouleau. L’augmentation de température a diminué la biomasse et devancé la phénologie, ce qui pourrait éventuellement affecter la synchronisation entre les besoins nutritifs du caribou et la disponibilité du bouleau. L’augmentation en nutriments dans le sol tendait à augmenter la biomasse et la qualité des feuilles, alors que le broutement élevé combiné à l’augmentation des nutriments dans le sol a diminué la concentration en fibres, améliorant la qualité des feuilles de bouleau. La dominance des effets du broutement pourrait n’être que temporaire, puisque les effets du réchauffement climatique sont directionnels et qu’ils pourraient graduellement surpasser les effets du broutement sur les arbustes. / The main objective of this study was to assess the effects of climate warming (increased temperature and soil nutrient availability) and simulated browsing on the availability and quality of dwarf birch (Betula glandulosa), one of the summer resources of migratory caribou (Rangifer tarandus). Moderate and heavy browsing reduced the biomass of dwarf birch leaves. Warming reduced the biomass and advanced the phenology of dwarf birch, potentially affecting the synchrony between caribou nutritive requirements and the availability of dwarf birch. An increase in soil nutrients availability tended to increase the biomass and quality of dwarf birch leaves, while heavy browsing combined to indirect effects of climate warming reduced fiber concentration, improving the quality of dwarf birch leaves. Dominance of browsing effects could be temporary, since climate warming effects are directional and could gradually surpass browsing effects on shrubs.

La performance à long terme du réseau routier au Québec est fortement influencée par le climat et les conditions météorologiques (Doré et Zubeck, 2008). Parmi ces facteurs, des niveaux élevés de saturation dans les sols et les matériaux de chaussée peuvent être une cause importante de la détérioration des routes. Selon des scénarios de changements climatiques établis par Ouranos, le sud du Québec subira des augmentations de précipitations moyennes mensuelles de -0,1 % à 8,45 % sur un horizon de 2010 à 2039. Le but de ce projet est de quantifier l’effet de ces futures augmentations sur le comportement mécanique des structures de chaussées et des sols. Basé sur les données collectées dans des sections de routes instrumentées, une relation entre les précipitations et l’augmentation du niveau de saturation des couches de chaussée est proposée. Afin de déterminer la relation entre les propriétés mécaniques et la teneur en eau, les comportements en module réversible et en déformation permanente (quatre différents sols) ont été déterminés avec des essais triaxiaux et ont été validés avec un simulateur de charge routier. Une analyse d’endommagement est effectuée afin de quantifier la diminution de la durée de vie utile causée par les changements climatiques. Il a été montré que les augmentations de précipitations, prévues au Québec, auront un impact significatif sur la performance des chaussées. Ainsi, certaines techniques de mitigation seront proposées pour pallier la diminution de la vie utile des chaussées. / The long-term performance of the road network of the province of Quebec (Canada) is strongly influenced by the climate and the weather conditions. Amongst other factors, high levels of saturation in soils and pavement materials are believed to be an important cause of pavement deterioration. According to the climate change scenarios established by Ouranos (2010), the North and South of Quebec will undergo an average precipitation increase from -0.1% to 8.45% in a future horizon of 2010-2039. The purpose of this project is to quantify the effect of these expected precipitation increases on the mechanical behavior of road structures, materials and soils. Based on literature and on data collected on instrumented road sections, a relationship between precipitation increase and saturation level of pavement layers is proposed. The resilient modulus and permanent deformation behavior for various water contents and four different subgrade soils was determined using triaxial tests, which were validated using small-scale heavy vehicle simulator, in order to determine the existing relationship between mechanical properties and moisture contents. Using the precipitation increase scenarios and the preset models, a damage analysis is performed to quantify the decrease of pavements service life caused by climate change. It is found that climate change, and more precisely the increase of precipitation expected in the Province of Quebec, will have a significant impact on pavement performance and that adapted pavement structures and materials, such as improved drainage, increased structural capacity or materials with reduced sensitivity to water, are possible options to reduce the loss of pavement life associated with climate change.

Le métabolisme microbien du dimethylsulfoniopropionate (DMSP), précurseur du gaz climatiquement actif diméthylsulfure (DMS), a été caractérisé dans l‟Arctique canadien en septembre 2008 à l‟aide du traceur 35S-DMSP et de la méthode MAR-CARD-FISH. Le DMSP contribuait en moyenne à 7 et 51% des besoins bactériens en carbone et en soufre, respectivement. Entre 40 et 65% de la communauté bactérienne assimilait le DMSP, avec une forte contribution des groupes Gammaproteobacteria et Alphaproteobacteria (autre que Roseobacter). Les différents assemblages bactériens étaient associés à deux masses d‟eau de surface : celle de l‟Arctique (AS) et celle de la Baie de Baffin (BBS), plus chaude et plus salée. Dans BBS, qui contenait plus de bactéries, la prise de DMSP était plus élevée que dans AS. Ces résultats montrent que les altérations des masses d‟eau anticipées par les modèles climatiques actuels pourraient se répercuter sur la dynamique du DMS(P) et le flux résultant de DMS vers l‟atmosphère. / The bacterial metabolism of dimethylsulfoniopropionate (DMSP), the precursor of the climate-active gas dimethylsulfide (DMS), was characterized in Arctic waters in September 2008, using the radioisotopic tracer 35S-DMSP and the MAR-CARD-FISH method. DMSP contributed on average 7 and 51% of the bacterial requirements in carbon and sulfur, respectively. Between 40 and 65% of the bacterial community assimilated DMSP, with a stronger contribution from the groups Gammaproteobacteria and Alphaproteobacteria (other than Roseobacter). The different bacterial assemblages were associated with two surface water masses: the Arctic Surface (AS) water, and the warmer and saltier Baffin Bay Surface (BBS) water. In BBS, which showed a higher bacterial abundance, the DMSP uptake was higher than in AS. Those results show that the alteration of water mass characteristics anticipated by current climate models could affect the DMS(P) dynamics and the resulting DMS flux to the atmosphere.

L’analyse de photographies de la région de la rivière Boniface (Nunavik) a montré une forte expansion du couvert arbustif entre 1957 et 2008. L’objectif principal de ce projet était de déterminer si cette expansion est observable dans les assemblages polliniques récents des sédiments lacustres de cette région. Pour ce faire, une analyse pollinique à fine résolution temporelle a été effectuée sur les sédiments de deux lacs. Les pourcentages polliniques du bouleau glanduleux, la principale espèce responsable de la densification du couvert arbustif, sont plus élevés dans les sédiments récents de ces lacs. De plus, les influx polliniques calculés à l’aide de la datation au 210Pb montrent une augmentation des genres Betula et Alnus depuis environ 1975. Ces influx polliniques doivent toutefois être interprétés avec prudence en raison de difficultés liées à l’élaboration de la chronologie des sédiments. L’analyse pollinique à fine résolution semble néanmoins permettre de détecter l’expansion récente des arbustes.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2008-2009 / Les flux de carbone organIque particulaire (POC) sont des déterminants clés de la régulation du climat sur Terre. Les rapides changements climatiques qui affectent les hautes latitudes nordiques sont susceptibles de modifier abruptement la dynamique des flux de POC dans l'océan Arctique. Les implications de telles modifications pour le réseau alimentaire marin et le climat sont inconnues. Afin d'examiner le contrôle physique et biologique des flux verticaux de POC dans l'océan Arctique, des pièges séquentiels à particules ont été mouillés dans la région du plateau canadien de la mer de Beaufort. Une . première étude sur le talus continental a révélé le fonctionnement horizontal partiel des flux de POC dans les mers englacées, puisqu'environ la moitié du flux vertical récolté à cet endroit a été relié à l'advection latérale de particules par la convection thermohaliQe. Ainsi, advenant un accroissement de la quantité de matériel organique ·produite dans un océan Arctique à couvert de glace estival absent, la séquestration de POC serait augmentée à cause de la formation d'un plus grand volume de glace saisonnière. L'investigation subséquente du cycle annuel d'export de POC sur le plateau a permis d'élaborer un modèle conceptuel de cet écosystème qui varie selon l'importance relative des apports de POC labile ou réfractaire. Ce modèle suggère que si les changements climatiques induisent une augmentation des apports terrigènes et de détritus sur les plateaux arctiques, le développement d'un écosystème basé sur le recyclage du POC serait favorisé. Finalement, l'examen de la variabilité interannuelle des flux de POC en marge du plateau ~ révélé que les flux verticaux les plus élevés étaient détectés en concomitance avec un retrait tardif des glaces, une forte floraison pulsée et un mismatch apparent avec les consommateurs pélagiques. Ce résultat supporte qu'une débâcle précoce et une plus longue saison de production primaire dans les mers arctiques profiterait davantage au transfert trophique de POC qu'à son export vertical. La conclusion générale de la thèse illustre la complexité de prédire la nature et la magnitude des flux verticaux de POC dans un océan Arctique affecté par de multiples changements environnementaux. Le maintien à long terme d'observatoires marins (ex. lignes de mouillage) est donc essentiel à une meilleure synthèse des liens qui unissent le climat et le cycle du carbone organique dans les hautes latitudes nordiques.

Cette étude porte sur la mesure des conséquences du réchauffement climatique sur la croissance et le développement des graines de l'épinette noire en zone boréale. Comme la densité maximale atteinte par les forêts d'épinette noire et la limite nordique et altitudinale de l'espèce sont sous contrôle thermique, dans le contexte des changements climatiques attendus, il nous est apparu important de préciser le régime des températures qui affecte la répartition de l'épinette. Cette étude porte plus précisément sur l'interface entre la toundra forestière et la limite septentrionale de la forêt boréale fermée dominée par l'épinette noire. L'objectif de cette thèse est de déterminer comment la température influence la croissance des arbres et la maturation des semences de l'épinette noire à la limite nord de la forêt boréale fermée. La thèse est divisée en trois parties. La première porte sur les caractéristiques thermiques de la limite nord de la forêt boréale fermée au cours de la période de 1970 à 2002. La deuxième porte sur la relation entre la température et la production de graines par l'épinette noire. Enfin, la troisième partie porte sur les effets climatiques de la mise en eau d'un vaste réservoir hydroélectrique sur la croissance et la production de graines par l'épinette noire des milieux néoinsulaires. Entre 1970 et 2002, une augmentation significative de la température moyenne annuelle, de la somme thermique et de la durée de la période sans gel est survenue particulièrement dans la zone de la forêt boréale ouverte et de la toundra forestière. La limite septentrionale de la forêt boréale fermée coïncide avec la position du seuil thermique de 800 à 950 degrés-jours, un intervalle qui correspond à la chaleur nécessaire pour la maturation des graines de l'épinette noire. De plus le seuil de 800 degrés-jours s'est déplacé d'un degré de latitude vers le nord entre 1970 et 2002, passant de 53° 12' N à 54° 16' N. Compte tenu du réchauffement climatique attendu, le seuil de maturation des graines d'épinette noire pourrait se déplacer vers le nord de façon importante et permettre ainsi la maturation des graines d'épinette dans une zone où les conditions thermiques nécessaires ne sont pas actuellement rencontrées. Si le processus physiologique de la maturation des graines est confirmé comme le facteur déterminant de la position de la limite nord de la forêt boréale, les conséquences pourraient être un déplacement vers le nord de la limite septentrionale de la forêt boréale fermée. Sur l'ensemble de son aire de répartition, Pépinette noire requiert une moyenne de 750 degrés-jours pour produire des graines quasi-matures (stade 22 sur une totalité de 23). Une variabilité existe selon les populations et leur position géographique. Néanmoins, le pourcentage de germination peut être prédit par une fonction sigmoïde des degrés-jours cumulés le long d'un gradient -latitudinal des températures. Les conditions locales peuvent toutefois produire des différences. Ainsi, les enclaves plus froides que constituent les réservoirs hydroélectriques ont une production de graines équivalente à celles des sites de un à trois degrés de latitude plus au nord. Enfin, les conditions climatiques des sites néoinsulaires du réservoir Rober Bourassa diffèrent des autres sites de la région, principalement au printemps. En effet, c'est pendant cette période que le déficit de chaleur est le plus marqué. Les conséquences sur la croissance radiale des arbres matures ne sont pas notables à l'échelle des peuplements. Les arbres répondent plus à un changement régional du climat qu'à l'effet local du réservoir. Les effets climatiques du plan d'eau (une augmentation du bois de compression dû à une déstabilisation par le vent et une diminution de la fréquence des cernes de gel) se limitent aux sites les plus ouverts et les plus exposés. Enfin, la croissance apicale de la strate de régénération montre que le réservoir peut avoir des conséquences positives sur la croissance des plantules dans les secteurs protégés par la végétation (arbres ou arbustes) et à la périphérie des îles.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2016-2017 / On observe présentement une augmentation du couvert arbustif à l'échelle circumpolaire, mais les déterminants de cette arbustation à l'échelle du paysage sont encore mal compris. Dans cette étude, nous avons adopté une approche de modélisation des changements de couvert afin d'identifier les variables déterminant l'augmentation du couvert arbustif à proximité d'Umiujaq, Québec. Nous avons modélisé les changements survenus entre 1990/1994 et 2010 et avons utilisé ces modèles afin de générer des prédictions sur l'augmentation du couvert arbustif. Nous avons trouvé que les arbustes tendent à coloniser plus facilement les sites à basse altitude ainsi que les pentes modérées, alors que leur couvert est plus susceptible d'augmenter à proximité de peuplements arbustifs déjà existants. Les probabilités d'arbustation prédites dans la région étaient cohérentes avec les patrons de changement inférés par les observations sur le terrain, mais pas par les changements récents de NDVI (normalized difference vegetation index, ou indice de végétation par différence normalisée) dans la région. Nos résultats indiquent que davantage de recherche est requise sur les impacts à grande échelle de l'arbustation sur les écosystèmes nordiques et sur le lien entre les changements de couvert végétal et les variations dans les indices dérivés de données de télédétection. / A circumpolar increase in shrub cover has been underway in northern ecosystems in recent decades, but the landscape-scale drivers of this shrubification are not well understood. In this study, we adopted a land cover change modelling approach to identify variables driving shrub cover increase near Umiujaq, Québec. We modelled land cover change in this area over the period 1990/1994-2010 and used these models to generate spatially explicit predictions of shrub cover increase. We found that shrubs tended to colonize low-elevation areas and moderate slopes more easily, whereas their cover was more likely to increase in the vicinity of existing shrub stands. The predicted probabilities of shrub cover increase in the region were consistent with patterns of change inferred from field observations, but not from recent local changes in NDVI (normalized difference vegetation index). Our findings warrant further research on the large-scale impacts of shrubification on ecosystem function and the link between land cover changes and variations in indices derived from remote sensing.