Modeling methane production and emission from thawing sub&#8208;sea permafrost on the warming Arctic Shelf
Authors: Ridolfi Emilia , Wilkenskjeld Stiig , Miesner Frederieke , Brovkin Victor , Overduin Paul , Arndt Sandra 
Published: March 2021

&lt;p&gt;The Arctic shelf hosts a large, yet poorly quantified reservoir of relic permafrost. It has been suggested that global warming, which is amplified in polar regions, will accelerate the thawing of this subsea permafrost, thus potentially unlocking large stocks of comparably reactive organic matter (OM). The microbial degradation of OM in the thawing and generally anoxic permafrost layer has the potential of producing and, ultimately, releasing important fluxes of&amp;#160; CH&lt;sub&gt;4&lt;/sub&gt; to the atmosphere. Because CH&lt;sub&gt;4&lt;/sub&gt; is a potent greenhouse gas, such a release would further intensify global warming. However, the potential role of subsea permafrost thaw on microbial CH&lt;sub&gt;4&lt;/sub&gt; production and CH&lt;sub&gt;4&lt;/sub&gt; emissions from Arctic sediments currently remains unconstrained.&lt;br&gt;Here, we use a nested model approach to address this critical knowledge gap. We developed a pseudo-three-dimensional reaction-transport model for permafrost bearing sediments on the Arctic shelf to estimate the production, consumption, and, efflux of CH&lt;sub&gt;4&lt;/sub&gt; on the Arctic shelf in response to projected subsea permafrost thaw. The model accounts for the most pertinent biogeochemical processes affecting methane and sulfur cycling in permafrost bearing marine sediments. &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160;&lt;br&gt;It is initialized based on a published submarine permafrost map (SuPerMap, [1]) and forced by a range of projected thawing rate scenarios derived from the Max Planck Institute Earth System Model (MPI-ESM) simulation results for the period 1850-2100. Critical model parameters, such as permafrost OM content and its apparent reactivity are chosen based on a comprehensive analysis of published experimental data. Here, we present the output of this environmental scenario ensemble. &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160;&lt;br&gt;Simulation results reveal that CH&lt;sub&gt;4&lt;/sub&gt; production rates are highly sensitive to changes in the apparent reactivity of permafrost OM. Although simulated CH&lt;sub&gt;4&lt;/sub&gt; production rates vary over a large range (0.001-130 PgC produced over 250 years), they generally highlight the potential for producing and, thus releasing large amounts of methane from thawing subsea permafrost on the warming Arctic Shelf.&lt;/p&gt;&lt;p&gt;[1] Overduin, P. P., Schneider von Deimling, T., Miesner, F., Grigoriev, M. N., Ruppel, C. D., Vasiliev, A., et al. (2019).&lt;br&gt;&amp;#160;Submarine permafrost map in the Arctic modeled using 1&amp;#8208;D transient heat flux (SuPerMAP). Journal of Geophysical Research: Oceans, 124, 3490&amp;#8211;3507. https://doi.org/10.1029/2018JC014675&lt;br&gt;&lt;br&gt;&lt;/p&gt;