Arctic permafrost thaw holds the potential to drastically alter the Earth’s surface in Northern high latitudes. We utilize high-resolution Large Eddy Simulations to investigate the impact of the changing surfaces onto the neutrally stratified Atmospheric Boundary Layer (ABL). A stochastic surface model based on Gaussian Random Fields modeling typical permafrost landscapes is established in terms of two land cover classes: grass land and open water bodies, which exhibit different surface roughness length and surface sensible heat flux. A set of experiments is conducted where two parameters, the lake areal fraction and the surface correlation length, are varied to study the sensitivity of the boundary layer with respect to surface heterogeneity. Our key findings from the simulations are the following: The lake areal fraction has a substantial impact on the aggregated sensible heat flux at the blending height. The larger the lake areal fraction, the smaller the sensible heat flux. This result gives rise to a potential feedback mechanism. When the Arctic dries due to climate heating, the interaction with the ABL may accelerate permafrost thaw. Furthermore, the blending height shows significant dependency on the correlation length of the surface features. A longer surface correlation length causes an increased blending height. This finding is of relevance for land surface models concerned with Arctic permafrost as they usually do not consider a heterogeneity metric comparable to the surface correlation length.

Fire is a crucial factor in terrestrial ecosystems playing a role in disturbance for vegetation dynamics. Process-based fire models quantify fire disturbance effects in stand-alone dynamic global vegetation models (DGVMs) and their advances have incorporated both descriptions of natural processes and anthropogenic drivers. Nevertheless, these models show limited skill in modeling fire events at the global scale, due to stochastic characteristics of fire occurrence and behavior as well as the limits in empirical parameterizations in process-based models. As an alternative, machine learning has shown the capability of providing robust diagnostics of fire regimes. Here, we develop a deep-learning-based fire model (DL-fire) to estimate daily burnt area fraction at the global scale and couple it within JSBACH4, the land surface model used in the ICON ESM. The stand-alone DL-fire model forced with meteorological, terrestrial and socio-economic variables is able to simulate global total burnt area, showing 0.8 of monthly correlation (rm) with GFED4 during the evaluation period (2011-15). The performance remains similar with the hybrid modeling approach JSB4-DL-fire (rm=0.79) outperforming the currently used uncalibrated standard fire model in JSBACH4 (rm=-0.07). We further quantify the importance of each predictor by applying layer-wise relevance propagation (LRP). Overall, land properties, such as fuel amount and water content in soil layers, stand out as the major factors determining burnt fraction in DL-fire, paralleled by meteorological conditions over tropical and high latitude regions. Our study demonstrates the potential of hybrid modeling in advancing fire prediction in ESMs by integrating deep learning approaches in physics-based dynamical models.