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).

Histogram SAR tomography, or SAR histomography, is an imaging technique that utilizes single-baseline interferometric synthetic aperture radar (InSAR) histograms to infer the vertical structure of vegetation. InSAR histograms are obtained by binning SAR image pixel values, such as backscatter, according to their interferometric phase values. InSAR histograms are a proxy for vegetation structure, from which surface topography and tree height can be extracted. This paper reviews the histomography algorithm and proposes a model for evaluating the algorithm performance. New histomogram results obtained at tropical forests imaged in previous NASA/JPL UAVSAR airborne campaigns are reported to support the method.

Several studies have found that the vertical Fourier transform of lidar, interferometric Synthetic Aperture Radar (SAR), and stereo photogrammetric profiles at empirically-determined spatial frequencies enables high-performance forest aboveground biomass (AGB) estimation. Linear combinations of real and imaginary parts of Fourier transforms of Tomographic (multi-baseline) SAR (TomoSAR) profiles, from Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne data, generate similar to 20%-precision estimates of AGB in the Saskatchewan area of Canada. We found that this 20% precision can be improved to similar to 15%, a factor of 30% improvement in root mean square error (RMSE) if, in addition to using Fourier transforms of the profile itself, we use Fourier transforms of the spatial, vertical derivative of the profile. The formulation of this "derivative" algorithm is the subject of this paper.

Vegetation information is essential for analyzing aboveground biomass and understanding subsurface characteristics, such as root biomass, soil organic matter, and soil moisture conditions. In this study, we mapped boreal forest canopy height (FCH) and forest species (FS) distributions in the Delta Junction region of interior Alaska, by integrating multi-source remote sensing observations within a machine learning framework based on the extreme gradient boosting technique. Model inputs included multi-frequency (C-/L-/P-band) SAR observations from Sentinel-1, UAVSAR (Uninhabited Aerial Vehicle SAR) and AirMOSS (Airborne Microwave Observatory of Subcanopy and Subsurface), and Sentinel-2 optical reflectance data. LVIS (Land Vegetation and Ice Sensor) LiDAR measurements (RH98) and Tanana Valley State Forest timber inventory data were used as respective canopy height and species ground truth data. The combination of multi-source datasets produced the best model performance (RMSE 1.62 m for FCH, and 84.27% overall FS classification accuracy) over other models developed from single source observations. The resulting FCH and FS maps using multi-source datasets were derived at 30 m spatial resolution and showed favorable agreement with plot level field measurements from the Forest Inventory and Analysis record. The model results also captured characteristic differences in stand structure between dominant species and from post-fire vegetation succession. Our results show the potential of multi-source remote sensing observations, including low frequency microwave sensors, for monitoring boreal forest complexity and changes due to global warming.

Abstract Seasonal 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.

Abstract. Permafrost-affected ecosystems of the Arctic-boreal zone in northwestern North America are undergoing profound transformation as a result of 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, and 2019 ABoVE airborne campaigns. We provide a detailed description of the ~80 SAR flight lines and how each fits into the ABoVE experimental design. Extensive maps, tables, and hyperlinks give direct access to every flight plan as well as individual flight lines. We illustrate the powerful 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 BERMS site in northern Saskatchewan; 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 SAR data.

The GLISTIN-A instrument was flown on the NASA/AFRC C-20 (Gulfstream III) in December 2022 to observe the Mauna Loa eruption event in Hawaii, USA. As the volcano was actively erupting, several of the swaths were repeated on both the same and successive days to observe changes in the lava flow thicknesses and lava fronts. These repeated swaths provided a unique opportunity to re-address our understanding of the calibration of GLISTIN-A and its ability to precisely and accurately compute the topography of significantly sloped terrains.After calibration with localized troposphere estimates from nearby GNSS sites, adjustments to the roll from post-processed Applanix data, and applying an empirical temperature model, the average slope difference is 4.3 millidegrees and the average RMS difference is 1.7 meters between overlapping swaths. The slope and RMS differences were computed for look angles from 11 to 52 degrees.

Knowledge of the spatial and temporal distribution of active layer thickness (ALT) throughout northern Alaska would help to understand the effects of climate change in the region, as well as to quantify how much the permafrost degradation manifestly in progress there is contributing to the accumulation of greenhouse gases in the atmosphere. For this reason, we are developing extensive high-resolution maps of ALT in northern Alaska. We use machine learning along with an extensive set of spatial data layers to upscale ALT from thousands of training pixels taken from high resolution swaths of estimated ALT derived from airborne polarimetric P-band synthetic aperture radar (SAR). The resulting maps of up-scaled ALT have been compared to thousands of validation samples set aside from the PolSAR-derived swaths and to in situ ALT measurements. The maps have achieved root-mean-square errors (RMSEs) of 5–7 cm relative to validation samples, and RMSEs of approximately 10–12 cm relative to in situ ALT measurements.

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.

Organic matter (OM) content and a shallow water table are two key variables that govern the physical properties of the subsurface within the active layer of arctic soils underlain by permafrost, where the majority of biogeochemical activities take place. A detailed understanding of the soil moisture and OM profile behavior over short vertical distances through the active layer is needed to adequately model the subsurface physical processes. To observe and characterize the profiles of soil properties in the active layer, we conducted detailed soil sampling at five sites along Dalton Highway on Alaska’s North Slope. These data were used to derive a generalized logistics function to characterize the total OM and water saturation fraction behavior through the profile. Furthermore, a new pedotransfer function was developed to estimate the soil bulk density and porosity—information that is largely missing from existing soil datasets—within each layer, solely from the soil texture (organic and mineral properties). Given the currently sparse soil database of the Alaskan Arctic, these profile models can be highly beneficial for radar remote sensing models to study active layer dynamics.

Accurate vegetation information is essential for analyzing above-ground biomass and understanding subsurface characteristics, such as root biomasss, soil or ganic ma tter and soil moisture profiles. This paper investigates novel mappings of forest species and canopy height in interior Alaska. We employ Random Forests to train a regression model for canopy height mapping and a classification model for forest species mapping utilizing L-band and P-band Uninhabited Aerial Vehicle Synthetic Aperture Radar(UAVSAR). For canopy height, canopy height model (CHM) data derived from Goddard's LiDAR, Hyperspectral, and Thermal Imager (G-LiHT) are treated as ground truth. For forest species prediction, Tanana Valley State Forest (TVSF) Timber Inventory and Forest Inventory and Analysis (FIA) data are used as reference. The experimental results show the proposed method yields a root-mean-square error of 1.90 m for forest height estimation and overall accuracy of 79.54% for forest species classification. They also demonstrate the feasibility of obtaining precise vegetation information by data-driven methods, which can be further used to enhance forest radar scattering forward models.

Surface organic carbon content and soil moisture (SM) represent first‐order controls on permafrost thaw and vulnerability, yet remain challenging to map accurately. Here we explored the links between surface organic soil properties and SM dynamics in the Alaska North Slope through data analysis and process‐based modeling. Our analysis, based on in situ SM and brightness temperature data from the Soil Moisture Active Passive (SMAP) mission, indicated that the SM drydown process in Arctic tundra is closely related to surface soil organic carbon (SOC) properties. More rapid drydown was generally observed in areas with high SOC concentration (SOCC) or low bulk density. The drydown timescale derived from the SMAP polarization ratio (PR) was significantly correlated with SoilGrids surface (0–5 cm) SOCC data ( R = −0.54 ∼ −0.68, p < 0.01) at regional scale. To understand the process, we used a coupled permafrost hydrology and microwave emission model to simulate changes in the L‐band PR during the thaw season. The model accounts for the variations in organic soil hydraulic and dielectric properties with SOC content and decomposition state. Model sensitivity runs showed larger L‐band PR decreases during the early thaw season in soils with higher SOCC consistent with the above analysis, whereby highly organic soils (SOCC > 34.8%) drain water more easily with a larger amount of water discharged or lost (through evapotranspiration) relative to soils with less carbon concentration (SOCC < 17.4%). Our findings indicate that satellite L‐band observations are sensitive to tundra SM and carbon properties, and may provide critical constraints on predictions of Arctic permafrost thaw and vulnerability.

Interferometric synthetic aperture radar (InSAR) has been used to quantify a range of surface and near surface physical properties in permafrost landscapes. Most previous InSAR studies have utilized spaceborne InSAR platforms, but InSAR datasets over permafrost landscapes collected from airborne platforms have been steadily growing in recent years. Most existing algorithms dedicated toward retrieval of permafrost physical properties were originally developed for spaceborne InSAR platforms. In this study, which is the first in a two part series, we introduce a series of calibration techniques developed to apply a novel joint retrieval algorithm for permafrost active layer thickness retrieval to an airborne InSAR dataset acquired in 2017 by NASA's Uninhabited Aerial Vehicle Synthetic Aperture Radar over Alaska and Western Canada. We demonstrate how InSAR measurement uncertainties are mitigated by these calibration methods and quantify remaining measurement uncertainties with a novel method of modeling interferometric phase uncertainty using a Gaussian mixture model. Finally, we discuss the impact of native SAR resolution on InSAR measurements, the limitation of using few interferograms per retrieval, and the implications of our findings for cross-comparison of airborne and spaceborne InSAR datasets acquired over Arctic regions underlain by permafrost.

In permafrost regions, active layer thickness (ALT) observations measure the effects of climate change and predict hydrologic and elemental cycling. Often, ALT is measured through direct ground-based measurements. Recently, synthetic aperture radar (SAR) measurements from airborne platforms have emerged as a method for observing seasonal thaw subsidence, soil moisture, and ALT in permafrost regions. This study validates airborne SAR-derived ALT estimates in three regions of Alaska, USA using calibrated ground penetrating radar (GPR) geophysical data. The remotely sensed ALT estimates matched the field observations within uncertainty for 79% of locations. The average uncertainty for the GPR-derived ALT validation dataset was 0.14 m while the average uncertainty for the SAR-derived ALT in pixels coincident with GPR data was 0.19 m. In the region near Utqiaġvik, the remotely sensed ALT appeared slightly larger than field observations while in the Yukon-Kuskokwim Delta region, the remotely sensed ALT appeared slightly smaller than field observations. In the northern foothills of the Brooks Range, near Toolik Lake, there was minimal bias between the field data and remotely sensed estimates. These findings suggest that airborne SAR-derived ALT estimates compare well with in situ probing and GPR, making SAR an effective tool to monitor permafrost measurements.

Active layer thickness (ALT) is a critical metric for monitoring permafrost. How soil moisture influences ALT depends on two competing hypotheses: (a) increased soil moisture increases the latent heat of fusion for thaw, resulting in shallower active layers, and (b) increased soil moisture increases soil thermal conductivity, resulting in deeper active layers. To investigate their relative influence on thaw depth, we analyzed the Field Measurements of Soil Moisture and Active Layer Thickness (SMALT) in Alaska and Canada dataset, consisting of thousands of measurements of thaw depth and soil moisture collected at dozens of sites across Alaska and Canada as part of NASA’s Arctic Boreal Vulnerability Experiment (ABoVE). As bulk volumetric water content (VWC) integrated over the entire active layer increases, ALT decreases, supporting the latent heat hypothesis. However, as VWC in the top 12 cm of soil increases, ALT increases, supporting the thermal conductivity hypothesis. Regional temperature variations determine the baseline thaw depth while precipitation may influence the sensitivity of ALT to changes in VWC. Soil latent heat dominates over thermal conductivity in determining ALT, and the effect of bulk VWC on ALT appears consistent across sites.

Detailed information on the spatial and temporal distribution of active layer thickness (ALT) throughout the North Slope of Alaska, were it available, could offer valuable insights into the effects of climate change throughout the region and facilitate the estimation of greenhouse gas emissions resulting from permafrost degradation. We are, therefore, developing extensive high-resolution maps of ALT on the North Slope of Alaska. To do this, we use a machine learning algorithm to extrapolate ALT from high resolution strips of estimated ALT derived from airborne P-band synthetic aperture radar (SAR) acquired over two sets of flights in each of three different years. Our results indicate upscaling root-mean-square error (RMSE) of about 4 cm relative to thousands of randomly-selected SAR-derived ALT validation samples, and RMSE of approximately 10 cm relative to a small number of in-situ ALT measurements.

We used the full polarimetric L-band and P-band SAR data collected from recent NASA Arctic Boreal Vulnerability Experiment (ABoVE) airborne campaign to understand the sensitivity of longwave radar backscatter intensity and phase to the post-fire recovery process of Arctic tundra. The 2007 Anaktuvuk River fire was used as a case study. At 10-years post-fire, we observed a strong increase (>∼4 dB) in both the P- and L-band radar backscatter in the severely burned areas, in contrast to limited backscatter differences (VV, VH) between burned and unburned areas at C-band. The polarimetric target decomposition analysis indicated a general trend towards more random surface scattering, and strong increases of the double-bounce and volumetric scattering power at both P- and L-band in the burned areas. Large differences were also observed in the Pauli phase angle and the dominant-scattering-type Touzi phase angle between burned and adjacent unburned areas. The above changes are likely caused by increasing surface roughness and microtopography due to thermokarst development and ice degradation, and increasing subsurface scattering due to an overall drier and deeper active layer in the burned areas.