Biogeochemical and Ecological Responses to Warming Climate in High Arctic Polar Deserts
Authors: Ota Mitsuaki , Muller Amanda , Dhilon Gurbir , Siciliano Steven 
Published: March 2022

&lt;p&gt;High Arctic polar deserts cover 26% of the Arctic and are predicted to transform dramatically with rapidly rising temperatures. Previous studies found that polar deserts store larger amounts of soil organic carbon (SOC) in the permafrost than previously expected and can emit greenhouse gases (GHGs) at rate comparable to mesic Arctic ecosystems. However, the mechanism of the GHG production is not clear, which contributes to a great source of uncertainty regarding ecological feedbacks to the warming climate. Extreme climate conditions thaw the uppermost part of the permafrost, and the accumulated soil nutrients are ejected into the overlying soil layers where the subsurface nutrient patches (diapirs) form to increase carbon and nitrogen (N) contents by 7% and 20%, respectively. Previous mechanical models suggest that the ejection is facilitated by the increase in soil viscosity in the overlying soil layer. We previously found that diapirs developed about 30% of sorted circles in our study site and that the dominant vascular plant (&lt;em&gt;Salix arctica&lt;/em&gt;) increased root biomass and nitrogen uptake from diapirs. To understand a GHG-feedback to the warming climate, we collected 40 soil samples with diapirs and 40 without diapirs during July and August 2013 to investigate gross N transformation rates and GHG emissions associated with diapirs in laboratory. Our study site encompasses two Canadian High Arctic polar deserts and is located near Alexandra Fjord (78&amp;#176;51&amp;#8242;N, 75&amp;#176;54&amp;#8242;W), Ellesmere Island, Nunavut, Canada. To deal with small amounts of nitrous oxide (N&lt;sub&gt;2&lt;/sub&gt;O) emissions near or below the detection limit, we employed the hurdle models including (1) a Bernoulli component that models whether the data cross the detection limit based on covariates and (2) generalized linear model component that models the data above the detection limit. Our results showed that diapirs decreased gross N mineralization up to 48% and slowed carbon dioxide and methane emissions. Consistently, we found that diapirs contained more recalcitrant SOC using attenuated total reflectance Fourier transformed mid-infrared (ATR-FTIR) spectroscopy. ATR-FTIR also showed higher amounts of polysaccharides known to raise soil viscosity. The hurdle model approach showed that diapirs increased the estimated N&lt;sub&gt;2&lt;/sub&gt;O emissions by up to 49% under wet conditions and suggested that the increase links to the increase in the probability of N&lt;sub&gt;2&lt;/sub&gt;O emissions. On the other hand, under dry conditions, the hurdle models suggested that the increase in the estimated N&lt;sub&gt;2&lt;/sub&gt;O emissions from diapirs links to the increase in the magnitude of the N&lt;sub&gt;2&lt;/sub&gt;O emissions. The higher abundance of polysaccharides and recalcitrant SOC may indicate that biological factors are involved in forming diapirs and that diapirs supply vascular plants with nutrients as a result of a mutualistic relationship. Our study showed that diapirs altered GHG emissions and suggest that future research should include plant-microbe relationship in diapirs and other factors such as occlusion in soil aggregates for a more robust evaluation of diaper-GHG production. Furthermore, we suggest that the hurdle model may be a useful tool for evaluating N&lt;sub&gt;2&lt;/sub&gt;O emissions that are locally small but could be critical in total in the Arctic.&lt;/p&gt;