Deforestación y Pérdida de hábitat en Bosques de montaña en la Cuenca alta del Río Zamora (Loja, Ecuador)

Los bosques húmedos de montaña se encuentran reconocidos como uno de los ecosistemas más amenazados en el mundo, llegando inclusive a ser considerado como un “hotspot” por su alta diversidad y endemismo. La acelerada pérdida de cobertura vegetal de estos bosques ha ocasionado que, en la actualidad, se encuentren restringidos a una pequeña fracción de su área de distribución histórica. Pese a esto, los estudios realizados sobre cual es efecto de la deforestación, fragmentación, cambios de uso de suelo y su efecto en las comunidades de plantas presentes en este tipo de vegetación aún son muy escuetos, en comparación a los realizados con sus similares amazónicos. En este trabajo, el cual se encuentra dividido en seis capítulos, abordaremos los siguientes objetivos: a) Comprender cuál es la dinámica que han seguido los diferentes tipos de bosques montanos andinos de la cuenca del Rio Zamora, Sur de Ecuador durante entre 1976 y 2002. b) Proveer de evidencia de las tasas de deforestación y fragmentación de todos los tipos diferentes de bosques montanos andinos presentes en la cuenca del Rio Zamora, Sur de Ecuador entre 1976 y 2002. c) Determinar qué factores inducen a la fragmentación de bosques de montaña en la cuenca alta del río Zamora entre 1976 y 2002. d) Determinar cuáles son y cómo afectan los factores ambientales y socioeconómicos a la dinámica de la deforestación y regeneración (pérdida y recuperación del hábitat) sufrida por los bosques de montaña dentro de la zona de estudio y e) Determinar si la deforestación y fragmentación actúan sobre la diversidad y estructura de las comunidades de tres tipos de organismos (comunidades de árboles, comunidades de líquenes epífitos y comunidades de hepáticas epífitas). Este estudio se centró en el cuenca alta del río Zamora, localizada al sur de Ecuador entre las coordenadas 3º 00´ 53” a 4º 20´ 24.65” de latitud sur y 79º 49´58” a 78º 35´ 38” de longitud oeste, que cubre alrededor de 4300 km2 de territorio situado entre las capitales de las provincias de Loja y Zamora-Chinchipe. Con objeto de predecir la dinámica futura de la deforestación en la región de Loja y cómo se verán afectados los diferentes tipos de hábitat, así como para detectar los factores que más influyen en dicha dinámica, se han construido modelos basados en la historia de la deforestación derivados de fotografías aéreas e imágenes satelitales de tres fechas (1976, 1989 y 2002). La cuantificación de la deforestación se realizó mediante la tasa de interés compuesto y para la caracterización de la configuración espacial de los fragmentos de bosque nativo se calcularon índices de paisaje los cuales fueron calculados utilizando el programa Fragstats 3.3. Se ha clasificado el recubrimiento del terreno en forestal y no forestal y se ha modelado su evolución temporal con Modelos Lineales Generalizados Mixtos (GLMM), empleando como variables explicativas tanto variables ambientales espacialmente explícitas (altitud, orientación, pendiente, etc) como antrópicas (distancia a zonas urbanizadas, deforestadas, caminos, entre otras). Para medir el efecto de la deforestación sobre las comunidades modelo (de árboles, líquenes y hepáticas) se monitorearon 11 fragmentos de vegetación de distinto tamaño: dos fragmentos de más de cien hectáreas, tres fragmentos de entre diez y noventa ha y seis fragmentos de menos de diez hectáreas. En ellos se instalaron un total de 38 transectos y 113 cuadrantes de 20 x 20 m a distancias que se alejaban progresivamente del borde en 10, 40 y 80 m. Nuestros resultados muestran una tasa media anual de deforestación del 1,16% para todo el período de estudio, que el tipo de vegetación que más alta tasa de destrucción ha sufrido, es el páramo herbáceo, con un 2,45% anual. El análisis de los patrones de fragmentación determinó un aumento en 2002 de más del doble de fragmentos presentes en 1976, lo cual se repite en el análisis del índice de densidad promedio. El índice de proximidad media entre fragmentos muestra una reducción progresiva de la continuidad de las áreas forestadas. Si bien las formas de los fragmentos se han mantenido bastante similares a lo largo del período de estudio, la conectividad entre estos ha disminuido en un 84%. Por otro lado, de nuestros análisis se desprende que las zonas con mayor probabilidad de deforestarse son aquellas que están cercanas a zonas previamente deforestadas; la cercanía a las vías también influye significativamente en la deforestación, causando un efecto directo en la composición y estructura de las comunidades estudiadas, que en el caso de los árboles viene mediado por el tamaño del fragmento y en el caso del componente epífito (hepáticas y líquenes), viene mediado tanto por el tamaño del fragmento como por la distancia al borde del mismo. Se concluye la posibilidad de que, de mantenerse esta tendencia, este tipo de bosques desaparecerá en corto tiempo y los servicios ecosistémicos que prestan, se verán seriamente comprometidos. ABSTRACT Mountain rainforests are recognized as one of the most threatened ecosystems in the world, and have even come to be considered as a “hotspot” due to their high degree of diversity and endemism. The accelerated loss of plant cover of these forests has caused them to be restricted today to a small fraction of their area of historic distribution. In spite of this, studies done on the effect of deforestation, fragmentation, changes in soil use and their effect on the plant communities present in this type of vegetation are very brief compared to those done on their analogues in the Amazon region. In this study, which is divided into six chapters, we will address the following objectives: a) To understand what the dynamic followed by the different types of Andean mountain forests in the Zamora River watershed of southern Ecuador has been between 1976 and 2002. b) To provide evidence of the rates of deforestation and fragmentation of all the different types of Andean mountain forests existing in the upper watershed of the Zamora River between 1976 and 2002. c) To determine the factors that induces fragmentation of all different types of Andean mountain forests existing in the upper watershed of the Zamora River between 1976 and 2002. d) To determine what the environmental and anthropogenic factors are driving the dynamic of deforestation and regeneration (loss and recuperation of the habitat) suffered by the mountain forests in the area of the study and e) To determine if the deforestation and fragmentation act upon the diversity and structure of three model communities: trees, epiphytic lichens and epiphytic liverworts. This study is centered on the upper Zamora River watershed, located in southern Ecuador between 3º 00´ 53” and 4º 20´ 24.65 south latitude and 79º 49´ 58” to 78º 35´ 38” west longitude, and covers around 4,300 km2 of territory located between Loja and Zamora-Chinchipe provinces. For the purpose of predicting the future dynamic of deforestation in the Loja region and how different types of habitats will be affected, as well as detecting the environmental and socioeconomic factors that influence landscape dynamics, models were constructed based on deforestation history, derived from aerial photographs and satellite images for three dates (1976, 1989 and 2002). Quantifying the deforestation was done using the compound interest rate; to characterize the spatial configuration of fragments of native forest, landscape indices were calculated with Fragstats 3.3 program. Land cover was classified as forested and not forested and its evolution over time was modeled with Generalized Linear Mixed Models (GLMM), using spatially explicit environmental variables (altitude, orientation, slope, etc.) as well as anthropic variables (distance to urbanized, deforested areas and roads, among others) as explanatory variables. To measure the effects of fragmentation on three types of model communities (forest trees and epiphytic lichen and liverworts), 11 vegetation fragments of different sizes were monitored: two fragments of more than one hundred hectares, three fragments of between ten and ninety ha and six fragments of fewer than ten hectares . In these fragments, a total of 38 transects and 113 20 x 20 m quadrats were installed at distances that progressively moved away from the edge of the fragment by 10, 40 and 80 m. Our results show an average annual rate of deforestation of 1.16% for the entire period of the study, and that the type of vegetation that suffered the highest rate of destruction was grassy paramo, with an annual rate of 2.45%. The analysis of fragmentation patterns determined the number of fragments in 2002 more than doubled the number of fragments present in 1976, and the same occurred for the average density index. The variation of the average proximity index among fragments showed a progressive reduction of the continuity of forested areas. Although fragment shapes have remained quite similar over the period of the study, connectivity among them has diminished by 84%. On the other hand, it emerged from our analysis that the areas of greatest probability of deforestation were those that are close to previously deforested areas; proximity to roads also significantly favored the deforestation causing a direct effect on the composition of our model communities, that in the case of forest trees is determined by the size of the fragment, and in the case of the epiphyte communities (liverworts and lichens), is determined, by the size of the fragment as well as the distance to edge. A subject under discussion is the possibility that if this tendency continues, this type of forest will disappear in a short time, and the ecological services it provides, will be seriously endangered.

Os impactos ambientais no meio físico: erosão e assoreamento na bacia hidrográfica do rio Taquari, MS, em decorrência da pecuária

O aumento da inundação em áreas do baixo curso do rio Taquari, no Pantanal do estado do Mato Grosso do Sul, tem transformado a pecuária desta região numa atividade com baixa rentabilidade, à medida que extensas áreas de campo passaram a ser inundadas vários meses durante o ano a partir da década de 70. A pecuária realizada em campos naturais de regiões úmidas do Pantanal indica que há necessidade de se investigar metodologias apropriadas para avaliação de impacto ambiental, que abordem impactos diretos, indiretos, cumulativos e processos do meio físico que alteram, de maneira prejudicial, o meio ambiente. Supõe-se que a inundação na planície do rio Taquari esteja relacionada com a ocupação antrópica nas áreas de planalto da bacia do rio Taquari. O presente trabalho tem por objetivo avaliar os impactos ambientais na planície de inundação do baixo curso do Taquari, decorrentes da ocupação antrópica da bacia hidrográfica do rio Taquari em sua totalidade, considerando os impactos ambientais causados pela pecuária à medida que se configura como principal atividade econômica da bacia bem como os processos erosivos e de assoreamento no quadro atual do regime de inundações. As etapas de caracterização da área, de análise dos impactos e as propostas de ações mitigadoras, previstas num Estudo de Impacto Ambiental, foram aqui analisadas. Foram utilizadas informações sobre as características do meio físico, biótico e socioeconômico, selecionadas a partir do levantamento dos dados existentes com recorte efetuado para a bacia hidrográfica do rio Taquari. Na maior parte dos temas, este foi um processo de levantamento, ordenamento e recuperação de informações, na escala original de 1:250.000, do Plano de Conservação da Bacia do Alto Paraguai-PCBAP, gerenciado no SPRING. Foram também realizadas viagens de campo para a complementação dos dados e para o levantamento de atividades antrópicas com verificações \"in loco\" da ocorrência de impacto ambiental. A maioria dos dados socioeconômicos compilados para o presente trabalho teve por base os censos agropecuários e demográficos realizados pelo IBGE. Os resultados obtidos demonstram que os impactos ambientais decorrentes da pecuária no planalto interferem no regime de inundação na planície da bacia, o que só foi possível de ser identificado a partir de análises integradas em toda a bacia hidrográfica do rio Taquari. Verificou-se que os métodos de EIA são adequados para identificar os impactos diretos decorrentes da pecuária, mas não são adequados para identificar os processos e seus efeitos cumulativos na extensão da bacia hidrográfica do rio Taquari. Além disto, a abordagem da avaliação ambiental estratégica, como procedimento para análise ambiental em políticas, planos e programas, mostra-se adequada para as análises na BHRT à medida que está centralizada nos efeitos do ambiente sobre as necessidades e oportunidades de desenvolvimento. Contudo, somente a recuperação de danos ambientais, o controle das origens dos impactos no ambiente e um sistema de gestão consciente de seus compromissos podem levar, juntamente com a melhora dos procedimentos técnicos e administrativos para análises ambientais, à uma maior proximidade da sustentabilidade ambiental na BHRT.

Características de ecossistema e sustentabilidade: estudo sobre a vegetação original e remanescente da bacia hidrográfica do Rio Canoas

Com o objetivo de comparar em termos de paisagem, a vegetação florestal nativa com os fragmentos da vegetação remanescente na bacia hidrográfica do Rio Canoas, situada na divisa dos Estados de São Paulo (à Nordeste) e Minas Gerais (à Centro-Oeste), busca-se dar subsídios aos atores locais para trabalhar políticas públicas com vistas à sustentabilidade do uso dos recursos naturais na região. Foram mapeados sete tipos de vegetação, com destaque para o cerrado e a floresta estacional semidecidual. O método do estudo é exploratório, realizado em fontes secundárias, pesquisa bibliográfica e documental, baseado no estudo de mapas e suas simulações, dados do IBGE, Embrapa, Ministério do Meio Ambiente e trabalhos de autores, como: Cogo (1987), Coutinho (2000) e Bozini (2005). Observou-se que a bacia hidrográfica objeto em estudo é degradada historicamente pelos diversos usos e ocupações do solo e sua maior parte ainda está sujeita à um elevado nível de perturbação, seja pela pecuária, ou mesmo pelas plantações de laranja, milho, eucalipto, cana-de-açúcar e, principalmente, pela cultura do café.

Eastern Africa showing expedition route of Speke and Grant, 1863 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Map of the route explored by Captns. Speke & Grant from Zanzibar to Egypt : showing the outfall of the Nile from the Victoria Nyanza (Lake) and the various Negro territories discovered by them. It was published by Edward Stanford in 1863. Scale [ca. 1:5,800,000]. Covers portions of north and eastern Africa including parts of Sudan, Eritrea, Ethiopia, Uganda, Kenya, Rwanda, Burundi, and Tanzania. The image inside the map neatline is georeferenced to the surface of the earth and projected to the 'World Mercator' projection. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, expedition routes of John Speke and James Grant, cities and other human settlements, tribe and territorial boundaries, and more. Relief is shown by hachures. Includes location map and text. This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection and the Harvard University Library as part of the Open Collections Program at Harvard University project: Organizing Our World: Sponsored Exploration and Scientific Discovery in the Modern Age. Maps selected for the project correspond to various expeditions and represent a range of regions, originators, ground condition dates, scales, and purposes.

Remote Sensing of Vegetation Change Across a Latitudinal Gradient in the Canadian Arctic

Global air surface temperatures and precipitation have increased over the last several decades resulting in a trend of greening across the Circumpolar Arctic. The spatial variability of warming and the inherent effects on plant communities has not proven to be uniform or homogeneous on global or local scales. We can apply remote sensing vegetation indices such as the Normalized Difference Vegetation Index (NDVI) to map and monitor vegetation change (e.g., phenology, greening, percent cover, and biomass) over time. It is important to document how Arctic vegetation is changing, as it will have large implications related to global carbon and surface energy budgets. The research reported here examined vegetation greening across different spatial and temporal scales at two disparate Arctic sites: Apex River Watershed (ARW), Baffin Island, and Cape Bounty Arctic Watershed Observatory (CBAWO), Melville Island, NU. To characterize the vegetation in the ARW, high spatial resolution WorldView-2 data were processed to create a supervised land-cover classification and model percent vegetation cover (PVC) (a similar process had been completed in a previous study for the CBAWO). Meanwhile, NDVI data spanning the past 30 years were derived from intermediate resolution Landsat data at the two Arctic sites. The land-cover classifications at both sites were used to examine the Landsat NDVI time series by vegetation class. Climate variables (i.e., temperature, precipitation and growing season length (GSL) were examined to explore the potential relationships of NDVI to climate warming. PVC was successfully modeled using high resolution data in the ARW. PVC and plant communities appear to reside along a moisture and altitudinal gradient. The NDVI time series demonstrated an overall significant increase in greening at the CBAWO (High Arctic site), specifically in the dry and mesic vegetation type. However, similar overall greening was not observed for the ARW (Low Arctic site). The overall increase in NDVI at the CBAWO was attributed to a significant increase in July temperatures, precipitation and GSL.

(Table 2) Stable carbon and oxygen isotope record of benthic foraminifera of eastern Mediterranean Sea sediment cores

The paleoproductivity, paleo-oxygenation, and paleohydrographic configuration of the southeastern Mediterranean during the late Holocene was reconstructed on the basis of the isotopic composition of the epibenthic Heterolepa floridana, shallow-endobenthic Uvigerina mediterranea, and the deeper endobenthic Bulimina inflata from two high-resolution cores GA-112 (470 m) and GA-110 (670 m). The Delta d13C between H. floridana and U. mediterranea reveals four intervals of enhanced productivity, from 3.3-2.6, 2.3-1.9, 1.5-1.1, and 0.8-0.4 kyr BP, coinciding with increased nutrient supply by the Nile River. The entire basin was well aerated, with oxygen consumption varying between 1.0 and 3.5 mL O2/L. Oxygen consumption increases toward present day, probably because of higher accumulation of total organic carbon at 1.7 kyr BP, coinciding with the appearance of the mesotropic benthic species. The hydrographic configuration of the basin has changed during the course of the last 3.75 kyr. The Levantine Intermediate Water (LIW) deepens below 470 m between 3.3 and 2.0 kyr, and especially between 2.5 and 2.0 kyr. During the last 1.5 kyr, the LIW becomes shallower than 470 m, similar to the present day. The change in the hydrographic configuration reflects changes in evaporation/precipitation ratio and in temperature.

Sr and Nd composition, stable oxygen isotope record and lithogenic flux in eastern Mediterranean sediments

Isotopic ratios of Sr and Nd from lithogenic components of three isochronous core sections recovered from an east-west transect in the Eastern Mediterranean Sea (EMS) have been analyzed. The data are used for a quantitative estimate of the temporal and spatial variation of detrital flux to the EMS, assuming Saharan dust and Aegean/Nile particulate matter as dominant end members. It was established that the carbonate-free Saharan dust flux during deposition of the nonsapropel layers of marine oxygen isotope stage 5.4 (MIS 5.4) was similar to the present flux. During the deposition of sapropels S5 and S6, however, the Saharan dust input was drastically reduced and was not balanced by a change in the riverine influx at this time. Denser vegetation cover during more humid conditions may have reduced physical erosion and sediment removal in the source area. During marine oxygen isotope stage 6.2 (MIS 6.2) a pronounced increase of Saharan dust and detrital influx from the Aegean region is evident and implies more arid conditions in the southern and northern catchment areas. During this period, intersite variations are interpreted in terms of their geographic location relative to the seaways connecting the Aegean Sea and EMS. The width of the straits and hence the amount of sediment entering the eastern basins may have been affected by a low sea level that impeded interbasin sediment dispersal.

Organic and inorganic geochemical measurements of sediment core GeoB7702-3

It has long been known that extreme changes in North African hydroclimate occurred during the late Pleistocene yet many discrepancies exist between sites regarding the timing, duration and abruptness of events such as Heinrich Stadial (HS) 1 and the African Humid Period (AHP). The hydroclimate history of the Nile River is of particular interest due to its lengthy human occupation history yet there are presently few continuous archives from the Nile River corridor, and pre-Holocene studies are rare. Here we present new organic and inorganic geochemical records of Nile Basin hydroclimate from an eastern Mediterranean (EM) Sea sediment core spanning the past 28 ka BP. Our multi-proxy records reflect the fluctuating inputs of Blue Nile versus White Nile material to the EM Sea in response to gradual changes in local insolation and also capture abrupt hydroclimate events driven by remote climate forcings, such as HS1. We find strong evidence for extreme aridity within the Nile Basin evolving in two distinct phases during HS1, from 17.5 to 16 ka BP and from 16 to 14.5 ka BP, whereas peak wet conditions during the AHP are observed from 9 to 7 ka BP. We find that zonal movements of the Congo Air Boundary (CAB), and associated shifts in the dominant moisture source (Atlantic versus Indian Ocean moisture) to the Nile Basin, likely contributed to abrupt hydroclimate variability in northern East Africa during HS1 and the AHP as well as to non-linear behavior of hydroclimate proxies. We note that different proxies show variable gradual and abrupt responses to individual hydroclimate events, and thus might have different inherent sensitivities, which may be a factor contributing to the controversy surrounding the abruptness of past events such as the AHP. During the Late Pleistocene the Nile Basin experienced extreme hydroclimate fluctuations, which presumably impacted Paleolithic cultures residing along the Nile corridor.

Physical Characteristics and Estimated Hydraulic Conductivity from Grain-Size Distribution for Vashon Till, South-Central Puget Lowland, Washington

Population growth, urban development, and increased commercial and industrial activity in the south-central Puget Lowlands of Washington State has led to an increased demand for groundwater. The Vashon till is a glacially consolidated, low-permeability unit comprising unstratified clay, silt, cobbles and boulders with ubiquitous coarse-grained lenses and is an extensive surficial unit throughout the south-central Puget Lowland. Thus, understanding the physical and hydrological characteristics – specifically, the hydraulic conductivity – of this unit is a necessary component of a groundwater model. This study provides (1) a record of the physical characteristics of Vashon till deposits within the study area; and (2) an estimate of the highest, lowest, and average value of saturated hydraulic conductivity based on the grain-size distribution of Vashon till samples collected from six field sites in the Puyallup River Watershed. Analysis shows that the average moisture content ranges between about 1 and 6%, average dry bulk density is about 2.20 g/cm3, and average porosity is about 17%. Grain-size distributions show that half of the samples analyzed are well graded, while the other half is poorly graded. Grain-size distributions also show an average d10 value of about 0.20 mm, and average ff values ≤ 16%, which are key values in estimating the saturated hydraulic conductivity of over-consolidated glacial deposits. Based on these observed values, the estimates of hydraulic conductivity range from a minimum of 0.02 m/d to a maximum of 1.38 m/d in within the general Vashon till.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Radiocarbon dating and alkenone-based proxies of sediment core GeoB7702-3

In this study we utilize two organic geochemical proxies, the Uk'37 index and TEX86, to examine past sea surface temperatures (SST) from a site located near the Nile River Delta in the eastern Mediterranean (EM) Sea. The Uk'37 and TEX86 records generally are in agreement and indicate SST ranges of 14°C-26°C and 14°C-28°C, respectively, during the last 27 cal ka. During the Holocene, TEX86-based SST estimates are usually higher than Uk'37-based SST estimates, which is likely due to seasonal differences between the timing of the haptophyte and crenarchaeota blooms in the EM and is related to the onset of the modern flow regime of the Nile River. Both records show that SST varied on centennial to millennial timescales in response to global climate events, i.e., cooling during the Last Glacial Maximum (LGM), Heinrich event 1 (H1), and the Younger Dryas (YD) and warming during the Bølling-Allerød and in the early Holocene during deposition of sapropel S1. The H1 cooling was particularly severe and is marked by a drop in SST of ~4.5°C in comparison to pre-H1 SST, with temperatures >1°C cooler than during the LGM. In contrast to high-latitude and western Mediterranean records, which indicate both an abrupt onset and termination of the YD event, the transition from the YD to the Holocene was much more gradual in the EM.

Remote Sensing of Vegetation Change Across a Latitudinal Gradient in the Canadian Arctic

Global air surface temperatures and precipitation have increased over the last several decades resulting in a trend of greening across the Circumpolar Arctic. The spatial variability of warming and the inherent effects on plant communities has not proven to be uniform or homogeneous on global or local scales. We can apply remote sensing vegetation indices such as the Normalized Difference Vegetation Index (NDVI) to map and monitor vegetation change (e.g., phenology, greening, percent cover, and biomass) over time. It is important to document how Arctic vegetation is changing, as it will have large implications related to global carbon and surface energy budgets. The research reported here examined vegetation greening across different spatial and temporal scales at two disparate Arctic sites: Apex River Watershed (ARW), Baffin Island, and Cape Bounty Arctic Watershed Observatory (CBAWO), Melville Island, NU. To characterize the vegetation in the ARW, high spatial resolution WorldView-2 data were processed to create a supervised land-cover classification and model percent vegetation cover (PVC) (a similar process had been completed in a previous study for the CBAWO). Meanwhile, NDVI data spanning the past 30 years were derived from intermediate resolution Landsat data at the two Arctic sites. The land-cover classifications at both sites were used to examine the Landsat NDVI time series by vegetation class. Climate variables (i.e., temperature, precipitation and growing season length (GSL) were examined to explore the potential relationships of NDVI to climate warming. PVC was successfully modeled using high resolution data in the ARW. PVC and plant communities appear to reside along a moisture and altitudinal gradient. The NDVI time series demonstrated an overall significant increase in greening at the CBAWO (High Arctic site), specifically in the dry and mesic vegetation type. However, similar overall greening was not observed for the ARW (Low Arctic site). The overall increase in NDVI at the CBAWO was attributed to a significant increase in July temperatures, precipitation and GSL.