Joel Carey Rowland

Thawing permafrost has the potential to dramatically alter the physical and ecological structure of northern landscapes. Warming of the Arctic and subsequent degradation of permafrost have created a need to assess the stability and movement of soils on hillslopes and the potential impacts on ecosystem structure. In this work, we explore the relationships among vegetation heterogeneity, soil temperature, and soil surface displacements observed from 2019 - 2022 in a watershed in the discontinuous permafrost region on the Seward Peninsula of Alaska. Vegetation heterogeneity was measured as the standard deviation (SD) of the normalized vegetation difference index (NDVI) from 3 m PlanetScope satellite imagery around each soil temperature and active layer thickness observation. Locations of observations were clustered into three soil thermal groups, warm, intermediate, and cold, based on soil temperature and active layer thickness. Average annual horizontal surface displacements were significantly lower for soils within the warm thermal group (median = 0.033 m yr-1) compared to soils within the cold thermal group (median = 0.090 m yr-1, p = < 0.001). Conversely, vegetation heterogeneity was significantly higher in the warm (median = 0.014 SD NDVI; p = 0.002) and intermediate (median = 0.015 SD NDVI; p = 0.002) groups compared with the cold thermal group (median = 0.012 SD NDVI), suggesting a warming-induced shift in vegetation community complexity. Because of the observed associations of ground surface displacement rates and vegetation heterogeneity with soil thermal state, we hypothesize that warming soil conditions induce changes in the rates and patterns of hillslope erosion due to an increase in surface movement as near-surface permafrost thaws, followed by a decrease as the permafrost table deepens and excess ice content diminishes. The transition to warm soils promotes surface ecosystem transformation, shifting the dominant vegetation at the site, given the warming climatic conditions of the region. We integrated our observations of soil temperature, vegetation heterogeneity, and soil surface displacements into a conceptual model that describes the co-evolution of hillslopes and vegetation in warming permafrost environments, which is currently unrepresented in earth system models.

Permafrost influences 25% of land in the Northern Hemisphere, where it stabilizes the ground beneath communities and infrastructure and sequesters carbon. However, the coevolution of permafrost, river dynamics, and vegetation in Arctic environments remains poorly understood. As rivers meander, they erode the floodplain at cutbanks and build new land through bar deposition, creating sequences of landforms with distinct formation ages. Here we mapped these sequences along the Koyukuk River floodplain, Alaska, analyzing permafrost occurrence, and landform and vegetation types. We used radiocarbon and optically stimulated luminescence (OSL) dating to develop a floodplain age map. Deposit ages ranged from modern to 10 ka, with more younger deposits near the modern channel. Permafrost rapidly reached 50% areal extent in all deposits older than 200 years then gradually increased up to ∼85% extent for deposits greater than 4 Kyr old. Permafrost extent correlated with increases in black spruce and wetland abundance, as well as increases in permafrost extent within wetland, and shrub and scrub vegetation classes. We developed an inverse model to constrain permafrost formation rate as a function of air temperature. Permafrost extent initially increased by ∼25% per century, in pace with vegetation succession, before decelerating to <10% per millennia as insulating overbank mud and moss slowly accumulated. Modern permafrost extent on the Koyukuk floodplain therefore reflects a dynamic balance between widespread, time‐varying permafrost formation and rapid, localized degradation due to cutbank erosion that might trigger a rapid loss of permafrost with climatic warming.

Permafrost soils are a critical component of the global carbon cycle and are locally important because they regulate the hydrologic flux from uplands to rivers. Furthermore, degradation of permafrost soils causes land surface subsidence, damaging infrastructure that is crucial for local communities. Regional and hemispherical maps of permafrost are too coarse to resolve distributions at a scale relevant to assessments of infrastructure stability or to illuminate geomorphic impacts of permafrost thaw. Here we train machine learning models to generate meter‐scale maps of near‐surface permafrost for three watersheds in the discontinuous permafrost region. The models were trained using ground truth determinations of near‐surface permafrost presence from measurements of soil temperature and electrical resistivity. We trained three classifiers: extremely randomized trees (ERTr), support vector machines (SVM), and an artificial neural network (ANN). Model uncertainty was determined using k‐fold cross validation, and the modeled extents of near‐surface permafrost were compared to the observed extents at each site. At‐a‐site near‐surface permafrost distributions predicted by the ERTr produced the highest accuracy (70%–90%). However, the transferability of the ERTr to the sites outside of the training data set was poor, with accuracies ranging from 50% to 77%. The SVM and ANN models had lower accuracies for at‐a‐site prediction (70%–83%), yet they had greater accuracy when transferred to the non‐training site (62%–78%). These models demonstrate the potential for integrating high‐resolution spatial data and machine learning models to develop maps of near‐surface permafrost extent at resolutions fine enough to assess infrastructure vulnerability and landscape morphology influenced by permafrost thaw.

Permafrost degradation is altering biogeochemical processes throughout the Arctic. Thaw‐induced changes in organic matter transformations and mineral weathering reactions are impacting fluxes of inorganic carbon (IC) and alkalinity (ALK) in Arctic rivers. However, the net impact of these changing fluxes on the concentration of carbon dioxide in the atmosphere (pCO₂) is relatively unconstrained. Resolving this uncertainty is important as thaw‐driven changes in the fluxes of IC and ALK could produce feedbacks in the global carbon cycle. Enhanced production of sulfuric acid through sulfide oxidation is particularly poorly quantified despite its potential to remove ALK from the ocean‐atmosphere system and increase pCO₂, producing a positive feedback leading to more warming and permafrost degradation. In this work, we quantified weathering in the Koyukuk River, a major tributary of the Yukon River draining discontinuous permafrost in central Alaska, based on water and sediment samples collected near the village of Huslia in summer 2018. Using measurements of major ion abundances and sulfate () sulfur (³⁴S/³²S) and oxygen (¹⁸O/¹⁶O) isotope ratios, we employed the MEANDIR inversion model to quantify the relative importance of a suite of weathering processes and their net impact on pCO₂. Calculations found that approximately 80% of in mainstem samples derived from sulfide oxidation with the remainder from evaporite dissolution. Moreover, ³⁴S/³²S ratios, ¹³C/¹²C ratios of dissolved IC, and sulfur X‐ray absorption spectra of mainstem, secondary channel, and floodplain pore fluid and sediment samples revealed modest degrees of microbial sulfate reduction within the floodplain. Weathering fluxes of ALK and IC result in lower values of pCO₂ over timescales shorter than carbonate compensation (∼10⁴ yr) and, for mainstem samples, higher values of pCO₂ over timescales longer than carbonate compensation but shorter than the residence time of marine (∼10⁷ yr). Furthermore, the absolute concentrations of and Mg²⁺ in the Koyukuk River, as well as the ratios of and Mg²⁺ to other dissolved weathering products, have increased over the past 50 years. Through analogy to similar trends in the Yukon River, we interpret these changes as reflecting enhanced sulfide oxidation due to ongoing exposure of previously frozen sediment and changes in the contributions of shallow and deep flow paths to the active channel. Overall, these findings confirm that sulfide oxidation is a substantial outcome of permafrost degradation and that the sulfur cycle responds to permafrost thaw with a timescale‐dependent feedback on warming.

Arctic coastal environments are eroding and rapidly changing. A lack of pan-Arctic observations limits our ability to understand controls on coastal erosion rates across the entire Arctic region. Here, we capitalize on an abundance of geospatial and remotely sensed data, in addition to model output, from the North Slope of Alaska to identify relationships between historical erosion rates and landscape characteristics to guide future modeling and observational efforts across the Arctic. Using existing datasets from the Alaska Beaufort Sea coast and a hierarchical clustering algorithm, we developed a set of 16 coastal typologies that captures the defining characteristics of environments susceptible to coastal erosion. Relationships between landscape characteristics and historical erosion rates show that no single variable alone is a good predictor of erosion rates. Variability in erosion rate decreases with increasing coastal elevation, but erosion rate magnitudes are highest for intermediate elevations. Areas along the Alaskan Beaufort Sea coast (ABSC) protected by barrier islands showed a three times lower erosion rate on average, suggesting that barrier islands are critical to maintaining mainland shore position. Finally, typologies with the highest erosion rates are not broadly representative of the ABSC and are generally associated with low elevation, north- to northeast-facing shorelines, a peaty pebbly silty lithology, and glaciomarine deposits with high ice content. All else being equal, warmer permafrost is also associated with higher erosion rates, suggesting that warming permafrost temperatures may contribute to higher future erosion rates on permafrost coasts. The suite of typologies can be used to guide future modeling and observational efforts by quantifying the distribution of coastlines with specific landscape characteristics and erosion rates.

Rapid Arctic warming may increase erosion and stream channel formation, which alters the flux of sediments, carbon, and nutrients in these sensitive ecosystems. Yet, understanding landscape change is hampered by a lack of predictive tools applicable to permafrost settings.

Assessing the lateral and vertical extent of permafrost is critical to understanding the fate of Arctic ecosystems under climate change. Yet, direct measurements of permafrost distribution and temperature are often limited to a small number of borehole locations. Here, we assess the use of co‐located shallow temperature and electrical resistivity tomography (ERT) measurements to estimate at high‐resolution the distribution of permafrost in three watersheds underlain by discontinuous permafrost. Synthetic modeling shows that co‐located temperature and ERT measurements allow for supervised classification schemes that provide 60% higher accuracy compared to unsupervised methods. Linking resistivity and size of the identified permafrost bodies to surface observations, we show that tall vegetation (>0.5 m) and gentle slopes (<15°) are related to warmer and smaller permafrost bodies, and a more frequent occurrence of taliks.

The Arctic is experiencing rapid climate change, and the effect on hydrologic processes and resulting geomorphic changes to hillslopes and channels is unclear because we lack quantitative models and theory for rapid changes resulting from thawing permafrost. The presence of permafrost modulates water flow and the stability of soil‐mantled slopes, implying that there should be a signature of permafrost processes, including warming‐driven disturbance, in channel network extent. To inform understanding of hillslope‐channel dynamics under changing climates, we examined soil‐mantled hillslopes within a ∼300 km² area of the Seward Peninsula, western Alaska, where discontinuous permafrost is particularly susceptible to thaw and rapid landscape change. In this study, we pair high‐resolution topographic and satellite data to multi‐annual observations of InSAR‐derived surface displacement over a 5‐year period to quantify spatial variations in topographic change across an upland landscape. We find that neither the basin slope nor the presence of knickzones controls the magnitude of recent surface displacements within the study basin, as may be expected under conceptual models of temperate hillslope evolution. Rather, the highest displacement magnitudes tended to occur at the broad hillslope‐channel transition zone. In this study area, this zone is occupied by water tracks, which are zero‐order ecogeomorphic features that concentrate surface and subsurface flow paths. Our results suggest that water tracks, which appear to occupy hillslope positions between saturation and incision thresholds, are vulnerable to warming‐induced subsidence and incision. We hypothesize that gullying within water tracks will outpace infilling by hillslope processes, resulting in the growth of the channel network under future warming.

Whether permafrost systematically alters the rate of riverbank erosion is a fundamental geomorphic question with significant importance to infrastructure, water quality, and biogeochemistry of high‐latitude watersheds. For over four decades, this question has remained unanswered due to a lack of data. Using remotely sensed imagery, we addressed this knowledge gap by quantifying riverbank erosion rates across the Arctic and subarctic. To compare these rates to non‐permafrost rivers, we assembled a global data set of published riverbank erosion rates. We found that erosion rates in rivers influenced by permafrost are on average nine times lower than non‐permafrost systems; erosion rate differences increase up to 40 times for the largest rivers. To test alternative hypotheses for the observed erosion rate difference, we examined differences in total water yield and erosional efficiency between these rivers and non‐permafrost rivers. Neither of these factors nor differences in river sediment loads provided compelling alternative explanations, leading us to conclude that permafrost limits riverbank erosion rates. This conclusion was supported by field investigations of rates and patterns of erosion along three rivers flowing through discontinuous permafrost in Alaska. Our results show that permafrost limits maximum bank erosion rates on rivers with stream powers greater than 900 Wm⁻¹. On smaller rivers, however, hydrology rather than thaw rate may be the dominant control on bank erosion. Our findings suggest that Arctic warming and hydrological changes should increase bank erosion rates on large rivers but may reduce rates on rivers with drainage areas less than a few thousand km².

Rapid climate and socioeconomic changes are transforming Arctic human–earth systems. An integral part of these systems is mobility, which encompasses the transport of humans and goods into, out of, and between Arctic regions. Impacts of climate and socioeconomic drivers on Arctic mobility are heterogenous. Methodologies are needed to quantify these impacts in measures that can be linked with broader socioeconomic systems. This article reviews existing methods and organizes them into a conceptual framework to understand trends and gaps in the literature. We found methods quantifying impacts of a range of climate drivers on most transportation modes present in the Arctic, but few methods focused on socioeconomic drivers. In addition, underrepresented were methods explicitly considering adaptive capacity of transportation systems. We provide insight into the data and relationships relevant to understanding impacts of Arctic change on transportation systems, laying a foundation for future work that investigates how these impacts fit into broader human–arth systems.