Search

Whether the presence of 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 dataset of published riverbank erosion rates. We found that erosion rates in rivers influenced by permafrost are on average six 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 W/m-1. On smaller rivers, however, hydrology rather thaw rate may be 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 km2.

Coastal cliff erosion is alongshore-variable and episodic, with retreat rates that depend upon sediment as either tools of abrasion or protective cover. However, the feedbacks between coastal cliff planform morphology, retreat rate, and sediment cover are poorly quantified. This study investigates Sargent Beach, Texas, USA, at the annual to interannual scale to explore (1) the relationship between temporal and spatial variability in cliff retreat rate, roughness, and sinuosity and (2) the response of retreat rate and roughness to changes in sand and shell hash cover of the underlying mud substrate as well as the impact of major storms using field measurements of sediment cover, erosion, and aerial images to measure shore platform morphology and retreat. A storm event in 2009 increased the planform roughness and sinuosity of the coastal cliff at Sargent Beach. Following the storm, aerial-image-derived shorelines with annual resolution show a decrease in average alongshore erosion rates from 12 to 4 m yr−1, coincident with a decrease in shoreline roughness and sinuosity (smoothing). Like the previous storm, a storm event in 2017 increased the planform roughness and sinuosity of the cliff. Over shorter timescales, monthly retreat of the sea cliff occurred only when the platform was sparsely covered with sediment cover on the shore platform, indicating that the tools and cover effects can significantly affect short-term erosion rates. The timescale to return to a smooth shoreline following a storm or roughening event, given a steady-state erosion rate, is approximately 24 years, with the long-term rate suggesting a maximum of ∼107 years until Sargent Beach breaches, compromising the Gulf Intracoastal Waterway (GIWW) under current conditions and assuming no future storms or intervention. The observed retreat rate varies, both spatially and temporally, with cliff face morphology, demonstrating the importance of multi-scale measurements and analysis for interpretation of coastal processes and patterns of cliff retreat.

Arctic river systems erode permafrost in their banks and mobilize particulate organic carbon (OC). Meandering rivers can entrain particulate OC from permafrost many meters below the depth of annual thaw, potentially enabling OC oxidation and the production of greenhouse gases. However, the amount and fate of permafrost OC that is mobilized by river erosion is uncertain. To constrain OC fluxes due to riverbank erosion and deposition, we collected riverbank and floodplain sediment samples along the Koyukuk River, which meanders through discontinuous permafrost in central Alaska. We measured sediment total OC (TOC), radiocarbon content, water content, bulk density, grain size, and floodplain stratigraphy. Radiocarbon abundance and TOC were higher in samples dominated by silt as compared to sand, which we used to map OC content onto floodplain stratigraphy and estimate carbon fluxes due to river meandering. Results showed that sediment being eroded from cutbanks and deposited as point bars had similar OC stocks (mean ± 1SD of 125.3 ± 13.1 kgOC m−2 in cutbanks versus 114.0 ± 15.7 kgOC m−2 in point bars) whether or not the banks contained permafrost. We also observed radiocarbon-depleted biospheric OC in both cutbanks and permafrost-free point bars. These results indicate that a significant fraction of aged biospheric OC that is liberated from floodplains by bank erosion is subsequently re-deposited in point bars, rather than being oxidized. The process of aging, erosion, and re-deposition of floodplain organic material may be intrinsic to river-floodplain dynamics, regardless of permafrost content.

Permafrost soils store approximately twice the amount of carbon currently present in Earth’s atmosphere and are acutely impacted by climate change due to the polar amplification of increasing global temperature. Many organic-rich permafrost sediments are located on large river floodplains, where river channel migration periodically erodes and redeposits the upper tens of meters of sediment. Channel migration exerts a first-order control on the geographic distribution of permafrost and floodplain stratigraphy and thus may affect microbial habitats. To examine how river channel migration in discontinuous permafrost environments affects microbial community composition, we used amplicon sequencing of the 16S rRNA gene on sediment samples from floodplain cores and exposed riverbanks along the Koyukuk River, a large tributary of the Yukon River in west-central Alaska. Microbial communities are sensitive to permafrost thaw: communities found in deep samples thawed by the river closely resembled near-surface active-layer communities in nonmetric multidimensional scaling analyses but did not resemble floodplain permafrost communities at the same depth. Microbial communities also displayed lower diversity and evenness in permafrost than in both the active layer and permafrost-free point bars recently deposited by river channel migration. Taxonomic assignments based on 16S and quantitative PCR for the methyl coenzyme M reductase functional gene demonstrated that methanogens and methanotrophs are abundant in older permafrost-bearing deposits but not in younger, nonpermafrost point bar deposits. The results suggested that river migration, which regulates the distribution of permafrost, also modulates the distribution of microbes potentially capable of producing and consuming methane on the Koyukuk River floodplain. IMPORTANCE Arctic lowlands contain large quantities of soil organic carbon that is currently sequestered in permafrost. With rising temperatures, permafrost thaw may allow this carbon to be consumed by microbial communities and released to the atmosphere as carbon dioxide or methane. We used gene sequencing to determine the microbial communities present in the floodplain of a river running through discontinuous permafrost. We found that the river’s lateral movement across its floodplain influences the occurrence of certain microbial communities—in particular, methane-cycling microbes were present on the older, permafrost-bearing eroding riverbank but absent on the newly deposited river bars. Riverbank sediment had microbial communities more similar to those of the floodplain active-layer samples than permafrost samples from the same depth. Therefore, spatial patterns of river migration influence the distribution of microbial taxa relevant to the warming Arctic climate.

Coastal cliff erosion is alongshore-variable and episodic, with retreat rates that depend upon sediment as either tools of abrasion or protective cover. However, the feedbacks between coastal cliff planform morphology, retreat rate, and sediment cover are poorly quantified. This study investigates Sargent Beach, Texas, USA at the annual to interannual scale to explore (1) the relationship between temporal and spatial variability in both cliff retreat rate and roughness and (2) the response of retreat rate and roughness to changes in sediment cover of the underlying mud substrate and the impact of major storms, using the low-lying mudstone cliff as a rapidly evolving model of a larger cliff system. A storm event in 2009 increased the planform roughness and sinuosity of the coastal cliff at Sargent Beach, TX. Following the storm, satellite image-derived shorelines with annual resolution show a decrease in average alongshore erosion rates from 4 to 12 m yr−1, coincident with a decrease in shoreline roughness and sinuosity (smoothing). A storm event in 2017 again increased the planform roughness and sinuosity of the cliff. The occurrence of storms and the presence of sediment to laterally erode the cliff influence the planform morphology and subsequent retreat. Over shorter timescales, monthly retreat of the sea cliff occurred only when the platform was sparsely covered with sediment cover on the wave cut platform, indicating that the tools and cover effects can significantly affect short-term erosion rates. The timescale to return to a smooth shoreline with a long-term steady-state erosion rate following a storm or roughening event is approximately five years, with the long-term rate suggesting a minimum of ~38 years until Sargent Beach breaches, compromising the Gulf Intracoastal Waterway (GIWW) under current conditions and assuming no future storms or intervention. The observed retreat rate varies, both spatially and temporally, with cliff face morphology, demonstrating the importance of multi-scale measurements and analysis for interpretation of coastal processes and patterns of cliff retreat.

17 USC 105 interim-entered record; under review. The article of record as published may be found at https://doi.org/10.3390/atmos11101090 Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent. This research was supported by the Regional and Global Model Analysis (RGMA) component of the Earth and Environmental System Modeling (EESM) program, within the U.S. Department of Energy’s Office of Science. It is a direct contribution to the HiLAT-RASM project. Funds were also provided by the closely related RUBISCO biogeochemistry validation effort, home to ILAMB and IOMB (International Land and Ocean Model Benchmarking). This research was funded by U.S. Department of Energy: OBER RGMA.

The abundant lakes dotting arctic deltas are hotspots of methane emissions and biogeochemical activity, but seasonal variability in lake extents introduces uncertainty in estimates of lacustrine carbon emissions, typically performed at annual or longer time scales. To characterize variability in lake extents, we analyzed summertime lake area loss (i.e., shrinkage) on two deltas over the past 20 years, using Landsat‐derived water masks. We find that monthly shrinkage rates have a pronounced structured variability around the channel network with the shrinkage rate systematically decreasing farther away from the channels. This pattern of shrinkage is predominantly attributed to a deeper active layer enhancing near‐surface connectivity and storage and greater vegetation density closer to the channels leading to increased evapotranspiration rates. This shrinkage signal, easily extracted from remote sensing observations, may offer the means to constrain estimates of lacustrine methane emissions and to develop process‐based estimates of depth to permafrost on arctic deltas., Key Points Summertime lake shrinkage in arctic deltas is highly structured around the channel networkHigher shrinkage rates closer to channels are attributed to enhanced near‐surface flow and storage as well as greater vegetation densityLake methane emission estimates on arctic deltas should account for spatially and seasonally variable lake extents to reduce uncertainty

Vegetation is an important component of constructional landscapes, as plants enhance deposition and provide organic sediment that can increase aggradation rates to combat land loss. We conducted two sets of laboratory experiments using alfalfa (Medicago sativa) to determine the effects of plants on channel organization and large-scale delta dynamics. In the first set, we found that rapid vegetation colonization enhanced deposition but inhibited channelization via increased form drag that reduced the shear stress available for sediment entrainment and transport. A second set of experiments used discharge fluctuations between flood and base flow (or interflood). Interfloods were critical for reworking the topset via channel incision and lateral migration to create channel relief and prevent rapid plant colonization. These low flow periods also greatly reduced the topset slope in the absence of vegetation by removing topset sediment and delivering it to the shoreline. Floods decreased relief by filling channels with sediment, resulting in periods of rapid progradation and enhanced aggradation over the topset surface, which was amplified by vegetation. The combination of discharge fluctuations and vegetation thus provided a balance of vertical aggradation and lateral progradation. We conclude that plants can inhibit channelization in depositional systems, and that discharge fluctuations encourage channel network organization to naturally balance against aggradation. Thus, variations in discharge are an important aspect of understanding the ecomorphodynamics of aggrading surfaces and modeling vegetated deltaic systems, and the combined influences of plants and discharge variations can act to balance vertical and lateral delta growth.

Lenticular sand bodies in the Upper Jurassic Norphlet Sandstone in the northeastern Gulf of Mexico basin have been interpreted as preserved large linear dunes that subsided into the underlying Louann Salt layer. Due to imaging difficulties, the thickness of the underlying salt is largely unknown. We performed physical experiments and created a mathematical model to understand dune and salt deformation and to determine the controls on the final preserved dune topography. Nonuniform dune geometry results in unequal loading by sand that drives locally varying subsidence into the salt substrate, resulting in varying degrees of dune flattening and dune relief. Our results indicate that initial salt thickness relative to dune size governs the initial rate of subsidence, thus affecting the final dune topography. However, interdune spacing and dune geometry are also important because they set the pressure gradient and determine the maximum possible subsidence. We conclude that preserved dunes with high relief indicate areas of relatively thin salt thickness and/or close dune spacing, and that dunes with low relief indicate areas of relatively thick salt thickness and/or wide dune spacing. For the case of the Norphlet Formation, dune spacing and salt thickness likely played the major roles in controlling dune subsidence.

[1] River mouths building into standing bodies of water have strikingly varied growth habits. This presents a compelling pattern formation problem that is also of great practical relevance for subsurface prediction and managing coastal wetlands. Here we present a generalized 2.5-dimensional potential vorticity (PV) theory that explains sedimentation patterns of a sediment-laden stationary jet by coupling an understanding of vorticity with suspended sediment concentration fields. We explore the physical meaning of this new sediment-PV definition, and its impact on outflow depositional patterns, by analyzing data from a shallow wall-bounded plane jet experiment and by discussing new theoretical insights. A key result is that lateral advection and diffusion of suspended sediment are directly proportional to jet vorticity, a feature that reveals the mechanistic process that forms elongated channels by focused levee deposition. The new PV theory constitutes a more generalized mathematical framework that expands the Rouse theory for the equilibrium of suspended sediment.