Arctic-Boreal Vulnerability Experiment Research Articles

The ABoVE L-band and P-band airborne synthetic aperture radar surveys

Abstract. Permafrost-affected ecosystems of the Arctic–boreal zone in northwestern North America are undergoing profound transformation due to rapid climate change. NASA's Arctic Boreal Vulnerability Experiment (ABoVE) is investigating characteristics that make these ecosystems vulnerable or resilient to this change. ABoVE employs airborne synthetic aperture radar (SAR) as a powerful tool to characterize tundra, taiga, peatlands, and fens. Here, we present an annotated guide to the L-band and P-band airborne SAR data acquired during the 2017, 2018, 2019, and 2022 ABoVE airborne campaigns. We summarize the ∼80 SAR flight lines and how they fit into the ABoVE experimental design (Miller et al., 2023; https://doi.org/10.3334/ORNLDAAC/2150). The Supplement provides hyperlinks to extensive maps, tables, and every flight plan as well as individual flight lines. We illustrate the interdisciplinary nature of airborne SAR data with examples of preliminary results from ABoVE studies including boreal forest canopy structure from TomoSAR data over Delta Junction, AK, and the Boreal Ecosystem Research and Monitoring Sites (BERMS) area in northern Saskatchewan and active layer thickness and soil moisture data product validation. This paper is presented as a guide to enable interested readers to fully explore the ABoVE L- and P-band airborne SAR data (https://uavsar.jpl.nasa.gov/cgi-bin/data.pl).

Assessing contemporary Arctic habitat availability for a woolly mammoth proxy

Interest continues to grow in Arctic megafaunal ecological engineering, but, since the mass extinction of megafauna ~ 12–15 ka, key physiographic variables and available forage continue to change. Here we sought to assess the extent to which contemporary Arctic ecosystems are conducive to the rewilding of megaherbivores, using a woolly mammoth (M. primigenius) proxy as a model species. We first perform a literature review on woolly mammoth dietary habits. We then leverage Oak Ridge National Laboratories Distributive Active Archive Center Global Aboveground and Belowground Biomass Carbon Density Maps to generate aboveground biomass carbon density estimates in plant functional types consumed by the woolly mammoth at 300 m resolution on Alaska’s North Slope. We supplement these analyses with a NASA Arctic Boreal Vulnerability Experiment dataset to downgrade overall biomass estimates to digestible levels. We further downgrade available forage by using a conversion factor representing the relationship between total biomass and net primary productivity (NPP) for arctic vegetation types. Integrating these estimates with the forage needs of woolly mammoths, we conservatively estimate Alaska’s North Slope could support densities of 0.0–0.38 woolly mammoth km−2 (mean 0.13) across a variety of habitats. These results may inform innovative rewilding strategies.

Species Abundance Modelling of Arctic-Boreal Zone Ducks Informed by Satellite Remote Sensing

The Arctic-Boreal zone (ABZ) covers over 26 million km2 and is home to numerous duck species; however, understanding the spatiotemporal distribution of their populations across this vast landscape is challenging, in part due to extent and data scarcity. Species abundance models for ducks in the ABZ commonly use static (time invariant) habitat covariates to inform predictions, such as wetland type and extent maps. For the first time in this region, we developed species abundance models using high-resolution, time-varying wetland inundation data produced using satellite remote sensing methods. This data captured metrics of surface water extent and inundated vegetation in the Peace Athabasca Delta, Canada, which is within the NASA Arctic Boreal Vulnerability Experiment core domain. We used generalized additive mixed models to demonstrate the improved predictive value of this novel data set over time-invariant data. Our findings highlight both the potential complementarity and efficacy of dynamic wetland inundation information for improving estimation of duck abundance and distribution at high latitudes. Further, these data can be an asset to spatial targeting of biodiversity conservation efforts and developing model-based metrics of their success under rapidly changing climatic conditions.

Combining a climate-permafrost model with fine resolution remote sensor products to quantify active-layer thickness at local scales

Quantification of active-layer thickness (ALT) over seasonally frozen terrains is critical to understand the impacts of climate warming on permafrost ecosystems in cold regions. Current large-scale process-based models cannot characterize the heterogeneous response of local landscapes to homogeneous climatic forcing. Here we linked a climate-permafrost model with a machine learning solution to indirectly quantify soil conditions reflected in the edaphic factor using high resolution remote sensor products, and then effectively estimated ALT across space and time down to local scales. Our nine-year field measurements during 2014–2022 and coincident high resolution airborne hyperspectral, lidar, and spaceborne sensor products provided a unique opportunity to test the developed protocol across two permafrost experiment stations in lowland terrains of Interior Alaska. Our developed model could explain over 60% of the variance of the field measured ALT for estimating the shallowest and deepest ALT in 2015 and 2019, suggesting the potential of the designed procedure for projecting local varying terrain response to long-term climate warming scenarios. This work will enhance the National Aeronautics and Space Administration’s Arctic-Boreal Vulnerability Experiment’s mission of combining field, airborne, and spaceborne sensor products to understand the coupling of permafrost ecosystems and climate change.

A New Methodology for Estimating Surface Albedo in Heterogeneous Areas from Satellite Imagery

Precise mapping and continuous monitoring of fine-scale surface albedo are indispensable for assessing and optimizing renewable energy sources. Understanding the variations in surface reflectivity is crucial in capturing the full potential of renewable technologies, as it directly impacts the efficiency of harnessing solar energy for sustainable power generation. Satellite remote sensing stands out as the sole practical approach for estimating surface albedo at both regional and global scales. Although there are different methods to calculate albedo from satellite data, most satellite products result in low spatial resolution for large heterogeneous areas, such as urban and peri-urban areas. This paper evaluates and compares several methodologies to calculate surface albedo from Landsat 8 imagery. As a result, a new methodology for estimating surface albedo for heterogeneous areas has been proposed. The new methodology has been compared with direct and indirect albedo measurements, improving the original methodologies of Baldinelli and Silva with respect to the Arctic-Boreal Vulnerability Experiment (ABOVE) albedo images by reducing the RMSE by 85% and 52%, respectively.

Uncertain Spatial Pattern of Future Land Use and Land Cover Change and Its Impacts on Terrestrial Carbon Cycle Over the Arctic–Boreal Region of North America

AbstractLand use and land cover change (LULCC) represents a key process of human‐Earth system interaction and has profound impacts on terrestrial ecosystem carbon cycling. As a key input for ecosystem models, future gridded LULCC data is typically spatially downscaled from regional LULCC projections by integrated assessment models, such as the Global Change Analysis Model (GCAM). The uncertainty associated with the different spatial downscaling methods and its impacts on the subsequent model projections have been historically ignored and rarely examined. This study investigated this problem using two representative spatial downscaling methods and focused on their impacts on the carbon cycle over the Arctic‐Boreal Vulnerability Experiment (ABoVE) domain, where extensive LULCC is expected. Specifically, we used the Future Land Use Simulation model (FLUS) and the Demeter model to generate 0.25° gridded LULCC data (i.e., LULCCFLUS and LULCCDemeter, respectively) with the same input of regional LULCC projections from GCAM, under both the low (i.e., SSP126) and high (i.e., SSP585) greenhouse gas emission scenarios. The two sets of downscaled LULCC were used to drive the Community Land Model version 5 and prognostically simulate the terrestrial carbon cycle dynamics over the 21st century. The results suggest large spatial‐temporal differences between LULCCFLUS and LULCCDemeter, and the spatial distributions of the needleleaf evergreen boreal tree, broadleaf deciduous boreal tree, broadleaf deciduous boreal shrub, and C3 arctic grass are particularly different under both SSP126 and SSP585. Additionally, the spatiotemporal differences are larger under SSP126 than SSP585, due to more intensive LULCC under SSP126 than SSP585 from GCAM projection. The differences in LULCC further lead to large discrepancies in the spatial patterns of projected gross primary productivity, ecosystem respiration, and net ecosystem exchange, which represent more than 79% of the contributions of future LULCC in 2100. Additionally, the difference in carbon flux under SSP126 is generally larger than those under SSP585. This study highlights the importance of considering the uncertainties induced by the spatial downscaling process in future LULCC projections and carbon cycle simulations.

Simulating spatio-temporal dynamics of surface PM2.5 emitted from Alaskan wildfires

Wildfire is a major disturbance agent in Arctic boreal and tundra ecosystems that emits large quantities of atmospheric pollutants, including PM2.5. Under the substantial Arctic warming which is two to three times of global average, wildfire regimes in the high northern latitude regions are expected to intensify. This imposes a considerable threat to the health of the people residing in the Arctic regions. Alaska, as the northernmost state of the US, has a sizable rural population whose access to healthcare is greatly limited by a lack of transportation and telecommunication infrastructure and low accessibility. Unfortunately, there are only a few air quality monitoring stations across the state of Alaska, and the air quality of most remote Alaskan communities is not being systematically monitored, which hinders our understanding of the relationship between wildfire emissions and human health within these communities. Models simulating the dispersion of pollutants emitted by wildfires can be extremely valuable for providing spatially comprehensive air quality estimates in areas such as Alaska where the monitoring station network is sparse. In this study, we established a methodological framework that is based on an integration of the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model, the Wildland Fire Emissions Inventory System (WFEIS), and the Arctic-Boreal Vulnerability Experiment (ABoVE) Wildfire Date of Burning (WDoB) dataset, an Arctic-oriented fire product. Through our framework, daily gridded surface-level PM2.5 concentrations for the entire state of Alaska between 2001 and 2015 for which wildfires are responsible can be estimated. This product reveals the spatio-temporal patterns of the impacts of wildfires on the regional air quality in Alaska, which, in turn, offers a direct line of evidence indicating that wildfire is the dominant driver of PM2.5 concentrations over Alaska during the fire season. Additionally, it provides critical data inputs for research on understanding how wildfires affect human health which creates the basis for the development of effective and efficient mitigation efforts.

Burned area and carbon emissions across northwestern boreal North America from 2001–2019

Abstract. Fire is the dominant disturbance agent in Alaskan and Canadian boreal ecosystems and releases large amounts of carbon into the atmosphere. Burned area and carbon emissions have been increasing with climate change, which have the potential to alter the carbon balance and shift the region from a historic sink to a source. It is therefore critically important to track the spatiotemporal changes in burned area and fire carbon emissions over time. Here we developed a new burned-area detection algorithm between 2001–2019 across Alaska and Canada at 500 m (meters) resolution that utilizes finer-scale 30 m Landsat imagery to account for land cover unsuitable for burning. This method strictly balances omission and commission errors at 500 m to derive accurate landscape- and regional-scale burned-area estimates. Using this new burned-area product, we developed statistical models to predict burn depth and carbon combustion for the same period within the NASA Arctic–Boreal Vulnerability Experiment (ABoVE) core and extended domain. Statistical models were constrained using a database of field observations across the domain and were related to a variety of response variables including remotely sensed indicators of fire severity, fire weather indices, local climate, soils, and topographic indicators. The burn depth and aboveground combustion models performed best, with poorer performance for belowground combustion. We estimate 2.37×106 ha (2.37 Mha) burned annually between 2001–2019 over the ABoVE domain (2.87 Mha across all of Alaska and Canada), emitting 79.3 ± 27.96 Tg (±1 standard deviation) of carbon (C) per year, with a mean combustion rate of 3.13 ± 1.17 kg C m−2. Mean combustion and burn depth displayed a general gradient of higher severity in the northwestern portion of the domain to lower severity in the south and east. We also found larger-fire years and later-season burning were generally associated with greater mean combustion. Our estimates are generally consistent with previous efforts to quantify burned area, fire carbon emissions, and their drivers in regions within boreal North America; however, we generally estimate higher burned area and carbon emissions due to our use of Landsat imagery, greater availability of field observations, and improvements in modeling. The burned area and combustion datasets described here (the ABoVE Fire Emissions Database, or ABoVE-FED) can be used for local- to continental-scale applications of boreal fire science.

Patterns of regional site index across a North American boreal forest gradient

Forest structure—the height, cover, vertical complexity, and spatial patterns of trees—is a key indicator of productivity variation across forested extents. During the 2017 and 2019 growing seasons, NASA’s Arctic-Boreal Vulnerability Experiment collected full-waveform airborne LiDAR using the land, vegetation and imaging sensor, sampling boreal and tundra landscapes across a variety of ecological regions from central Canada westward through Alaska. Here, we compile and archive a geo-referenced gridded suite of these data that include vertical structure estimates and novel horizontal cover estimates of vegetation canopy cover derived from vegetation’s vertical LiDAR profile. We validate these gridded estimates with small footprint airborne LiDAR, and link >36 million of them with stand age estimates from a Landsat time-series of tree-canopy cover that we confirm with plot-level disturbance year data. We quantify the regional magnitude and variability in site index, the age-dependent rates of forest growth, across 15 boreal ecoregions in North America. With this open archive suite of forest structure data linked to stand age, we bound current forest productivity estimates across a boreal structure gradient whose response to key bioclimatic drivers may change with stand age. These results, derived from a reduction of a large archive of airborne LiDAR and a Landsat time series, quantify forest productivity bounds for input into forest and ecosystem growth models, to update forecasts of changes in North America’s boreal forests by improving the regional parametrization of forest growth rates.

Oil and natural gas wells across the NASA ABoVE domain: fugitive methane emissions and broader environmental impacts

Arctic-boreal regions are experiencing major anthropogenic disturbances in addition to intensifying natural disturbance regimes as a consequence of climate change. Oil and natural gas (OG) activities are extensive in the Arctic-boreal region of western North America, a large portion of which is underlain by permafrost. The total number and distribution of OG wells and their potential fate remain unclear. Consequently, the collective impacts of OG wells on natural and cultural resources, human health and emissions of methane (CH4), are poorly understood. Using public OG well databases, we analysed the distribution of OG wells drilled between 1984 and 2018 across the Core Domain of the NASA Arctic-Boreal Vulnerability Experiment (‘ABoVE domain’). We identified 242 007 OG wells drilled as of 2018 in the ABoVE domain, of which almost two thirds are now inactive or abandoned OG wells. We found that annual drilling has increased from 269 to 8599 OG wells from 1984 to 2014 with around 1000, 700 and 1800 OG wells drilled annually in evergreen forest, deciduous forest and herbaceous land cover types, respectively. 65 588 OG well sites were underlain by permafrost in 2012. Fugitive CH4 emissions from active and abandoned OG wells drilled in the Canadian portion of the ABoVE domain accounted for approximately 13% of the total anthropogenic CH4 emissions in Canada in 2018. Our analysis identified OG wells as an anthropogenic disturbance in the ABoVE domain with potentially non-negligible consequences to local populations, ecosystems, and the climate system.

Tipping point in North American Arctic-Boreal carbon sink persists in new generation Earth system models despite reduced uncertainty

Estimating the impacts of climate change on the global carbon cycle relies on projections from Earth system models (ESMs). While ESMs currently project large warming in the high northern latitudes, the magnitude and sign of the future carbon balance of Arctic-Boreal ecosystems are highly uncertain. The new generation of increased complexity ESMs in the Intergovernmental Panel on Climate Change Sixth Assessment Report (IPCC AR6) is intended to improve future climate projections. Here, we benchmark the Coupled Model Intercomparison Project (CMIP) 5 and 6 (8 CMIP5 members and 12 CMIP6 members) with the International Land Model Benchmarking (ILAMB) tool over the region of NASA’s Arctic-Boreal vulnerability experiment (ABoVE) in North America. We show that the projected average net biome production (NBP) in 2100 from CMIP6 is higher than that from CMIP5 in the ABoVE domain, despite the model spread being slightly narrower. Overall, CMIP6 shows better agreement with contemporary observed carbon cycle variables (photosynthesis, respiration, biomass) than CMIP5, except for soil carbon and turnover time. Although both CMIP ensemble members project the ABoVE domain will remain a carbon sink by the end of the 21st century, the sink strength in CMIP6 increases with CO2 emissions. CMIP5 and CMIP6 ensembles indicate a tipping point defined here as a negative inflection point in the NBP curve by 2050–2080 independently of the shared socioeconomic pathway (SSP) for CMIP6 or representative concentration pathway (RCP) for CMIP5. The model ensembles therefore suggest that, if the carbon sink strength keeps declining throughout the 21st century, the Arctic-Boreal ecosystems in North America may become a carbon source over the next century.

Permafrost Dynamics Observatory (PDO): 2. Joint Retrieval of Permafrost Active Layer Thickness and Soil Moisture From L‐Band InSAR and P‐Band PolSAR

AbstractSeasonal subsidence induced by ground ice melt can be measured by interferometric synthetic aperture radar (InSAR) techniques to infer active layer thickness (ALT) in permafrost regions. The magnitude of subsidence depends on both how deep the soil thawed and how much ice/water content existed in the active layer soil. To provide the later, P‐band polarimetric synthetic aperture radar (PolSAR) backscatter is used due to its sensitivity to subsurface soil moisture and freeze/thaw conditions. In this study, which is the second in a two‐part series of Permafrost Dynamics Observatory (PDO), we exploit L‐band InSAR subsidence and P‐band PolSAR backscatter in a joint retrieval scheme to simultaneously estimate ALT and soil moisture profile of permafrost active layer. Both subsidence and backscatter are explicitly characterized by physics‐based models and share a common set of soil parameters including porosity and water saturation profiles. The PDO joint retrieval has been applied to the L‐ and P‐band SAR data acquired by National Aeronautics and Space Administration/Jet Propulsion Laboratory's Uninhabited Aerial Vehicle Synthetic Aperture Radar over Alaska and western Canada during the 2017 Arctic‐Boreal Vulnerability Experiment (ABoVE) airborne campaign. This high‐resolution (30 m) regional estimates of ALT and soil moisture profile spanning over the ABoVE study domain can help link the ground‐based field surveys with satellite observations to further understand the permafrost and active layer soil process dynamics to disturbances and climate change occurring across the northern circumpolar region.

Evaluating photosynthetic activity across Arctic-Boreal land cover types using solar-induced fluorescence

Photosynthesis of terrestrial ecosystems in the Arctic-Boreal region is a critical part of the global carbon cycle. Solar-induced chlorophyll Fluorescence (SIF), a promising proxy for photosynthesis with physiological insight, has been used to track gross primary production (GPP) at regional scales. Recent studies have constructed empirical relationships between SIF and eddy covariance-derived GPP as a first step to predicting global GPP. However, high latitudes pose two specific challenges: (a) Unique plant species and land cover types in the Arctic–Boreal region are not included in the generalized SIF-GPP relationship from lower latitudes, and (b) the complex terrain and sub-pixel land cover further complicate the interpretation of the SIF-GPP relationship. In this study, we focused on the Arctic-Boreal vulnerability experiment (ABoVE) domain and evaluated the empirical relationships between SIF for high latitudes from the TROPOspheric Monitoring Instrument (TROPOMI) and a state-of-the-art machine learning GPP product (FluxCom). For the first time, we report the regression slope, linear correlation coefficient, and the goodness of the fit of SIF-GPP relationships for Arctic-Boreal land cover types with extensive spatial coverage. We found several potential issues specific to the Arctic-Boreal region that should be considered: (a) unrealistically high FluxCom GPP due to the presence of snow and water at the subpixel scale; (b) changing biomass distribution and SIF-GPP relationship along elevational gradients, and (c) limited perspective and misrepresentation of heterogeneous land cover across spatial resolutions. Taken together, our results will help improve the estimation of GPP using SIF in terrestrial biosphere models and cope with model-data uncertainties in the Arctic-Boreal region.

Tracking transient boreal wetland inundation with Sentinel-1 SAR: Peace-Athabasca Delta, Alberta and Yukon Flats, Alaska

ABSTRACT Accurate and frequent mapping of transient wetland inundation in the boreal region is critical for monitoring the ecological and societal functions of wetlands. Satellite Synthetic Aperture Radar (SAR) has long been used to map wetlands due to its sensitivity to surface inundation and ability to penetrate clouds, darkness, and certain vegetation canopies. Here, we track boreal wetland inundation by developing a two-step modified decision-tree algorithm implemented in Google Earth Engine using Sentinel-1 C-band SAR and Sentinel-2 Multispectral Instrument (MSI) time-series data as inputs. This approach incorporates temporal as well as spatial characteristics of SAR backscatter and is evaluated for the Peace-Athabasca Delta, Alberta (PAD), and Yukon Flats, Alaska (YF) from May 2017 to October 2019. Within these two boreal study areas, we map spatiotemporal patterns in wetland inundation classes of Open Water (OW), Floating Plants (FP), Emergent Plants (EP), and Flooded Vegetation (FV). Temporal variability, frequency, and maximum extents of transient wetland inundation are quantified. Retrieved inundation estimates are compared with in-situ field mapping obtained during the NASA Arctic-Boreal Vulnerability Experiment (ABoVE), and a multi-temporal Landsat-derived surface water map. Over the 2017–2019 study period, we find that fractional inundation area ranged from 18.0% to 19.0% in the PAD, and from 10.7% to 12.1% in the YF. Transient wetland inundation covered ~595 km2 of the PAD, comprising ~9.1% of its landscape, and ~102 km2 of the YF, comprising ~3.6%. The implications of these findings for wetland function monitoring, and estimating landscape-scale methane emissions are discussed, together with limitations and uncertainties of our approach. We conclude that time series of Sentinel-1 C-band SAR backscatter, screened with Sentinel-2 MSI optical imagery and validated by field measurements, offer a valuable tool for tracking transient boreal wetland inundation.

Attenuated atmospheric backscatter profiles measured by the CO2 Sounder lidar in the 2017 ASCENDS/ABoVE airborne campaign

Abstract. A series of attenuated atmospheric backscatter profiles were measured at 1572 nm by the CO2 Sounder lidar during the eight flights of the 2017 ASCENDS/ABoVE (Active Sensing of CO2 Emission over Nights, Days, and Seasons, Arctic Boreal Vulnerability Experiment) airborne campaign. In addition to measuring the column average CO2 mixing ratio from the laser signals reflected by the ground, the CO2 Sounder lidar also recorded the height-resolved attenuated atmospheric backscatter profiles beneath the aircraft. We have recently processed these vertical profiles with a 15 m vertical spacing and 1 s integration time along the flight path (∼ 200 m) for all the 2017 flights and have posted the results at NASA Distributed Active Archive Center (DAAC) for Biogeochemical Dynamics https://doi.org/10.3334/ORNLDAAC/2051 (Sun et al., 2022). This paper describes the measurement principles, the data processing technique, and the signal to noise ratios.

Using atmospheric trace gas vertical profiles to evaluate model fluxes: a case study of Arctic-CAP observations and GEOS simulations for the ABoVE domain

Abstract. Accurate estimates of carbon–climate feedbacks require an independent means for evaluating surface flux models at regional scales. The altitude-integrated enhancement (AIE) derived from the Arctic Carbon Atmospheric Profiles (Arctic-CAP) project demonstrates the utility of this bulk quantity for surface flux model evaluation. This bulk quantity leverages background mole fraction values from the middle free troposphere, is agnostic to uncertainties in boundary layer height, and can be derived from model estimates of mole fractions and vertical gradients. To demonstrate the utility of the bulk quantity, six airborne profiling surveys of atmospheric carbon dioxide (CO2), methane (CH4), and carbon monoxide (CO) throughout Alaska and northwestern Canada between April and November 2017 were completed as part of NASA's Arctic–Boreal Vulnerability Experiment (ABoVE). The Arctic-CAP sampling strategy involved acquiring vertical profiles of CO2, CH4, and CO from the surface to 5 km altitude at 25 sites around the ABoVE domain on a 4- to 6-week time interval. All Arctic-CAP measurements were compared to a global simulation using the Goddard Earth Observing System (GEOS) modeling system. Comparisons of the AIE bulk quantity from aircraft observations and GEOS simulations of atmospheric CO2, CH4, and CO highlight the fidelity of the modeled surface fluxes. The model–data comparison over the ABoVE domain reveals that while current state-of-the-art models and flux estimates are able to capture broad-scale spatial and temporal patterns in near-surface CO2 and CH4 concentrations, more work is needed to resolve fine-scale flux features that are captured in CO observations.

Mapping causal agents of disturbance in boreal and arctic ecosystems of North America using time series of Landsat data

The arctic and boreal biomes are changing as temperatures increase, including changes in the type, frequency, intensity, and seasonality of disturbances. However, our understanding of the frequency, extent, and causes of disturbance events remains incomplete. Disturbances such as fire, forest harvest, drought, wind, flooding, and insects and pathogens occur at different frequencies and severities, posing challenges to characterize and assess them under a single framework. We used the Continuous Change Detection and Classification (CCDC) algorithm on all available Landsat observations from 1984 to 2014 to detect land cover and land condition change. We mapped the following causes of disturbances annually across the study domain of NASA's Arctic Boreal Vulnerability Experiment (ABoVE): fire, logging, and pest damage. Differences between Landsat Tasseled Cap (TC) values pre- and post-disturbance were used in a random forest classifier to map causal agents. For forested ecosystems, we mapped causal agents including fire, insect, and logging. In areas that were not forest before disturbance, only the fire class was mapped. The result shows that multidimensional spectral-temporal change information is useful for mapping the causes of disturbance in arctic and boreal biomes. We employed two rounds of post-processing and used the information obtained from the comparison between the map and reference data to improve the final map. The user's and producer's accuracies of an aggregated disturbance map were 94.6% (± 2.37%) and 89.3% (± 21.78%) (95% confidence intervals in parenthesis). When evaluating the causal agents, insect damage was found the most challenging to map and validate. We estimated that 10.8% of the ABoVE core domain was disturbed between 1987 and 2012, with a margin of error of 0.5% at the 95% confidence level. Rates of disturbance due to logging remained constant over time, while fires were more episodic, and insect damage was highest between 2005 and 2010. Overall, fires affected 8.8% of the study area, while logging was 1.4% and insect damage 0.6%. Our maps indicate that pest damage became a significant issue after 2000, but it was more severe for forest ecosystems in Western Canada than in Alaska.

Sensitivity of Multifrequency Polarimetric SAR Data to Postfire Permafrost Changes and Recovery Processes in Arctic Tundra

We used full-polarimetric L-band and P-band synthetic aperture radar (SAR) data collected from the recent NASA Arctic Boreal Vulnerability Experiment (ABoVE) airborne campaign and Sentinel-1 C-band dual-polarization data to understand the sensitivity of radar backscatter intensity and phase to fire-induced changes in the surface and subsurface soil processes in Arctic tundra underlain by permafrost. The 2007 Anaktuvuk River fire on the Alaska North Slope was used as a case study. At ~10-year postfire, we observed a strong increase (>~3–4 dB) in the low-frequency radar backscatter in severely burned areas during the thaw season, in contrast to limited (< ~0.5 dB) C-band backscatter differences (VV and VH) between burned and unburned areas. However, C-band winter backscatter is generally higher (>1 dB) in burned areas than the adjacent unburned areas. Polarimetric decomposition analysis indicated a general trend toward more random surface scattering, and strong increases in double-bounce scattering and volume scattering power at both P- and L-band in the burned areas. The ice-rich yedoma region shows the largest backscatter increases in burned areas and the highest correlation with burn severity and microtopography changes. The above backscatter changes are attributed to increasing surface roughness and microtopography due to ice-wedge degradation and thermokarst development and increasing subsurface scattering due to an overall drier and deeper active layer in burned areas. Among all frequencies, P-band shows consistently larger contrast in backscatter power and phase between burned and unburned areas, which makes it potentially more useful to study fire–permafrost interactions in the Arctic over decadal time scales.

Geomorphological patterns of remotely sensed methane hot spots in the Mackenzie Delta, Canada

We studied geomorphological controls on methane (CH4) hotspots in the Mackenzie Delta region in northern Canada using airborne imaging spectroscopy collected as part of the Arctic Boreal Vulnerability Experiment. Methane emissions hotspots were retrieved at ∼25 m2 spatial resolution from a ∼10 000 km2 NASA’s Next Generation Airborne Visible/Infrared Imaging Spectrometer survey of the Mackenzie Delta acquired 31 July–3 August 2017. Separating the region into the permafrost plateau and the lowland delta, we refined the domain wide power law of CH4 enhancements detected as a function of distance to standing water in different ecoregions. We further studied the spatial decay of the distance to water relationship as a function of land cover across the Delta. We show that geomorphology exerts a strong control on the spatial patterns of emissions at regional to sub-regional scales: compared to methane hotspots detected in the upland, we find that methane hotspots detected in the lowland have a more gradual power law curve indicating a weaker spatial decay with respect to distance from water. Spatial decay of CH4 hotspots in uplands is more than 2.5 times stronger than in lowlands, which is due to differences in topography and geomorphological influence on hydrology. We demonstrate that while the observed spatial distributions of CH4 follow expected trends in lowlands and uplands, these quantitatively complement knowledge from conventional wetland and freshwater CH4 mapping and modeling.

Multi-Dimensional Remote Sensing Analysis Documents Beaver-Induced Permafrost Degradation, Seward Peninsula, Alaska

Beavers have established themselves as a key component of low arctic ecosystems over the past several decades. Beavers are widely recognized as ecosystem engineers, but their effects on permafrost-dominated landscapes in the Arctic remain unclear. In this study, we document the occurrence, reconstruct the timing, and highlight the effects of beaver activity on a small creek valley confined by ice-rich permafrost on the Seward Peninsula, Alaska using multi-dimensional remote sensing analysis of satellite (Landsat-8, Sentinel-2, Planet CubeSat, and DigitalGlobe Inc./MAXAR) and unmanned aircraft systems (UAS) imagery. Beaver activity along the study reach of Swan Lake Creek appeared between 2006 and 2011 with the construction of three dams. Between 2011 and 2017, beaver dam numbers increased, with the peak occurring in 2017 (n = 9). Between 2017 and 2019, the number of dams decreased (n = 6), while the average length of the dams increased from 20 to 33 m. Between 4 and 20 August 2019, following a nine-day period of record rainfall (&gt;125 mm), the well-established dam system failed, triggering the formation of a beaver-induced permafrost degradation feature. During the decade of beaver occupation between 2011 and 2021, the creek valley widened from 33 to 180 m (~450% increase) and the length of the stream channel network increased from ~0.6 km to more than 1.9 km (220% increase) as a result of beaver engineering and beaver-induced permafrost degradation. Comparing vegetation (NDVI) and snow (NDSI) derived indices from Sentinel-2 time-series data acquired between 2017 and 2021 for the beaver-induced permafrost degradation feature and a nearby unaffected control site, showed that peak growing season NDVI was lowered by 23% and that it extended the length of the snow-cover period by 19 days following the permafrost disturbance. Our analysis of multi-dimensional remote sensing data highlights several unique aspects of beaver engineering impacts on ice-rich permafrost landscapes. Our detailed reconstruction of the beaver-induced permafrost degradation event may also prove useful for identifying degradation of ice-rich permafrost in optical time-series datasets across regional scales. Future field- and remote sensing-based observations of this site, and others like it, will provide valuable information for the NSF-funded Arctic Beaver Observation Network (A-BON) and the third phase of the NASA Arctic-Boreal Vulnerability Experiment (ABoVE) Field Campaign.