Understory vegetation mediates permafrost active layer dynamics and carbon dioxide fluxes in open-canopy larch forests of northeastern Siberia

PLOS ONE, Mar 2018

Arctic ecosystems are characterized by a broad range of plant functional types that are highly heterogeneous at small (~1–2 m) spatial scales. Climatic changes can impact vegetation distribution directly, and also indirectly via impacts on disturbance regimes. Consequent changes in vegetation structure and function have implications for surface energy dynamics that may alter permafrost thermal dynamics, and are therefore of interest in the context of permafrost related climate feedbacks. In this study we examine small-scale heterogeneity in soil thermal properties and ecosystem carbon and water fluxes associated with varying understory vegetation in open-canopy larch forests in northeastern Siberia. We found that lichen mats comprise 16% of understory vegetation cover on average in open canopy larch forests, and lichen abundance was inversely related to canopy cover. Relative to adjacent areas dominated by shrubs and moss, lichen mats had 2–3 times deeper permafrost thaw depths and surface soils warmer by 1–2°C in summer and less than 1°C in autumn. Despite deeper thaw depths, ecosystem respiration did not differ across vegetation types, indicating that autotrophic respiration likely dominates areas with shrubs and moss. Summertime net ecosystem exchange of CO2 was negative (i.e. net uptake) in areas with high shrub cover, while positive (i.e. net loss) in lichen mats and areas with less shrub cover. Our results highlight relationships between vegetation and soil thermal dynamics in permafrost ecosystems, and underscore the necessity of considering both vegetation and permafrost dynamics in shaping carbon cycling in permafrost ecosystems.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Understory vegetation mediates permafrost active layer dynamics and carbon dioxide fluxes in open-canopy larch forests of northeastern Siberia

March Understory vegetation mediates permafrost active layer dynamics and carbon dioxide fluxes in open-canopy larch forests of northeastern Siberia Michael M. Loranty 0 1 2 Logan T. Berner 0 2 Eric D. Taber 0 1 2 Heather Kropp 0 1 2 Susan M. Natali 0 2 Heather D. Alexander 0 2 Sergey P. Davydov 0 2 Nikita S. Zimov 0 2 0 Funding: Colgate University ( 1 Department of Geography, Colgate University , Hamilton , NY United States of America, 2 School of Informatics , Computing, and Cyber Systems , Northern Arizona University , Flagstaff , AZ United States of America, 3 Woods Hole Research Center, Falmouth, MA United States of America, 4 Department of Forestry, Forest and Wildlife Research Center, Mississippi State University, Mississippi State, MS United States of America, 5 Northeast Science Station, Pacific Institute for Geography , Far East Branch , Russian Academy of Sciences, Cherskii, Republic of Sakha , Yakutia , Russia 2 Editor: Riikka Rinnan, University of Copenhagen , DENMARK Arctic ecosystems are characterized by a broad range of plant functional types that are highly heterogeneous at small (~1±2 m) spatial scales. Climatic changes can impact vegetation distribution directly, and also indirectly via impacts on disturbance regimes. Consequent changes in vegetation structure and function have implications for surface energy dynamics that may alter permafrost thermal dynamics, and are therefore of interest in the context of permafrost related climate feedbacks. In this study we examine small-scale heterogeneity in soil thermal properties and ecosystem carbon and water fluxes associated with varying understory vegetation in open-canopy larch forests in northeastern Siberia. We found that lichen mats comprise 16% of understory vegetation cover on average in open canopy larch forests, and lichen abundance was inversely related to canopy cover. Relative to adjacent areas dominated by shrubs and moss, lichen mats had 2±3 times deeper permafrost thaw depths and surface soils warmer by 1±2ÊC in summer and less than 1ÊC in autumn. Despite deeper thaw depths, ecosystem respiration did not differ across vegetation types, indicating that autotrophic respiration likely dominates areas with shrubs and moss. Summertime net ecosystem exchange of CO2 was negative (i.e. net uptake) in areas with high shrub cover, while positive (i.e. net loss) in lichen mats and areas with less shrub cover. Our results highlight relationships between vegetation and soil thermal dynamics in permafrost ecosystems, and underscore the necessity of considering both vegetation and permafrost dynamics in shaping carbon cycling in permafrost ecosystems. - 0732944, OPP-1044610, and PLR-1417700 to SN; and PLR-1304040 to HA. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. thaw is likely underway and poised to continue with climate warming. Increased permafrost thaw is important for a number of reasons. Among the most prominent is the large amount of organic carbon stored in permafrost [ 2 ] that will become vulnerable to decomposition and transfer to the atmosphere as greenhouse gasses when thawed [ 3 ], constituting a potentially large climate feedback [ 4 ]. Permafrost thaw may also alter local hydrology and nutrient availability, leading to changes in vegetation composition [ 5,6 ] that will result in multiple climate feedbacks related to altered land surface albedo and evapotranspiration [7]. In permafrost soils, the depth of the seasonally thawed active layer is often used as a diagnostic measure of permafrost status. Annual thaw depth measurements can be used to monitor permafrost responses to temporal variation in climate [ 8,9 ] while spatially distributed measurements can help to elucidate processes that underlie climate responses at scales ranging from ecosystems to continents [10±12]. In addition to serving as an indicator of permafrost status, thaw depth measurements characterize the portion of the soil column that is available for biological activity, including root growth, plant acquisition of water and nutrients, and decomposition of organic matter. Deepening of the active layer is often associated with enhanced ecosystem respiration (RECO [13]). Understanding both drivers and consequences of active layer dynamics and how they co-vary within and between ecosystems is critical for predicting carbon cycle responses to continued climate warming in permafrost ecosystems. In the simplest terms, active layer warming and permafrost thaw is driven by ground heat flux, which is governed by ecosystem influences on surface energy partitioning [ 14 ] and soil properties that influence heat transfer [ 15 ]. At local scales, surface energy exchange may be influenced by the effects of vegetation cover on ground temperature and moisture via radiation interception [ 16 ] and partitioning of sensible and latent heat fluxes [17±20], the insulating effects of snow cover [ 21 ], and changes in thermal conductivity associated with soil water content [ 22 ]. A potential consequence of such variability is persistent subsurface heterogeneity in active layer depths that, in addition to controlling the amount of unfrozen carbon, may also influence RECO via impacts on the distribution of heat and water within the soil column [ 23 ]. In this study, we examine spatial variability in active layer depth and ecosystem carbon fluxes associated with understory vegetation composition and biomass in an open-canopy larch forest in northeastern Siberia. Larch forests comprise a large portion of the continuous permafrost zone and have low canopy cover [ 24 ], meaning that understory vegetation plays a crucial role in regional permafrost and carbon dynamics [ 25 ]. Specifically, we address two research questions: 1) do thaw depth and soil temperature vary between common understory vegetation types? and 2) does ecosystem respiration differ among key vegetation types, and if so are these differences related to active layer depth? 2. Methods 2.1 Study site Our study was conducted in two small watersheds (Y3 & Y4 ~20km2; Fig 1) [ 26 ] at the Northeast Science Station (NESS) near the town of Cherskii (68Ê47'N, 161Ê20'E) along the Kolyma River in the Sakha (Yakutia) Republic of Russia on private land with permissions arranged via NESS. For the 1986±2015 period the mean annual temperature was -10ÊC, with mean January and July temperatures of -32ÊC and 13ÊC, respectively and average annual precipitation is 218 mm with 85 mm occurring as rain, and 133 mm as snow [ 27 ]. Mean summer temperatures increased ~1ÊC from 1938 to 2009, though there was no systematic change in annual precipitation over this period [ 28 ]. Forests in the area are comprised entirely of larch (Larix cajanderi). Understory vegetation includes deciduous (Betula nana exilis, Betula nana divartica, Salix spp, Alnus fruticosa) and 2 / 17 Fig 1. Map of the study area showing the Northeast Science Station and measurement locations. The yellow circle denotes the flux and soil temperature timeseries measurement site, and gray points indicate locations of lichen and canopy cover measurements. Sites in the western half of the map, including the flux plots, are in the Y4 watershed while those along the road in the eastern half are in the Y3 watershed. Stands in the Y3 watershed are approximately 50 years old and of relatively highdensity, while the stands in Y4 are older (~150 years) and have generally lower forests density. Background image is a WorldView2 true color composite from 21 August 2012 provided by the Polar Geospatial Center. evergreen (Vaccinium vitis idaea, Empetrum nigrum, and Ledum decumbens) shrubs, mosses (including Aulacomnium turgidum and Polytrichum spp), and lichens (including Cladonia rangiferina, Cetraria cuculata, and Stereocaulon tomentosum). Larch density is a function of fire history and ranges from 0.05 to ~4 tree m-2 across early to late successional stands (20±200 years; [ 29 ]). For this study, we characterized differences between lichen mats and adjacent shrub-moss patches and measured the distribution of lichen mats in stands throughout the Y3 and Y4 watershed. Here we note that the spatial distribution of vegetation can be strongly influenced by microtpography related to geomorphic dynamics [ 30 ] or permafrost dynamics such as ice-wedge polygons [ 31 ]. This does not appear to be the case for our study sties because there are virtually no vertical differences between lichen mats and adjacent shrub-moss patches in our study sites. Instead, it is likely that the distribution of understory vegetation, and specifically lichen mats are a function of interactive influences of overstory forest cover, variability in post-fire soil conditions, and species interactions. In this study we examine how understory vegetation varies, and what conditions vary with vegetation, rather than seeking to determine specific causal mechanisms. 2.2 Data collection To understand the distribution of understory vegetation across the Y3 and Y4 watersheds, we measured percent cover of understory vegetation along with larch canopy cover in 35 stands (Fig 1). At each stand, we established three parallel 20 m long transects spaced 8±10 m apart. 3 / 17 At the end of each transect, we visually estimated understory percent cover in 1 m2 plots (6 per stand). The plots were divided into four quadrats and in each we estimated deciduous and evergreen shrub, forb, graminoid, moss, lichen, and non-vegetated percent cover. Time constraints prevented us from using the point intercept method, which can be more accurate [ 32 ]. At the center of each transect, we measured canopy cover using a hemispherical densiometer. Within one low-density (0.06 ± 0.03 trees m-2) late succession (~178 year old) larch stand in the Y4 watershed we examined relationships between understory vegetation cover, permafrost thaw depth (TD), ecosystem thermal characteristics, and fluxes of CO2 and H2O. At this particular stand, understory vegetation cover is generally dominated by Betula nana middendorffi, atop a moss layer, interspersed with mats of Cladonia spp. ranging from 0.25±2.0 m2 in size. To examine seasonal differences in soil temperature at this stand, we installed HOBO Pendant temperature sensors (Onset Corp, Bourne, MA) at a depth of 10 cm in the soil beneath 0.25 m2 lichen mats (n = 9) and shrub/moss patches (n = 9). The sensors recorded hourly soil temperature (Tsoil) from July 2012 until June 2014. Here we also collected soil samples from lichen mats (n = 6) and shrub/moss patches (n = 6) in order to characterize soil properties. We used a soil saw to collect 10cm by 10cm to a depth of approximately 20cm. For each sample we recorded the depth of the organic horizon. Subsamples from the organic and mineral horizons were oven dried at 60ÊC for 48 hours to determine gravimetric soil moisture, and then placed in a muffle furnace at 450ÊC for 5 hours to quantify organic matter content using the loss-on-ignition method. In summer 2013, we measured CO2 and H2O fluxes and thaw depth at separate 0.25 m2 plots dominated by lichen (L; n = 5), low-density shrubs with moss understory (SM; n = 5), and high-density shrubs with little moss in the understory (S; n = 5). These flux plots were located approximately 20m away from those instrumented for Tsoil measurements and where the soil samples were collected. A total of 174 flux measurements (87 light and dark) were made at these plots between July 17 and August 5, 2013 on seven days with low wind and no rain. On each day we randomly selected the order in which plots were measured. At each plot, we measured CO2 and H2O fluxes using a LI-COR 840 infrared gas analyzer (IRGA; LI-COR Biosciences, Lincoln, NE) with a manual closed chamber system. A transparent acrylic chamber (50x50x50 cm) was used to measure evapotranspiration (ET) and net ecosystem exchange of CO2 (NEE), and an opaque cover was used to measure fluxes under dark conditions (i.e., RECO). A plastic skirt and chain were used to establish an airtight seal between the chamber and the ground. Concentrations of CO2 and H2O were recorded every second for approximately two minutes with a tablet connected to the IRGA. To calculate gas fluxes, we fit a slope to each set of concentration measurements, using only the linear portion of curve, omitting data without significant linear relationships [ 33 ]. At each plot, with each set of fluxes, we also measured soil temperature from 0-5cm depth (Tsoil) and air temperature (Tair) using a thermocouple (Fisher Scientific, Waltham, MA), photosynthetically active radiation (QSO-PAR Decagon Devices, Pullman, WA), surface soil moisture (GS-3 Decagon Devices, Pullman, WA), and radiometric surface temperature (Tsurf; Appogee Instruments, Logan, UT). We also measured TD on the north and south sides of each plot by inserting a graduated metal rod into the ground until firm resistance was met. Twice at the beginning of the field season we measured soil thermal conductivity (KS) integrated over 0±5 cm depth (KD-2, Decagon Devices, Pullman, WA). Half hourly values of Tair and air pressure logged to a HOBO Microstation (Onset Corp, Bourne, MA) were matched to the corresponding fluxes according to time of measurement. Tair from the meteorological station were used to fill chamber-level gaps due to instrument failure. On three occasions, we also measured the normalized difference vegetation index (NDVI) for each plot (SRS-Nr, Decagon Devices, Pullman, WA). At each flux plot, we measured the basal diameter for Betula spp. and 4 / 17 Salix spp. shrubs and then calculated aboveground biomass using regional allometric equations [ 34 ]. We also visually estimated percent cover of the dominant plant functional types. At the same low-density stand where fluxes were measured we quantified subsurface heterogeneity related to surface vegetation using electrical resistivity imaging [ 35,36 ]. We measured ERI along ~10m linear transects centered on lichen mats approximately 1±2 m in diameter. ERI is a minimally invasive technique applied by injecting a direct current into the ground and measuring differences in potential along a two-dimensional electrode array. The result of a resistivity survey is a series of point measures that depend upon electrode geometry and an assumption of subsurface homogeneity. An inversion model was used to calculate true resistivity values for the heterogeneous subsurface and generate resistivity images. Roll-along resistivity surveys were conducted using a Syscal Kid resistivity meter (Iris Instruments, Orleans, France) with 24 electrodes spaced 20 cm apart. The 30cm graphite electrodes were inserted in the ground at 45-degree angles in order to maximize electrode contact and minimize vertical distortion, and data were sampled using a Schlumberger sampling array. At 20cm intervals along each ~10 m long transect, we measured TD with a graduated metal rod. All resistivity surveys were inverted using Res2DInv software [ 37 ]. The high resistivity values and contrasts present in permafrost ecosystems introduce challenges with respect to data inversion [ 38 ]. In order to avoid over-fitting that leads to extremely high resistivity values [ 39 ], we limited the number of model iterations to five. We used a robust inversion because it more adequately captures sharp gradients characteristic of permafrost ecosystems [ 35 ]. Additionally, prior research shows resistivity inversions in permafrost ecosystems to be highly sensitive to a model smoothness parameter (λ) [ 38 ] that controls the degree to which the model is constrained by individual data points. We used the default initial value of 0.15 in Res2Dinv, and then confirmed our results by generating additional inversions with higher λ values in order to confirm that modeled features were real and not inversion artifacts [ 40 ]. The resistivity inversion uses the finite element method for forward modeling, where the surface of interest is discretized into a series of cells, and we varied the number of cell divisions between points as well. All data analyses were performed in R version 3.2.3 [ 41 ]. Preliminary data analyses were conducted to ensure that all data satisfied assumptions of normality associated with each statistical test. For time series Tsoil and soils data collected beneath lichen mats and shrub/moss patches we tested for differences in soil characteristics among vegetation types using two-sample t-tests. We used a one-way Analysis of Variance (ANOVA) with vegetation type as the fixed effect to examine differences in biological and physical variables related to surface energy partitioning, and also to carbon dynamics measured at the flux plots. A post-hoc Tukey's Honest Significant Difference test was used to determine whether differences between each vegetation type were significant. Ordinary least squared regressions were performed across all days to examine variability in RECO responses to soil and air temperature across each surface cover type (e.g. shrub, lichen, and shrub-lichen). All data and code used for analyses are freely available at the following url https://github.com/mloranty/lichen_pft/tree/v2.0 [ 33 ]. 3. Results 3.1 Vegetation distribution and active layer properties Lichen occurred widely in the study area and had notable impacts on resistivity, thaw depth, and soil temperature. Areal percent cover of lichen mats was 16.5 ± 2.5% (±1 SD) with a range of 3±33% across the 35 forest stands sampled in this watersheds. Mean percent cover of deciduous shrubs and moss were (37.6 ± 2.9% and 39.2 ± 4.1%, respectively). The distribution of lichen mats was not uniform, but rather exhibited a significant (p < 0.05) inverse relationship with larch canopy cover (Fig 2). The ERI tomography and thaw depth measurements in mid5 / 17 Fig 2. Relationship between larch canopy cover and lichen percent cover. Areal percent cover of lichen declined non-linearly with increasing larch canopy cover for 35 stands in the study area. Lichen cover was estimated visaually and larch canopy cover was measured with a hemispherical densitometer. July revealed low resistivity beneath lichen mats relative to adjacent areas dominated by shrubs and moss, indicative of deeper thaw and/or saturated conditions beneath lichen mats (Fig 3). Concurrent measurements of thaw depth aligned well with resistivity patterns, showing deeper thaw beneath lichen mats (Fig 3). At the intensively measured low-density stand soil temperature measured at 10 cm depth from July 2012 ±June 2014 revealed elevated temperatures beneath lichen mats relative to adjacent vegetation patches dominated by shrubs and mosses (Fig 4). These differences were largest during the early to mid growing season (~2ÊC), and gradually decreased towards the end of the growing season. Differences were minimal at the onset of fall freeze-back, but then increased again as soils beneath lichen mats spent longer time in the `zero curtain' period, after which differences were minimal during the winter. Aggregated seasonal temperatures (Table 1) were significantly different during the growing season (June-August; p < 0.01) and autumn (September±November), but not during winter (December-February) or spring (March-May). The soil organic layer was twice as thick (Table 2; p = 0.00001) beneath shrub/moss plots (12.3 ± 1.5 cm) in comparison to lichen mats (6.0 ±1.4 cm). Gravimetric soil moisture and soil 6 / 17 Fig 3. Surface and subsurface characteristics of understory lichen mats. Photograph showing lichen mat and resistivity survey near Cherskii, Siberia. Two representative resistivity profiles where colors indicate resistivity, black lines represent thaw depth measurements, and green boxes labeled lichen above each profile denote the location of lichen mats. Note the deep areas of low resistivity (blue) and thaw depth beneath lichen. organic matter content in the organic soil horizon were significantly higher beneath shrub/ moss patches relative to lichen mats (Table 2; p = 0.00008 and p = 0.0002 respectively); however, soil moisture and organic matter content within the mineral soil did not differ significantly between vegetation types (Table 2). 3.2 Vegetation influences on carbon and water dynamics Plots utilized for flux measurements exhibited clear differences in biological and physical characteristics among vegetation types (Table 3). Shrub plots had an average of 210 ± 46 g shrub aboveground biomass, which was significantly higher than the 69 ± 15 g in shrub-moss plots, while only one lichen mat had a shrub (0.9 ± 0.9 g). Percent moss cover did not differ significantly between shrub-moss plots (52 ± 11%) and shrub plots (34 ± 9%). Despite significant 7 / 17 Fig 4. Soil temperatures beneath lichen and shrub understory vegetation. (A) Time series of daily mean Tsoil at 10 cm depth beneath lichen mats (n = 9) and shrub patches (n = 9) from July 2012 ±June 2014. Dotted lines indicate one standard deviation. (B) The difference between daily mean lichen and shrub temperatures (lichen-shrub). differences in shrub biomass, NDVI did not differ significantly between shrub-moss (0.61 ± 0.01) and shrub plots (0.63 ± 0.01), though both were significantly higher (ANOVA, p < 0.05) than lichen mats (0.40 ± 0.01). Thermal conductivity was significantly higher (ANOVA, p < 0.05) in lichen mats (0.40 ± 0.06 W m-1 K-1) in comparison to shrub-moss (0.12 ± 0.02 W m-1 K-1) and shrub plots (0.07 ± 0.01 W m-1 K-1). Thaw depth below lichen mats was significantly higher (ANOVA, p < 0.05) at the beginning (July 19; 72 ± 4 cm) and end (Aug 3; 79 ± 3 cm) of the study period than either shrub-moss (26 ± 4 cm and 36 ± 6 cm) or shrub plots (34 ± 9 cm and 48 ± 8 cm). Differences in thermal conductivity and thaw depth between shrub-moss and shrub plots were not statistically significant. Lichen2 0.33 (0.11) -6.34 (1.46) -5.61 (0.49) 3.17 (0.48) 0.15 (0.13) -3.57 (0.85) -4.61 (0.52) Shrub/Moss 0.01 (0.23) -6.62 (1.56) -5.58 (0.59) 1.68 (0.47) -0.16 (0.28) -3.96 (0.85) -4.52 (0.53) 8 / 17 Organic Soil Layer Soil Moisture (%) 52.7 (5) 70.9 (3.6) Examination of variables related to surface energy partitioning may help elucidate drivers of observed differences in TD. To accomplish this, we present plot-level differences in Tsoil, Tair, Tsurf, PAR, and ET (latent heat flux) for two days with contrasting environmental conditions (Table 4). On 1 August 2013 all fifteen plots were sampled between 16:00 and 18:00 when temperature was moderate and PAR was low, and both were invariant across the sample period. During this time, the difference between Tsurf and Tair (Tdif) was slightly positive (warmer) on lichen mats while the shrub patches were approximately 1ÊC cooler than the air, but differences in Tsurf and Tdif between vegetation type were not statistically significant. On 3 August 2013 all plots were sampled between 11:30 and 13:30 when Tair and PAR were relatively high compared to 1 August. On this occasion, Tsurf was substantially higher than Tair across all vegetation types; with lichen Tsurf being on average 12.0 ± 2.6ÊC warmer than Tair whereas Tsurf was 6.3 ± 1.5ÊC and 4.8 ± 1.6ÊC warmer than Tair for shrub-moss and shrubs, respectively (Table 4). Differences in ET between vegetation types were not significantly different on either day. However, on August 3 when vegetation surfaces were substantially warmer than the air, lichen had the lowest ET and the highest Tsurf, while on August 1 there were no appreciable patterns in ET and small differences in Tsurf. On both days there were significant differences in Tsoil (ANOVA, p < 0.05) that were consistent with observed patterns of TD and time series observations of Tsoil at the site. Deeper thaw depths and warmer soil temperatures associated with lichen did not lead to higher RECO (Table 5). Across the study period, average RECO was not significantly different across vegetation types; 2.42±0.27 μmol m-2 sec-1, 2.92±0.27 μmol m-2 sec-1, and 2.91±0.36 μmol m-2 sec-1 for lichen, shrub-moss, and shrub plots, respectively. Average NEE across the study period was negative in shrub plots (-0.54±0.67 μmol m-2 sec-1), indicating net carbon uptake. This was significantly lower than NEE values of 2.01±0.58 μmol m-2 sec-1 and 2.46±0.30 μmol 1Letters indicate significant differences in measured variables between plot type 2Understory precent cover estimate includes moss cover beneath shrub canopies. 3Means are an average of three observations taken during the study, no temporal patterns were observed 4Lichen n = 5; Shrub-Moss n = 5; Shrub n = 5 9 / 17 1Letters indicate significant differences in measured variables between vegetation type. Variables without letters have no significant differences at p < 0.05 2Photosynthetically Active Radiation. 3 For each day Lichen n = 5; Shrub-Moss n = 5; Shrub n = 5 m-2 sec-1 in shrub-moss and lichen plots, respectively (ANOVA; p < 0.05). Individual measurements of RECO were linearly related to air temperature (adj r2 = 0.19, p < 0.01) but not 0±5 cm soil temperature (Fig 5). We did not observe significant relationships between RECO and air or soil temperatures within individual vegetation types. 4. Discussion 4.1 Vegetation and active layer dynamics Our results illustrate the influence of understory plant functional type on permafrost thaw depth, likely due to differences in surface energy partitioning. Specifically, we found that thaw depths were 2±3 times greater under lichen mats than under shrubs or mosses. Low moisture content and limited evapotranspiration from lichens prevent dissipation of latent heat [ 18 ], potentially allowing for deeper thaw. On the other hand mosses have a low thermal conductivity and high insulating capacity when dry, due to a high air volume [ 42,43 ]. Mosses can thus prevent thawing of permafrost by reducing the transfer of solar heat into the soil [ 20,44,45 ]. When damp, mosses have high rates of evapotranspiration that minimize sensible and ground heat fluxes, thus buffering against temperature variations [ 22 ] and reducing soil temperature [ 20,46 ]. Similarly, a combination of ground shading and relatively high latent heat fluxes may minimize permafrost thaw beneath shrub canopies [16]. The absence of a root system in mosses and lichens limits soil water loss, and may lead to the development of saturated soil [ 44 ]. A study by Stoy et al [ 23 ] found higher soil temperatures (1 cm and 5 cm depth) beneath 1Letters indicate significant differences in measured variables between vegetation type at p < 0.05. Variables without letters have no significant differences. 2Photosynthetically Active Radiation 3Net Ecosystem Exchange of CO2. 4Lichen n = 29; Shrub-Moss n = 29; Shrub n = 29 10 / 17 Fig 5. Variation in ecosystem respiration with air and soil temperatures. (A) RECO was not significantly related to Tsoil. (B) RECO was was positively related to air temperature when observations were polled across plant communities and sampling periods (intercept = 1.16, slope = 0.09, adj r2 = 0.19, p < 0.01). lichen and feather moss relative to nearby sphagnum, however, modeled ground heat flux for lichen was relatively low in comparison to mosses. Our study is broadly inline with these previous results; we found that lichen mats had both thinner organic layers with less organic matter and correspondingly higher thermal conductivity, and also higher surface temperatures associated with lower latent heat dissipation. Both of these sets of factors are likely to contribute to warmer soils and deeper active layers beneath lichen mats. It is also important to consider that variability in micortopography and soil moisture can also influence vegetation distribution [ 30,31,47 ]. This is commonly observed in areas with icewedge poloygons, where permafrost dynamics create micortopography with low wet areas adjacent to higher drier areas. Elevation differences in these cases are typically on the order of 10s of cm and can be the dominant controls on moisture [ 48 ], which in turn impacts vegetation community distribution. This does not seem to be the case at our sites as there is no microtopography of this nature. However it is still possible that, in addition to overstory forest cover, lichen distribution is partially determined by post-fire soil conditions that have persisted to the present. In this case the differences in thermal properties that we observed may actually reinforce differences in active layer properties partially responsible for vegetation distribution. Fire typically alters active layer depths for up to fifty years post-fire [ 49,50 ], and ecosystem recovery typically promotes recovery of the permafrost table [ 51,52 ]. In our study sites the high-density stands are among the youngest (~50yrs), whereas the lower density sites with higher lichen abundance are upwards of 150 years old, indicating that our estimates of active layer heterogeneity are not likely the results of fire legacy effects. In any case, improved understanding of active layer heterogeneity in Siberian larch forests requires comprehensive patchscale observations of understory vegetation. 11 / 17 4.2 Thaw depth and ecosystem respiration Previous studies suggest that increased thaw depth can lead to greater carbon loss through increased heterotrophic respiration if soil carbon and quality with depth [17±20,53±55]; however, we found no relationship between thaw depth and RECO. Autotrophic respiration can comprise a substantial portion of RECO in permafrost ecosystems [ 13 ]. It is likely that higher autotrophic respiration compensated for lower heterotrophic respiration in shrub plots, which had higher aboveground biomass and shallower permafrost thaw depths than other vegetation types. The observed dependence of RECO on air temperature rather than soil temperature supports this idea. Additionally, differences in the quantity and quality of soil carbon beneath different vegetation types may also contribute to variability in RECO, and our observations of thicker organic soil layers with higher organic matter content suggest soil carbon differences may have been a factor in this study. For example, the low biomass, slow turnover and low litter input of lichens could limit the amount of soil carbon available for decomposition in lichen patches relative to areas dominated by mosses or shrubs [ 56 ]. Thus, decreased soil organic matter may negate carbon flux from increased thaw depth. Conversely, high soil organic matter content of soil beneath shrub patches may enhance RECO from moderate thaw depths. Alternatively, labile carbon may have already decomposed beneath lichen patches with deeper thaw depths. Spatial variation in the amount and lability of soil C are also likely influenced by fire, which is the dominant disturbance in the region [ 29,57 ]. Soil moisture may also influence relationships between thaw depth and RECO via moisture limitations on decomposition [58±60] that can be regulated by vegetation [ 44,58,61 ]. Several studies in Alaska have identified accumulation of water at the base of the active layer as a potential explanation for lack of observed relationships between thaw depth and RECO [ 21,53 ], and this could be plausible at our study site as well. Similarly, Zona et al [62] observed higher RECO at micro sites with shallower thaw depth and lower water tables. In these cases, suppression of RECO by soil moisture may be offset by CH4 efflux associated with anaerobic respiration [ 63,64 ]. Our observations of delayed autumn freeze-back suggest a higher amount of soil moisture beneath lichen mats relative to adjacent areas dominated by shrubs and mosses; however, this does not necessarily correspond to saturated soil at the base of the active layer and could also result simply from deeper thaw depths. Other studies have observed spatial variability in CH4 across different vegetation types [65], but it is unclear if CH4 efflux varies with thaw depth at our study site. Our measurements of RECO are snapshots taken during the peak of the growing season; however, non-growing season fluxes often determine sign and magnitude of the annual carbon balance for an ecosystem [ 53,66 ]. Thus year round measurements of CO2 and CH4 are required to understand the effects of vegetation-mediated active layer dynamics on the understory carbon fluxes in Siberian larch forests. 4.3 Implications for ecosystem change The relationship between plant functional type and thaw depth supports the assertion that vegetation-mediated variability in surface energy partitioning may alter soil thermal and biogeochemical dynamics in permafrost ecosystems [ 23 ]. Lichen constitutes upwards of 16% of understory aboveground biomass in low-density larch stands in northeastern Siberia [ 29 ] and 8±32% of aboveground biomass among a network of tundra sites in western Siberia [ 67 ]. Accounting for this variability will be necessary for accurate estimates of carbon dynamics in these ecosystems. The same is likely to be true for other permafrost ecosystems, assuming the occurrence of similar relationships between thaw depth and plant functional type, particularly where changes in herbivore browsing patterns (e.g. Rangifer tarandus; [ 68 ]), fire [ 69 ] or replacement by other forms of vegetation [ 70,71 ] lead to rapid changes in lichen distribution. 12 / 17 In this context, it will be particularly important to determine whether depressed RECO associated with lichen cover is a function of saturated soils, low rates of heterotrophic respiration, or some combination of both. At the landscape scale, the type of subsurface heterogeneity that we observed may translate to bidirectional responses of RECO to interannual climate variability e.g. [ 72 ], potentially leading to hotspots for CO2 and CH4 emissions [ 64 ] that vary spatially from year to year. This sort of space-time variation may present challenges in modeling future greenhouse gas emissions from permafrost ecosystems. On the other hand, high-resolution spectral data combined with observed relationships between surface vegetation and thaw depth in these ecosystems may lead to improved understanding of ecosystem-scale variability in active layer dynamics [73±76]. 5. Conclusions In open canopy larch forests in northeastern Siberia, variations in understory vegetation over short distances (e.g. < 1 m) correspond to large differences in thaw depth. Our results illustrate a strong interactions between vegetation and active layer dynamics in these ecosystems, where lower latent heat fluxes and higher thermal conductivity in lichen mats lead to deeper thaw depths. These vegetation types have clear spectral differences, and offer the possibility for remote detection of active layer heterogeneity. Despite increased thaw depth beneath lichen mats, we did not observe elevated ecosystem respiration. Here, soil moisture, along with lower soil carbon content may help to explain the absence of elevated ecosystem respiration, relative to adjacent areas with shallower thaw depths. A better understanding of this subsurface variability, particularly thaw depth controls on soil moisture redistribution, will be necessary to accurately quantify the effects of permafrost thaw on ecosystem carbon cycling. To this end, our results indicate the usefulness of electrical resistivity imaging for visualizing active layer heterogeneity. In Siberian larch forests, fire controls on tree density and understory vegetation succession will exert strong controls on variation in active layer and ecosystem carbon dynamics under current and future climates. Acknowledgments We thank Ellen Squires, numerous other Polaris Project student researchers, and NESS personnel for assistance with field work. Author Contributions Conceptualization: Michael M. Loranty, Logan T. Berner, Eric D. Taber, Susan M. Natali, Heather D. Alexander, Sergey P. Davydov. Data curation: Michael M. Loranty. Formal analysis: Michael M. Loranty, Heather Kropp. Funding acquisition: Michael M. Loranty, Susan M. Natali, Heather D. Alexander. Investigation: Michael M. Loranty, Logan T. Berner, Eric D. Taber, Susan M. Natali, Heather D. Alexander, Sergey P. Davydov, Nikita S. Zimov. Resources: Michael M. Loranty, Sergey P. Davydov, Nikita S. Zimov. Writing ± original draft: Michael M. Loranty. Writing ± review & editing: Logan T. Berner, Eric D. Taber, Heather Kropp, Susan M. Natali, Heather D. Alexander, Nikita S. Zimov. 13 / 17 14 / 17 15 / 17 16 / 17 1. Romanovsky VE , Smith SL , Christiansen HH . Permafrost thermal state in the polar Northern Hemisphere during the international polar year 2007 ± 2009 : a synthesis . Permafrost Periglac Process . 2010 ; 21 : 106 ± 116 . https://doi.org/10.1002/ppp.689 2. Hugelius G , Strauss J , Zubrzycki S , Harden JW , Schuur EAG , Ping CL , et al. Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps . Biogeosciences . 2014 ; 11 : 6573 ± 6593 . https://doi.org/10.5194/bg-11- 6573 -2014 3. Harden JW , Koven CD , Ping C-L , Hugelius G , David McGuire A , Camill P , et al. Field information links permafrost carbon to physical vulnerabilities of thawing . Geophys Res Lett . 2012 ; 39: n/a±n/a . https:// doi.org/10.1029/2012GL051958 4. Schuur EAG , Mcguire AD , SchaÈdel C , Grosse G , Harden JW , Hayes DJ , et al. Climate change and the permafrost carbon feedback . Nature . 2015 ; 520 : 171 ± 179 . https://doi.org/10.1038/nature14338 PMID: 25855454 5. Schuur EAG , Crummer KG , Vogel JG , Mack MC . Plant Species Composition and Productivity following Permafrost Thaw and Thermokarst in Alaskan Tundra . Ecosystems. 2007 ; 10 : 280 ± 292 . https://doi.org/ 10.1007/s10021-007-9024-0 6. Jorgenson MT , HARDEN J , Kanevskiy M. Reorganization of vegetation, hydrology and soil carbon after permafrost degradation across heterogeneous boreal landscapes . Env Res Lett . 2013 . https://doi.org/ 10.1088/ 1748 -9326/8/3/035017 7. Pearson RG , Phillips SJ , Loranty MM , Beck PSA , Damoulas T , Knight SJ , et al. Shifts in Arctic vegetation and associated feedbacks under climate change . Nature Climate Change . 2013 ; 3 : 673 ± 677 . https://doi.org/10.1038/nclimate1858 8. Shiklomanov NI , Streletskiy DA , Nelson FE , Hollister RD , Romanovsky VE , Tweedie CE , et al. Decadal variations of active-layer thickness in moisture-controlled landscapes , Barrow, Alaska. J Geophys Res . 2010 ; 115 : G00I04. https://doi.org/10.1029/2009JG001248 9. Shiklomanov NI , Streletskiy DA , Neloson FE . Northern Hemisphere component of the Global Circumpolar Active Layer Monitoring (CALM) Program . 2012 . 10. Streletskiy DA , Sherstiukov AB , Frauenfeld OW , Nelson FE . Changes in the 1963±2013 shallow ground thermal regime in Russian permafrost regions . Environ Res Lett. IOP Publishing ; 2015 ; 10 : 1± 10 . https://doi.org/10.1088/ 1748 -9326/10/12/125005 11. Nelson F , Hinkel K , Shiklomanov N , Mueller G , Miller LL , Walker DA . Active-layer thickness in north central Alaska: Systematic sampling, scale, and spatial autocorrelation . J Geophys Res . 1998 ; 103 : 28963 ± 28973 . 12. Hinkel KM , Nelson FE . Spatial and temporal patterns of active layer thickness at Circumpolar Active Layer Monitoring (CALM) sites in northern Alaska , 1995 ± 2000 . J Geophys Res . 2003 ; 108 . https://doi. org/10.1029/2001jd000927 13. Pries CEH , Schuur EAG , Crummer KG . Thawing permafrost increases old soil and autotrophic respiration in tundra: Partitioning ecosystem respiration using δ13C and Δ14C . Global Change Biol . 2013 ; 19 : 649 ± 661 . https://doi.org/10.1111/gcb.12058 PMID: 23504799 14. Beringer J , Chapin FS , Thompson CC , Mcguire AD . Surface energy exchanges along a tundra-forest transition and feedbacks to climate. Agr Forest Meteorol . Elsevier; 2005 ; 131 : 143 ± 161 . https://doi.org/ 10.1016/j.agrformet. 2005 . 05 .006 15. O 'Donnell JA , Romanovsky VE , Harden JW , Mcguire AD . The Effect of Moisture Content on the Thermal Conductivity of Moss and Organic Soil Horizons From Black Spruce Ecosystems in Interior Alaska . Soil Science . 2009 ; 174 : 646 ± 651 . https://doi.org/10.1097/ss.0b013e3181c4a7f8 16. Blok D , Heijmans M , Schaepman-Strub G , Kononov A , Maximov T , Berendse F . Shrub expansion may reduce summer permafrost thaw in Siberian tundra . Global Change Biol . Wiley Online Library; 2010 ; 16 : 1296 ± 1305 . 17. Heijmans MMPD , Arp WJ , Chapin FS III. Carbon dioxide and water vapour exchange from understory species in boreal forest . Agr Forest Meteorol . 2004 ; 123 : 135 ± 147 . https://doi.org/10.1016/j.agrformet. 2003 . 12 .006 18. Lafleur PM , Schreader CP . Water loss from the floor of a subarctic forest . Arctic and Alpine Research . 1994 ; 26 : 152 . https://doi.org/10.2307/1551778 19. Lantz TC , Gergel SE , Kokelj SV . Spatial Heterogeneity in the Shrub Tundra Ecotone in the Mackenzie Delta Region, Northwest Territories: Implications for Arctic Environmental Change . Ecosystems . 2010 ; : 1 ± 11 . https://doi.org/10.1007/s10021-009-9310-0 20. Blok D , Heijmans MMPD , Schaepman-Strub G , Ruijven J , Parmentier FJW , Maximov TC , et al. The Cooling Capacity of Mosses: Controls on Water and Energy Fluxes in a Siberian Tundra Site . Ecosystems. 2011 ; 14 : 1055 ± 1065 . https://doi.org/10.1007/s10021-011-9463-5 21. Natali SM , Schuur EAG , Trucco C , Hicks Pries CE , Crummer KG , Baron Lopez AF . Effects of experimental warming of air, soil and permafrost on carbon balance in Alaskan tundra . Global Change Biol . 2011 ; 17 : 1394 ± 1407 . https://doi.org/10.1111/j.1365- 2486 . 2010 . 02303 .x 22. Romanovsky VE , Osterkamp TE . Effects of unfrozen water on heat and mass transport processes in the active layer and permafrost . Permafrost Periglac Process . Wiley Online Library; 2000 ; 11 : 219 ± 239 . 23. Stoy PC , Street LE , Johnson AV , Prieto-Blanco A , Ewing SA . Temperature, Heat Flux, and Reflectance of Common Subarctic Mosses and Lichens under Field Conditions: Might Changes to Community Composition Impact Climate-Relevant Surface Fluxes? Arctic, Antarctic, and Alpine Research. 2012 ; 44 : 500 ± 508 . https://doi.org/10.1657/1938-4246- 44 .4. 500 24. Loranty MM , Liberman-Cribbin W , Berner LT , Natali SM , Goetz SJ , Alexander HD , et al. Spatial variation in vegetation productivity trends, fire disturbance, and soil carbon across arctic-boreal permafrost ecosystems . Environ Res Lett. IOP Publishing ; 2016 ; 11 : 1± 13 . https://doi.org/10.1088/ 1748 -9326/11/ 9/095008 25. Xue B-L , Kumagai T , Iida S , Nakai T , Matsumoto K , Komatsu H , et al. Influences of canopy structure and physiological traits on flux partitioning between understory and overstory in an eastern Siberian boreal larch forest . Ecol Model . Elsevier B.V; 2011 ; 222 : 1479 ± 1490 . https://doi.org/10.1016/j. ecolmodel. 2011 . 01 .021 26. Webb EE , Heard K , Natali SM , Bunn AG , Alexander HD , Berner LT , et al. Variability in Above and Belowground Carbon Stocks in a Siberian Larch Watershed . Biogeosciences Discuss . 2017 ; : 1 ± 39 . https://doi.org/10.5194/bg-2017 -88 27. Menne MJ , Durre I , Vose RS , Gleason BE , Houston TG . An Overview of the Global Historical Climatology Network-Daily Database . Journal of Atmospheric and Oceanic Technology . 2012 ; 29 : 897 ± 910 . https://doi.org/10.1175/ JTECH-D- 11- 00103 . 1 28. Berner LT , Beck PSA , Bunn AG , Goetz SJ . Plant response to climate change along the forest-tundra ecotone in northeastern Siberia . Global Change Biol . 2013 ; 19 : 3449 ± 3462 . https://doi.org/10.1111/ gcb.12304 PMID: 23813896 29. Alexander HD , Mack MC , Goetz SJ , Loranty MM , Beck PS , Earl K , et al. Carbon Accumulation Patterns During Post-Fire Succession in Cajander Larch (Larix cajanderi) Forests of Siberia . Ecosystems. Springer; 2012 ; 15 : 1065 ± 1082 . https://doi.org/10.1007/s10021-012-9567-6 30. Peterson KM , Billings WD . Tundra Vegetational Patterns and Succession in Relation to Microtopography near Atkasook, Alaska . Arctic and Alpine Research . 1980 ; 12 : 473 . https://doi.org/10.2307/ 1550495 31. Gamon JA , Kershaw GP , Williamson S , Hik DS . Microtopographic patterns in an arctic baydjarakh field: do fine-grain patterns enforce landscape stability? Environ Res Lett . 2012 ; 7 : 015502 . https://doi.org/ 10.1088/ 1748 -9326/7/1/015502 32. GodÂõnez-Alvarez H , Herrick JE , Mattocks M , Toledo D , Van Zee J. Comparison of three vegetation monitoring methods: Their relative utility for ecological assessment and monitoring . Ecological Indicators . 2009 ; 9 : 1001 ± 1008 . https://doi.org/10.1016/j.ecolind. 2008 . 11 .011 33. Loranty MM . Larch understory vegetation and code . GitHub Repository . 2nd ed. 2018 . https://doi.org/ 10.5281/zenodo.1146037 34. Berner LT , Alexander HD , Loranty MM , Ganzlin P , Mack MC , Davydov SP , et al. Biomass allometry for alder, dwarf birch, and willow in boreal forest and tundra ecosystems of far northeastern Siberia and north-central Alaska. Forest Ecol Manag . Elsevier B.V; 2015 ; 337 : 110 ± 118 . https://doi.org/10.1016/j. foreco. 2014 . 10 .027 35. Lewkowicz AG , EtzelmuÈller B , Smith SL . Characteristics of Discontinuous Permafrost based on Ground Temperature Measurements and Electrical Resistivity Tomography , Southern Yukon, Canada. Permafrost Periglac Process. 2011 ; 22 : 320 ± 342 . https://doi.org/10.1002/ppp.703 36. Dafflon B , Oktem R , Peterson J , Ulrich C , Tran AP , Romanovsky V , et al. Coincident aboveground and belowground autonomous monitoring to quantify covariability in permafrost, soil, and vegetation properties in Arctic tundra . Journal of Geophysical Research Biogeosciences . 2017 ; 122 : 1321 ± 1342 . https:// doi.org/10.1002/2016JG003724 37. Loke MH , Barker RD . Rapid least-squares inversion of apparent resistivity pseudosections by a quasiNewton method1 . Geophys Prospect . 1996 ; 44 : 131 ± 152 . https://doi.org/10.1111/j.1365- 2478 . 1996 . tb00142.x 38. Hauck C , MuÈhll DV. Inversion and interpretation of two-dimensional geoelectrical measurements for detecting permafrost in mountainous regions . Permafrost Periglac Process . 2003 ; 14 : 305 ± 318 . https:// doi.org/10.1002/ppp.462 39. Kneisel C , Hauck C. Electrical Methods . In: Hauck C , Kneisel C , editors. Applied Geophysics in Periglacial Environments . Cambridge: Cambridge University Press; 2008 . pp. 1 ± 24 . 40. Hauck C , Vonder MuÈhll D , Maurer H. Using DC resistivity tomography to detect and characterize mountain permafrost . Geophys Prospect . Wiley Online Library; 2003 ; 51 : 273 ± 284 . 41. R Development Core Team. R: A language and environment for statistical computing . R Foundation for Statistical Computing . Vienna, Austria: R Foundation for Statistical Computing; 2014 . Available: http:// www.R-project.org/ 42. Beringer J , Lynch A , Chapin F III , Mack M , Bonan G. The representation of arctic soils in the land surface model: the importance of mosses . Journal of Climate . 2001 ; 14 : 3324 ± 3335 . 43. O 'Donnell JA , Romanovsky VE , Harden JW , Mcguire AD . The Effect of Moisture Content on the Thermal Conductivity of Moss and Organic Soil Horizons From Black Spruce Ecosystems in Interior Alaska . Soil Science . 2009 ; 174 : 646 ± 651 . https://doi.org/10.1097/SS.0b013e3181c4a7f8 44. Zimov SA , Chuprynin VI , Oreshko AP , Chapin FS III, Reynolds JF , Chapin MC . Steppe-tundra transition: a herbivore-driven biome shift at the end of the Pleistocene. American Naturalist . JSTOR; 1995 ;: 765 ± 794 . 45. Gornall JL , Jonsdottir IS , Woodin SJ , Van Der Wal R . Arctic mosses govern below-ground environment and ecosystem processes . Oecologia . 2007 ; 153 : 931 ± 941 . https://doi.org/10.1007/s00442-007 -0785- 0 PMID: 17618466 46. Heijmans MM , Arp WJ , Chapin FS . Controls on moss evaporation in a boreal black spruce forest . Wiley Online Library; 2004 ; 18. https://doi.org/10.1029/2003GB002128 47. Minke M , Donner N , Karpov N , de Klerk P , Joosten H. Patterns in vegetation composition, surface height and thaw depth in polygon mires in the Yakutian Arctic (NE Siberia): a microtopographical characterisation of the active layer . Permafrost Periglac Process . 2009 ; 20 : 357 ± 368 . https://doi.org/10. 1002/ppp.663 48. Engstrom R , Hope A , Kwon H , Stow D , Zamolodchikov D. Spatial distribution of near surface soil moisture and its relationship to microtopography in the Alaskan Arctic coastal plain . Hydrol Res . 2005 ; 36 : 219 ± 234 . 49. Burn CR . The response ( 1958 ± 1997 ) of permafrost and near-surface ground temperatures to forest fire, Takhini River valley, southern Yukon Territory . Can J Earth Sci . 1998 ; 35 : 184 ± 199 . https://doi.org/ 10.1139/e97- 105 50. Yoshikawa K , Bolton WR , Romanovsky VE , Fukuda M , Hinzman LD . Impacts of wildfire on the permafrost in the boreal forests of Interior Alaska . J Geophys Res . 2003 ; 108 : 8148 . https://doi.org/10.1029/ 2001JD000438 51. Jafarov EE , Romanovsky VE , Genet H , David McGuire A , Marchenko SS . The effects of fire on the thermal stability of permafrost in lowland and upland black spruce forests of interior Alaska in a changing climate . Environ Res Lett . 2013 ; 8 : 035030 . https://doi.org/10.1088/ 1748 -9326/8/3/035030 52. Smith SL , Riseborough DW , Bonnaventure PP. Eighteen Year Record of Forest Fire Effects on Ground Thermal Regimes and Permafrost in the Central Mackenzie Valley , NWT , Canada. Permafrost Periglac Process. 2015 ; : n/a±n/a . https://doi.org/10.1002/ppp.1849 53. Trucco C , Schuur EAG , Natali SM , Belshe EF , Bracho R , Vogel J . Seven-year trends of CO 2exchange in a tundra ecosystem affected by long-term permafrost thaw . J Geophys Res . 2012 ; 117 : G02031. https://doi.org/10.1029/2011JG001907 54. Schuur EAG , Vogel JG , Crummer KG , Lee H , Sickman JO , Osterkamp TE . The effect of permafrost thaw on old carbon release and net carbon exchange from tundra . Nature . 2009 ; 459 : 556 ± 559 . https:// doi.org/10.1038/nature08031 PMID: 19478781 55. Vogel J , Schuur EAG , Trucco C , Lee H . Response of CO 2exchange in a tussock tundra ecosystem to permafrost thaw and thermokarst development . J Geophys Res . 2009 ; 114 : G04018. https://doi.org/10. 1029/2008JG000901 56. Shaver GR , Chapin FS III. Production: biomass relationships and element cycling in contrasting arctic vegetation types . Ecological Monographs. Eco Soc America ; 1991 ; 61 : 1± 31 . 57. Berner LT , Beck PSA , Loranty MM , Alexander HD , Mack MC , Goetz SJ . Cajander larch (Larix cajanderi) biomass distribution, fire regime and post-fire recovery in northeastern Siberia . Biogeosciences . 2012 ; 9 : 3943 ± 3959 . https://doi.org/10.5194/bg-9- 3943 -2012 58. Howard DM , Howard P . Relationships between CO2 evolution, moisture content and temperature for a range of soil types . Soil Biology and Biochemistry . 1993 ; 25 : 1537 . 59. Shaver GR , Giblin AE , Nadelhoffer KJ , Thieler KK , Downs MR , Laundre JA , et al. Carbon turnover in Alaskan tundra soils: effects of organic matter quality, temperature, moisture and fertilizer . J Ecology . Wiley Online Library; 2006 ; 94 : 740 ± 753 . https://doi.org/10.1111/j.1365- 2745 . 2006 . 01139 .x 60. Oberbauer SF , Tweedie CE , Welker JM , Fahnestock JT , Henry GH , Webber PJ , et al. Tundra CO2 fluxes in response to experimental warming across latitudinal and moisture gradients . Ecological Monographs. Eco Soc America ; 2007 ; 77 : 221 ± 238 . 61. Elberling B , Nordstrøm C , Grøndahl L , Søgaard H , Friborg T , Christensen T , et al. High-Arctic Soil CO2 and CH4 Production Controlled by Temperature, Water, Freezing and Snow . Advances in Ecological Research . 2008 ; 40 : 441 ± 472 . 62. Zona D , Lipson DA , Zulueta RC , Oberbauer SF , Oechel WC . Microtopographic controls on ecosystem functioning in the Arctic Coastal Plain . J Geophys Res . 2011 ; 116 : G00I08. https://doi.org/10.1029/ 2009JG001241 63. Zona D , Oechel WC , Kochendorfer J , Paw U KT , Salyuk AN , Olivas PC , et al. Methane fluxes during the initiation of a large-scale water table manipulation experiment in the Alaskan Arctic tundra . 2009 ; 23 : 1± 11 . https://doi.org/10.1029/2009GB003487 64. Natali SM , Schuur EAG , Mauritz M , Schade JD , Celis G , Crummer KG , et al. Permafrost thaw and soil moisture driving CO 2 and CH 4 release from upland tundra . Journal of Geophysical Research Biogeosciences . 2015 ; 120 : 525 ± 537 . https://doi.org/10.1002/2014jg002872 65. Davidson SJ , Sloan VL , Phoenix GK , Wagner R , Fisher JP , Oechel WC , et al. Vegetation Type Dominates the Spatial Variability in CH4 Emissions Across Multiple Arctic Tundra Landscapes . Ecosystems. Springer US; 2016 ;: 1 ± 17 . https://doi.org/10.1007/s10021-016-9991-0 66. Commane R , Lindaas J , Benmergui J , Luus KA , Chang RYW , Daube BC , et al. Carbon dioxide sources from Alaska driven by increasing early winter respiration from Arctic tundra . P Natl Acad Sci Usa . 2017 ; 114 : 5361 ± 5366 . https://doi.org/10.1073/pnas.1618567114 PMID: 28484001 67. Walker DA , Epstein HE , Raynolds MK , Kuss P , Kopecky MA , Frost GV , et al. Environment, vegetation and greenness (NDVI) along the North America and Eurasia Arctic transects . Environ Res Lett . 2012 ; 7 : 015504 ± 18 . https://doi.org/10.1088/ 1748 -9326/7/1/015504 68. Joly K , Jandt RR , Klein DR . Decrease of lichens in Arctic ecosystems: the role of wildfire, caribou, reindeer, competition and climate in north-western Alaska . Polar Research . 2009 ; 28 : 433 ± 442 . https://doi. org/10.1111/j.1751- 8369 . 2009 . 00113 .x 69. Jandt R , Joly K , Randy Meyers C , Racine C. Slow Recovery of Lichen on Burned Caribou Winter Range in Alaska Tundra: Potential Influences of Climate Warming and Other Disturbance Factors . Arctic, Antarctic, and Alpine Research. 2008 ; 40 : 89 ± 95 . https://doi.org/10.1657/ 1523 - 0430 ( 06 ±122) [JANDT] 2.0 .CO; 2 70. Walker M , Wahren C , Hollister R , Henry G , Ahlquist LE , Alatalo JM , et al. Plant community responses to experimental warming across the tundra biome . Proceedings of the National Academy of Sciences . 2006 ; 103 : 1342 ± 1346 . 71. Fraser RH , Lantz TC , Olthof I , Kokelj SV , Sims RA . Warming-Induced Shrub Expansion and Lichen Decline in the Western Canadian Arctic . Ecosystems. 2014 ; 17 : 1151 ± 1168 . https://doi.org/10.1007/ s10021-014-9783-3 72. Riveros-Iregui DA , Mcglynn BL , Emanuel RE , Epstein HE . Complex terrain leads to bidirectional responses of soil respiration to inter-annual water availability . Global Change Biol . 2011 ; 18 : 749 ± 756 . https://doi.org/10.1111/j.1365- 2486 . 2011 . 02556 .x 73. TheÂau J , Peddle DR , Duguay CR . Mapping lichen in a caribou habitat of Northern Quebec, Canada, using an enhancement_classification method and spectral mixture analysis . Remote Sens Environ . Elsevier; 2005 ; 94 : 232 ± 243 . https://doi.org/10.1016/j.rse. 2004 . 10 .008 74. Korpela IS . Mapping of understory lichens with airborne discrete-return LiDAR data . Remote Sens Environ . Elsevier; 2008 ; 112 : 3891 ± 3897 . 75. Regmi P , Grosse G , Jones MC , Jones BM , Anthony KW . Characterizing post-drainage succession in thermokarst lake basins on the Seward Peninsula, Alaska with TerraSAR-X backscatter and Landsatbased NDVI data . Remote Sensing . Multidisciplinary Digital Publishing Institute; 2012 ; 4 : 3741 ± 3765 . https://doi.org/10.3390/rs4123741 76. Dafflon B , Oktem R , Peterson J , Ulrich C , Tran AP , Romanovsky V , et al. Coincident aboveground and belowground autonomous monitoring to quantify covariability in permafrost, soil, and vegetation properties in Arctic tundra . Journal of Geophysical Research Biogeosciences . 2017 ; 164 : 322 ± 22 . https://doi. org/10.1002/2016JG003724

This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0194014&type=printable

Michael M. Loranty, Logan T. Berner, Eric D. Taber, Heather Kropp, Susan M. Natali, Heather D. Alexander, Sergey P. Davydov, Nikita S. Zimov. Understory vegetation mediates permafrost active layer dynamics and carbon dioxide fluxes in open-canopy larch forests of northeastern Siberia, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0194014