Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Circumpolar dataset of Soil Organic Carbon north of treeline derived from ENVISAT ASAR GM, link to GeoTIFF

A new approach for the estimation of soil organic carbon (SOC) pools north of the tree line has been developed based on synthetic aperture radar (SAR; ENVISAT Advanced SAR Global Monitoring mode) data. SOC values are directly determined from backscatter values instead of upscaling using land cover or soil classes. The multi-mode capability of SAR allows application across scales. It can be shown that measurements in C band under frozen conditions represent vegetation and surface structure properties which relate to soil properties, specifically SOC. It is estimated that at least 29 Pg C is stored in the upper 30 cm of soils north of the tree line. This is approximately 25 % less than stocks derived from the soil-map-based Northern Circumpolar Soil Carbon Database (NCSCD). The total stored carbon is underestimated since the established empirical relationship is not valid for peatlands or strongly cryoturbated soils. The approach does, however, provide the first spatially consistent account of soil organic carbon across the Arctic. Furthermore, it could be shown that values obtained from 1 km resolution SAR correspond to accounts based on a high spatial resolution (2 m) land cover map over a study area of about 7 × 7 km in NE Siberia. The approach can be also potentially transferred to medium-resolution C-band SAR data such as ENVISAT ASAR Wide Swath with ~120 m resolution but it is in general limited to regions without woody vegetation. Global Monitoring-mode-derived SOC increases with unfrozen period length. This indicates the importance of this parameter for modelling of the spatial distribution of soil organic carbon storage.

(Table 2) Survivorship of transplanted seedlings of Picea glauca inside and outside three tree islands in 1994, 1995, 2009 and after 15 years

The northern boundary of boreal forest and the ranges of tree species are expected to shift northward in response to climate warming, which will result in a decrease in the albedo of areas currently covered by tundra vegetation, an increase in terrestrial carbon sequestration, and an alteration of biodiversity in the current Low Arctic. Central to the prediction of forest expansion is an increase in the reproductive capacity and establishment of individual trees. We assessed cone production, seed viability, and transplanted seedling success of Picea glauca (Moench.) Voss. (white spruce) in the early 1990s and again in the late 2000s at four forest stand sites and eight tree island sites (clonal populations beyond present treeline) in the Mackenzie Delta region of the Northwest Territories, Canada. Over the past 20 years, average temperatures in this region have increased by 0.9 °C. This area has the northernmost forest-tundra ecotone in North America and is one of the few circumpolar regions where the northern limit of conifer trees reaches the Arctic Ocean. We found that cone production and seed viability did not change between the two periods of examination and that both variables decreased northward across the forest-tundra ecotone. Nevertheless, white spruce individuals at the northern limit of the forest-tundra ecotone produced viable seeds. Furthermore, transplanted seedlings were able to survive in the northernmost sites for 15 years, but there were no signs of natural regeneration. These results indicate that if climatic conditions continue to ameliorate, reproductive output will likely increase, but seedling establishment and forest expansion within the forest-tundra of this region is unlikely to occur without the availability of suitable recruitment sites. Processes that affect the availability of recruitment sites are likely to be important elsewhere in the circumpolar ecotone, and should be incorporated into models and predictions of climate change and its effects on the northern forest-tundra ecotone.

Establishment, expansion and regeneration of aspen (Populus tremula) in the Torneträsk area, subarctic Sweden

In subarctic Sweden, recent decadal colonization and expansion of aspen (Populus tremula L.) were recorded. Over the past 100 years, aspen became c. 16 times more abundant, mainly as a result of increased sexual regeneration. Moreover, aspen now reach tree-size (>2 m) at the alpine treeline, an ecotone that has been dominated by mountain birch (Betula pubescens ssp. czerepanovii) for at least the past 4000 years. We found that sexual regeneration in aspen probably occurred seven times or more within the last century. Whereas sexual regeneration occurred during moist years following a year with an exceptionally high June-July temperature, asexual regeneration was favored by warm and dry summers. Disturbance to the birch forest by cyclic moth population outbreaks was critical in aspen establishment in the subalpine area. At the treeline, aspen colonization was less determined by these moth outbreaks, and was mainly restricted by summer temperature. If summer warming persists, aspen spread may continue in subarctic Sweden, particularly at the treeline. However, changing disturbance regimes, future herbivore population dynamics and the responses of aspen's competitors birch and pine to a changing climate may result in different outcomes.

Arctic Siberia: Modern pollen dataset and pollen spectra changes over the last 7,000 years from a lacustrine sediment core (southern Taymyr region)

Palaeoecological investigations in the larch forest-tundra ecotone in northern Siberia have the potential to reveal Holocene environmental variations, which likely have consequences for global climate change because of the strong high-latitude feedback mechanisms. A sediment core, collected from a small lake (radius ~100 m), was used to reconstruct the development of the lake and its catchment as well as vegetation and summer temperatures over the last 7100 calibrated years. A multi-proxy approach was taken including pollen and sedimentological analyses. Our data indicate a gradual replacement of open larch forests by tundra with scattered single trees as found today in the vicinity of the lake. An overall trend of cooling summer temperature from a ~2 °C warmer-than-present mid-Holocene summer temperatures until the establishment of modern conditions around 3000 years ago is reconstructed based on a regional pollen-climate transfer function. The inference of regional vegetation changes was compared to local changes in the lake's catchment. An initial small water depression occurred from 7100 to 6500 cal years BP. Afterwards, a small lake formed and deepened, probably due to thermokarst processes. Although the general trends of local and regional environmental change match, the lake catchment changes show higher variability. Furthermore, changes in the lake catchment slightly precede those in the regional vegetation. Both proxies highlight that marked environmental changes occurred in the Siberian forest-tundra ecotone over the course of the Holocene.

Genetic data from Staurosira (Fragilariophyceae) sedimentary DNA

Genetic investigations on eukaryotic plankton confirmed the existence of modern biogeographic patterns, but analyses of palaeoecological data exploring the temporal variability of these patterns have rarely been presented. Ancient sedimentary DNA proved suitable for investigations of past assemblage turnover in the course of environmental change, but genetic relatedness of the identified lineages has not yet been undertaken. Here, we investigate the relatedness of diatom lineages in Siberian lakes along environmental gradients (i.e. across treeline transects), over geographic distance and through time (i.e. the last 7000 years) using modern and ancient sedimentary DNA. Our results indicate that closely-related Staurosira lineages occur in similar environments and less-related lineages in dissimilar environments, in our case different vegetation and co-varying climatic and limnic variables across treeline transects. Thus our study reveals that environmental conditions rather than geographic distance is reflected by diatom-relatedness patterns in space and time. We tentatively speculate that the detected relatedness pattern in Staurosira across the treeline could be a result of adaptation to diverse environmental conditions across the arctic boreal treeline, however, a geographically-driven divergence and subsequent repopulation of ecologically different habitats might also be a potential explanation for the observed pattern.

Pollen, DNA and Vegetation records from surface sediments of lacustrine lakes along a forest tundra - tundra transect at the Taymyr peninsula, northern Siberia

Reliable information of past vegetation changes are important to project future changes, especially for areas undergoing rapid transitioning such as the boreal treeline. The application of detailed sedDNA records has the potential to enhance our understanding of vegetation changes gained mainly from pollen studies of lake sediments. This study investigates sedDNA and pollen records from 31 lakes along a gradient of increasing larch forest cover in northern Siberia (Taymyr Peninsula) and compares them with vegetation field surveys within the lake's catchment. With respect to vegetation richness, sedDNA recorded 114 taxa, about half of them to species level, while pollen analyses identified 43 pollen taxa. Both approaches exceed the 31 taxa revealed by vegetation field surveys of 400 m**2 plots. From north to south, Larix percentages increase, as is consistently recorded by all three methods. Furthermore, tundra sites are separated from forested sites in the plots of the principal component analyses. Comparison of ordination results by Procrustes and Protest analyses yields a significant fit among all compared pairs of records. Despite the overall comparability of sedDNA and pollen analyses certain idiosyncrasies in the compositional signal are observed, such as high percentages of Alnus and Betula in all pollen spectra and high percentages of Salix in all sedDNA spectra. In conclusion, our results from the treeline show that sedDNA analyses perform better than pollen in recording site-specific richness (i.e. presence/absence of certain vegetation taxa in the direct vicinity of the lake) and perform as good as pollen in tracing regional vegetation composition.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.