An Arctic delta reduced-complexity model and its reproduction of key geomorphological structures

Abstract

Arctic river deltas define the interface between the terrestrial Arctic and the Arctic Ocean. They are the site of sediment, nutrient, and soil organic carbon discharge to the Arctic Ocean. Arctic deltas are unique globally because they are underlain by permafrost and acted on by river and sea ice, and many are surrounded by a broad shallow ramp. Such ramps may buffer the delta from waves, but as the climate warms and permafrost thaws, the evolution of Arctic deltas will likely take a different course, with implications for both the local scale and the wider Arctic Ocean. One important way to understand and predict the evolution of Arctic deltas is through numerical models. Here we present ArcDelRCM.jl, an improved reduced-complexity model (RCM) of arctic delta evolution based on the DeltaRCM-Arctic model (Lauzon et al., 2019), which we have reconstructed in Julia language using published information. Unlike previous models, ArcDelRCM.jl is able to replicate the ramp around the delta. We have found that the delayed breakup of the so-called “bottom-fast ice” (i.e. ice that is in direct contact with the bed of the channel or the sea, also known as “bed-fast ice”) on and around the deltas is ultimately responsible for the appearance of the ramp feature in our models. However, changes made to the modelling of permafrost erosion and the protective effects of bottom-fast ice are also important contributors. Graph analyses of the delta network performed on ensemble runs show that deltas produced by ArcDelRCM.jl have more interconnected channels and contain less abandoned subnetworks. This may suggest a more even feeding of sediments to all sections of the delta shoreline, supporting ramp growth. Moreover, we showed that the morphodynamic processes during the summer months remain active enough to contribute significant sediment input to the growth and evolution of Arctic deltas and thus should not be neglected in simulations gauging the multi-year evolution of delta features. Finally, we tested a strong climate-warming scenario on the simulated deltas of ArcDelRCM.jl, with temperature, discharge, and ice conditions consistent with RCP7–8.5. We found that the ramp features degrade on the timescale of centuries and effectively disappear in under 1 millennium. Ocean processes, which are not included in these models, may further shorten the timescale. With the degradation of the ramps, any dissipative effects on wave energy they offered would also decrease. This could expose the sub-aerial parts of the deltas to increased coastal erosion, thus impacting permafrost degradation, nutrients, and carbon releases.