Are vegetation influences on Arctic–boreal snow melt rates detectable across the Northern Hemisphere?

The timing and rate of northern high latitude spring snowmelt plays a critical role in surface albedo, hydrology, and soil carbon cycling. Ongoing changes in the abundance and distribution of trees and shrubs in tundra and boreal ecosystems can alter snowmelt via canopy impacts on surface energy partitioning. It is unclear whether vegetation-related processes observed at the ecosystem scale influence snowmelt patterns at regional or continental scales. We examined the influence of vegetation cover on snowmelt across the boreal and Arctic region across a ten-year reference period (2000–2009) using a blended snow water equivalent (SWE) data product and gridded estimates of surface temperature, tree cover, and land cover characterized by the dominant plant functional type. Snow melt rates were highest in locations with a late onset of melt, higher temperatures during the melt period, and higher maximum SWE before the onset of melt. After controlling for temperature, melt onset, and the maximum SWE, we found snow melt rates were highest in evergreen needleleaf forest, mixed boreal forest, and herbaceous tundra compared to deciduous needleleaf forest and deciduous shrub tundra. Tree canopy cover had little effect on snowmelt rate within each land cover type. While accounting for the influence of vegetative land cover type is necessary for predictive understanding of snowmelt rate variability across the Arctic–Boreal region. The relationships differed from observations at the ecosystem and catchment scales in other studies. Thus highlighting the importance of spatial scale in identifying snow-vegetation relationships

Vegetation type is an important predictor of the arctic summer land surface energy budget

Despite the importance of high-latitude surface energy budgets (SEBs) for land-climate interactions in the rapidly changing Arctic, uncertainties in their prediction persist. Here, we harmonize SEB observations across a network of vegetated and glaciated sites at circumpolar scale (1994-2021). Our variance-partitioning analysis identifies vegetation type as an important predictor for SEB-components during Arctic summer (June-August), compared to other SEB-drivers including climate, latitude and permafrost characteristics. Differences among vegetation types can be of similar magnitude as between vegetation and glacier surfaces and are especially high for summer sensible and latent heat fluxes. The timing of SEB-flux summer-regimes (when daily mean values exceed 0 Wm(-2)) relative to snow-free and -onset dates varies substantially depending on vegetation type, implying vegetation controls on snow-cover and SEB-flux seasonality. Our results indicate complex shifts in surface energy fluxes with land-cover transitions and a lengthening summer season, and highlight the potential for improving future Earth system models via a refined representation of Arctic vegetation types.An international team of researchers finds high potential for improving climate projections by a more comprehensive treatment of largely ignored Arctic vegetation types, underscoring the importance of Arctic energy exchange measuring stations.Peer reviewe

Global transpiration data from sap flow measurements : the SAPFLUXNET database

Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land-atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80 % of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50 % of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56 % of the datasets. Many datasets contain data for species that make up 90 % or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (https://doi.org/10.5281/zenodo.3971689; Poyatos et al., 2020a). The "sapfluxnetr" R package - designed to access, visualize, and process SAPFLUXNET data - is available from CRAN.Peer reviewe

Vegetation type is an important predictor of the arctic summer land surface energy budget

Despite the importance of high-latitude surface energy budgets (SEBs) for land-climate interactions in the rapidly changing Arctic, uncertainties in their prediction persist. Here, we harmonize SEB observations across a network of vegetated and glaciated sites at circumpolar scale (1994–2021). Our variance-partitioning analysis identifies vegetation type as an important predictor for SEB-components during Arctic summer (June-August), compared to other SEB-drivers including climate, latitude and permafrost characteristics. Differences among vegetation types can be of similar magnitude as between vegetation and glacier surfaces and are especially high for summer sensible and latent heat fluxes. The timing of SEB-flux summer-regimes (when daily mean values exceed 0 Wm−2) relative to snow-free and -onset dates varies substantially depending on vegetation type, implying vegetation controls on snow-cover and SEB-flux seasonality. Our results indicate complex shifts in surface energy fluxes with land-cover transitions and a lengthening summer season, and highlight the potential for improving future Earth system models via a refined representation of Arctic vegetation types

GEOG 338

This course will focus the ways in which terrestrial ecosystems influence exchanges of carbon dioxide, water, and energy between the land surface and the atmosphere. Through a combination of lecture, discussion, and laboratory exercises we will learn about the physical and biological processes that regulate these exchanges and how they vary in space and time. Using examples from specific ecosystems we will explore how these processes are likely to respond to climatic change

Modelling soil CO2 production and transport with dynamic source and diffusion terms: Testing the steady-state assumption using DETECT v1.0

The flux of CO2 from the soil to the atmosphere (soil respiration, Rsoil) is a major component of the global carbon cycle. Methods to measure and model Rsoil, or partition it into different components, often rely on the assumption that soil CO2 concentrations and fluxes are in steady state, implying that Rsoil is equal to the rate at which CO2 is produced by soil microbial and root respiration. Recent research, however, questions the validity of this assumption. Thus, the aim of this work was two-fold: (1) to describe a non-steady state (NSS) soil CO2 transport and production model, DETECT, and (2) to use this model to evaluate the environmental conditions under which Rsoil and CO2 production are likely in NSS. The backbone of DETECT is a non-homogeneous, partial differential equation (PDE) that describes production and transport of soil CO2, which we solve numerically at fine spatial and temporal resolution (e.g., 0.01 m increments to 1 m, every 6 hours). Production of soil CO2 is simulated for every depth and time increment as the sum of root respiration and microbial decomposition of soil organic matter, both of which can be driven by current and antecedent soil water content and temperature, which can also vary by time and depth. We also analytically solved the ordinary differential equation (ODE) corresponding to the steady-state (SS) solution to the PDE model. We applied the DETECT NSS and SS models to the 6-month growing season period representative of a native grassland in Wyoming. Simulation experiments were conducted with both model versions to evaluate factors that could affect departure from SS: (1) varying soil texture; (2) shifting the timing or frequency of precipitation; and (3) with and without the environmental antecedent drivers. For a coarse-textured soil, Rsoil from the SS model closely matched that of the NSS model. However, in a fine-textured (clay) soil, growing season Rsoil was ~3% higher under the assumption of NSS (versus SS). These differences were exaggerated in clay soil at daily time-scales whereby Rsoil under the SS assumption deviated from NSS by up to ~20% in the 10 days following a major precipitation event. Moreover, incorporation of antecedent drivers increased the magnitude of Rsoil by 15% to 37% for coarse- and fine-textured soils, respectively. However, the responses of Rsoil to the timing of precipitation and antecedent drivers did not differ between SS and NSS assumptions. In summary, the assumption of SS conditions can be violated depending on soil type and soil moisture status, as affected by precipitation inputs, and the DETECT model provides a framework for accommodating NSS conditions to better predict Rsoil and associated soil carbon cycling processes

Vegetation Indices Do Not Capture Forest Cover Variation in Upland Siberian Larch Forests

Boreal forests are changing in response to climate, with potentially important feedbacks to regional and global climate through altered carbon cycle and albedo dynamics. These feedback processes will be affected by vegetation changes, and feedback strengths will largely rely on the spatial extent and timing of vegetation change. Satellite remote sensing is widely used to monitor vegetation dynamics, and vegetation indices (VIs) are frequently used to characterize spatial and temporal trends in vegetation productivity. In this study we combine field observations of larch forest cover across a 25 km2 upland landscape in northeastern Siberia with high-resolution satellite observations to determine how the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) are related to forest cover. Across 46 forest stands ranging from 0% to 90% larch canopy cover, we find either no change, or declines in NDVI and EVI derived from PlanetScope CubeSat and Landsat data with increasing forest cover. In conjunction with field observations of NDVI, these results indicate that understory vegetation likely exerts a strong influence on vegetation indices in these ecosystems. This suggests that positive decadal trends in NDVI in Siberian larch forests may correspond primarily to increases in understory productivity, or even to declines in forest cover. Consequently, positive NDVI trends may be associated with declines in terrestrial carbon storage and increases in albedo, rather than increases in carbon storage and decreases in albedo that are commonly assumed. Moreover, it is also likely that important ecological changes such as large changes in forest density or variable forest regrowth after fire are not captured by long-term NDVI trends

Modeling soil CO2 production and transport with dynamic source and diffusion terms: testing the steady-state assumption using DETECT v1.0

The flux of CO2 from the soil to the atmosphere (soil respiration, R-soil) is a major component of the global carbon (C) cycle. Methods to measure and model R-soil, or partition it into different components, often rely on the assumption that soil CO2 concentrations and fluxes are in steady state, implying that R-soil is equal to the rate at which CO2 is produced by soil microbial and root respiration. Recent research, however, questions the validity of this assumption. Thus, the aim of this work was two-fold: (1) to describe a non-steady state (NSS) soil CO2 transport and production model, DETECT, and (2) to use this model to evaluate the environmental conditions under which R-soil and CO2 production are likely in NSS. The backbone of DETECT is a non-homogeneous, partial differential equation (PDE) that describes production and transport of soil CO2, which we solve numerically at fine spatial and temporal resolution (e.g., 0.01m increments down to 1 m, every 6 h). Production of soil CO2 is simulated for every depth and time increment as the sum of root respiration and microbial decomposition of soil organic matter. Both of these factors can be driven by current and antecedent soil water content and temperature, which can also vary by time and depth. We also analytically solved the ordinary differential equation (ODE) corresponding to the steady-state (SS) solution to the PDE model. We applied the DETECT NSS and SS models to the six-month growing season period representative of a native grassland in Wyoming. Simulation experiments were conducted with both model versions to evaluate factors that could affect departure from SS, such as (1) varying soil texture; (2) shifting the timing or frequency of precipitation; and (3) with and without the environmental antecedent drivers. For a coarse-textured soil, R-soil from the SS model closely matched that of the NSS model. However, in a fine-textured (clay) soil, growing season R-soil was similar to 3% higher under the assumption of NSS (versus SS). These differences were exaggerated in clay soil at daily time scales whereby R-soil under the SS assumption deviated from NSS by up to 35% on average in the 10 days following a major precipitation event. Incorporation of antecedent drivers increased the magnitude of R-soil by 15 to 37% for coarse-and fine-textured soils, respectively. However, the responses of R-soil to the timing of precipitation and antecedent drivers did not differ between SS and NSS assumptions. In summary, the assumption of SS conditions can be violated depending on soil type and soil moisture status, as affected by precipitation inputs. The DETECT model provides a framework for accommodating NSS conditions to better predict R soil and associated soil carbon cycling processes

Evaluating Post-Fire Vegetation Recovery in Cajander Larch Forests in Northeastern Siberia Using UAV Derived Vegetation Indices

The ability to monitor post-fire ecological responses and associated vegetation cover change is crucial to understanding how boreal forests respond to wildfire under changing climate conditions. Uncrewed aerial vehicles (UAVs) offer an affordable means of monitoring post-fire vegetation recovery for boreal ecosystems where field campaigns are spatially limited, and available satellite data are reduced by short growing seasons and frequent cloud cover. UAV data could be particularly useful across data-limited regions like the Cajander larch (Larix cajanderi Mayr.) forests of northeastern Siberia that are susceptible to amplified climate warming. Cajander larch forests require fire for regeneration but are also slow to accumulate biomass post-fire; thus, tall shrubs and other understory vegetation including grasses, mosses, and lichens dominate for several decades post-fire. Here we aim to evaluate the ability of two vegetation indices, one based on the visible spectrum (GCC; Green Chromatic Coordinate) and one using multispectral data (NDVI; Normalized Difference Vegetation Index), to predict field-based vegetation measures collected across post-fire landscapes of high-latitude Cajander larch forests. GCC and NDVI showed stronger linkages with each other at coarser spatial resolutions e.g., pixel aggregated means with 3-m, 5-m and 10-m radii compared to finer resolutions (e.g., 1-m or less). NDVI was a stronger predictor of aboveground carbon biomass and tree basal area than GCC. NDVI showed a stronger decline with increasing distance from the unburned edge into the burned forest. Our results show NDVI tended to be a stronger predictor of some field-based measures and while GCC showed similar relationships with the data, it was generally a weaker predictor of field-based measures for this region. Our findings show distinguishable edge effects and differentiation between burned and unburned forests several decades post-fire, which corresponds to the relatively slow accumulation of biomass for this ecosystem post-fire. These findings show the utility of UAV data for NDVI in this region as a tool for quantifying and monitoring the post-fire vegetation dynamics in Cajander larch forests

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

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