A multi-temporal approach to model endangered species distribution in Europe: the case of the Eurasian otter in Italy

We analysed the relationship between changes in land cover patterns and the Eurasian otter occurrence over the course of about 20 years (1985-2006) using multi-temporal Species Distribution Models (SDMs). The study area includes five river catchments covering most of the otter's Italian range. Land cover and topographic data were used as proxies of the ecological requirements of the otter within a 300-m buffer around river courses. We used species presence, pseudo-absence data, and environmental predictors to build past (1985) and current (2006) SDMs by applying an ensemble procedure through the BIOMOD modelling package. The performance of each model was evaluated by measuring the area under the curve (AUC) of the receiver-operating characteristic (ROC). Multi-temporal analyses of species distribution and land cover maps were performed by comparing the maps produced for 1985 and 2006. The ensemble procedure provided a good overall modelling accuracy, revealing that elevation and slope affected the otter's distribution in the past; in contrast, land cover predictors, such as cultivations and forests, were more important in the present period. During the transition period, 20.5% of the area became suitable, with 76% of the new otter presence data being located in these newly available areas. The multi-temporal analysis suggested that the quality of otter habitat improved in the last 20 years owing to the expansion of forests and to the reduction of cultivated fields in riparian belts. The evidence presented here stresses the great potential of riverine habitat restoration and environmental management for the future expansion of the otter in Italy

Local and Regional Scale Determinants of Biodiversity Patterns in Boreal Agricultural Landscapes

This thesis examines the local and regional scale determinants of biodiversity patterns using existing species and environmental data. The research focuses on agricultural environments that have experienced rapid declines of biodiversity during past decades. Existing digital databases provide vast opportunities for habitat mapping, predictive mapping of species occurrences and richness and understanding the speciesenvironment relationships. The applicability of these databases depends on the required accuracy and quality of the data needed to answer the landscape ecological and biogeographical questions in hand. Patterns of biodiversity arise from confounded effects of different factors, such as climate, land cover and geographical location. Complementary statistical approaches that can show the relative effects of different factors are needed in biodiversity analyses in addition to classical multivariate models. Better understanding of the key factors underlying the variation in diversity requires the analyses of multiple taxonomic groups from different perspectives, such as richness, occurrence, threat status and population trends. The geographical coincidence of species richness of different taxonomic groups can be rather limited. This implies that multiple geographical regions should be taken into account in order to preserve various groups of species. Boreal agricultural biodiversity and in particular, distribution and richness of threatened species is strongly associated with various grasslands. Further, heterogeneous agricultural landscapes characterized by moderate field size, forest patches and non-crop agricultural habitats enhance the biodiversity of rural environments. From the landscape ecological perspective, the major threats to Finnish agricultural biodiversity are the decline of connected grassland habitat networks, and general homogenization of landscape structure resulting from both intensification and marginalization of agriculture. The maintenance of key habitats, such as meadows and pastures is an essential task in conservation of agricultural biodiversity. Furthermore, a larger landscape context should be incorporated in conservation planning and decision making processes in order to respond to the needs of different species and to maintain heterogeneous rural landscapes and viable agricultural diversity in the future.

Facing the broader dimensions of biological invasions

Invasive species are an excellent opportunity to think about the nature society desires, particularly in the face of global changes. Nature and human views of nature are rapidly evolving; our approach to bio- logical invasions through biosecurity institutions and land management policies must evolve in tandem with these changes. We review three dimensions that are insufficiently addressed. First, biological inva- sions are culturally shaped and interpreted. Humans play a major role in the movement and nurturing of alien life, and esthetics, perception, and emotion are deeply implicated in the management of invasive species. What people fear or regret with invasive species are not their effects on nature per se, but their effects on a particular desired nature, and policymaking must reflect this. Second, biological invasions are not restricted to negative impacts. Invasions take place in landscapes where many natural condi- tions have been altered, so policy tools must recognize that invasive species are a functional, structural, and compositional part of transformed ecosystems. In some cases, native species benefit from changes in resource availability caused by invasions or from protections provided by an invasive plant. Finally, invasive species can help ecosystems and people to adapt to global change by maintaining ecosystem processes such as productivity, carbon storage, and nutrient cycling in a context of climate change or land cover transformations. While recognition is growing among ecologists that novel, invaded ecosystems have value, and while the on-the-ground application of biosecurity policies has of necessity adjusted to local contexts and other agendas, invasion biology could aid policymaking by better addressing the three complexities inherent in the three dimensions highlighted above.

Water consumption and soil moisture distribution in melon crop with mulching and in a protected environment

Mulching has become an important technique for land cover, but there are some technical procedures which should be adjusted for these new modified conditions to establish optimum total water depth. It is also important to observe the soil-water relations as soil water distribution and wetted volume dimensions. The objective of the present study was to estimate melon evapotranspiration under mulching in a protected environment and to verify the water spatial distribution around the melon root system in two soil classes. Mulching provided 27 mm water saving by reducing water evaporation. In terms of volume each plant received, on average, the amount of 175.2 L of water in 84 days of cultivation without mulching, while when was used mulching the water requirement was 160.2 L per plant. The use of mulching reduced the soil moisture variability throughout the crop cycle and allowed a greater distribution of soil water that was more intense in the clay soil. The clayey soil provided on average 43 mm more water depth retention in 0.50 m soil deep relative to the sandy loam soil, and reduced 5.6 mm the crop cycle soil moisture variation compared to sandy loam soil.

Predicting the future effectiveness of protected areas for bird conservation in Mediterranean ecosystems under climate change and novel fire regime scenarios

AimGlobal environmental changes challenge traditional conservation approaches based on the selection of static protected areas due to their limited ability to deal with the dynamic nature of driving forces relevant to biodiversity. The Natura 2000 network (N2000) constitutes a major milestone in biodiversity conservation in Europe, but the degree to which this static network will be able to reach its long-term conservation objectives raises concern. We assessed the changes in the effectiveness of N2000 in a Mediterranean ecosystem between 2000 and 2050 under different combinations of climate and land cover change scenarios. LocationCatalonia, Spain. MethodsPotential distribution changes of several terrestrial bird species of conservation interest included in the European Union's Birds Directive were predicted within an ensemble-forecasting framework that hierarchically integrated climate change and land cover change scenarios. Land cover changes were simulated using a spatially explicit fire-succession model that integrates fire management strategies and vegetation encroachment after the abandonment of cultivated areas as the main drivers of landscape dynamics in Mediterranean ecosystems. ResultsOur results suggest that the amount of suitable habitats for the target species will strongly decrease both inside and outside N2000. However, the effectiveness of N2000 is expected to increase in the next decades because the amount of suitable habitats is predicted to decrease less inside than outside this network. Main conclusionsSuch predictions shed light on the key role that the current N2000may play in the near future and emphasize the need for an integrative conservation perspective wherein agricultural, forest and fire management policies should be considered to effectively preserve key habitats for threatened birds in fire-prone, highly dynamic Mediterranean ecosystems. Results also show the importance of considering landscape dynamics and the synergies between different driving forces when assessing the long-term effectiveness of protected areas for biodiversity conservation.

Flash flood evolution in North-Western Mediterranean

The present paper shows an in-depth analysis of the evolution of floods and precipitation in Catalonia for the period 1981-2010. In order to have homogeneous information, and having in mind that not gauge data was available for all the events, neither for all the rivers and stream flows, daily press from a specific newspaper has been systematically analysed for this period. Furthermore a comparison with a longer period starting in 1900 has been done. 219 flood events (mainly flash flood events) have been identified for the period of 30 years (375 starting in 1900), 79 of them were ordinary, 117 of them were extraordinary and 23 of them were catastrophic, being autumn and summer the seasons with the maxima values. 19% of the events caused a total of 110 casualties. 60% of them died when they tried to cross the street or the stream. Factors like the evolution of precipitation, population density and other socio-economical aspects have been considered. The trend analysis shows an increase of 1 flood/decade that probably has been mainly due to inter-annual and intra-annual changes in population density and in land-use and land-cover.

Distribution Patterns of Soil Properties in a Rural Mediterranean Area in Northeastern Spain

Soil properties on the Cap de Creus Peninsula, NE Spain depend primarily on scarce agricultural practices and early abandonment. In the study area, 90% of which is mainly covered by Cistus shrubs, 8 environments representing variations in land use/land cover and soil properties at different depths were identified. In each environment variously vegetated areas were selected and sampled. The soils, collected at different depths, were classified as Lithic Xerorthents according to the United States Department of Agriculture system of soil classification (USDA-NRCS 1975). Differences in soil properties were largely found according to the evolution of the plant canopy and the land use history. To identify underlying patterns in soil properties related to environmental evolution, factor analysis was performed and factor scores were used to determine how the factor patterns varied between soil variables, soil depths and selected environments. The three-factor model always accounted for 80% of the total variation in the data at the different soil depths. Organic matter was the more relevant soil property at 0–2 cm depth, whereas active minerals (silt and clay) were found to be the most relevant soil parameters controlling soil dynamics at the other depths investigated. Results showed that vineyards and olive tree soils are poorly developed and present worse conditions for mineral and organic compounds. Analysis of factor scores allowed independent assessment of soils, depth and plant cover and demonstrated that soils present the best physico-chemical characteristics under Erica arborea and meadows. In contrast, soils under Cistus monspeliensis were less nutrient rich and less well structured

Automated, adaptive methods for forest inventory

Forest inventories are used to estimate forest characteristics and the condition of forest for many different applications: operational tree logging for forest industry, forest health state estimation, carbon balance estimation, land-cover and land use analysis in order to avoid forest degradation etc. Recent inventory methods are strongly based on remote sensing data combined with field sample measurements, which are used to define estimates covering the whole area of interest. Remote sensing data from satellites, aerial photographs or aerial laser scannings are used, depending on the scale of inventory. To be applicable in operational use, forest inventory methods need to be easily adjusted to local conditions of the study area at hand. All the data handling and parameter tuning should be objective and automated as much as possible. The methods also need to be robust when applied to different forest types. Since there generally are no extensive direct physical models connecting the remote sensing data from different sources to the forest parameters that are estimated, mathematical estimation models are of "black-box" type, connecting the independent auxiliary data to dependent response data with linear or nonlinear arbitrary models. To avoid redundant complexity and over-fitting of the model, which is based on up to hundreds of possibly collinear variables extracted from the auxiliary data, variable selection is needed. To connect the auxiliary data to the inventory parameters that are estimated, field work must be performed. In larger study areas with dense forests, field work is expensive, and should therefore be minimized. To get cost-efficient inventories, field work could partly be replaced with information from formerly measured sites, databases. The work in this thesis is devoted to the development of automated, adaptive computation methods for aerial forest inventory. The mathematical model parameter definition steps are automated, and the cost-efficiency is improved by setting up a procedure that utilizes databases in the estimation of new area characteristics.

Measurement and simulation of solar radiation availability in relation to the growth of coffee plants in an agroforestry system with rubber trees

Solar radiation is an important factor for plant growth, being its availability to understory crops strongly modified by trees in an Agroforestry System (AFS). Coffee trees (Coffea arabica - cv. Obatã IAC 1669-20) were planted at a 3.4 x 0.9 m spacing inside and aside rows of monocrops of 12 year-old rubber trees (Hevea spp.), in Piracicaba-SP, Brazil (22º42'30" S, 47º38'00" W - altitude: 546m). One-year-old coffee plants exposed to 25; 30; 35; 40; 45; 80; 90; 95 and 100% of the total solar radiation were evaluated according to its biophysical parameters of solar radiation interception and capture. The Goudriaan (1977) adapted by Bernardes et al. (1998) model for radiation attenuation fit well to the measured data. Coffee plants tolerate a decrease in solar radiation availability to 50% without undergoing a reduction on growth and LAI, which was approximately 2m².m-2 under this condition. Further reductions on the availability of solar radiation caused a reduction in LAI (1.5m².m-2), thus poor land cover and solar radiation interception, resulting in growth reduction.

Landscape structure in the northern coast of Paraná state, a hotspot for the brazilian Atlantic Forest conservation

The "Serra do Mar" region comprises the largest remnant of the Brazilian Atlantic Forest. The coast of the Paraná State is part of the core area of the "Serra do Mar" corridor and where actions for biodiversity conservation must be planned. In this study we aimed at characterizing the landscape structure in the APA-Guaraqueçaba, the largest protected area in this region, in order to assist environmental policies of this region. Based on a supervised classification of a mosaic of LANDSAT-5-TM satellite images (from March 2009), we developed a map (1:75,000 scale) with seven classes of land use and land cover and analyzed the relative quantities of forests and modified areas in slopes and lowlands. The APA-Guaraqueçaba is comprised mainly by the Dense Ombrophilous Forest (68.6% of total area) and secondary forests (9.1%), indicating a forested landscape matrix; anthropogenic and bare soil areas (0.8%) and the Pasture/Grasslands class (4.2%) were less representative. Slopes were less fragmented and more preserved (96.3% of Dense Ombrophilous Forest and secondary forest) than lowlands (71.3%), suggesting that restoration initiatives in the lowlands must be stimulated in this region. We concluded that most of the region sustains well-conserved ecosystems, highlighting the importance of Paraná northern coast for the biodiversity maintenance of the Atlantic Forest.

Georeferenced database generation with the purpose of hydrologic molding in reservoirs of the hydrographic basin of Jaguaribe river in the state of Ceará, Brazil

The edafoclimatic conditions of the Brazilian semiarid region favor the water loss by surface runoff. The state of Ceará, almost completely covered by semiarid, has developed public policies for the construction of dams in order to attend the varied water demand. Several hydrological models were developed to support decisive processes in the complex management of reservoirs. This study aimed to establish a methodology for obtaining a georeferenced database suitable for use as input data in hydrological modeling in the semiarid of Ceará. It was used images of Landsat satellite and SRTM Mission, and soil maps of the state of Ceará. The Landsat images allowed the determination of the land cover and the SRTM Mission images, the automatic delineation of hydrographic basins. The soil type was obtained through the soil map. The database was obtained for Jaguaribe River hydrographic basin, in the state of Ceará, and is applicable to hydrological modeling based on the Curve Number method for estimating the surface runoff.

Applicabilité de la texture couleur à la différentiation des classes d’occupation du territoire sur des images satellitales multispectrales

La texture est un élément clé pour l’interprétation des images de télédétection à fine résolution spatiale. L’intégration de l’information texturale dans un processus de classification automatisée des images se fait habituellement via des images de texture, souvent créées par le calcul de matrices de co-occurrences (MCO) des niveaux de gris. Une MCO est un histogramme des fréquences d’occurrence des paires de valeurs de pixels présentes dans les fenêtres locales, associées à tous les pixels de l’image utilisée; une paire de pixels étant définie selon un pas et une orientation donnés. Les MCO permettent le calcul de plus d’une dizaine de paramètres décrivant, de diverses manières, la distribution des fréquences, créant ainsi autant d’images texturales distinctes. L’approche de mesure des textures par MCO a été appliquée principalement sur des images de télédétection monochromes (ex. images panchromatiques, images radar monofréquence et monopolarisation). En imagerie multispectrale, une unique bande spectrale, parmi celles disponibles, est habituellement choisie pour générer des images de texture. La question que nous avons posée dans cette recherche concerne justement cette utilisation restreinte de l’information texturale dans le cas des images multispectrales. En fait, l’effet visuel d’une texture est créé, non seulement par l’agencement particulier d’objets/pixels de brillance différente, mais aussi de couleur différente. Plusieurs façons sont proposées dans la littérature pour introduire cette idée de la texture à plusieurs dimensions. Parmi celles-ci, deux en particulier nous ont intéressés dans cette recherche. La première façon fait appel aux MCO calculées bande par bande spectrale et la seconde utilise les MCO généralisées impliquant deux bandes spectrales à la fois. Dans ce dernier cas, le procédé consiste en le calcul des fréquences d’occurrence des paires de valeurs dans deux bandes spectrales différentes. Cela permet, en un seul traitement, la prise en compte dans une large mesure de la « couleur » des éléments de texture. Ces deux approches font partie des techniques dites intégratives. Pour les distinguer, nous les avons appelées dans cet ouvrage respectivement « textures grises » et « textures couleurs ». Notre recherche se présente donc comme une analyse comparative des possibilités offertes par l’application de ces deux types de signatures texturales dans le cas spécifique d’une cartographie automatisée des occupations de sol à partir d’une image multispectrale. Une signature texturale d’un objet ou d’une classe d’objets, par analogie aux signatures spectrales, est constituée d’une série de paramètres de texture mesurés sur une bande spectrale à la fois (textures grises) ou une paire de bandes spectrales à la fois (textures couleurs). Cette recherche visait non seulement à comparer les deux approches intégratives, mais aussi à identifier la composition des signatures texturales des classes d’occupation du sol favorisant leur différentiation : type de paramètres de texture / taille de la fenêtre de calcul / bandes spectrales ou combinaisons de bandes spectrales. Pour ce faire, nous avons choisi un site à l’intérieur du territoire de la Communauté Métropolitaine de Montréal (Longueuil) composé d’une mosaïque d’occupations du sol, caractéristique d’une zone semi urbaine (résidentiel, industriel/commercial, boisés, agriculture, plans d’eau…). Une image du satellite SPOT-5 (4 bandes spectrales) de 10 m de résolution spatiale a été utilisée dans cette recherche. Puisqu’une infinité d’images de texture peuvent être créées en faisant varier les paramètres de calcul des MCO et afin de mieux circonscrire notre problème nous avons décidé, en tenant compte des études publiées dans ce domaine : a) de faire varier la fenêtre de calcul de 3*3 pixels à 21*21 pixels tout en fixant le pas et l’orientation pour former les paires de pixels à (1,1), c'est-à-dire à un pas d’un pixel et une orientation de 135°; b) de limiter les analyses des MCO à huit paramètres de texture (contraste, corrélation, écart-type, énergie, entropie, homogénéité, moyenne, probabilité maximale), qui sont tous calculables par la méthode rapide de Unser, une approximation des matrices de co-occurrences, c) de former les deux signatures texturales par le même nombre d’éléments choisis d’après une analyse de la séparabilité (distance de Bhattacharya) des classes d’occupation du sol; et d) d’analyser les résultats de classification (matrices de confusion, exactitudes, coefficients Kappa) par maximum de vraisemblance pour conclure sur le potentiel des deux approches intégratives; les classes d’occupation du sol à reconnaître étaient : résidentielle basse et haute densité, commerciale/industrielle, agricole, boisés, surfaces gazonnées (incluant les golfs) et plans d’eau. Nos principales conclusions sont les suivantes a) à l’exception de la probabilité maximale, tous les autres paramètres de texture sont utiles dans la formation des signatures texturales; moyenne et écart type sont les plus utiles dans la formation des textures grises tandis que contraste et corrélation, dans le cas des textures couleurs, b) l’exactitude globale de la classification atteint un score acceptable (85%) seulement dans le cas des signatures texturales couleurs; c’est une amélioration importante par rapport aux classifications basées uniquement sur les signatures spectrales des classes d’occupation du sol dont le score est souvent situé aux alentours de 75%; ce score est atteint avec des fenêtres de calcul aux alentours de11*11 à 15*15 pixels; c) Les signatures texturales couleurs offrant des scores supérieurs à ceux obtenus avec les signatures grises de 5% à 10%; et ce avec des petites fenêtres de calcul (5*5, 7*7 et occasionnellement 9*9) d) Pour plusieurs classes d’occupation du sol prises individuellement, l’exactitude dépasse les 90% pour les deux types de signatures texturales; e) une seule classe est mieux séparable du reste par les textures grises, celle de l’agricole; f) les classes créant beaucoup de confusions, ce qui explique en grande partie le score global de la classification de 85%, sont les deux classes du résidentiel (haute et basse densité). En conclusion, nous pouvons dire que l’approche intégrative par textures couleurs d’une image multispectrale de 10 m de résolution spatiale offre un plus grand potentiel pour la cartographie des occupations du sol que l’approche intégrative par textures grises. Pour plusieurs classes d’occupations du sol un gain appréciable en temps de calcul des paramètres de texture peut être obtenu par l’utilisation des petites fenêtres de traitement. Des améliorations importantes sont escomptées pour atteindre des exactitudes de classification de 90% et plus par l’utilisation des fenêtres de calcul de taille variable adaptées à chaque type d’occupation du sol. Une méthode de classification hiérarchique pourrait être alors utilisée afin de séparer les classes recherchées une à la fois par rapport au reste au lieu d’une classification globale où l’intégration des paramètres calculés avec des fenêtres de taille variable conduirait inévitablement à des confusions entre classes.

Suivi des changements des utilisations/occupations du sol en milieu urbain par imagerie satellitale de résolution spatiale moyenne : le cas de la région métropolitaine de Montréal

De nos jours les cartes d’utilisation/occupation du sol (USOS) à une échelle régionale sont habituellement générées à partir d’images satellitales de résolution modérée (entre 10 m et 30 m). Le National Land Cover Database aux États-Unis et le programme CORINE (Coordination of information on the environment) Land Cover en Europe, tous deux fondés sur les images LANDSAT, en sont des exemples représentatifs. Cependant ces cartes deviennent rapidement obsolètes, spécialement en environnement dynamique comme les megacités et les territoires métropolitains. Pour nombre d’applications, une mise à jour de ces cartes sur une base annuelle est requise. Depuis 2007, le USGS donne accès gratuitement à des images LANDSAT ortho-rectifiées. Des images archivées (depuis 1984) et des images acquises récemment sont disponibles. Sans aucun doute, une telle disponibilité d’images stimulera la recherche sur des méthodes et techniques rapides et efficaces pour un monitoring continue des changements des USOS à partir d’images à résolution moyenne. Cette recherche visait à évaluer le potentiel de telles images satellitales de résolution moyenne pour obtenir de l’information sur les changements des USOS à une échelle régionale dans le cas de la Communauté Métropolitaine de Montréal (CMM), une métropole nord-américaine typique. Les études précédentes ont démontré que les résultats de détection automatique des changements dépendent de plusieurs facteurs tels : 1) les caractéristiques des images (résolution spatiale, bandes spectrales, etc.); 2) la méthode même utilisée pour la détection automatique des changements; et 3) la complexité du milieu étudié. Dans le cas du milieu étudié, à l’exception du centre-ville et des artères commerciales, les utilisations du sol (industriel, commercial, résidentiel, etc.) sont bien délimitées. Ainsi cette étude s’est concentrée aux autres facteurs pouvant affecter les résultats, nommément, les caractéristiques des images et les méthodes de détection des changements. Nous avons utilisé des images TM/ETM+ de LANDSAT à 30 m de résolution spatiale et avec six bandes spectrales ainsi que des images VNIR-ASTER à 15 m de résolution spatiale et avec trois bandes spectrales afin d’évaluer l’impact des caractéristiques des images sur les résultats de détection des changements. En ce qui a trait à la méthode de détection des changements, nous avons décidé de comparer deux types de techniques automatiques : (1) techniques fournissant des informations principalement sur la localisation des changements et (2)techniques fournissant des informations à la fois sur la localisation des changements et sur les types de changement (classes « de-à »). Les principales conclusions de cette recherche sont les suivantes : Les techniques de détection de changement telles les différences d’image ou l’analyse des vecteurs de changements appliqués aux images multi-temporelles LANDSAT fournissent une image exacte des lieux où un changement est survenu d’une façon rapide et efficace. Elles peuvent donc être intégrées dans un système de monitoring continu à des fins d’évaluation rapide du volume des changements. Les cartes des changements peuvent aussi servir de guide pour l’acquisition d’images de haute résolution spatiale si l’identification détaillée du type de changement est nécessaire. Les techniques de détection de changement telles l’analyse en composantes principales et la comparaison post-classification appliquées aux images multi-temporelles LANDSAT fournissent une image relativement exacte de classes “de-à” mais à un niveau thématique très général (par exemple, bâti à espace vert et vice-versa, boisés à sol nu et vice-versa, etc.). Les images ASTER-VNIR avec une meilleure résolution spatiale mais avec moins de bandes spectrales que LANDSAT n’offrent pas un niveau thématique plus détaillé (par exemple, boisés à espace commercial ou industriel). Les résultats indiquent que la recherche future sur la détection des changements en milieu urbain devrait se concentrer aux changements du couvert végétal puisque les images à résolution moyenne sont très sensibles aux changements de ce type de couvert. Les cartes indiquant la localisation et le type des changements du couvert végétal sont en soi très utiles pour des applications comme le monitoring environnemental ou l’hydrologie urbaine. Elles peuvent aussi servir comme des indicateurs des changements de l’utilisation du sol. De techniques telles l’analyse des vecteurs de changement ou les indices de végétation son employées à cette fin.

Régionalisation et synthèse des patrons de la végétation du Québec : utilisation d'indices de patrons à l'échelle provinciale.

Le Québec est une immense province à l’intérieur de laquelle existe une grande diversité de conditions bioclimatiques et où les perturbations anthropiques et naturelles du couvert végétal sont nombreuses. À l’échelle provinciale, ces multiples facteurs interagissent pour sculpter la composition et la distribution des paysages. Les objectifs généraux de cette recherche visaient à explorer et comprendre la distribution spatiale des patrons des paysages du Québec, de même qu’à caractériser les patrons observés à partir d’images satellitaires. Pour ce faire, les patrons des paysages ont été quantifiés avec un ensemble complet d’indices calculés à partir d’une cartographie de la couverture végétale. Plusieurs approches ont été développées et appliquées pour interpréter les valeurs d’indices sur de vastes étendues et pour cartographier la distribution des patrons des paysages québécois. Les résultats ont révélé que les patrons de la végétation prédits par le Ministère des Ressources naturelles du Québec divergent des patrons de la couverture végétale observée. Ce mémoire dresse un portrait des paysages québécois et les synthétise de manière innovatrice, en plus de démontrer le potentiel d’utilisation des indices comme attributs biogéographiques à l’échelle nationale.