Tundra Regions Research Articles

Divergent future trends in global land greening and carbon sequestration under climate change scenarios

Abstract Over the past four decades, global land greening has promoted carbon (C) storage in terrestrial ecosystems. However, whether the future trajectory of this positive greening effect on ecosystem C sequestration is sustainable under various climate scenarios remains uncertain. Here, using projections from ten Earth system models, we found divergent trends in the relationship between global land greening and C storage among three distinct shared socioeconomic pathways (SSPs). We identified global transition times for their positive relationships, which will occur in 2034, 2038, and 2048 for SSP1-2.6, SSP2-4.5, and SSP5-8.5, respectively. We found a widespread decoupling trend of vegetation greening and ecosystem C storage under medium to high emission scenarios, particularly in the tundra, boreal forest, and grassland regions, by delving deeper into six transition modes. These findings underscore the highly uncertain role of vegetation greening in land C sequestration under prospective climate change scenarios.

Slow post-fire carbon balance recovery despite increased net uptake rates in Alaskan tundra

Abstract Increasing wildfire occurrence and intensity have immediate effects on northern ecosystems due to combustion of aboveground vegetation and belowground soil organic matter. These immediate impacts have indirect and longer term effects, including deepening of the active layer, changes in soil decomposition rates, and shifts in plant community composition. Despite the increasing fire impacts across the tundra region, the implications of wildfire on ecosystem carbon balance are not well understood. Using paired eddy covariance towers in unburned and burned tundra, we examined the effects of a 2015 wildfire on carbon dioxide and methane fluxes in a wetland tundra ecosystem in the Yukon-Kuskokwim Delta, Alaska, from 2020 to 2022. Wildfire increased the amplitude and variability of carbon uptake and release on seasonal and annual timescales and increased the temperature sensitivity of soil respiration. Seven years post fire, there was annual net uptake in both unburned and burned tundra based on net ecosystem exchange, with the sink strength of burned tundra exceeding that of the unburned tundra by 1.18 to 1.64 times. However, when considering emissions, it would take approximately 86 years to recover the carbon lost from the wildfire itself. Soil moisture was a dominant driver of fluxes and positively associated with higher rates of carbon dioxide uptake and release and methane release. This study underscores the importance of understanding the effects of wildfire-induced shifts on tundra carbon cycling, allowing better predictions of long-term landscape-scale climate feedbacks as the climate continues to warm.

Spatial variability in Arctic–boreal fire regimes influenced by environmental and human factors

Wildfire activity in Arctic and boreal regions is rapidly increasing, with severe consequences for climate and human health. Regional long-term variations in fire frequency and intensity characterize fire regimes. The spatial variability in Arctic–boreal fire regimes and their environmental and anthropogenic drivers, however, remain poorly understood. Here we present a fire tracking system to map the sub-daily evolution of all circumpolar Arctic–boreal fires between 2012 and 2023 using 375 m Visible Infrared Imaging Radiometer Suite active fire detections and the resulting dataset of the ignition time, location, size, duration, spread and intensity of individual fires. We use this dataset to classify the Arctic–boreal biomes into seven distinct ‘pyroregions’ with unique climatic and geographic environments. We find that these pyroregions exhibit varying responses to environmental drivers, with boreal North America, eastern Siberia and northern tundra regions showing the highest sensitivity to climate and lightning density. In addition, anthropogenic factors play an important role in influencing fire number and size, interacting with other factors. Understanding the spatial variability of fire regimes and its interconnected drivers in the Arctic–boreal domain is important for improving future predictions of fire activity and identifying areas at risk for extreme events.

A synthesized field survey database of vegetation and active-layer properties for the Alaskan tundra (1972–2020)

Abstract. Studies in recent decades have shown strong evidence of physical and biological changes in the Arctic tundra, largely in response to rapid rates of warming. Given the important implications of these changes for ecosystem services, hydrology, surface energy balance, carbon budgets, and climate feedbacks, research on the trends and patterns of these changes is becoming increasingly important and can help better constrain estimates of local, regional, and global impacts as well as inform mitigation and adaptation strategies. Despite this great need, scientific understanding of tundra ecology and change remains limited, largely due to the inaccessibility of this region and less intensive studies compared to other terrestrial biomes. A synthesis of existing datasets from past field studies can make field data more accessible and open up possibilities for collaborative research as well as for investigating and informing future studies. Here, we synthesize field datasets of vegetation and active-layer properties from the Alaskan tundra, one of the most well-studied tundra regions. Given the potentially increasing intensive fire regimes in the tundra, fire history and severity attributes have been added to data points where available. The resulting database is a resource that future investigators can employ to analyze spatial and temporal patterns in soil, vegetation, and fire disturbance-related environmental variables across the Alaskan tundra. This database, titled the Synthesized Alaskan Tundra Field Database (SATFiD), can be accessed at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC) for Biogeochemical Dynamics (Chen et al., 2023: https://doi.org/10.3334/ORNLDAAC/2177).

Spatiotemporal coverage capability of moon-based SAR for water dynamic monitoring in global vegetation tipping elements

ABSTRACT Changes in water content dynamics can serve as predictive indicators for catastrophic transitions in vegetation tipping elements. Vegetation water content is a rapidly changing parameter with sub-daily fluctuations, and current space-borne sensors fall short of providing adequate spatiotemporal coverage for effective analysis. The Moon-based Earth observation, i.e. Synthetic Aperture Radar has recently showcased prolonged durations, rapid revisiting rates, and extensive monitoring capabilities with a decent spatial resolution (10-100 m). This opens up new possibilities for vegetation water content monitoring which has not yet been discussed. The study has assessed its monitoring capability in four critical vegetation tipping elements and considered both overall temporal observation availability fraction and diurnal variation acquisition. The results indicate that Amazon and Sahelian regions can be observed for ∼30% of a nutation period (∼18.6 years). Diurnal variations in water content can be captured with a 1-hour interval every 16–23 days for the Amazon and Sahelian regions. For polar tundra and boreal regions, continuous observation spanning 7–12 h can be achieved. These outstanding characteristics offer a deeper understanding of water content dynamics in vegetation tipping elements, which will enhance our prediction for their critical transitions.

Review article: Terrestrial dissolved organic carbon in northern permafrost

Abstract. As the permafrost region warms and permafrost soils thaw, vast stores of soil organic carbon (C) become vulnerable to enhanced microbial decomposition and lateral transport into aquatic ecosystems as dissolved organic carbon (DOC). The mobilization of permafrost soil C can drastically alter the net northern permafrost C budget. DOC entering aquatic ecosystems becomes biologically available for degradation as well as other types of aquatic processing. However, it currently remains unclear which landscape characteristics are most relevant to consider in terms of predicting DOC concentrations entering aquatic systems from permafrost regions. Here, we conducted a systematic review of 111 studies relating to, or including, concentrations of DOC in terrestrial permafrost ecosystems in the northern circumpolar region published between 2000 and 2022. We present a new permafrost DOC dataset consisting of 2845 DOC concentrations, collected from the top 3 m in permafrost soils across the northern circumpolar region. Concentrations of DOC ranged from 0.1 to 500 mg L−1 (median = 41 mg L−1) across all permafrost zones, ecoregions, soil types, and thermal horizons. Across the permafrost zones, the highest median DOC concentrations were in the sporadic permafrost zone (101 mg L−1), while lower concentrations were found in the discontinuous (60 mg L−1) and continuous (59 mg L−1) permafrost zones. However, median DOC concentrations varied in these zones across ecosystem type, with the highest median DOC concentrations in each ecosystem type of 66 and 63 mg L−1 found in coastal tundra and permafrost bog ecosystems, respectively. Coastal tundra (130 mg L−1), permafrost bogs (78 mg L−1), and permafrost wetlands (57 mg L−1) had the highest median DOC concentrations in the permafrost lens, representing a potentially long-term store of DOC. Other than in Yedoma ecosystems, DOC concentrations were found to increase following permafrost thaw and were highly constrained by total dissolved nitrogen concentrations. This systematic review highlights how DOC concentrations differ between organic- or mineral-rich deposits across the circumpolar permafrost region and identifies coastal tundra regions as areas of potentially important DOC mobilization. The quantity of permafrost-derived DOC exported laterally to aquatic ecosystems is an important step for predicting its vulnerability to decomposition.

Analysis of taiga and tundra lake browning trends from 2002 to 2021 using MODIS data

Lakes in taiga and tundra regions may be silently undergoing changes due to global warming. One of those changes is browning in lake color. The browning interacts with the carbon cycle, ecosystem dynamics, and water quality in freshwater systems. However, spatiotemporal variabilities of browning in these regions have not been well documented. Using MODIS remote sensing reflectance at near ultraviolet wavelengths from 2002 to 2021 on the Google Earth Engine platform, we quantified long-term browning trends across 7616 lakes (larger than 10 km2) in taiga and tundra biomes. These lakes showed an overall decreased trend in browning (Theil-Sen Slope = 0.00015), with ∼36% of these lakes showing browning trends, and ∼1% of these lakes showing statistically significant (p-value <0.05) browning trends. The browning trends more likely occurred in small lakes in high latitude, low ground ice content regions, where air temperature increased and precipitation decreased. While temperature is projected to increase in response to climate change, our results provide one means to understand how biogeochemical cycles and ecological dynamics respond to climate change.

Observed links between heatwaves and wildfires across Northern high latitudes

Data on Arctic and Sub-Arctic summer heat events are limited due to the sparse network of surface observation stations. Here, we analyze large heat events within 60°–80°N using land surface temperature (LST) observations from the moderate resolution imaging spectroradiometer (MODIS) sensor aboard the Terra satellite. Our heatwave (HW) detection method uses exceedances of the climatological 90th percentile of LST across summer months, and a spatio-temporal density-based clustering algorithm to distinguish space-time coherent events across Northern Hemispheric high latitudes for the summers of 2000–2022. We find a close link between HW duration, spatial extent and amplitude across the study region (correlations ranging from 0.63 to 0.73). MODIS-derived burned area data show that wildfire seasons are significantly correlated to summer HW activity, particularly in Siberia (r = 0.87 at p < 0.05) and Alaska and NW Canada (r = 0.45 at p < 0.05), and are also spatially co-located. MODIS active fire data also show substantial increase during larger heat events. For the strongest HWs in Siberia, the peak in daily fire count (from the MODIS active fire archive) lags behind peak HW activity. We conclude that there is a close link between intense fire summers and extensive HWs over boreal and shrub tundra regions (Interior Alaska, the Canadian Prairies and Eastern Siberia).

Wing lengths of three Arctic butterfly species decrease as summers warm in Alaska

Climate warming can cause arthropods to express plastic and/or evolved changes in morphology. Previous studies have demonstrated that body sizes of Arctic butterflies are influenced by the temperatures experienced as larvae. To investigate whether this was occurring among Alaskan butterflies, we analyzed temporal trends in the wing sizes of three Holarctic species, Colias hecla, Boloria chariclea and Boloria freija, using museum specimens collected in Arctic tundra regions of Alaska between 1971 and 1995. Wing length was compared to accumulated growing degree days (GDD) during both the spring of the year collected and the previous year's summer during the normal period of larval development. We used mixed‐effects models to test if spring and summer temperatures affected adult morphology. Results show that for every 1°C increase in average seasonal temperature, wingspans decreased between 0.7 and 5 mm, with B. freija the most strongly affected. Our results suggest that the morphological sensitivity of Arctic butterflies to warming is the outcome of interactions between life‐history traits and regional climate, with all species sensitive to warming the summer before the flight year as well as warming the spring of the flight year. Boloria freija, which overwinters as late instar larvae that do not feed before pupation the following spring, was particularly strongly affected by summer warming.

Variability and drivers of winter near-surface temperatures over boreal and tundra landscapes

Abstract. Winter near-surface air temperatures have important implications for ecosystem functioning such as vegetation dynamics and carbon cycling. In cold environments, the persistence of seasonal snow cover can exert a strong control on the near-surface temperatures. However, the lack of in situ measurements of both snow cover duration and surface temperatures over high latitudes has made it difficult to estimate the spatio-temporal variability in this relationship. Here, we quantified the fine-scale variability in winter near-surface air temperatures (+2 cm) and snow cover duration (calculated from temperature time series) using a total of 441 microclimate loggers in seven study areas across boreal and tundra landscapes in Finland during 2019–2021. We further examined the drivers behind this variation using a structural equation model and the extent to which near-surface air temperatures are buffered from free-air temperatures during winter. Our results show that while average winter near-surface temperatures stay close to 0 ∘C across the study domain, there are large differences in their fine-scale variability among the study areas. Areas with large topographical variation, as well as areas with shallow snowpacks, showed the greatest variation in near-surface temperatures and in snow cover duration. In the tundra, for example, differences in minimum near-surface temperatures between study sites were close to 30 ∘C and topography was shown to be an important driver of this variability. In contrast, flat topography and long snow cover duration led to little spatial variation, as well as long periods of decoupling between near-surface and air temperatures. Quantifying and understanding the landscape-wide variation in winter microclimates improves our ability to predict the local effects of climate change in the rapidly warming boreal and tundra regions.

Traditional nature management as a way to prevent the loss of wildlife species in a changing environment

Wildlife species of tundra and sub-arctic boreal regions are facing an immediate threat to its existence owning to environmental changes. This paper highlights how traditional nature management in northwestern Russia and northeastern Canada is able to reduce threats to the environment, focusing on vulnerability of biodiversity to changing climatic conditions. The goal of this study is to identify key areas of wildlife species loss during a changing climate by exploring the ability of traditional nature management to support environmentally sustainable habitats for the existence of the most typical biomes of tundra and sub-arctic boreal landscapes. The differentiating biodiversity method was used to determinate presence of rare species as a criterion of non-disturbed areas. This research is based on statistical data on biodiversity dynamics, meteorological data, reports on environmental conditions, cartographic materials, satellite images collected from open sources, and fieldworks. The author indicates non-disturbed sites in terms of biological resources protection in the studied regions. Obtained results confirm that on territories where traditional nature management is carried out, the reduction of biodiversity is much lower than on areas located in equal environmental conditions.

Global temperature homogenization can obliterate temporal isolation in migratory animals with potential loss of population structure.

Climate change is expected to increase the spatial autocorrelation of temperature, resulting in greater synchronization of climate variables worldwide. Possibly such 'homogenization of the world' leads to elevated risks of extinction and loss of biodiversity. In this study, we develop an empirical example on how increasing synchrony of global temperatures can affect population structure in migratory animals. We studied two subspecies of bar-tailed godwits Limosa lapponica breeding in tundra regions in Siberia: yamalensis in the west and taymyrensis further east and north. These subspecies share pre- and post-breeding stopover areas, thus being partially sympatric, but exhibiting temporal segregation. The latter is believed to facilitate reproductive isolation. Using satellite tracking data, we show that migration timing of both subspecies is correlated with the date of snowmelt in their respective breeding sites (later at the taymyrensis breeding range). Snow-cover satellite images demonstrate that the breeding ranges are on different climate trajectories and become more synchronized over time: between 1997 and 2020, the date of snowmelt advanced on average by 0.5 days/year in the taymyrensis breeding range, while it remained stable in the yamalensis breeding range. Previous findings showed how taymyrensis responded to earlier snowmelt by advancing arrival and clutch initiation. In the predicted absence of such advancements in yamalensis, we expect that the two populations will be synchronized by 2036-2040. Since bar-tailed godwits are social migrants, this raises the possibility of population exchange and prompts the question whether the two subspecies can maintain their geographic and morphological differences and population-specific migratory routines. The proposed scenario may apply to a wide range of (social) migrants as temporal segregation is crucial for promoting and maintaining reproductive isolation in many (partially sympatric) migratory populations. Homogenization of previously isolated populations could be an important consequence of increasing synchronized environments and hence climate change.

Upscaling Wetland Methane Emissions From the FLUXNET‐CH4 Eddy Covariance Network (UpCH4 v1.0): Model Development, Network Assessment, and Budget Comparison

AbstractWetlands are responsible for 20%–31% of global methane (CH4) emissions and account for a large source of uncertainty in the global CH4 budget. Data‐driven upscaling of CH4 fluxes from eddy covariance measurements can provide new and independent bottom‐up estimates of wetland CH4 emissions. Here, we develop a six‐predictor random forest upscaling model (UpCH4), trained on 119 site‐years of eddy covariance CH4 flux data from 43 freshwater wetland sites in the FLUXNET‐CH4 Community Product. Network patterns in site‐level annual means and mean seasonal cycles of CH4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash‐Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH4 emissions of 146 ± 43 TgCH4 y−1 for 2001–2018 which agrees closely with current bottom‐up land surface models (102–181 TgCH4 y−1) and overlaps with top‐down atmospheric inversion models (155–200 TgCH4 y−1). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH4 fluxes has the potential to produce realistic extra‐tropical wetland CH4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid‐to‐arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/2253).

Tundra Browning in the Indigirka Lowlands (North‐Eastern Siberia) Explained by Drought, Floods and Small‐Scale Vegetation Shifts

AbstractContrary to the general “greening of the Arctic”, the Siberian Indigirka Lowlands show strong “browning” (a decrease in the Normalized Difference Vegetation Index or “NDVI”) in various recent satellite records. Since greening and browning are generally indicative of increases and losses in photosynthetically active biomass, this browning trend may have implications for the carbon balance and vegetation of this Arctic tundra region. To explore potential mechanisms responsible for this trend break from general Arctic greening, we studied timeseries of Landsat summer maximum NDVI, weather data, and high‐resolution maps of vegetation compositional change, topography, geomorphology and hydrology. We find that a significant proportion of browning (lower summer NDVI) is explained by moisture dynamics, with high snow depths and resulting floods as well as summer drought coinciding with low NDVI. Relations between seasonal weather variables and NDVI are spatially heterogeneous, with floodplains, drained thaw lake basins and Yedoma ridges showing different patterns of association with weather variables. Low summer NDVI after high snowfall was particularly evident in floodplains, likely explained by early summer floods. Local small‐scale vegetation changes explained only small amounts of variance in browning rates in Landsat NDVI. Local expansion of Sphagnum vegetation in particular may have contributed to recent browning of our study site, but higher resolution NDVI timeseries are necessary to accurately constrain the role of small‐scale vegetation shifts. Overall, associations identified in this study suggest that future increases in Arctic precipitation variability and extremes may limit tundra greening, but to different extents even across comparatively small topographical contrasts.

Trees on the tundra: warmer climate might not favor prostrate Larix tree but Betula nana shrub growth in Siberian tundra (Lena River Delta)

Tundra is primarily a habitat for shrub growth, not trees, but growth of prostrate forms of trees has been reported occasionally from the subarctic tundra region. In the light of on-going climate change, climate sensitivity studies of these unique trees are essential to predict vegetation dynamics and potential northward expansion of boreal forest tree species into tundra. Here we studied one of the northernmost Larix Mill. trees and Betula nana L. shrubs (72°N) from the Siberian tundra for the common period 1980-2017. We took advantage of the discovery of a single cohort of prostrate Larix trees within a tundra ecosystem, i.e., ca. 60 km northwards from the northern treeline, and compared climate-growth relationships of the two species. Both woody plants were sensitive to the July temperature, however this relationship was stable across the entire study period (1980-2017) only for Betula nana chronology. Additionally, radial growth of Larix trees became negatively correlated to temperatures during the previous summer. In recent period moisture sensitivity between Larix trees and Betula nana shrubs was contrasting, with generally wetter soil conditions favoring Larix trees growth and dryer conditions promoting Betula nana growth. Our study indicates that Larix trees radial growth in recent years is more sensitive to moisture than to summer air temperatures, whereas temperature sensitivity of Betula nana shrub is stable over time. We provide first detailed insight into the annual resolution on Larix tree growth sensitivity to climate in the heart of the tundra. The potentially higher Betula nana shrub resistance to warmer and drier climate versus Larix trees on a tundra revealed in our study needs to be further examined across habitats of various soil, moisture and permafrost status.

Observation of spatial and temporal patterns of seasonal ground deformation in central Yakutia using time series InSAR data in the freezing season

The seasonal freeze/thaw variations of the active layer cause seasonal ground deformation, such as thaw subsidence in summer and frost heave in winter due to the cyclic phase changes of soil water. The Synthetic Aperture Radar (SAR) remote sensing, especially interferometric SAR (InSAR) techniques, has been proven to be a very useful tool for observing surface displacements of permafrost environments. Most previous studies on seasonal InSAR displacement have been concentrated in several study sites in the Arctic tundra regions. The objective of this study is to extend the InSAR analysis of seasonal ground deformations to boreal forest regions, where monitoring ground deformation is particularly important in terms of better understanding active layer dynamics under various ecological and geological conditions. In this study, we presented a strategy to analyze the time series SAR data in the freezing season to estimate the differential phase for seasonal ground deformation while overcoming the problem of coherence losses due to the rapid growth of vegetation during the short growing season. We applied the SBAS-based InSAR processing method to Sentinel-1 SAR in the central Yakutian lowlands, where boreal forests and thermokarst landforms are widely distributed, and retrieved frost-related displacement signals successfully. In addition to the magnitude of ground deformation in the freezing season, we further examined the temporal deformation pattern related to the progression of the freezing front in the soil during active freezing. Based on the Stefan solution, we proposed a piecewise linear model for the time series deformation. The proposed model divides the time series displacement according to the progression of freezing into three linear segments, and the slope parameter of each linear model can be interpreted as indicating the hydrological characteristics of the upper, middle, and lower parts of the soil. The results for the Yakutian lowlands illustrated the potential of SAR remote sensing to observe and quantify spatial details of the thermal and hydrological properties inside the active layer soil even in highly vegetated subarctic boreal forest areas.

FIRST DATA ON THE DYTISCIDAE (COLEOPTERA) FROM THE COASTAL TUNDRA OF YUGORSKY PENINSULA, POLAR RUSSIA

In July 2018, diving beetles were collected in the most typical water bodies and other habitats in the vicinity of Amderma, Kara Sea coast, northern Yugorsky Peninsula. The fauna comprises 15 dytiscid species from 7 genera. The most diverse genera were Agabus and Hydroporus. The record of Dytiscus lapponicus was the northernmost for this species, while the reports of Hydroporus cf. fuscipennis and Agabus pallens were the northernmost for the Palaearctic parts of their distribution ranges. Original photographs illustrate peculiar specimens of Dytiscus lapponicus with dark head, pronotum and scutellum. Most of the species show vast Holarctic Arcto-Boreal or Arcto-Boreal-Montane distributions. The species composition of the Dytiscidae is quite similar to those of the northernmost mainland regions of the European Northeast (Pakhancheskaya Bay, Kara Tundra and Pamal), as well as the Vaigatch and Dolgyi islands, being less similar to those of the south tundra regions (Kanin Peninsula, Bolvanskaya Bay and Antipayuta village, Gydan Peninsula). Among the study habitats, the most diverse beetle assemblage inhabited thermokarst lakes (up to 14 species, the maximum density noted was 1150 ind./100 trap-days). On the contrary, on salt marshes, only 2–3 beetle species were recorded and the total abundance was significantly lower (1.8–3.6 ind./100 trap-days).

Uncertainty propagation in a global biogeochemical model driven by leaf area data

Satellite-observed leaf area index (LAI) is often used to depict vegetation canopy structure and photosynthesis processes in terrestrial biogeochemical models. However, it remains unclear how the uncertainty of LAI among different satellite products propagates to the modeling of carbon (C), nitrogen (N), and phosphorus (P) cycles. Here, we separately drive a global biogeochemical model by three satellite-derived LAI products (i.e., GIMMS LAI3g, GLASS, and GLOBMAP) from 1982 to 2011. Using a traceability analysis, we explored the propagation of LAI-driven uncertainty to modeled C, N, and P storage among different biomes. The results showed that the data uncertainty of LAI was more considerable in the tropics than in non-tropical regions, whereas the modeling uncertainty of C, N, and P stocks showed a contrasting biogeographic pattern. The spread of simulated C, N, and P storage derived by different LAI datasets resulted from assimilation rates of elements in shrubland and C3 grassland but from the element residence time (τ) in deciduous needle leaf forest and tundra regions. Moreover, the assimilation rates of elements are the main contributing factor, with 67.6, 93.2, and 93% of vegetated grids for the modeled uncertainty of C, N, and P storage among the three simulations. We further traced the variations in τ to baseline residence times of different elements and the environmental scalars. These findings indicate that the data uncertainty of plant leaf traits can propagate to ecosystem processes in global biogeochemical models, especially in non-tropical forests.

Moss and Liverwort Covers Structure Soil Bacterial and Fungal Communities Differently in the Icelandic Highlands.

Cryptogamic covers extend over vast polar tundra regions and their main components, e.g., bryophytes and lichens, are frequently the first visible colonizers of deglaciated areas. To understand their role in polar soil development, we analyzed how cryptogamic covers dominated by different bryophyte lineages (mosses and liverworts) influence the diversity and composition of edaphic bacterial and fungal communities as well as the abiotic attributes of underlying soils in the southern part of the Highlands of Iceland. For comparison, the same traits were examined in soils devoid of bryophyte covers. We measured an increase in soil C, N, and organic matter contents coupled with a lower pH in association with bryophyte cover establishment. However, liverwort covers showed noticeably higher C and N contents than moss covers. Significant changes in diversity and composition of bacterial and fungal communities were revealed between (a) bare and bryophyte-covered soils, (b) bryophyte covers and the underlying soils, and (c) moss and liverworts covers. These differences were more obvious for fungi than bacteria, and involved different lineages of saprotrophic and symbiotic fungi, which suggests a certain specificity of microbial taxa to particular bryophyte groups. In addition, differences observed in the spatial structure of the two bryophyte covers may be also responsible for the detected differences in microbial community diversity and composition. Altogether, our findings indicate that soil microbial communities and abiotic attributes are ultimately affected by the composition of the most conspicuous elements of cryptogamic covers in polar regions, which is of great value to predict the biotic responses of these ecosystems to future climate change.

First Evidence of Microplastic Occurrence in the Marine and Freshwater Environments in a Remote Polar Region of the Kola Peninsula and a Correlation with Human Presence

Microplastics (MPs) have even been detected in remote environments, including high-latitude regions, where human activities are restricted or strongly limited. We investigated the surface water of the bays of the Barents Sea and the freshwater lakes that are located close to and several kilometers from a year-round resident field station in the remote tundra region of the Kola Peninsula. The microplastics' presence in aquatic environments in this region has not been indicated yet. Microplastics were detected in all samples collected from the Barents Sea (<4800 items·m-3) and the lakes (<3900 items·m-3). Fibers made from polyethylene terephthalate (PET)-the most common thermoplastic polymer of the polyester family-and semi-synthetic cellulosic rayon were the most dominant. This indicated that the source of fiber contamination may come from protective clothes, ropes, ship equipment, and fishing nets. Small microplastics can spread through current and atmospheric transport. The Norwegian Current is likely responsible for the lack of correlations found between MP contamination and the distance from the field station between the studied bays of the Barents Sea. On the contrary, a significant correlation with human presence was observed in the concentration of microfibers in the water of the tundra lakes. The number of MP fibers decreased with an increase in the distance from the field station. This is the first study, to the best of our knowledge, that reports such a correlation in a remote region. We also discuss implications for animals. Our results show that even the most isolated ecosystems are not free from microplastic pollution.