Daily Field Observations of Retrogressive Thaw Slump Dynamics in the Canadian High Arctic

With observed increases in retrogressive thaw slump (RTS) number, rates, and size in recent decades, there is a need to understand these highly dynamic landforms as they impact surrounding ecosystems and infrastructure. There is a general lack of detailed (e.g., daily) field observations of change in RTSs; we help fill this gap by monitoring three RTSs for much of the 2017 thaw period by setting up and tracking survey transects on a near daily basis. We correlated mean daily and cumulative retreat to mean daily air temperature (MDAT), total daily precipitation (TDP), and cumulative thawing degree days (TDD) using various polynomial regressions and Pearson correlation techniques. Our results show that July retreat was highly variable, and periods of increased RTS retreat did not always align with periods of increased air temperature. Also, multiple periods of increased retreat largely driven by sediment distribution in the RTS floor could occur within a single period of increased air temperature. Retreat rates decreased suddenly in early August, indicating a threshold of either air temperature, solar radiation or a combination of both must be reached for increased retreat rates. A statistically significant correlation was found between daily mean and mean cumulative retreat with MDAT (p > 0.001) and TDD (p > 0.001 and > 0.0001) but not with total daily precipitation. Correlating mean cumulative retreat and TDD using polynomial regression generated R2 values greater than 0.99 for all three sites. Both cumulative retreat and TDD account for past and current conditions, as well as lag responses, within the monitoring period. The high R2 values for the correlation of mean cumulative retreat and TDD suggest the potential for accurately modelling RTS retreat with minimal field data (air temperature and headwall position), however modelling is currently restricted to individual RTSs and only within short time scales. Monitoring RTSs on a daily scale allows RTS behaviour and trends to be identified that may be obscured at annual time scales and highlights the importance of all system inputs when considering RTS retreat dynamics.

Thaw-driven mass wasting couples slopes with downstream systems, and effects propagate through Arctic drainage networks

The intensification of thaw-driven mass wasting is transforming glacially conditioned permafrost terrain, coupling slopes with aquatic systems, and triggering a cascade of downstream effects. Within the context of recent, rapidly evolving climate controls on the geomorphology of permafrost terrain, we (A) quantify three-dimensional retrogressive thaw slump enlargement and describe the processes and thresholds coupling slopes to downstream systems, (B) investigate catchment-scale patterns of slope thermokarst impacts and the geomorphic implications, and (C) map the propagation of effects through hydrological networks draining permafrost terrain of northwestern Canada. Power-law relationships between retrogressive thaw slump area and volume (R2=0.90), as well as the thickness of permafrost thawed (R2=0.63), combined with the multi-decadal (1986–2018) increase in the areal extent of thaw slump disturbance, show a 2 order of magnitude increase in catchment-scale geomorphic activity and the coupling of slope and hydrological systems. Predominant effects are to first- and second-order streams where sediment delivery, often indicated by formation of recent debris tongue deposits, commonly exceeds the transport capacity of headwater streams by orders of magnitude, signaling centennial- to millennial-scale perturbation of downstream systems. Assessment of hydrological networks indicates that thaw-driven mass wasting directly affects over 5538 km of stream segments, 889 km of coastline, and 1379 lakes in the 994 860 km2 study area. Downstream propagation of slope thermokarst indicates a potential increase in the number of affected lakes by at least a factor of 4 (n>5692) and impacted stream length by a factor of 8 (>44 343 km), and it defines several major impact zones on lakes, deltas, and coastal areas. Prince of Wales Strait is the receiving marine environment for greatly increased sediment and geochemical fluxes from numerous slump-impacted hydrological networks draining Banks Island and Victoria Island. The Peel and Mackenzie rivers are globally significant conveyors of the slope thermokarst cascade, delivering effects to North America's largest Arctic delta and the Beaufort Sea. Climate-driven erosion of ice-rich slopes in permafrost-preserved glaciated terrain has triggered a time-transient cascade of downstream effects that signal the rejuvenation of post-glacial landscape evolution. Glacial legacy, ground-ice conditions, and continental drainage patterns dictate that terrestrial, freshwater, coastal, and marine environments of western Arctic Canada will be an interconnected hotspot of thaw-driven change through the coming millennia.

Mobilization of particulate organic matter and minerals in Zackenberg valley, Greenland

<p>Ongoing warming of the Northern high latitudes has intensified abrupt thaw processes throughout the permafrost zone. The resulting terrain disturbances are prone to release large amounts of particulate organic matter (OM) from deeper permafrost soils with thus far poorly constrained decay kinetics. Organo-mineral interactions may inhibit OM decomposition, thereby mediating the release of carbon to the atmosphere. Yet how these interactions evolve upon release and during transport along the fluvial continuum is still insufficiently understood. Here we investigate the mobilization of particulate OM from disturbed permafrost soils to the aquatic environment in the Zackenberg watershed in Northeastern Greenland. We collected soil samples in a thermo-erosion gully and a retrogressive thaw slump, as well as suspended solids and stream sediments along the glacio-nival Zackenberg River, including its tributaries, and a small headwater stream (Grænselv) affected by abrupt permafrost thaw. To evaluate the organic and mineral material transported, we compare mineral element and organic carbon (OC) concentrations, bulk carbon isotopes (<sup>13</sup>C and <sup>14</sup>C), together with source-specific molecular biomarkers (plant-wax lipids and branched glycerol dialkyl glycerol tetraethers, brGDGTs) for the suspended load with their soil and sediment counterparts.</p><p>Preliminary results show large contrasts in OC concentrations as well as Δ<sup>14</sup>C between the glacio-nival river and the headwater stream, as well as between the different thaw features. The retrogressive thaw slump mobilizes relatively OC-poor material with very low Δ<sup>14</sup>C signatures suggesting a petrogenic contribution, while soil samples from the thermo-erosion gully had higher OC concentrations and Δ<sup>14</sup>C values. For Grænselv, Δ<sup>14</sup>C values of the particulate OC were lower close to the eroding stream bank, whereas the Zackenberg main stem displayed fairly constant Δ<sup>14</sup>C values, with some of the Zackenberg tributaries delivering relatively organic-rich particles low in Δ<sup>14</sup>C.</p><p>Molecular biomarker analyses will provide additional information on specific OM sources, while X-ray Diffraction (XRD) and X-ray Fluorescence (XRF) analyses on the soils, sediments and suspended mineral load will give more detailed insights into the composition of the mineral matrices. By combining these analytical methods, we aim to improve our understanding of the interactions between minerals and OM and thereby help to constrain the fate of mobilized OM upon permafrost thaw.</p>

Detailed Characterization and Monitoring of a Retrogressive Thaw Slump from Remotely Piloted Aircraft Systems and Identifying Associated Influence on Carbon and Nitrogen Export

Ice-rich permafrost landscapes are sensitive to ongoing changes in climate. Permafrost retrogressive thaw slumps (RTSs) represent one of the more abrupt and prolonged disturbances, which occur along Arctic river and lake shorelines. These features impact local travel and infrastructure, and there are many questions regarding associated impacts on biogeochemical cycling. Predicting the duration and magnitude of impacts requires that we enhance our knowledge of RTS geomorphological drivers and rates of change. Here we demonstrate the utility of remotely piloted aircraft systems (RPAS) for documenting the volumetric change, associated drivers and potential impacts of the largest active RTS along the Old Crow River in Old Crow Flats, Yukon, Canada. RPAS surveys revealed that 29,174 m3 of sediment was exported during the initial evacuation in June 2016 and an additional 18,845 m3 continued to be exported until June 2019. More sediment export occurred during the warmer 2017 summer that experienced less cumulative rainfall than summer 2018. However, several rain events during 2017 were of higher intensity than during 2018. Overall mean soil organic carbon (SOC) and total nitrogen (TN) within sampled thaw slump sediment was 1.36% and 0.11%, respectively. A combination of multispectral, thermal and irradiance (derived from the RPAS digital surface model) data provided detailed classification of thaw slump floor terrain types including raised dry clay lobes, shaded and relatively stable, and low-lying evacuation-prone sediments. Notably, the path of evacuation-prone sediments extended to a series of ice wedges in the northern headwall, where total irradiance was highest. Using thaw slump floor mean SOC and TN values in conjunction with sediment bulk density and thaw slump fill volume, we estimated that 713 t SOC and 58 t TN were exported to the Old Crow River during the three-year study. Findings showcase the utility of high-resolution RPAS datasets for refining our knowledge of thaw slump geomorphology and associated impacts.

Thermokarst amplifies fluvial inorganic carbon cycling and export across watershed scales on the Peel Plateau, Canada

As climate warming and precipitation increase at high latitudes, permafrost terrains across the circumpolar north are poised for intensified geomorphic activity and sediment mobilization that are expected to persist for millennia. In previously glaciated permafrost terrain, ice-rich deposits are associated with large stores of reactive mineral substrate. Over geological timescales, chemical weathering moderates atmospheric CO2 levels, raising the prospect that mass wasting driven by terrain consolidation following thaw (thermokarst) may enhance weathering of permafrost sediments and thus climate feedbacks. The nature of these feedbacks depends upon the mineral composition of sediments (weathering sources) and the balance between atmospheric exchange of CO2 vs. fluvial export of carbonate alkalinity (Σ[HCO3-, CO32-]). Working in the fluvially incised, ice-rich glacial deposits of the Peel Plateau in northwestern Canada, we determine the effects of slope thermokarst in the form of retrogressive thaw slump (RTS) activity on mineral weathering sources, CO2 dynamics, and carbonate alkalinity export and how these effects integrate across watershed scales (∼ 2 to 1000 km2). We worked along three transects in nested watersheds with varying connectivity to RTS activity: a 550 m transect along a first-order thaw stream within a large RTS, a 14 km transect along a stream which directly received inputs from several RTSs, and a 70 km transect along a larger stream with headwaters that lay outside of RTS influence. In undisturbed headwaters, stream chemistry reflected CO2 from soil respiration processes and atmospheric exchange. Within the RTS, rapid sulfuric acid carbonate weathering, prompted by the exposure of sulfide- and carbonate-bearing tills, appeared to increase fluvial CO2 efflux to the atmosphere and propagate carbonate alkalinity across watershed scales. Despite covering less than 1 % of the landscape, RTS activity drove carbonate alkalinity to increase by 2 orders of magnitude along the largest transect. Amplified export of carbonate alkalinity together with isotopic signals of shifting DIC and CO2 sources along the downstream transects highlights the dynamic nature of carbon cycling that may typify glaciated permafrost watersheds subject to intensification of hillslope thermokarst. The balance between CO2 drawdown in regions where carbonic acid weathering predominates and CO2 release in regions where sulfides are more prevalent will determine the biogeochemical legacy of thermokarst and enhanced weathering in northern permafrost terrains. Effects of RTSs on carbon cycling can be expected to persist for millennia, indicating a need for their integration into predictions of weathering–carbon–climate feedbacks among thermokarst terrains.

Thermokarst amplifies fluvial inorganic carbon cycling and export across watershed scales on the Peel Plateau, Canada

The chemical weathering of minerals is a primary control on atmospheric CO2 levels and Earth's climate over geological timescales. As climate warming and precipitation intensify at high latitudes, glaciated terrains across the circumpolar north are poised for rapid geomorphic change and associated changes in mineral weathering dynamics. Here, we determine how the effects of permafrost thaw on mineral weathering sources and inorganic carbon cycling and export integrate across watershed scales (from ~ 2 to 1000 km2) in a permafrost terrain within a former glacial margin and dominated by relatively inorganic sediments (Peel Plateau, Canada). Our work was conducted along three nested transects with varying intensities of retrogressive thaw slump (RTS) thermokarst activity: a 550 m transect along a first-order thaw stream within a RTS; a 14 km transect along a stream which directly received RTS inputs; and a 70 km transect along a larger stream which received inputs from RTS-affected tributaries. In the thaw stream, rapid sulfuric acid weathering of carbonate tills appeared to amplify CO2 efflux to the atmosphere and HCO3− export downstream, where DIC and CO2 stable isotopes revealed a shift to an abiotic-inorganic driven aquatic carbon cycle. Along the intermediate transect, DIC concentrations were ten times higher in the RTS-affected reach than in the undisturbed headwaters, and decreased downstream with decreasing RTS area. Along the largest transect, HCO3− concentrations increased by two orders of magnitude in association with RTS activity, despite RTSs covering only ~ 0.5 % of the landscape. Statistical modeling of hydrochemical measurements and geospatial landscape data showed that RTSs were a primary landscape driver of HCO3− export across watershed scales. Constraining sources and rates of mineral weathering across diverse permafrost terrains will help to understand future changes in Arctic aquatic carbon cycling, as our results suggest that abiotic-inorganic processes may become prevalent.

Using ArcticDEM to identify and quantify pan-Arctic retrogressive thaw slump activity

<p>With the increased availability and coverage of high resolution satellite imagery, characterizing processes at the pan-Arctic scale is now possible. This baseline pan-Arctic product will enable us to highlight areas for future research efforts and to standardize observations that are currently locally or regionally focused. The ArcticDEM project (www.arcticdem.org) has released a large collection of 2 meter resolution Digital Elevation Models (DEMs) for all land areas above 60 °N. These DEMs are created using high resolution (~0.5 m) stereo paired satellite images (by DigitalGlobe and include Worldview- 1 (launched 2007), 2 (2009), 3 (2014) and GeoEye-1 (2008) satellites). Using repeat DEMs, we are developing algorithms for automated detection to identify and quantify land surface topographic changes from Arctic volcano eruptions and mass wasting events to create a pan-Arctic mass wasting inventory, including retrogressive thaw slumps. Currently, retreat rates reported for retrogressive thaw slumping activity differ between studies, and our dataset will enable rates to be standardized for slump activity after 2007. Furthermore, our mass wasting inventory will enable us to investigate the triggers of mass wasting events and to analyze the linkages to the contributing factors including climate, topography, and geology. We will be presenting preliminary results focusing specifically on retrogressive thaw slumps, including time series analysis for topographic change detection and using field observations for validation. We welcome collaborators who can share the field or remote sensing observations to aid in our validation efforts.</p>

Multi-methodological investigation of a retrogressive thaw slump in the Richardson Mountains, Northwest Territories, Canada

<p>The Mackenzie-Delta Region is known for strong morphological activity in context of global warming and permafrost degradation, which reveals in a large number of retrogressive thaw slumps. These are frequently found along the shorelines of inland lakes and the coast; however, this geomorphological phenomenon also occurs at inland ​​streams and creeks of the Peel Plateau and the Richardson Mountains, located in the southwest of the delta. Here several active retrogressive thaw slumps are found of which some have reached an extent of several hectares, e.g. the mega slump at the Dempster Creek.</p><p>In this study we investigated a recent retrogressive thaw slump at the edge of the Richardson Mountains close to the Dempster Highway to determine the subsurface properties using non-invasive geophysical methods. We performed three-dimensional Ground Penetrating Radar (GPR) surveys, as well as quasi-three-dimensional Electrical Resistivity Tomography (ERT) surveys in order to investigate the subsurface characteristics adjacent to the retreating headwall of the slump. These measurements provide information on the topography of the permafrost table, ice content and/or water pathways on top, within or under the permafrost layer. Additionally, we performed manual measurements of the active layer thickness for validation of the geophysical models. The approach was complemented by the analysis of high-resolution photogrammetric digital elevation models (DEM) that were generated using in situ drone acquisitions.</p><p>The measured active layer depths show a strong influence of the relief and especially of small creeks on the permafrost table topography. Likely, this influence also is the primary trigger for the initial slump activity. In addition, the ERT measurements show strong variations of the electrical resistivity values in the upper few meters, which are indicative for heterogeneities, also within the ice-rich permafrost body. Especially noticeable is a layer of low resistivity values in an area adjacent to the slump headwall. This layer is found at depths between 4m to 7m, which approximately corresponds to the base of the headwall. Here, the low resistivity values could be indicative for an unfrozen or water-rich layer below the ice-rich permafrost. Consequently, this layer may have contributed to the initial formation of the slump and is important for the spatial extension of the slump.</p><p>These results present new insights into the subsurface of an area adjacent to an active retrogressive thaw slump and may contribute to a better understanding of slump development.</p>

A landsystems approach to understanding the evolution of ice-cored topography and distribution of retrogressive thaw slumps, western Canadian Arctic

<p>The landscape of the Tuktoyaktuk Coastlands, western Canadian Arctic is dominated by glacial and geocryological processes that have modified, imprinted and sculpted the surface, depositing surficial materials upon underlying bedrock. Climate warming continues in this region at a rate that is twice the global average, and retrogressive thaw slump (RTS) activity is increasing. Recently, RTS distribution was associated with glacial limits reached by the Laurentide Ice Sheet and corresponding morainal deposits, but RTS are common in other local terrain units. In this glacial-marginal region, permafrost existed pre-glacially, and non-glacial geomorphic processes occurred throughout the Late Quaternary. Superimposed on these conditions are the effects of thermokarst during the Holocene climatic optimum, followed by a period of cooling. Collectively, these processes and associated forms and deposits have contributed variously to preservation, development, or degradation of permafrost and ground ice. The multifaceted Late Quaternary history in this region has impeded understanding of the distributions of ice-cored topography and RTS. For example, rather than glaciogenic ice, the long reigning regional model for ice-cored topography is according to post-glacial development of intrasedimental segregation-intrusion ice. Toward better understanding the evolution of the whole landscape and the distribution of climate-sensitive terrain, we use a landsystems approach as a means to understand how the ice-cored topography developed where RTS form, through analysing the cryostratigraphy. To this end, we identify 6 RTS representing a suite of ice-cored topographic settings, including: (i) preserved basal glacial ice facies within clayey diamict that has been thrusted and folded by glacial push representing morainal deposits of the Sitidgi Stade; (ii) ice contact outwash sediments associated with the Sitidgi Stade, overlying a thermo-erosional contact with underlying basal glacial icy diamict of the Toker Point Stade; (iii) deformed basal glacial ice, eroded down by meltwater-deposited outwash sands some time between the Toker Point and Sitidgi Stades (could be ca. 12.9 kyr BP); (iv) massive, undeformed segregation-intrusion basal ice, likely formed subglacially by freezing of intrasedimental water in pre-existing Pleistocene sands into the base of the glacier, overlain by glacial diamicton;  (v) deformed basal ice facies of intermediate Toker Point – Sitidgi Stades, with an upper layer that may be supra-glacial melt-out till into which segregated ice formed; and (vi) segregation ice that formed as permafrost aggraded into glaciolacustrine clays deposited in proglacial or glacially dammed basins, that was subsequently eroded down by glaciofluvial outwash (Sitidgi Stade). To summarize, the distribution of RTS reflects primarily the distribution of icy basal glacial diamict preserved in moraines, but also basal ice and icy basal diamict that are preserved beneath glaciofluvial deposits, segregation ice in glaciolacustrine deposits, and massive segregation-intrusion ice in Pleistocene sands beneath a till plain.</p>