Downstream Evolution of Particulate Organic Matter Composition From Permafrost Thaw Slumps

Permafrost soils, which store almost half of the global belowground organic carbon (OC), are susceptible to thaw upon climate warming. On the Peel Plateau of northwestern Canada, the number and size of retrogressive thaw slumps (RTS) has increased in recent decades due to rising temperatures and higher precipitation. These RTS features caused by the rapid thaw of ice-rich permafrost release organic matter dominantly as particulate organic carbon (POC) to the stream network. In this study, we sampled POC and streambank sediments along a fluvial transect (∼12 km) downstream from two RTS features and assessed the composition and degradation status of the mobilized permafrost OC. We found that RTS features add old, Pleistocene-aged permafrost POC to the stream system that is traceable kilometers downstream. The POC released consists mainly of recalcitrant compounds that persists within stream networks, whereas labile compounds originate from the active layer and appear to largely degrade within the scar zone of the RTS feature. Thermokarst on the Peel Plateau is likely to intensify in the future, but our data suggest that most of the permafrost OC released is not readily degradable within the stream system and thus may have little potential for atmospheric evasion. Possibilities for the recalcitrant OC to degrade over decadal to millennial time scales while being transported via larger river networks, and within the marine environment, do however, still exist. These findings add to our understanding of the vulnerable Arctic landscapes and how they may interact with the global climate.

Stable isotope (carbon, hydrogen) variation in plants and lake surface sediments from northwestern Canada

<p>Hydrogen isotope ratios of leaf waxes are used to reconstruct past hydroclimate because they are correlated to meteoric/growth water hydrogen isotopes. The interpretation of these signatures from ancient sedimentary archives relies on a thorough understanding of the drivers of modern isotope variability. Studies in the high latitudes, regions that are particularly valuable in light of their vulnerability to rapid climate change, are scarce. We studied modern vegetation (22 plants) in two areas in the Northwestern Territories (Canada): Herschel Island and Peel Plateau, to understand the stable isotope variability found in plants of Arctic regions. Bulk biomass stable carbon isotope and radiocarbon composition have been measured as well as fatty acids (wax lipids coating the plant leaves) stable carbon and hydrogen isotopes. Furthermore, lake surface sediments and river bank sediments from the Mackenzie River Delta (surrounded by the same plants) have been similarly studied. Bulk carbon isotope composition of the plants show strong difference between plant type, i.e. herbs, shrubs, lichen and moss, as shown in previous studies. Whereas the commonly used average chain length (ACL) is not useful to differentiate the plants. In term of compound-specific isotope ratios, herbs are generally <sup>2</sup>H-enriched in comparison to shrubs as shown in other regions of the world, and the C<sub>28</sub> fatty acid present the most differences amongst plant type (from ~ -207‰ for herbs to ~ – 240‰ for shrubs). No major difference between the areas is noted indicating that the ~ 250 km (Herschel Island 69.5⁰N and Peel Plateau 67.3⁰N) have no impact on the hydrogen isotope composition of the fatty acids. As such we decided to compare the plant with the lake surface sediments (from the Mackenzie Delta, located between Herschel Island and the Peel Plateau). Short-chain fatty acids, sourced from organisms growing in the lake, from isolated lakes shows <sup>2</sup>H-enriched isotopic values indicating the effect of increased evaporation in the lake during summer plant growth. Whereas long-chain fatty acids do not show any differences and are enriched compared to the shorter-chain (~ -260‰ for long-chain vs ~ -260‰ to - 280‰ for short-chain). In conclusion, differences between plant fatty acids seems to be best represented by the C<sub>28</sub> fatty acid, indicating the potential to reconstruct past vegetation and hydrological conditions in the region using lacustrine archives.</p>

Resolving thermal histories via multi-kinetic apatite fission track data: A case study from Phanerozoic strata within the Peel Plateau NWT, Canada

<p>Resolving the thermal history of sedimentary basins through geological time is essential when evaluating the maturity of source rocks within petroleum systems. Traditional methods used to estimate maximum burial temperatures in prospective sedimentary basin such as and vitrinite reflectance (%Ro) are unable to constrain the timing and duration of thermal events. In comparison, low-temperature thermochronology methods, such as apatite fission track thermochronology (AFT), can resolve detailed thermal histories within a temperature range corresponding to oil and gas generation. In the Peel Plateau of the Northwest Territories, Canada, Phanerozoic sedimentary strata exhibit oil-stained outcrops, gas seeps, and bitumen occurrences. Presently, the timing of hydrocarbon maturation events are poorly constrained, as a regional unconformity at the base of Cretaceous foreland basin strata indicates that underlying Devonian source rocks may have undergone a burial and unroofing event prior to the Cretaceous. Published organic thermal maturity values from wells within the study area range from 1.59 and 2.46 %Ro for Devonian strata and 0.54 and 1.83 %Ro within Lower Cretaceous strata. Herein, we have resolved the thermal history of the Peel Plateau through multi-kinetic AFT thermochronology. Three samples from Upper Devonian, Lower Cretaceous and Upper Cretaceous strata have pooled AFT ages of 61.0 ± 5.1 Ma, 59.5 ± 5.2 and 101.6 ± 6.7 Ma, respectively, and corresponding U-Pb ages of 497.4 ± 17.5 Ma (MSWD: 7.4), 353.5 ± 13.5 Ma (MSWD: 3.1) and 261.2 ± 8.5 Ma (MSWD: 5.9). All AFT data fail the χ<sup>2</sup> test, suggesting AFT ages do not comprise a single statistically significant population, whereas U-Pb ages reflect the pre-depositional history of the samples and are likely from various provenances. Apatite chemistry is known to control the temperature and rates at which fission tracks undergo thermal annealing. The r<sub>mro</sub> parameter uses grain specific chemistry to predict apatite’s kinetic behaviour and is used to identify kinetic populations within samples. Grain chemistry was measured via electron microprobe analysis to derive r<sub>mro</sub> values and each sample was separated into two kinetic populations that pass the χ<sup>2</sup> test: a less retentive population with ages ranging from 49.3 ± 9.3 Ma to 36.4 ± 4.7 Ma, and a more retentive population with ages ranging from 157.7 ± 19 Ma to 103.3 ± 11.8 Ma, with r<sub>mr0</sub> benchmarks ranging from 0.79 and 0.82. Thermal history models reveal Devonian strata reached maximum burial temperatures (~165°C-185°C) prior to late Paleozoic to Mesozoic unroofing, and reheated to lower temperatures (~75°C-110°C) in the Late Cretaceous to Paleogene. Both Cretaceous samples record maximum burial temperatures (75°C-95°C) also during the Late Cretaceous to Paleogene. These new data indicate that Devonian source rocks matured prior to deposition of Cretaceous strata and that subsequent burial and heating during the Cretaceous to Paleogene was limited to the low-temperature threshold of the oil window. Integrating multi-kinetic AFT data with traditional methods in petroleum geosciences can help unravel complex thermal histories of sedimentary basins. Applying these methods elsewhere can improve the characterisation of petroleum systems.</p>

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.

Characterization of mobilized sediments and organic matter in retrogressive thaw slumps on the Peel Plateau, NWT, Canada

<p>The Peel Plateau in northwestern Canada hosts some of the fastest growing “mega slumps”, retrogressive thaw slumps exceeding 2000 m<sup>2</sup> in area. The region is located at the former margin of the Laurentide ice sheet and its landscape is dominated by ice-rich hummocky moraines. Rapid permafrost thaw resulting from enhanced warming and increases in summer precipitation has been identified as a major driver of sediment mobilization in the area, with some of the largest slumps relocating up to 10<sup>6</sup> m<sup>3</sup> of previously frozen sediments into fluvial networks. The biogeochemical transformation of this thawed substrate within fluvial networks may represent a source of CO<sub>2</sub> to the atmosphere and have a large impact on downstream ecosystems, yet its fate is currently unclear. Concentrations of dissolved organic matter are lowered in slump-impacted streams, while the particle loads increase. Here, we aim to characterize the mobilized material and its sources by analyzing active layer, Holocene and Pleistocene permafrost, debris (recently thawed, still at the headwall) and slump outflow samples from four different slumps on the Peel Plateau. We use sediment properties (mineral surface area, grain size distribution), carbon isotopes (<sup>13</sup>C, <sup>14</sup>C) and molecular markers (solvent-extractable lipids, lignin phenols, cutin acids, non-extractable compound classes analyzed by pyrolysis-GCMS) in order to assess the composition and quality of the mobilized sediment and organic matter and thereby improve our understanding of their fate and downstream effects. Preliminary results show that organic matter content and radiocarbon age in debris and outflow from all four slumps are dominantly derived from Holocene and Pleistocene permafrost soils with a smaller influence of the organic-rich active layer. Degradation proxies based on extractable lipid and lignin biomarkers suggest Holocene and Pleistocene permafrost organic matter to be more matured than the fresh plant material found in the active layer, while debris and outflow samples show a mixed signal. For the non-extractable organic matter, aromatics and phenols make up the largest fraction of all samples. Lignin markers are almost exclusively found in the active layer samples, which also contain a larger proportion of polysaccharides, while N-containing compounds and alkanes make up the remaining 2-25 % with no obvious patterns. Active layer soils also have the highest median grain sizes, whereas Pleistocene permafrost soils consist of much finer mineral grains. Samples collected at the slump outflow are significantly more homogeneous (i.e., showing a narrower grain size distribution) than any of the other samples. We thus infer that both organic matter degradation and hydrodynamic sorting during transport play a role within these slump features; determining their relative magnitudes will be crucial to better assess potential feedbacks of these increasingly abundant “mega slumps” to changing climate.</p>

Downstream persistence of particulate organic carbon released from thaw slumps on the Peel Plateau, NT, Canada

<p>Underlain by ice-rich permafrost, the Peel Plateau in western Canada is highly susceptible to rapid permafrost degradation in the form of retrogressive thaw slumps and has experienced a non-linear intensification in the area, volume, and thickness of permafrost thawed since 2002. These slumps tend to occur along stream networks, which flow directly into the Peel River, through the Mackenzie Delta, and into the Beaufort Sea. Thus, lateral transport of previously sequestered organic carbon from these features has the potential to propagate far downstream. Upstream-downstream comparisons have shown that thaw slumps mobilize material to stream systems primarily in the form of particulate organic carbon (POC), increasing organic carbon yields by orders of magnitude, and switching stream networks to particle-dominated systems. Furthermore, the bulk POC released from slumps can be upwards of 10,000 <sup>14</sup>C years old, and base-extracted fluorescence measurements suggest material is more reworked since terrestrial production compared to upstream material.</p><p>To determine how far this effect propagates downstream we measured particulate and dissolved organic carbon (DOC) fluxes across stream transects extending 0.4 to 1 km downstream of thaw slumps in 1<sup>st</sup> to 2<sup>nd</sup> order streams and found no consistent decrease in TSS or POC fluxes with transit downstream.  In addition, we measured the composition (%POC, C:N, fluorescence, D<sup>14</sup>C) and flux of DOC and POC within the ~1100 km<sup>2</sup> Stony Creek watershed, examining tributary streams representing different vegetative, slump-density, and geological units in addition to the Stony Creek mainstem, to determine contributions to downstream flux. We found organic carbon fluxes were dominated by slump-mobilized POC at all points downstream of disturbance, and that these organic carbon fluxes were greater than any non-disturbed tributary stream. The <sup>14</sup>C age of POC along the Stony Creek mainstem increased by thousands of years with the introduction of slump inputs and remained similarly depleted in <sup>14</sup>C at the watershed outlet. Using historical suspended sediment, POC, and discharge data for the 75,000 km<sup>2</sup> Peel River drainage basin containing the Stony Creek watershed, we will examine whether there have been increases in instantaneous sediment and POC fluxes during the thaw season to track the trends of intensifying slump activity that have been documented on the Peel Plateau. Constraining the downstream effect of these abrupt, localized disturbances may improve detection and prediction of change that will likely cascade through the region over the coming decades.</p>