Acta geographica Slovenica, 62-1, 2022, 105–124 REMOTE SENSING ANALYSIS TO MAP INTER-REGIONAL SPATIO-TEMPORAL VARIATIONS OF THE VEGETATION IN ICELAND DURING 2001–2018 Haraldur Olafsson, Iman Rousta Figure: Single shrubs growing in grasslands in Hveragerði-Iceland. IM A N R O U S TA 62-1_acta49-1.qxd 20.6.2022 10:43 Page 105 106 DOI: https://doi.org/10.3986/AGS.10390 UDC: 581.9:551.583(491.1)»2001/2018« COBISS: 1.01 Haraldur Olafsson1,2,3, Iman Rousta2,3,4 Remote sensing analysis to map inter-regional spatio-temporal variations of the vegetation in Iceland during 2001–2018 ABSTRACT: Changes in the vegetation of the Arctic and sub-Arctic regions have been used as indicators of the impact and seriousness of climate change. In this study, 342 MODIS NDVI images were used to monitor and assess the variability and long-term changes in the vegetation in Iceland in the period 2001–2018. An insignificant trend in the changes of the vegetation coverage (R = 0.16, p-value = 0.05) was obtained, however, it also resulted that the area with the low values of the NDVI (< 0.6) is decreasing, whereas the area with higher values of the NDVI (> 0.6, mostly forests) is increasing. The NDVI index during the study period rose for the area of about 3260 km2, while it declined for 1635 km2. The results of this study can be used for organizing the strategies preventing climate change and global warming. KEY WORDS: Iceland, vegetation dynamics, MODIS, NDVI, anomaly analysis Kartiranje spreminjanja vegetacije v prostoru in času v različnih regijah na Islandiji med letoma 2001 in 2018 s pomočjo daljinskega zaznavanja Spremembe v vegetaciji arktičnih in subarktičnih pokrajin so pokazatelj vpliva in pomembnosti podneb- nih sprememb. V pričujoči raziskavi smo 342 posnetkov NDVI senzorja MODIS uporabili za spremljanje in oceno spremenljivosti vegetacije in njenih dolgoročnih sprememb na območje Islandije v obdobje 2001–2018. Ugotovili smo neznačilen trend spreminjanja pokrovnosti vegetacije (R = 0,16, p = 0,05), polega tega pa smo opazili, da se površina z nizkimi vrednostmi NDVI (< 0,6) zmanjšuje, površina z visokimi vrednostmi NDVI (> 0,6, predvsem gozdovi) pa povečuje. Indeks NDVI je porasel v opazovanem obdobju na območju 3260km2, zmanjšal pa se je na 1635km2. Rezultati študije so lahko uporabni za pripravo strategij preprečevanja podnebnih sprememb in globalnega segrevanja. KLJUČNE BESEDE: Islandija, dinamika spreminjanja vegetacije, MODIS, NDVI, analiza anomalij The article was submitted for publication on September 5th, 2021. Uredništvo je prejelo prispevek 5. septembra 2021. 1 University of Iceland, Department of Physics, Reykjavik, Iceland haraldur@vedur.is (https://orcid.org/0000-0002-4181-0988) 2 University of Iceland, Institute for Atmospheric Sciences-Weather and Climate, Reykjavik, Iceland irousta@yazd.ac.ir (https://orcid.org/0000-0002-3694-6936) 3 Icelandic Meteorological Office (IMO), Reykjavik, Iceland 4 Yazd University, Department of Geography, Yazd, Iran 62-1_acta49-1.qxd 20.6.2022 10:43 Page 106 1 Introduction Vegetation is an indispensable element of the Earth that links soil, air, water, and other components of the environment (Foley et al. 2000; Cui et al. 2009; Rousta et al. 2018; Rousta et al. 2020a; Mansourmoghaddam et al. 2021; Rousta et al. 2021). Environmental conditions and variability can define varaiety of future land cover (Zhang et al. 2019b; Shen et al. 2020; Wang et al. 2020; Chao et al. 2021). Vegetation dynamics carry valuable information about the impacts of global warming (Pettorelli et al. 2005; Olafsson and Rousta 2021), land degradation (Metternicht et al. 2010; Rousta et al. 2020c), and desertification (Symeonakis and Drake 2004; Rousta et al. 2020b). The Intergovernmental Panel on Climate Change (IPCC) highlights that the high northern latitudes are warming faster than other regions on the planet (Masson-Delmotte et al. 2018). As some researches point out, this is due to Polar Amplification (PA), a phenomenon in which effects such as decreasing sea ice and lower albedo due to reduced snow cover cause the Arctic temperatures to rise disproportionately under increased greenhouse gas emissions (Polyakov et al. 2003). The Arctic and sub- Arctic vegetation are highly sensitive to climate change, and the changes in the vegetation of this region have been used as indicators of the global impact of climate change. With the advent of the space-borne era, using remote sensing and satellite images for assessment and monitoring of vegetation dynamics is common in research (Miao et al. 2018; Zhao et al. 2021a). Remotely sensed studies showed that since the 1980s the vegetation throughout northern latitudes has been increas- ing (Slayback et al. 2003; Bokhorst et al. 2009; Liu et al. 2015; Raynolds et al. 2015; Merrington 2019). One of the most widely used remotely-sensed indicators of vegetation is the Normalized Difference Vegetation Index (NDVI), which is highly sensitive to ecosystem conditions (Ollinger 2011; Lemenkova 2020a). NDVI is the normalized difference between the near-infrared (NIR) and visible red band’s reflectance (Rouse et al. 1974; Tucker 1979) and has great potential for vegetation monitoring. The greater the NDVI value is, the higher is the photosynthetic capacity (Tucker 1979; Gao and Goetzt 1995; Chen and Brutsaert 1998). It can be used for detecting changes in vegetation development, such as downtrends or uptrends (Alcaraz‐Segura et al. 2010; Ghafarian Malamiri et al. 2020). Over the last few decades, many scientists have established the effective role of NDVI on vegetation variability, or drought dynamics assessment, and monitoring of them (Kogan 1991; Kogan 1995; Yang et al. 1998; McVicar and Bierwirth 2001; Ji and Peters 2003; Wan, Wang and Li 2004; Zhao et al. 2021b). The AVHRR (Advanced Very High-Resolution Radiometer) derived NDVI is available since 1981, which creates an opportunity to monitor vegetation dynamics (VD) on var- ious scales, from global or country to small regions, using time series data for different periods. The increased NDVI was found in peak parts of the Arctic, with the growth of the summer surface temperature believed to be the reason for this. In-situ measurements performed at various locations across the tundra biome dur- ing 1980–2010 have also connected the increases in sampled vegetation to summer warming (Elmendorf et al. 2012). The highest NDVI trends for the AVHRR record were observed for Iceland in the period 1982–2010 (0.008 NDVI units per year) (Epstein et al. 2012). Greenland, which was mentioned as the sec- ond with the highest increases, also had an increase of 0.005 NDVI units per year, and the rest of the northern countries experienced increases equal to 0.003 per year or less (Epstein et al. 2012). Additionally, the NDVI has decreased in certain areas, in which reduction of air temperature occurred (Bhatt et al. 2013). For example, in parts of Scandinavia, recent decreases in NDVI have been linked to winter warming events that reduced the protective snow layer (Bokhorst et al. 2009) and to unfavorable summer conditions and insect outbreaks (Bjerke et al. 2014). Soil is an important factor that influences the processes between the land surface and atmosphere (Zhang et al. 2019a; Zhao et al. 2020; Zuo et al. 2020; Li et al. 2021; Li et al. 2022; Miao et al. 2022). In terms of world ecosystems, Iceland is unique, mostly because of its soils. The island’s geology is overwhelming- ly influenced by relatively recent extrusive volcanism, which forms the main parent material for its soils. On a geologic timescale, soil formation occurred about 10,000 years ago, which coincides with the end of the last glacial period (Arnalds et al. 1995). The basaltic tephra, which is the term for material of any size ejected by volcanoes, is one of the parent materials of soil types in Iceland (James, Chester and Duncan 2000; Arnalds 2008). The land-cover of Iceland is very dynamic. Erupting volcanoes produce lava flows and ash deposits, which dramatically transform landscapes. Despite the low population density (3.1 inhabitants/km2 in 2010), it’s the human influence that affects most of the countryside (Jóhannesson 2010) and its water surrounding (Mansourmoghaddam et al. 2022). Within an Icelandic context, understanding Acta geographica Slovenica, 62-1, 2022 107 62-1_acta49-1.qxd 20.6.2022 10:43 Page 107 biodiversity is of special importance, as Iceland has a landscape history of which has been characterized by frequent, energetic geological processes (Thorarinsson 1967) and pressure on its ecosystems due to the land-use changes imposed after Iceland was settled by humans (Dugmore et al. 2009; Sigurmundsson et al. 2014; Bates et al. 2021). This has led to relatively low biodiversity in Iceland compared to the lands at similar latitudes (Jóhannesdóttir et al. 2017). The reduction in vegetation cover in Iceland in historical times has resulted in feedback loop. With less vegetation present to bind the soil, it becomes vulnerable to aeo- lian processes (Arnalds et al. 2001). This wind erosion removes useful surface nutrients, making vegetation colonization problematic (Óskarsson et al. 2004). It is reported that these processes have left large areas of Iceland barren (Arnalds et al. 2001). In the sense of environmental sustainability, Iceland is subject to volcanic eruptions, local overgraz- ing, and soil erosion. Most of the vegetation is located in the coastal zone and the human settlement areas. On the other hand, the Icelandic vegetation, same as it is with the sub-arctic and arctic vegetation, is of high importance for wild and human life. Studying vegetation dynamic and their link to climate change is therefore of high importance for securing the ecological sustainability of Iceland. However, studies that focus on Iceland are rare and almost none of them focused specifically on the assessment of the vegeta- tion evolution for the whole Iceland during recent years. Therefore, the objective of this research is to assess the spatio-temporal NDVI variations in Iceland and its main sub-regions. 2 Material and methods 2.1 Study area and data collection Iceland is located between 63 and 67°N, and 13 and 25°W. It is the second-largest island in the Northern Atlantic, having an area of 103,000 km2 and a population of 360,000. More than two-thirds of the pop- ulation live in Reykjavík city and surrounding areas in the southwestern part of the country. Reykjavík is the capital and the largest city of Iceland. The country is volcanically and geologically active and con- sists of highlands, glaciers, and mountains. Iceland is a mountainous island, yet it does not have areas with high elevations. About 26% area of Iceland lies between 0–200 m, 36% between 200–600 m, 17% between 600–800 m, and 21% above 800 m. The highest mountain is 2119 m (Thorsteinsson, Olafsson and Van Dyne 1971). The climate of the country is temperate and the archipelago lies in a tundra veg- etation zone. In the present study, the country was divided into 8 main regions. The South and East regions with 9421 and 7786 km2 have the largest share of vegetation coverage and the South East Peninsula and Reykjavik area with 894.4 and 701.6 km2 have the smallest share of Iceland’s vegetation coverage (Table 1). The south coast is wetter, windier, and warmer than the north. The snowfall in winter is more common in the north than in the south. The climate of Iceland was described in detail by Einarsson (1984) and Ólafsson et al. (2007). While the country is close to the Arctic, the coastal area remains ice-free through the winter (Figure 1). The largest parts of Iceland are covered by grasslands (about 55%, 56,250 km2), barren land (about 28%, 28,958 km2), and permanent snow and ice cover (about 10%, 10,496 km2) (Table 2). 108 Table 1: Basic characteristics of the main regions of Iceland. Name Area (km2) Vegetated Nonvegetated Vegetated Max Average area (km2) area (km2) area (%) Elevation (m) Elevation (m) South (S) 24423 9421.7 15001.3 38.6 2015 513 Southwest Peninsula (SWP) 792 701.6 90.4 88.6 610 99 Reykjavik Area (RA) 1001 894.4 106.6 89.4 1050 219 West (W) 9346 6929.4 2416.6 74.1 1717 320 West Fjords (WF) 9038 4401.0 4637.0 48.7 1012 323 Northwest (NW) 12346 7148.9 5197.1 57.9 1796 494 Northeast (NE) 21712 6980.7 14731.3 32.2 2057 580 East (E) 22074 7786.8 14287.2 35.3 2076 606 62-1_acta49-1.qxd 20.6.2022 10:43 Page 108 Acta geographica Slovenica, 62-1, 2022 109 Iceland Greenland (Denmark) Jan Mayen (Norway) S E NE W NW WF RA SWP 5°0'0"W10°0'0"W15°0'0"W20°0'0"W25°0'0"W30°0'0"W35°0'0"W 71 °0 '0" N 70 °0 '0" N 70 °0 '0" N 69 °0 '0" N 69 °0 '0" N 68 °0 '0" N 68 °0 '0" N 67 °0 '0" N 67 °0 '0" N 66 °0 '0" N 66 °0 '0" N 65 °0 '0" N 65 °0 '0" N 64 °0 '0" N 64 °0 '0" N 63 °0 '0" N 63 °0 '0" N ® 0 260130 km –1–0 0–300 300–600 600–1200 1200–1800 1800–2076 Elevation (m) Content by: Iman Rousta Map by: Iman Rousta Source: Contact to irousta@yazd.ac.ir 71 °0 '0" N Figure 1: The map of the study area showing elevation and main regions. Table 2: Iceland land cover types derived from MODIS (MCD12Q1) – Yearly Land Cover Type 1: Annual International Geosphere-Biosphere Programme (IGBP) (Loveland et al. 1999; Didan et al. 2015). Land cover type Area (km2) Percent of the whole area (%) Evergreen Needleleaf Forests 9.7 0.009 Deciduous Needleleaf Forests 1.3 0.001 Deciduous Broadleaf Forests 4.8 0.005 Mixed Forests 6.1 0.006 Open Shrublands: 177.7 0.173 Woody Savannas 1.4 0.001 Savannas 2,953.5 2.867 Grasslands 56,255.2 54.617 Permanent Wetlands 1,823.0 1.770 Urban and Built-up Lands 47.9 0.046 Permanent Snow and Ice 10,496.5 10.191 Barren 28,960.8 28.117 Water Bodies 2,271.7 2.206 Total 103,000.00 100.000 62-1_acta49-1.qxd 20.6.2022 10:43 Page 109 2.2 Methods To explore the spatio-temporal variability of vegetation coverage in Iceland, a total of 342 NDVI images were downloaded from the Terra MODIS vegetation Indices (MOD13Q1) database from the Land Processes Distributed Active Archive Center (https://lpdaacsvc.cr.usgs.gov.appeears) (Didan 2015a). The MODIS analy- sis-ready 16-day composite product (MOD13Q1.006) with 250m spatial resolution was collected for the period from 1 January 2001 to 1 November 2018. ArcGIS 10.7 software and R environment were used to perform spatial and statistical analyses. The 1 arc-second (~30 m) spatial resolution SRTM (Shuttle Radar Topography Mission) elevation data was used to visualize the terrain profile of the study area (Zandbergen 2008). NDVI is the most commonly used vegetation index (VI) for detecting the greenness of vegetation and production patterns (Tarpley, Schneider and Money 1984; Gitelson et al. 2003; Thenkabail, Gamage and Smakhtin 2004; Dutta et al. 2015). Scientists from all over the world have been using VI for investigation of various aspects of vegetation dynamic, i.e. vegetation mapping, monitoring, phonological analysis, crop growth, yields, and many more (Running et al. 1995; Moulin et al. 1997; Dabrowska-Zielinska et al. 2002; Geerken, Zaitchik and Evans 2005; Martínez and Gilabert 2009; Moniruzzaman et al. 2021). NDVI is defined as: (1) The Bnir and Bred in Eq. 1 stands for the spectral reflectance in the near-infrared band and red band, respec- tively. NDVI ranges between –1 and +1, with the values from –1 to 0 indicating the absence of green leaves and the values from 0 to +1 indicating the greenest areas. Moderate NDVI values from 0.2 to 0.3 repre- sent shrub and grassland, while high NDVI values (0.6 to 0.8) indicate dense vegetation (Montandon and Small 2008; Atasoy 2018; Jovanović, Milanović and Zorn 2018; Rousta et al. 2020b). NDVI values close to 0 represent the bare ground, while negative NDVI values correspond to water bodies snow and ice (Dye and Tucker 2003; Gandhi et al. 2015). In the present research, the NDVI data were acquired from the Terra-MODIS Vegetation Indices MOD13Q1. To calculate the vegetation coverage, the number of pixels with vegetation, that is, having the NDVI m × 250 m = 0.0625 km2). The same procedure was used for calculating the yearly and seasonal veg- etation coverages. NDVI index was defined into 7 classes (0.2–0.3, 0.3–0.4, 0.4–0.5, 0.5–0.6, 0.6–0.7, 0.7–0.8 and > 0.8) for each of 8 selected regions separately, to calculate seasonal and annual vegetation coverages, along with inter-annual and inter-seasonal vegetation anomalies. The annual and seasonal NDVI anomalies were calculated for each pixel during 2001–2018 as: (2) (3) Where NDVIYA is a yearly anomaliy and NDVISA is a seasonal anomaliy for each pixel in each year, NDVIs is the average NDVI for each season and NDVIy is the average NDVI for each year, NDVIȳ is the average yearly NDVI for the whole study period 2001–2018, NDVIs- is the average seasonal NDVI for the whole study period 2001–2018, and is the standard deviation of NDVI for the whole study period 2001–2018, and is the standard deviation of NDVI for each season for the whole study period 2001–2018. The study used a linear regression to define the correlations significant at a level of 0.05, which were taken into account in further analyses. Linear regression is a statistical method to find the relationship between two variables by fitting a linear model, in which one variable acts as an explanatory variable and the other one is a dependent variable (Song et al., 2005). A linear regression model is defined as: (4) where a and b are the regression coefficients. The b coefficient can be obtained from the given pairs of (xi, yi). In the current study regression model was used for calculating the trend of vegetation (as the dependent variable) during the years (as the independent variable) in Iceland during 2001–2018. 110 NDVI = . Bnir + Bred Bnir – Bred yi = a + bxi , NDVIYA = ,Stdy NDVIy – NDVIȳ NDVISA = .Stds NDVIs – NDVIs- 62-1_acta49-1.qxd 20.6.2022 10:43 Page 110 Acta geographica Slovenica, 62-1, 2022 111 0 10,000 20,000 30,000 40,000 50,000 60,000 70,000 80,000 A re a (k m )2 Time 01 – N o v 16 – O c t 30 – S ep 14 – S ep 29 – A u g 13 – A u g 28 – J u l 12 – J u l 26 – J u n 10 – J u n 25 – M a y 09 – M a y 23 – A p r 07 – A p r 22 – M a r 06 – M a r 18 – F e b 02 – F e b 17 – J a n Average NDVI coverage 3 Results 3.1 NDVI variations In Figure 2 the average vegetation coverage (NDVI > 0.2) in Iceland for the study period is shown. The NDVI is steadily rising from the midst of February, however, it is very slow till the middle of March (12,948 km2). From the end of March, vegetation coverage starts to increase rapidly to reach 33,524 km2 only one month later (at the end of April). The uptrend is continued in the next months and the maximum coverage is reached in the period from the middle of July to late August, having an average value of 66,858 km2. After the end of August, the coverage of green vegetation decreases rapidly. Therefore, it can be stated that the average growing season (GS) in Iceland starts around 23rd March and ends at the end of August (Figure 2). Figure 3 shows the time series of the annual average vegetation coverage in Iceland for the study period. There is an insignificant decreasing trend in the NDVI coverage for the study area in 2001–2018, however sub- stantial interannual variability is also visible. In the years 2003, 2004, and 2017 the maxima in the annual average vegetation coverage were observed (being 43,229, 42,106, and 42,255km2, respectively), while in 2008, 2009, 2013, 2014, and 2015 the minima occurred (33,233, 33,962, 34,645, 34,051 and 34,398km2, respectively). In Figure 4 the similar time series as in Figure 3 are shown, but for individual classes of the NDVI. There is a decreasing trend in the NDVI for the ranges between 0.2–0.5, and an increasing trend in the NDVI from the range of 0.6–1. However, no long-term trend for NDVI for the range 0.5–0.6 was observed (Figure 4 and Table 3). 34,000 32,000 36,000 38,000 40,000 42,000 44,000 46,000 A re a (k m )2 Year 20 18 20 17 20 16 20 15 20 14 20 13 20 12 20 11 20 10 20 09 20 08 20 07 20 06 20 05 20 04 20 03 20 02 20 01 Average NDVI coverage Figure 3: The time series of average annual vegetation coverage (NDVI > 0.2) in Iceland. Figure 2: The average value of the NDVI coverage (> 2) in Iceland for the period 2001–2018. 62-1_acta49-1.qxd 20.6.2022 10:43 Page 111 112 A v e ra g e – 0 .3 N D V I 0 .2 4500 6500 4900 4700 5100 5500 5900 6100 6300 5300 5700 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 A ve ra ge N D V I 0 . 3 – 0 .4 5500 6500 7500 8500 6000 8000 9000 9500 7000 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 A ve ra ge N D V I 0 . 4 – 0 .5 7000 6500 7500 8500 9500 10,500 9 00 0 11, 00 0 6000 8000 10,000 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 A ve ra ge N D V I 0 . 5 – 0 .6 5500 6500 7500 8500 6000 7000 8000 9000 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 Figure 4: The time series of the average annual coverage of the vegetated areas with different ranges of NDVI values in Iceland. (p. 112–113) 62-1_acta49-1.qxd 20.6.2022 10:43 Page 112 Acta geographica Slovenica, 62-1, 2022 113 A v e ra g e N D V I 0 . 6 – 0 .7 4500 5500 6500 7500 5000 6000 7000 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 A ve ra ge N D V I 0 . 7 – 0 .8 2500 3500 4500 2 00 0 3 00 0 4 00 0 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 A ve ra ge N D V I > 0 .8 4 00 600 500 700 900 1100 1300 800 1000 1200 1400 20 10 2002 2003 2004 2005 2006 2007 2008 2009 2010 20 11 2012 2013 2014 2015 2016 2017 2018 Table 3: The trend of different NDVI classes in Iceland in the period 2001–2018. NDVI Categories Correlation Average NDVI 0.2–0.3 –0.46* Average NDVI 0.3–0.4 –0.59* Average NDVI 0.4–0.5 –0.55* Average NDVI 0.5–0.6 0.05 Average NDVI 0.6–0.7 0.43 Average NDVI 0.7–0.8 0.59* Average NDVI 0.8–0.9 0.70* Average NDVI 0.9–1 0.80* Note: * denotes significance at p= 0.05 62-1_acta49-1.qxd 20.6.2022 10:43 Page 113 114 3.2 Regional variability Table 4 presents both the seasonal and annual percentage of area covered by vegetation in the regions of Iceland. The Southwest Peninsuls (SWP) and Rekjavik Area (RA) regions have the highest vegetation cov- erage, with an average annual value of 77.1 and 67.3%, respectively. The Northeast (NE) and East (E) regions have the lowest values, 26.9%, and 31.3%, respectively, followed closely by the South region with an annu- al value of only 33.6%. In winter very low percentage of the area is covered by vegetation, except in the Southwest peninsula. (Table 4 and Figure 5). Table 4: The average seasonal and annual percentage of areas covered by vegetation (NDVI > 0.2) in the regions of Iceland in the period 2001–2018. Period S SWP RA W WF NW NE E Winter 2.2 36.8 4.5 2.8 0.0 0.1 0.2 1.6 Spring 31.5 85.1 78.2 60.5 23.6 40.2 16.7 17.9 Summer 53.6 94.0 95.5 92.1 77.4 80.9 51.3 57.4 Fall 47.1 92.5 90.8 83.5 61.5 66.7 39.3 48.2 Year 33.6 77.1 67.3 59.7 40.6 47.0 26.9 31.3 S E NE W NW WF RA SWP 15°0'0"W20°0'0"W25°0'0"W 66 °0 '0" N 64 °0 '0" N 0 330165 km< 0.2 0.2–0.3 0.3–0.4 0.4–0.5 0.5–0.6 0.6–0.7 0.7–0.8 > 0.8 Content by: Iman Rousta Map by: Iman Rousta Source: Contact to irousta@yazd.ac.ir S E NE W NW WF RA SWP 15°0'0"W20°0'0"W25°0'0"W 67 °0 '0" N 65 °0 '0" N 63 °0 '0" N S E NE W NW WF RA SWP 15°0'0"W20°0'0"W S E NE W NW WF RA SWP 15°0'0"W20°0'0"W25°0'0"W Winter Spring Summer Fall 63 °0 '0" N 66 °0 '0" N 63 °0 '0" N 66 °0 '0" N ® Figure 5: Average seasonal vegetation coverage in Iceland in the period 2001–2018. 62-1_acta49-1.qxd 20.6.2022 10:43 Page 114 In Table 5 information about the variability and trend of the NDVI during the study period for each of 8 regions of Iceland is provided. The highest average maximum NDVI values are found in the West, the South (0.66), and the East (0.63) regions. The lowest average maximum values are observed in the Northeast and Northwest (0.58) and the West Fjords (0.56). A small and insignificant downtrend in the maximum values in the East, and slightly higher, but also insignificant, uptrends in the Northeast, Reykjavik, Northwest, and West were found. A positive and statistically significant trend was found in the average maximum NDVI in the South. The average NDVI values represent mainly the proportion of mountains in the respective region, with the highest values observed for the Northeast, East, and West Fjords. The downtrends of average NDVI were observed for all the regions except the Southwest peninsula and the Northeast. The highest negative trend values were obtained for Reykjavik, the South, and the West regions, however, all of them were statistically non-significant. The standard deviation of the mean NDVI values is the smallest for Reykjavik and the Southwest Peninsula regions (about 0.1 std). The increasing trend of the standard deviation of the NDVI is significant in Reykjavik. Figure 6 shows the maximum annual NDVI for every region of Iceland for the whole study period (2001–2018). It was observed that the highest mean was in the South region (0.75), while the lowest was in the West Fjords (0.63). Table 5: The descriptive NDVI statistics for the regions of Iceland for the period 2001–2018. Name Avg. Max NDVI Avg. Max NDVI Avg. NDVI Avg. NDVI Avg. Std NDVI Avg. Std NDVI Trend Trend Trend South 0.66 0.55* 0.20 –0.26 0.19 0.01 Southwest Peninsula 0.60 –0.08 0.34 0.07 0.10 0.24 Reykjavik Area 0.62 0.19 0.34 –0.21 0.10 0.58* West 0.66 0.19 0.29 –0.22 0.14 0.50 West Fjords 0.56 –0.01 0.19 –0.13 0.14 0.00 Northwest 0.58 0.29 0.23 –0.14 0.15 0.14 Northeast 0.58 0.44 0.13 0.04 0.17 0.07 East 0.63 –0.14 0.15 –0.15 0.17 0.04 Note: * denotes significance at p= 0.05 Acta geographica Slovenica, 62-1, 2022 115 E NE NW WF W RA SWP S R eg io n n am e 0.55 0.60 0.65 0.70 0.75 NDVI Figure 6: Maximum annual NDVI for the regions of Iceland for the period 2001–2018. 62-1_acta49-1.qxd 20.6.2022 10:43 Page 115 116 Assessment of the vegetation coverage trends in Iceland (Table 6) indicates that a reduction in NDVI in all regions was occurring in the winter season, however, it was significant only in Reykjavik (–0.5). It is worth noting that rather high values were observed also for the Northwest, West, and West Fjords regions (–0.4), but they were insignificant. In the spring, a negative non-significant trend occurred in most regions. In the summer, a strong and significant average NDVI uptrend can be noticed in the Southwest penin- sula, while in other regions much weaker and non-significant trends occurred. In the fall, positive non-significant trends in the West, Northwest, and the West fjords can be noticed, while in other regions only minor changes occurred. Table 6: Trends of the seasonal average NDVI values in regions of Iceland in the period 2001–2018. Season S SWP RA W WF NW NE E Winter –0.3 –0.3 –0.5* –0.4 –0.4 –0.4 –0.2 –0.2 Spring –0.3 –0.1 –0.3 –0.2 –0.2 –0.2 0.0 –0.2 Summer 0.2 0.5* –0.3 –0.1 –0.1 0.0 0.2 –0.1 Fall 0.1 0.2 –0.1 0.2 0.3 0.0 0.1 0.1 Note: * denotes significance at p= 0.05 Table 7 shows the percentage of the area with decreasing (change from –0.12 to –0.007 NDVI/year), stable (from –0.007 to 0.005 NDVI/year), or increasing (from 0.005 to 0.08 NDVI/year) values of the NDVI index for each region in Iceland during the study period. In all regions, no constant change of vegetation was observed on more than 80% of the lands. The WF and E regions had the highest share of land with a decrease of the vegetation (2.8 and 1.8%, respectively), while the remaining regions had smaller fractions of the area, on which vegetation decreased. At the same time, the WF and E were also the regions having the largest share of the area with increasing vegetation (5.6 and 3.7%, respectively), while in the remain- ing regions vegetation index rose on about 2–3% of the area. In general, an increase in the NDVI index was observed on about 3260 km2, whereas a decrease on 1635 km2 (Table 7 and Figure 7). Table 7: The percentage of the area with decreasing, stable, or increasing vegetation of each region in Iceland during 2001–2018. Trend (NDVI/year) S SWP RA W WF NW NE E from –0.12 to –0.02 (Strong decrease) 0.2 0.0 0.1 0.1 0.2 0.1 0.1 0.2 from –0.02 to –0.007 (Decrease) 1.5 0.5 1.2 1.3 2.6 1.0 1.0 1.6 from –0.007 to –0.002 (No change) 8.1 3.3 7.3 6.6 12.3 5.6 6.9 9.4 from –0.002 to 0.005 (No change) 87.4 94.9 88.2 88.9 79.3 90.8 89.6 85.0 from 0.005 to 0.08 (Increase) 2.8 1.3 3.2 3.1 5.6 2.4 2.5 3.7 3.3 Location of the areas with the most extreme changes Several areas with relatively large changes in the annual average of the NDVI index have been detected. Such areas appear as small clusters of pixels with a uniform color. In Figure 7 three additional locations relative to pre-2013 locations indicated by Raynolds et al. (2015) were marked by the boxes. The first one is located along the shorelines of artificial lake Hálslón in eastern Iceland (Figure 7, box a). The lake began to be filled in the fall of 2006 and due to variability of the lake water level, a reduction in vegetation along its shoreline was observed (Aradóttir et al. 2013). A second location is the extended riverbed of Skaftá River (Figure 7, box b). The river flooded in the fall of 2015 and for the second time in the summer of 2018, causing a substantial level of destruction of vegetation (Óskarsdóttir 2016). The third one is a landslide in Hítardalur, in western Iceland, which is covering about 2 km2, including a small lake that was formed due to a dammed river (Dabiri et al. 2019). Although the landslide fell as late as at the end of the study period, in July 2018, it appeared in Figure 7 (box c) as a cluster of pixels with a downtrend. 62-1_acta49-1.qxd 20.6.2022 10:43 Page 116 4 Discussion Not many studies on the regional, seasonal and annual spatial vegetation dynamics of Iceland were under- taken up to now. Our study is the first study that has assessed the vegetation variation in the whole of Iceland and its sub-regions using remote sensing. In the recent works dealing with the vegetation dynamics over the studied area, only a part of the country was analyzed (Lemenkova 2020b; Bates et al. 2021). In the case of study made for the whole country, vegetation variation and the inter-regional differences were not studied, and additionally, they employed satellite images having a lower spatial resolution, and analyzes were done for a shorter period (Raynolds et al. 2015). The satellite images used by Raynolds et al. (2015) were not trustable, and the next version (V006), which was unavailable at that time, could solve the uncertainties. Raynolds et al. (2015, 9495) said: »Version 6 of the MODIS Vegetation Indices products will correct this problem, but is not yet available«. It occurred that the MODIS NDVI images are useful for monitoring vegetation dynam- ics in Iceland and the method used in the study can be used for other regions and countries. Acta geographica Slovenica, 62-1, 2022 117 (a) (b) (c) S E NE W NW WF RA SWP 15°0'0"W20°0'0"W25°0'0"W 67 °0 '0" N 67 °0 '0" N 66 °0 '0" N 66 °0 '0" N 65 °0 '0" N 65 °0 '0" N 64 °0 '0" N 64 °0 '0" N 63 °0 '0" N 63 °0 '0" N ® –0.12––0.02 (Very Decreasing) –0.02––0.007 (Decreasing) –0.007––0.002 (Normal) –0.002–0.005 (Normal) 0.005–0.08 (Increasing) Glacier 0 7537.5 km Content by: Iman Rousta Map by: Iman Rousta Source: Contact to irousta@yazd.ac.ir Figure 7: The map of the values of the NDVI index trends (NDVI/year) in Iceland during 2001–2018. (a) shorelines of artificial lake Hálslón, (b) Skaftá River, (c) landslide in Hítardalur. 62-1_acta49-1.qxd 20.6.2022 10:43 Page 117 The results of this study indicated that the maximum vegetation coverage in Iceland occurs from the middle of July to late August, with an average area of 66,858 km2 covered in vegetation (NDVI ≥ 0.2), which is about 65% of the country’s territory. Also, the results showed that the growing season starts on average from late March and lasts up to the end of August. Raynolds et al. (2015) in their study showed that a reduction in NDVI occurred in Iceland in the peri- od from 2002 to 2013. The present study dealt with the period 2001–2018 and the results are in line with the results of the Raynold’s team in the case of moderate decreasing trend and the high spatial variabili- ty of trends for the different regions of Iceland. In the study, significant downtrends for the lower values of the NDVI were found, whereas, for high- er values, positive uptrends were spotted. The reason for this effect is probably due to the growth of both planted, as well as natural forests, which again is in line with results presented in Raynolds et al. (2015). The total area of forests in Iceland has indeed increased from a total of 298 km2 in the year 2000 to ~514 km2 in 2020 (estimated value), whereas the area of other wooded lands increased from a total of 1344 km2 in 2000 to ~1495 km2 in 2020 (estimated value) (MacDicken 2015). A forest growing on grassland will most likely increase the NDVI for that area. As some of the new forests in Iceland do indeed grow on grass- lands, they may simultaneously contribute to a reduction of the areas with a low value of the NDVI and an increase of the areas with higher NDVI value. The high positive trend in the maximum annual NDVI index value in the southern part of Iceland is in line with the widespread afforestation of that region (Fries 2017). The other studies about northern latitude’s vegetation variations, indicate that there is a great increase in vegetation coverage between 45°N and 70°N (Myneni et al. 1997; Huang et al. 2017). Substantial interannual variability in all of the analyzed classes of the NDVI values was found. This observation requires to be explored further in connection with the actual state of the atmosphere, espe- cially with temperatures and precipitation. The percentage of land covered by vegetation is by far the highest in the Southwest peninsula and the Reykjavik area (these are the smallest regions of the 8 main regions of Iceland). This is primarily due to the fact that both regions are located mostly at terrain having lower elevation. On the other hand, field- work showed that the high-density vegetation of that area is established after reclamation through the use of inorganic fertilizers and birch trees, which is consistent with the result of the Nyirenda (Nyirenda 2020) as well as the Reykjavik area has the most dense site of reforestation in the whole country (Fries 2017). In three regions, the West fjords, the Northeast, and the South, anomalously low vegetation coverage was spot- ted. The West fjords had much less vegetation coverage than the West even though their mean elevation is similar. Similarly, the Northeast had less vegetation coverage than the East, and the South had less veg- etation than the Northwest, although their mean elevations are similar. The explanation for the West fjords having less vegetation than the West lies undoubtedly in a relatively cold climate in the West fjords (Einarsson 1984). As for the other regions with low-vegetation coverage, the explanation lies in the fact that the South and the Northeast are located in volcanically active zones, in which every few years eruptions occur, which in turn results in large areas with progressing erosion, even at relatively low altitudes (Arnalds, Ólafsson and Dagsson-Waldhauserova 2014; Lemenkova 2020b). Only minor changes in the vegetation coverage in Iceland during the study period were observed, with the trend of changes significantly smaller than the inter-annual variability, which in turn must be considered as substantial. The total vegetation coverage was almost 67,000 km2 (about 65% of the Iceland area), which is quite high for an island located in high latitudes, with a high elevation and a high share of land covered by glaciers and lakes. 5 Conclusion Climate change and global warming can affect all areas around the world. It is predicted that the impact of climate change will be more severe for regions with higher latitudes (King et al. 2018; Anderson, Bayer and Edwards 2020), such as Iceland. Vegetation variations are surface phenomena, which could be a good indi- cator of climate change. The present study attempted to identify and analyze the spatio-temporal variations of the NDVI in both Iceland and its 8 sub-regions using remotely-sensed satellite images. It was found that MODIS NDVI images can be useful for monitoring the vegetation variations in the studied area. 118 62-1_acta49-1.qxd 20.6.2022 10:43 Page 118 The NDVI retrieved from remote sensing for the period 2001–2018 showed considerable inter-annu- al variability and a minor, statistically insignificant, decrease in the vegetation coverage in Iceland. In the late summer, 65% of Iceland is covered with vegetation (NDVI > 0.2). The areas with low NDVI (NDVI < 0.6) had a decreasing trend, whereas for the areas with high NDVI (NDVI > 0.6) an increase in the vegetation coverage was observed. In general, an increase in the NDVI was observed in the area of about 3260 km2, whereas a decrease in the NDVI index was observed in the area of 1635 km2, which is in line with an increase in forested areas in Iceland. Vegetation coverage is influenced by such atmospheric parameters as temperature, precipitation, rel- ative humidity, etc. On the other hand, all of these atmospheric parameters are influenced by atmospheric patterns and teleconnections. Therefore, to have a more accurate assessment of inter-seasonal, inter-annu- al, and inter-regional vegetation variations in the studied area, the atmosphere dynamics should be studied simultaneously with the vegetation variations. Up to now, one paper, in which the relationships between the atmospheric circulations and teleconnections with vegetation dynamics of Iceland (Olafsson and Rousta 2021) has already been published. ACKNOWLEDGMENTS: Iman Rousta is deeply grateful to his supervisor (Haraldur Olafsson, Professor of Atmospheric Sciences, Institute for Atmospheric Sciences-Weather and Climate, and Department of Physics, University of Iceland, and Icelandic Meteorological Office (IMO)), for his great support, kind guid- ance, and encouragement. Author Contributions: H. O. and I. R. proposed the topic, commanded the data processing, analysis, and wrote the manuscript. Funding Statement: This work was supported by Vedurfelagid, Rannis and Rannsoknastofa i vedurfraedi. Conflicts of Interest: The authors declare no con- flict of interest. 6 References Alcaraz‐Segura, D., Chuvieco, E., Epstein, H. E., Kasischke, E. S., Trishchenko, A. 2010: Debating the greening vs. browning of the North American boreal forest: Differences between satellite datasets. Global Change Biology 16-2. DOI: https://doi.org/10.1111/j.1365-2486.2009.01956.x Anderson, R., Bayer, P. E., Edwards, D. 2020: Climate change and the need for agricultural adaptation. Current Opinion in Plant Biology 56. DOI: https://doi.org/10.1016/j.pbi.2019.12.006 Aradóttir, Á. L., Petursdottir, T., Halldorsson, G., Svavarsdottir, K., Arnalds, O. 2013: Drivers of ecologi- cal restoration: lessons from a  century of restoration in Iceland. Ecology and Society 18-4. DOI: https://doi.org/10.5751/ES-05946-180433 Arnalds, Ó. 2008: The soils of Iceland. World Soils Book Series. Dortrecht. DOI: https://doi.org/10.1007/978- 94-017-9621-7 Arnalds, Ó., Ellin Fjola, T., Sigmar, M., Asgeir, J., Arnor, A. 2001: Soil erosion in Iceland. Reykjavík. Arnalds, O., Hallmark, C. T., Wilding, L. P. 1995: Andisols from four different regions of Iceland. Soil Science Society of America Journal 59-1. DOI: https://doi.org/10.2136/sssaj1995.03615995005900010025x Arnalds, O., Ólafsson, H., Dagsson-Waldhauserova, P. 2014: Quantification of iron-rich volcanogenic dust emissions and deposition over the ocean from Icelandic dust sources. Biogeosciences 11. DOI: https://doi.org/10.5194/bg-11-6623-2014 Atasoy, M. 2018: Monitoring the urban green spaces and landscape fragmentation using remote sensing: A case study in Osmaniye, Turkey. Environmental Monitoring and Assessment 190. DOI: https://doi.org/ 10.1007/s10661-018-7109-1 Bates, R., Erlendsson, E., Eddudóttir, S. D., Möckel, S. C., Tinganelli, L., Gísladóttir, G. 2021: Landnám, land use and landscape change at Kagaðarhóll in Northwest Iceland. Environmental Archaeology 27-2. DOI: https://doi.org/10.1080/14614103.2021.1949680 Bhatt, U. S., Walker, D. A., Raynolds, M. K., Bieniek, P. A., Epstein, H. E., Comiso, J. C., Pinzon, J. E., et al. 2013: Recent declines in warming and vegetation greening trends over pan-Arctic tundra. Remote Sensing 5-9. DOI: https://doi.org/10.3390/rs5094229 Bjerke, J. W., Karlsen, S. R., Høgda, K. A., Malnes, E., Jepsen, J. U., Lovibond, S., Vikhamar-Schuler, D., et al. 2014: Record-low primary productivity and high plant damage in the Nordic Arctic Region in 2012 caused by multiple weather events and pest outbreaks. Environmental Research Letters 9-8. DOI: https://doi.org/10.1088/1748-9326/9/8/084006 Acta geographica Slovenica, 62-1, 2022 119 62-1_acta49-1.qxd 20.6.2022 10:43 Page 119 Bokhorst, S. F., Bjerke, J. W., Tømmervik, H., Callaghan, T. V., Phoenix, G. K. 2009: Winter warming events damage sub‐Arctic vegetation: consistent evidence from an experimental manipulation and a natural event. Journal of Ecology 97-6. DOI: https://doi.org/10.1111/j.1365-2745.2009.01554.x Chao, L., Zhang, K., Wang, J., Feng, J., Zhang, M. 2021: A comprehensive evaluation of five evapotran- spiration datasets based on ground and GRACE satellite observations: Implications for improvement of evapotranspiration retrieval algorithm. Remote Sensing 13-12. DOI: https://doi.org/10.3390/rs13122414 Chen, D., Brutsaert, W. 1998: Satellite-sensed distribution and spatial patterns of vegetation parameters over a tallgrass prairie. Journal of the Atmospheric Sciences 55-7. DOI: https://doi.org/10.1175/1520- 0469(1998)055<1225:SSDASP>2.0.CO;2 Cui, L., Shi, J., Yang, Y., Fan, W. 2009: Ten-day response of vegetation NDVI to the variations of temper- ature and precipitation in eastern China. Acta Geographica Sinica 64-7. DOI: https://doi.org/10.11821/ xb200907009 Dabiri, Z., Hölbling, D., Abad, L., Prasicek, G., Argentin, A.-L., Tsai, T.-T. 2019: An Object-Based Approach for Monitoring the Evolution of Landslide-Dammed Lakes and Detecting Triggering Landslides in Taiwan. ISPRS-International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences 4238. DOI: https://doi.org/10.5194/isprs-archives-XLII-3-W8-103-2019 Dabrowska-Zielinska, K., Kogan, F., Ciolkosz, A., Gruszczynska, M., Kowalik, W. 2002: Modelling of crop growth conditions and crop yield in Poland using AVHRR-based indices. International Journal of Remote Sensing 23-6. DOI: https://doi.org/10.1080/01431160110070744 Didan, K. 2015a: MOD13Q1 MODIS/Terra vegetation indices 16-day L3 global 250 m SIN grid V006. NASA EOSDIS Land Processes DAAC. DOI: https://doi.org/10.5067/MODIS/MOD13Q1.006 Didan, K., Munoz, A. B., Solano, R., Huete, A. 2015: MODIS vegetation index user’s guide (MOD13 series). Internet: https://vip.arizona.edu/documents/MODIS/MODIS_VI_UsersGuide_June_2015_C6.pdf (6. 5. 2022). Dugmore, A. J., Gísladóttir, G., Simpson, I. A., Newton, A. 2009: Conceptual models of 1200 years of Icelandic soil erosion reconstructed using tephrochronology. Journal of the North Atlantic 2-1. DOI: https://doi.org/ 10.3721/037.002.0103 Dutta, D., Kundu, A., Patel, N. R., Saha, S. K., Siddiqui, A. R. 2015: Assessment of agricultural drought in Rajasthan (India) using remote sensing derived Vegetation Condition Index (VCI) and Standardized Precipitation Index (SPI). The Egyptian Journal of Remote Sensing and Space Science 18-1. DOI: https://doi.org/10.1016/j.ejrs.2015.03.006 Dye, D. G., Tucker, C. J. 2003: Seasonality and trends of snow‐cover, vegetation index, and temperature in northern Eurasia. Geophysical Research Letters 30-7. DOI: https://doi.org/10.1029/2002GL016384 Einarsson, M. A., 1984: Climate of iceland. Climates of the Oceans. World Survey of Climatology. Amsterdam. Elmendorf, S. C., Henry, G. H., Hollister, R. D., Björk, R. G., Boulanger-Lapointe, N., Cooper, E. J., Cornelissen, J. H., et al. 2012: Plot-scale evidence of tundra vegetation change and links to recent summer warming. Nature Climate Change 2. DOI: https://doi.org/10.1038/nclimate1465 Epstein, H. E., Raynolds, M. K., Walker, D. A., Bhatt, U. S., Tucker, C. J., Pinzon, J. E. 2012: Dynamics of aboveground phytomass of the circumpolar Arctic tundra during the past three decades. Environmental Research Letters 7-1. DOI: https://doi.org/10.1088/1748-9326/7/1/015506 Foley, J. A., Levis, S., Costa, M. H., Cramer, W., Pollard, D. 2000: Incorporating dynamic vegetation cover within global climate models. Ecological Applications 10-6. DOI: https://doi.org/10.2307/2641227 Fries, E., 2017: Reforestation in the far north. Umeå. Gandhi, G. M., Parthiban, S., Thummalu, N., Christy, A. 2015: NDVI: Vegetation change detection using remote sensing and GIS–a  case study of Vellore District. Procedia Computer Science 57. DOI: https://doi.org/10.1016/j.procs.2015.07.415 Gao, B.-C., Goetzt, A. F. H. 1995: Retrieval of equivalent water thickness and information related to bio- chemical components of vegetation canopies from AVIRIS data. Remote Sensing of Environment 52-3. DOI: https://doi.org/10.1016/0034-4257(95)00039-4 Geerken, R., Zaitchik, B., Evans, J. 2005: Classifying rangeland vegetation type and coverage from NDVI time series using Fourier Filtered Cycle Similarity. International Journal of Remote Sensing 26-24. DOI: https://doi.org/10.1080/01431160500300297 120 62-1_acta49-1.qxd 20.6.2022 10:43 Page 120 Ghafarian Malamiri, H. R., Zare, H., Rousta, I., Olafsson, H., Izquierdo Verdiguier, E., Zhang, H., Mushore, T. D. 2020: Comparison of Harmonic Analysis of Time Series (HANTS) and Multi-Singular Spectrum Analysis (M-SSA) in reconstruction of long-gap missing data in NDVI time series. Remote Sensing 12-17. DOI: https://doi.org/10.3390/rs12172747 Gitelson, A. A., Viña, A., Arkebauer, T. J., Rundquist, D. C., Keydan, G., Leavitt, B. 2003: Remote estimation of leaf area index and green leaf biomass in maize canopies. Geophysical Research Letters 30. DOI: https://doi.org/10.1029/2002GL016450 Huang, M., Piao, S., Janssens, I. A., Zhu, Z., Wang, T., Wu, D., Ciais, P., et al. 2017: Velocity of change in vegetation productivity over northern high latitudes. Nature Ecology and Evolution 1. DOI: https://doi.org/10.1038/s41559-017-0328-y James, P., Chester, D., Duncan, A. 2000: Volcanic soils: Their nature and significance for archaeology. Geological Society, London, Special Publications, 171. DOI: https://doi.org/10.1144/GSL.SP.2000.171.01.23 Ji, L., Peters, A. J. 2003: Assessing vegetation response to drought in the northern Great Plains using veg- etation and drought indices. Remote Sensing of Environment 87-1. DOI: https://doi.org/10.1016/ S0034-4257(03)00174-3 Jóhannesdóttir, L., Alves, J. A., Gill, J. A., Gunnarsson, T. G. 2017: Reconciling biodiversity conservation and agricultural expansion in the sub-arctic environment of Iceland. Ecology and Society 22-1. DOI: https://doi.org/10.5751/ES-08956-220116 Jóhannesson, T. 2010: Agriculture in Iceland: Conditions and characteristics. Internet: https://rafhladan.is/ handle/10802/9353 (11. 5. 2022). Jovanović, M. M., Milanović, M. M., Zorn, M. 2018: The use of NDVI and CORINE Land Cover data- bases for forest management in Serbia. Acta geographica Slovenica 58-1. DOI: https://doi.org/10.3986/ AGS.818 King, M., Altdorff, D., Li, P., Galagedara, L., Holden, J., Unc, A. 2018: Northward shift of the agricultural climate zone under 21st-century global climate change. Scientific Reports 8. DOI: https://doi.org/ 10.1038/s41598-018-26321-8 Kogan, F. 1991: Observations of the 1990 US drought from the NOAA-11 polar-orbiting satellite. Drought Network News 3. Kogan, F. N. 1995: Droughts of the late 1980s in the United States as derived from NOAA polar-orbiting satellite data. Bulletin of the American Meteorological Society 76-5. DOI: https://doi.org/10.1175/1520- 0477(1995)076<0655:DOTLIT>2.0.CO;2 Lemenkova, P. 2020a: Hyperspectral vegetation indices calculated by Qgis using Landsat Tm image: A case study of Northern Iceland. Advanced Research in Life Sciences 4. DOI: https://doi.org/10.2478/ arls-2020-0021 Lemenkova, P. 2020b: SAGA GIS for information extraction on presence and conditions of vegetation of northern coast of Iceland based on the Landsat TM. Acta Biologica Marisiensis 3-2. DOI: https://doi.org/ 10.2478/abmj-2020-0007 Li, J., Charles, L. S., Yang, Z., Du, G., Fu, S. 2022: Differential mechanisms drive species loss under artifi- cial shade and fertilization in the alpine meadow of the Tibetan Plateau. Frontiers in Plant Science 13. DOI: https://doi.org/10.3389/fpls.2022.832473 Li, W., Shi, Y., Zhu, D., Wang, W., Liu, H., Li, J., Shi, N., et al. 2021: Fine root biomass and morphology in a temperate forest are influenced more by the nitrogen treatment approach than the rate. Ecological Indicators 130. DOI: https://doi.org/10.1016/j.ecolind.2021.108031 Liu, Y., Li, Y., Li, S., Motesharrei, S. 2015: Spatial and temporal patterns of global NDVI trends: Correlations with climate and human factors. Remote Sensing 7-10. DOI: https://doi.org/10.3390/rs71013233 Loveland, T. R., Zhu, Z., Ohlen, D. O., Brown, J. F., Reed, B. C., Yang, L. 1999: An analysis of the IGBP global land-cover characterization process. Photogrammetric Engineering and Remote Sensing 65-9. MacDicken, K. G. 2015: Global forest resources assessment 2015: what, why and how? Forest Ecology and Management 352. DOI: https://doi.org/10.1016/j.foreco.2015.02.006 Mansourmoghaddam, M., Ghafarian Malamiri, H. R., Rousta, I., Olafsson, H., Zhang, H. 2022: Assessment of Palm Jumeirah Island’s construction effects on the surrounding water quality and surface temperatures during 2001–2020. Water 14-4. DOI: https://doi.org/10.3390/w14040634 Acta geographica Slovenica, 62-1, 2022 121 62-1_acta49-1.qxd 20.6.2022 10:43 Page 121 Mansourmoghaddam, M., Rousta, I., Zamani, M., Mokhtari, M. H., Karimi Firozjaei, M., Alavipanah, S. K. 2021: Study and prediction of land surface temperature changes of Yazd city: Assessing the proximity and changes of land cover. Journal of RS and GIS for Natural Resources 12-4. Martínez, B., Gilabert, M. A. 2009: Vegetation dynamics from NDVI time series analysis using the wavelet transform. Remote sensing of environment 113-9. DOI: https://doi.org/10.1016/j.rse.2009.04.016 Masson-Delmotte, V., Zhai, P., Pörtner, H.-O., Roberts, D., Skea, J., Shukla, P. R., Pirani, A., et al. (eds.) 2018: Global Warming of 1.5 °C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. Internet: https://www.ipcc.ch/site/assets/uploads/sites/2/2019/06/SR15_Full_ Report_High_Res.pdf (6. 5. 2022). McVicar, T. R., Bierwirth, P. N. 2001: Rapidly assessing the 1997 drought in Papua New Guinea using com- posite AVHRR imagery. International Journal of Remote Sensing 22-11. DOI: https://doi.org/10.1080/ 01431160120728 Merrington, A. T., 2019: A time series analysis of vegetation succession on lava flow fields at Hekla Volcano: Assessing the utility of Landsat data. M.Sc. thesis, University of Iceland. Reykjavik. Metternicht, G., Zinck, J. A., Blanco, P. D., del Valle, H. F. 2010: Remote sensing of land degradation: Experiences from Latin America and the Caribbean. Journal of Environmental Quality 39-1. DOI: https://doi.org/10.2134/jeq2009.0127 Miao, R., Liu, Y., Wu, L., Wang, D., Liu, Y., Miao, Y., Yang, Z., et al. 2022: Effects of long-term grazing exclu- sion on plant and soil properties vary with position in dune systems in the Horqin Sandy Land. Catena 209-2. DOI: https://doi.org/10.1016/j.catena.2021.105860 Miao, R., Qiu, X., Guo, M., Musa, A., Jiang, D. 2018: Accuracy of space-for-time substitution for vegeta- tion state prediction following shrub restoration. Journal of Plant Ecology 11-2. DOI: https://doi.org/ 10.1093/jpe/rtw133 Moniruzzaman, M., Thakur, P. K., Kumar, P., Ashraful Alam, M., Garg, V., Rousta, I., Olafsson, H. 2021: Decadal urban land use/land cover changes and its impact on surface runoff potential for the Dhaka city and surroundings using remote sensing. Remote Sensing 13-1. DOI: https://doi.org/10.3390/ rs13010083 Montandon, L. M., Small, E. E. 2008: The impact of soil reflectance on the quantification of the green veg- etation fraction from NDVI. Remote Sensing of Environment 112-4. DOI: https://doi.org/10.1016/ j.rse.2007.09.007 Moulin, S., Kergoat, L., Viovy, N., Dedieu, G. 1997: Global-scale assessment of vegetation phenology using NOAA/AVHRR satellite measurements. Journal of Climate 10-6. DOI: https://doi.org/10.1175/1520- 0442(1997)010<1154:GSAOVP>2.0.CO;2 Myneni, R. B., Keeling, C. D., Tucker, C. J., Asrar, G., Nemani, R. R. 1997: Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 386. DOI: https://doi.org/10.1038/386698a0 Nyirenda, H. 2020: Soil carbon status after vegetation restoration in South West Iceland. Heliyon 6-10. DOI: https://doi.org/10.1016/j.heliyon.2020.e05254 Olafsson, H., Rousta, I. 2021: Influence of atmospheric patterns and North Atlantic Oscillation (NAO) on vegetation dynamics in Iceland using Remote Sensing. European Journal of Remote Sensing 54-1. DOI: https://doi.org/10.1080/22797254.2021.1931462 Ollinger, S. V. 2011: Sources of variability in canopy reflectance and the convergent properties of plants. New Phytologist 189-2. DOI: https://doi.org/10.1111/j.1469-8137.2010.03536.x Óskarsdóttir, G. 2016: Gróðurvöktun í Kringilsárrana: samanburður á samsetningu og þekju gróðurs árin 2006 og 2015. Internet: https://www.sjalfbaerni.is/static/research/files/lv-2016-064-grodurvoktun-i- kringilsarrana.pdf (6. 5. 2022). Óskarsson, H., Arnalds, Ó., Gudmundsson, J., Gudbergsson, G. 2004: Organic carbon in Icelandic Andosols: Geographical variation and impact of erosion. Catena 56-1,2,3. DOI: https://doi.org/10.1016/ j.catena.2003.10.013 Pettorelli, N., Vik, J. O., Mysterud, A., Gaillard, J.-M., Tucker, C. J., Stenseth, N. C. 2005: Using the satel- lite-derived NDVI to assess ecological responses to environmental change. Trends in ecology and evolution 20-9. DOI: https://doi.org/10.1016/j.tree.2005.05.011 122 62-1_acta49-1.qxd 20.6.2022 10:43 Page 122 Polyakov, I. V., Bekryaev, R. V., Alekseev, G. V., Bhatt, U. S., Colony, R. L., Johnson, M. A., Maskshtas, A. P., et al. 2003: Variability and trends of air temperature and pressure in the maritime Arctic, 1875–2000. Journal of Climate 16-12. DOI: https://doi.org/10.1175/1520-0442(2003)016<2067:VATOAT>2.0.CO;2 Raynolds, M., Magnússon, B., Metúsalemsson, S., Magnússon, S. H. 2015: Warming, sheep and volcanoes: Land cover changes in Iceland evident in satellite NDVI trends. Remote Sensing 7-8. DOI: https://doi.org/ 10.3390/rs70809492 Rouse, J. W., Haas, R. W., Schell, J. A., Deering, D. W., Harlan, J. C. 1974: Monitoring the vernal advance- ment and retrogradation (green wawe effect) of natural vegetation. NASA/GSFCT Type III Final report. Rousta, I., Javadizadeh, F., Dargahian, F., Olafsson, H., Shiri-Karimvandi, A., Vahedinejad, S. H., Doostkamian, M., et al. 2018: Investigation of vorticity during prevalent winter precipitation in Iran. Advances in Meteorology 2018. DOI: https://doi.org/10.1155/2018/6941501 Rousta, I., Karampour, M., Doostkamian, M., Olafsson, H., Zhang, H., Mushore, T. D., et al. 2020a: Synoptic- dynamic analysis of extreme precipitation in Karoun River Basin, Iran. Arabian Journal of Geosciences 13. DOI: https://doi.org/10.1007/s12517-020-5101-x Rousta, I., Olafsson, H., Moniruzzaman, M., Ardö, J., Zhang, H., Mushore, T. D., Shahin, S., et al. 2020b: The 2000–2017 drought risk assessment of the western and southwestern basins in Iran. Modeling Earth Systems and Environment 6. DOI: https://doi.org/10.1007/s40808-020-00751-8 Rousta, I., Olafsson, H., Moniruzzaman, M., Zhang, H., Liou, Y.-A., Mushore, T. D., Gupta, A. 2020c: Impacts of drought on vegetation assessed by vegetation indices and meteorological factors in Afghanistan. Remote Sensing 12-15. DOI: https://doi.org/10.3390/rs12152433 Rousta, I., Olafsson, H., Nasserzadeh, M. H., Zhang, H., Krzyszczak, J., Baranowski, P. 2021: Dynamics of daytime land surface temperature (LST) variabilities in the Middle East countries during 2001–2018. Pure and Applied Geophysics 178. DOI: https://doi.org/10.1007/s00024-021-02765-4 Running, S. W., Loveland, T. R., Pierce, L. L., Nemani, R. R., Hunt Jr, E. R. 1995: A remote sensing based vegetation classification logic for global land cover analysis. Remote Sensing of Environment 51-1. DOI: https://doi.org/10.1016/0034-4257(94)00063-S Shen, X., Liu, B., Jiang, M., Lu, X. 2020: Marshland loss warms local land surface temperature in China. Geophysical research letters 47-6. DOI: https://doi.org/10.1029/2020GL087648 Sigurmundsson, F. S., Gísladóttir, G., Óskarsson, H. 2014: Decline of birch woodland cover in Þjórsárdalur Iceland from 1587 to 1938. Human ecology 42. DOI: https://doi.org/10.1007/s10745-014-9670-8 Slayback, D. A., Pinzon, J. E., Los, S. O., Tucker, C. J. 2003: Northern hemisphere photosynthetic trends 1982–99. Global Change Biology 9-1. DOI: https://doi.org/10.1046/j.1365-2486.2003.00507.x Symeonakis, E., Drake, N. 2004: Monitoring desertification and land degradation over sub-Saharan Africa. International Journal of Remote Sensing 25-3. DOI: https://doi.org/10.1080/0143116031000095998 Tarpley, J. D., Schneider, S. R., Money, R. L. 1984: Global vegetation indices from the NOAA-7 meteoro- logical satellite. Journal of Climate and Applied Meteorology 23-3. Thenkabail, P. S., Gamage, M. S. D. N., Smakhtin, V. U. 2004: The use of remote sensing data for drought assessment and monitoring in Southwest Asia. Colombo. Thorarinsson, S. 1967: Some problems of volcanism in Iceland. Geologische Rundschau 57. DOI: https://doi.org/10.1007/BF01825713 Thorsteinsson, I., Olafsson, G., Van Dyne, G. M. 1971: Range resources of Iceland. Journal of Range Management 24-2. Tucker, C. J. 1979: Red and photographic infrared linear combinations for monitoring vegetation. Remote Sensing of Environment 8-2. DOI: https://doi.org/10.1016/0034-4257(79)90013-0 Wan, Z., Wang, P., Li, X. 2004: Using MODIS land surface temperature and normalized difference vege- tation index products for monitoring drought in the southern Great Plains, USA. International Journal of Remote Sensing 25-1. DOI: https://doi.org/10.1080/0143116031000115328 Wang, P., Wang, L., Leung, H., Zhang, G. 2020: Super-resolution mapping based on spatial–spectral correlation for spectral imagery. IEEE Transactions on Geoscience and Remote Sensing 59-3. DOI: https://doi.org/10.1109/TGRS.2020.3004353 Yang, L., Wylie, B. K., Tieszen, L. L., Reed, B. C. 1998: An analysis of relationships among climate forcing and time-integrated NDVI of grasslands over the US northern and central Great Plains. Remote Sensing of Environment 65-1. DOI: https://doi.org/10.1016/S0034-4257(98)00012-1 Acta geographica Slovenica, 62-1, 2022 123 62-1_acta49-1.qxd 20.6.2022 10:43 Page 123 Zandbergen, P. 2008: Applications of shuttle radar topography mission elevation data. Geography Compass 2-5. DOI: https://doi.org/10.1111/j.1749-8198.2008.00154.x Zhang, K., Ali, A., Antonarakis, A., Moghaddam, M., Saatchi, S., Tabatabaeenejad, A., Chen, R., et al. 2019a: The sensitivity of North American terrestrial carbon fluxes to spatial and temporal variation in soil moisture: An analysis using radar‐derived estimates of root‐zone soil moisture. Journal of Geophysical Research: Biogeosciences 124-11. DOI: https://doi.org/10.1029/2018JG004589 Zhang, K., Chao, L.-j., Wang, Q.-q., Huang, Y.-c., Liu, R.-h., Hong, Y., Tu, Y., et al. 2019b: Using multi- satellite microwave remote sensing observations for retrieval of daily surface soil moisture across China. Water Science and Engineering 12-2. DOI: https://doi.org/10.1016/j.wse.2019.06.001 Zhao, T., Shi, J., Entekhabi, D., Jackson, T. J., Hu, L., Peng, Z., Yao, P., et al. 2021a: Retrievals of soil mois- ture and vegetation optical depth using a multi-channel collaborative algorithm. Remote Sensing of Environment 257. DOI: https://doi.org/10.1016/j.rse.2021.112321 Zhao, T., Shi, J., Lv, L., Xu, H., Chen, D., Cui, Q., Jackson, T. J., et al. 2020: Soil moisture experiment in the Luan River supporting new satellite mission opportunities. Remote Sensing of Environment 240. DOI: https://doi.org/10.1016/j.rse.2020.111680 Zhao, X., Xia, H., Pan, L., Song, H., Niu, W., Wang, R., Li, R., et al. 2021b: Drought monitoring over Yellow River basin from 2003–2019 using reconstructed MODIS land surface temperature in Google Earth Engine. Remote Sensing 13-18. DOI: https://doi.org/10.3390/rs13183748 Zuo, Y., Jiang, S., Wu, S., Xu, W., Zhang, J., Feng, R., Yang, M., et al. 2020: Terrestrial heat flow and lithos- pheric thermal structure in the Chagan Depression of the Yingen‐Ejinaqi Basin, north central China. Basin Research 32-6. DOI: https://doi.org/10.1111/bre.12430 124 62-1_acta49-1.qxd 20.6.2022 10:43 Page 124