Climate, Earth surface processes and soil thermal hydrological conditions drive landscape development, ecosystem functioning and human activities in high latitude regions. These systems are at the focal point of concurrent global change studies as the ongoing shifts in climate regimes has already changed the dynamics of fragile and highly specialized environments across pan Arctic. This thesis aimed to 1) analyze and model extreme air temperatures, soil thermal and hydrological conditions, and the main Earth surface processes (ESP) (cryoturbation, solifluction, nivation and palsa mires) controlling the functioning of high latitude systems in current and future climate conditions; 2) identify the key environmental factors driving the spatial variation of the studied phenomena; and 3) develop methodology for producing novel high quality datasets. To accomplish these objectives, spatial analyses were conducted throughout geographical scales by utilizing multiple statistical modelling approaches, such as regression, machine learning techniques and ensemble forecasting. This thesis was based on unique datasets from the northern Fennoscandia; climate station records from Finland, Sweden and Norway, state of the art climate model simulations, fine scale field measurements collected in arctic alpine tundra and remotely sensed geospatial data. In paper I, accurate extreme air temperature maps were produced, which were notably improved after incorporating the influence of local factors such as topography and water bodies into the spatial models. In paper II, the results showed extreme variation in soil temperature and moisture over very short distances, while revealing the factors controlling the heterogeneity of ground thermal and hydrological conditions. Finally, the modelling outputs in papers III and IV provided new insights into the determination of geomorphic activity patterns across arctic alpine landscapes, while stressing the need for accurate climate data for predictive geomorphological distribution mapping. Importantly, Earth surface processes were found to be extremely climatic sensitivity, and drastic changes in geomorphic systems towards the end of 21st century can be expected. The increase of current temperature conditions by 2 ˚C was projected to cause a near complete loss of active ESPs in the high latitude study area. This thesis demonstrated the applicability of spatial modelling techniques as a useful framework in multiple key challenges of contemporary physical geography. Moreover, with the utilized model ensemble approach, the modelling uncertainty can be reduced while presenting the local trends in response variables more robustly. In future Earth system studies, it is essential to further assess the dynamics of arctic alpine landscapes under changing climatic conditions and identify potential tipping points of these sensitive systems.

Extreme temperatures are key drivers controlling both biotic and abiotic processes, and may be strongly modified by topography and land cover. We modelled mean and extreme temperatures in northern Fennoscandia by combining digital elevation and land cover data with climate observations from northern Finland, Norway and Sweden. Multivariate partitioning technique was utilized to investigate the relative importance of environmental variables for the variation of the three temperature parameters: mean annual absolute minima and maxima, and mean annual temperature. Generalized additive modeling showed good performance, explaining 84–95 % of the temperature variation. The inclusion of remotely sensed variables improved significantly the modelling of thermal extremes in this system. The water cover variables and topography were the most important drivers of minimum temperatures, whereas elevation was the most important factor controlling maximum temperatures. The spatial variability of mean temperatures was clearly driven by geographical location and the effects of topography. Partitioning technique gave novel insights into temperature-environment relationship at the meso-scale and thus proved to be useful tool for the study of the extreme temperatures in the high-latitude setting.

Methane (CH4) strongly contributes to observed global warming. As natural CH4 emissions mainly originate from wet ecosystems, it is important to unravel how climate change may affect these emissions. This is especially true for ebullition (bubble flux from sediments), a pathway that has long been underestimated but generally dominates emissions. Here we show a remarkably strong relationship between CH4 ebullition and temperature across a wide range of freshwater ecosystems on different continents using multi-seasonal CH4 ebullition data from the literature. As these temperature–ebullition relationships may have been affected by seasonal variation in organic matter availability, we also conducted a controlled year-round mesocosm experiment. Here 4 °C warming led to 51% higher total annual CH4 ebullition, while diffusion was not affected. Our combined findings suggest that global warming will strongly enhance freshwater CH4 emissions through a disproportional increase in ebullition (6–20% per 1 °C increase), contributing to global warming.

The collection of modern, spatially extensive pollen data is important for the interpretation of fossil pollen assemblages and the reconstruction of past vegetation communities in space and time. Modern datasets are readily available for percentage data but lacking for pollen accumulation rates (PARs). Filling this gap has been the motivation of the pollen monitoring network, whose contributors monitored pollen deposition in modified Tauber traps for several years or decades across Europe. Here we present this monitoring dataset consisting of 351 trap locations with a total of 2742 annual samples covering the period from 1981 to 2017. This dataset shows that total PAR is influenced by forest cover and climate parameters, which determine pollen productivity and correlate with latitude. Treeless vegetation produced PAR values of at least 140 grains cm⁻² yr⁻¹. Tree PAR increased by at least 400 grains cm⁻² yr⁻¹ with each 10 % increase in forest cover. Pollen traps situated beyond 200 km of the distribution of a given tree species still collect occasional pollen grains of that species. The threshold of this long-distance transport differs for individual species and is generally below 60 grains cm⁻² yr⁻¹. Comparisons between modern and fossil PAR from the same regions show similar values. For temperate taxa, modern analogues for fossil PARs are generally found downslope or southward of the fossil sites. While we do not find modern situations comparable to fossil PAR values of some taxa (e.g. Corylus), CO₂ fertilization and land use may cause high modern PARs that are not documented in the fossil record. The modern data are now publicly available in the Neotoma Paleoecology Database and aid interpretations of fossil PAR data.

There is broad consensus that, via changes in stomatal conductance, plants moderate the exchanges of water and carbon between the biosphere and atmosphere, playing a major role in global hydroclimate. Tree rings record atmospheric CO2 concentration (ca) and its isotopic composition (13C/12C)—mediated by stomatal and photosynthetic influences—that can be expressed in terms of intrinsic water-use efficiency (W). Here, we compile a global W dataset based on 422 tree-ring isotope series and report that W increased with ca over the twentieth century, but the rates of increase (dW/dca) declined by half. Angiosperms contributed more than gymnosperms to the slowdown, and in recent decades, dW/dca for angiosperms was close to zero. dW/dca varies widely across climatic regions and reflects pauses in emissions during the Great Depression and after World War II. There is strong spatial variability in climate forcing via an increasing W, which is weakening globally with time.

Species of Russula subsect. Xerampelinae are notoriously difficult to identify and name and have not been subject to molecular study. A group of species, referred to here as the R. clavipes complex, growing in association with Salix, Betula and Populus as well as coniferous tree species from temperate to arctic and alpine habitats, were examined. Analyses of the nuc rDNA internal transcribed spacer (ITS) region and a numerical analysis of morphological characters were used. The R. clavipes complex is a monophyletic group within Russula subsect. Xerampelinae, according to molecular results. The complex includes three species: R. nuoljae is a phylogenetically and morphologically well-supported species while the other two, R. clavipes and R. pascua, are similar based on ITS data and morphology but separate based on their ecology. Russula pseudoolivascens is conspecific with R. clavipes. Several combinations of characters traditionally used in the taxonomy of R. subsect. Xerampelinae are inappropriate for species delimitation in this group and the adequacy of the ITS for species identification in this group is discussed. Detailed microscopic observations on the type collection of R. nuoljae are presented and illustrated, along with a key to the European members of R. subsect. Xerampelinae.

1 Decomposition of plant litter, a key component of the global carbon budget, is hierarchically controlled by the triad: climate > litter quality > soil organisms. Given the sensitivity of decomposition to temperature, especially in cold biomes, it has been hypothesized that global warming will lead to increased litter decomposition rates, both through direct temperature effects and through indirect effects on litter quality and soil organisms. 2 A meta-analysis of experimental warming studies in cold biomes (34 site-species combinations) showed that warming resulted in slightly increased decomposition rates. However, this response was strongly dependent on the method used: open top chambers reduced decomposition rates, whereas heating lamps stimulated decomposition rates. The low responsiveness was mainly due to moisture-limited decomposition rates in the warming treatments, especially at mesic and xeric sites. This control of litter decomposition by both temperature and moisture was corroborated by natural gradient studies. 3 Interspecific differences in litter quality and decomposability are substantially larger than warming-induced phenotypic responses. Thus, the changes in the species composition and structure of plant communities that have been observed in medium-term warming studies in cold biomes will have a considerably greater impact on ecosystem litter decomposition than phenotypic responses. 4 Soil fauna communities in cold biomes are responsive to climate warming. Moreover, temperature-driven migration of the, hitherto absent, large comminuters to high-latitude sites may significantly increase decomposition rates. However, we do not know how far-reaching the consequences of changes in the species composition and structure of the soil community are for litter decomposition, as there is a lack of data on functional species redundancy and the species’ dispersal ability. 5 Global warming will lead to increased litter decomposition rates only if there is sufficient soil moisture. Hence, climate scenario and experimental studies should focus more on both factors and their interaction. As interspecific differences in potential decomposability and litter chemistry are substantially larger than phenotypic responses to warming, the focus of future research should be on the former. In addition, more light should be shed on the below-ground ‘darkness’ to evaluate the ecological significance of warming-induced soil fauna community changes for litter decomposition processes in cold biomes.

Climate scenarios for high-latitude areas predict not only increased summer temperatures, but also larger variation in snowfall and winter temperatures. By using open-top chambers, we experimentally manipulated both summer temperatures and winter and spring snow accumulations and temperatures independently in a blanket bog in subarctic Sweden, yielding six climate scenarios. We studied the effects of these scenarios on flowering phenology and flower production of Andromeda polifolia (woody evergreen) and Rubus chamaemorus (perennial herb) during 2 years. The second year of our study (2002) was characterized by unusually high spring and early summer temperatures. Our winter manipulations led to consistent increases in winter snow cover. As a result, average and minimum air and soil temperatures in the high snow cover treatments were higher than in the winter ambient treatments, whereas temperature fluctuations were smaller. Spring warming resulted in higher average, minimum, and maximum soil temperatures. Summer warming led to higher air and soil temperatures in mid-summer (June-July), but not in late summer (August-September). The unusually high temperatures in 2002 advanced the median flowering date by 2 weeks for both species in all treatments. Superimposed on this effect, we found that for both Andromeda and Rubus, all our climate treatments (except summer warming for Rubus) advanced flowering by 1-4 days. The total flower production of both species showed a more or less similar response: flower production in the warm year 2002 exceeded that in 2001 by far. However, in both species flower production was only stimulated by the spring-warming treatments. Our results show that the reproductive ecology of both species is very responsive to climate change but this response is very dependent on specific climate events, especially those that occur in winter and spring. This suggests that high-latitude climate change experiments should focus more on winter and spring events than has been the case so far.

Nutrient resorption from senescing plant tissues is an important determinant of the fitness of plant populations in nutrient-poor ecosystems, because it makes plants less dependent on current nutrient uptake. Moreover, it can have significant “afterlife” effects through its impact on litter chemistry and litter decomposability. Little is known about the effects of climate change on nutrient resorption. We studied the effects of climate change treatments (including winter snow addition, and spring and/or summer warming) on nutrient resorption of four dominant species in a nutrient-poor subarctic peatland. These species were Betula nana (woody deciduous), Vaccinium uliginosum (woody deciduous), Calamagrostis lapponica (graminoid) and Rubus chamaemorus (forb). After five years of treatments both mature and senesced leaf N concentrations showed a small but significant overall reduction in response to the climate treatments. However, the effects were species-specific. For example, in the controls the N concentration in senesced leaves of Calamagrostis (3.0 +/- 0.2 mg N g(-1)) was about four times lower than for Rubus (11.2 +/- 0.2 mg N g(-1)). There were no significant treatment effects on N resorption efficiency (% of the N pool in mature leaves that is resorbed during senescence). The nitrogen resorption efficiency of Calamagrostis (about 80%) was higher than in the other three species (about 60%). Thus, climate change has only a minor impact on nutrient resorption parameters. However, given the substantial interspecific differences in these parameters, substantial changes in plant-soil feedbacks may be expected as a result of the observed changes in the species composition of high-latitude vegetation. These changes are species-specific and thus difficult to predict.

Abstract Vast areas of (sub)arctic tundra are dominated by the ericoid dwarf shrub Empetrum hermaphroditum. Recent experimental and observational data have shown that Empetrum can be damaged heavily by recurrent extreme winter warming. In addition, summer warming leads to increased soil N availability in tundra ecosystems. In a 7-year experiment, I investigated the recovery of subarctic Empetrum-dominated tundra vegetation using a factorial combination of various degrees of aboveground Empetrum removal (simulating the damaging effects of extreme winter warming) and N addition (simulating one of the effects of summer warming). After 7 years no new species had established in the plots. The growth of planted Betula nana seedlings was stimulated by Empetrum removal and reduced by N addition. This Empetrum-dominated tundra ecosystem was resilient against severe disturbances. Only when Empetrum was 100% removed did it fail to recover, and only in combination with high N supply the subordinate species (notably Eriophorum vaginatum and Rubus chamaemorus, a graminoid and a forb) could benefit. In the 50% removal treatment Empetrum recovered in 7 years when no N was supplied and the cover of the subordinate species did not change. However, when N was added Empetrum recovered faster (in 4 years) and the subordinates decreased. When Empetrum was not removed and N was added, Empetrum even increased in abundance at the expense of the subordinate species. Thus, profound changes in tundra ecosystems can only be expected when Empetrum is very heavily damaged as a result of recurrent extreme winter warming and when soil N availability is increased as a result of summer warming. These changes in species composition upon extreme disturbance events may lead to a wide variety of ecosystem feedbacks and cascade processes as this tundra system is relatively species-poor, and can be hypothesized to have low functional redundancy.

Climate change in cold biomes not only involves higher summer temperatures, but also warmer springs and more winter precipitation. So far, little is known about species responses to these seasonal components of climate change.

 We experimentally manipulated spring and summer temperatures and winter snow accumulation and temperatures independently in a peatland in sub‐arctic Sweden. This yielded six climate scenarios and we studied the responses of the peat moss Sphagnum fuscum, the evergreen dwarf shrubs Empetrum hermaphroditum and Andromeda polifolia, the deciduous dwarf shrubs Betula nana and Vaccinium uliginosum, the grass Calamagrostis lapponica and the forb Rubus chamaemorus.                  

We found substantial interspecific differences in leaf nutrient and carbon exchange variables that reflect the response of tundra plants to climate change. S. fuscum had the lowest N and P concentrations, with increasing N and P concentrations (and decreasing C/N and C/P ratios) going from evergreen dwarf shrubs, to the grass, deciduous dwarf shrubs and the forb. Leaf N/P ratios varied between 10 and 14 which points to N‐limited plant growth.                  

The natural abundance of ¹⁵N varied very strongly among species and growth forms. These differences corresponded with the presence and type of mycorrhizal association in the plant roots. Leaf carbon isotope discrimination also differed strongly among species and growth forms, but the absolute differences were relatively small (< 5‰). The rank order was: forb and graminoid < moss and evergreen shrubs < deciduous shrubs.                  

After 4 years, the effects of realistic climate change manipulations on leaf nutrient and carbon exchange variables were idiosyncratic with respect to species and generally small compared to the differences among species and growth forms. At the phenotypic level, spring warming or winter snow addition effects occurred as frequently as summer warming effects.

This implies that the changes in the species composition and structure of plant communities that have been observed in medium‐term warming studies in cold biomes will have much more impact on plant‐mediated nutrient and carbon cycling pathways and rates than climate‐change induced phenotypic responses, irrespective of the seasonal timing of these climate changes.

Artic and subarctic environments are particularly sensitive to climate change with a faster warming compared to other latitudes. Vegetation is changing but its role on the biogeochemical cycling is poorly understood. In this study, we evaluated the distribution of trace elements in subarctic soils from different land covers at Abisko, northern Sweden: grassland, moor, broad-leaved forest, and peat bog. Using various multivariate analysis approaches, results indicated a spatial heterogeneity with a strong influence of soil horizon classes considered: lithogenic elements (e.g., Al, Cr, Ti) were accumulated in mineral horizon classes and surface process-influenced elements (e.g., Cd, Cu, Se) in organic horizon classes. Atmospheric influences included contamination by both local mines (e.g., Cu, Fe, Ni) and regional or long-range atmospheric transport (e.g., Cd, Pb, Zn). A non-negative matrix factorization was used to estimate, for each element, the contribution of various sources identified. For the first time, a comparison between geochemical and ecological data was performed to evaluate the influence of vegetation on element distribution. Apart from soil pH that could control dynamics of As, Cu, and Se, two vegetation classes were reported to be correlated to geochemical factors: forbs and shrubs/dwarf shrubs probably due to their annual vs. perennial activities, respectively. Since these are considered as the main vegetation classes that quickly evolve with climate change, we expect to see modifications in trace element biogeochemical cycling in the future.

Permafrost peatlands have the capacity to store significant amounts of carbon, and thus they act as important controllers of the climate. Approximately 14% of the world’s soil organic carbon pool is stored in permafrost peatlands, which are sensitive to climatic fluctuations due to their location in the high latitudes of the subarctic zone. Permafrost peatlands also act as a habitat for a large number of moisture-sensitive organisms, such as bryophytes and testate amoebae, which can be used to study how the hydrology of peatlands has changed and will continue to change throughout time, giving us an opportunity to predict the future of peatlands under a changing climate. In this Master’s Thesis I examined the testate amoebae composition and used these species as indicators to study hydrological fluctuations from three subarctic permafrost peatland cores extracted from Taavavuoma and Abisko in northern Sweden. The species compositions were combined with radiocarbon (14C) and lead (210Pb) dates to reconstruct the past water table levels for the late Holocene, spanning four climatic periods. The reconstructions were then compared to past studies on testate amoebae to understand how permafrost peatlands and their species assemblages respond to changes in the hydrology of the active layer of the peat. Out of the study sites only the Taavavuoma cores spanned the Dark Age Cold Period (DACP) and Medieval Climate Anomaly (MCA). Species compositions in both cores indicated fluctuating water tables during the DACP, but during the MCA the results began to contradict with one site showing a wetter, and the other a drier MCA. Two out of three study sites indicated a wetter Little Ice Age and a drier Post-Industrial Warming, supporting past studies indicating similar results, whereas one study site gave opposite results. The results indicated large variability in testate amoebae assemblages throughout time, indicating that the hydrology of peatlands can change very abruptly and vary considerably even on a local scale. Modelling is however complicated by the poorly known ecology of testate amoebae, which is why a multi-proxy approach is essential to reliably predict the future fate of permafrost peatlands.