Description: Permafrost is an Essential Climate Variable (ECV) within the Global Climate Observing System (GCOS), which is characterized by subsurface temperatures and the depth of the seasonal thaw layer. Complementing ground-based monitoring networks, the Permafrost CCI project funded by the European Space Agency (ESA) 2018-2021 will establish Earth Observation (EO) based products for the permafrost ECV spanning the last two decades. Since ground temperature and thaw depth cannot be directly observed from space-borne sensors, we will ingest a variety of satellite and reanalysis data in a ground thermal model, which allows to quantitatively characterize the changing permafrost systems in Arctic and High-Mountain areas. As recently demonstrated for the Lena River Delta in Northern Siberia, the algorithm uses remotely sensed data sets of Land Surface Temperature (LST), Snow Water Equivalent (SWE) and landcover to drive the transient permafrost model CryoGrid 2, which yields ground temperature at various depths, in addition to thaw depth. For the circumpolar CCI product, we aim for a spatial resolution between 10 and 1km, but ensemble runs will be performed for each pixel to represent the subgrid variability of snow and land cover. The performance of the transient algorithm crucially depends on the correct representation of ground properties, in particular ice and organic contents. Therefore, the project will compile a new subsurface stratigraphy product which also holds great potential for improving Earth System Model results in permafrost environments. We report on simulation runs for various permafrost regions and characterize the accuracy and ability to reproduce trends against ground-based data. Finally, we evaluate the feasibility of future “permafrost reanalysis” products, exploiting the information content of various satellite products to deliver the best possible estimate for the permafrost thermal state over a range of spatial scales.

Description: Permafrost is an Essential Climate Variable (ECV) within the Global Climate Observing System (GCOS), which is characterized by subsurface temperatures and the depth of the seasonal thaw layer. Complementing ground-based monitoring networks, the Permafrost CCI project funded by the European Space Agency (ESA) 2018-2021 will establish Earth Observation (EO) based products for the permafrost ECV spanning the last two decades. Since ground temperature and thaw depth cannot be directly observed from space-borne sensors, we will ingest a variety of satellite and reanalysis data in a ground thermal model, which allows to quantitatively characterize the changing permafrost systems in Arctic and High-Mountain areas. As recently demonstrated for the Lena River Delta in Northern Siberia, the algorithm uses remotely sensed data sets of Land Surface Temperature (LST), Snow Water Equivalent (SWE) and landcover to drive the transient permafrost model CryoGrid 2, which yields ground temperature at various depths, in addition to thaw depth. For the circumpolar CCI product, we aim for a spatial resolution of 1km, and ensemble runs will be performed for each pixel to represent the subgrid variability of snow and land cover. The performance of the transient algorithm crucially depends on the correct representation of ground properties, in particular ice and organic contents. Therefore, the project will compile a new subsurface stratigraphy product which also holds great potential for improving Earth System Model results in permafrost environments. We present simulation runs for various permafrost regions and characterize the accuracy and ability to reproduce trends against ground-based data. Finally, we evaluate the feasibility of future “permafrost reanalysis” products, exploiting the information content of various satellite products to deliver the best possible estimate for the permafrost thermal state over a range of spatial scales.

Description: With its high albedo, low thermal conductivity and large water storing capacity, snow strongly modulates the surface energy and water balance, which makes it a critical factor in mid- to high-latitude and mountain environments. However, estimating the snow water equivalent (SWE) is challenging in remote-sensing applications already at medium spatial resolutions of 1 km. We present an ensemble-based data assimilation framework that estimates the peak subgrid SWE distribution (SSD) at the 1 km scale by assimilating fractional snow-covered area (fSCA) satellite retrievals in a simple snow model forced by downscaled reanalysis data. The basic idea is to relate the timing of the snow cover depletion (accessible from satellite products) to the peak SSD. Peak subgrid SWE is assumed to be lognormally distributed, which can be translated to a modeled time series of fSCA through the snow model. Assimilation of satellite-derived fSCA facilitates the estimation of the peak SSD, while taking into account uncertainties in both the model and the assimilated data sets. As an extension to previous studies, our method makes use of the novel (to snow data assimilation) ensemble smoother with multiple data assimilation (ES-MDA) scheme combined with analytical Gaussian anamorphosis to assimilate time series of Moderate Resolution Imaging Spectroradiometer (MODIS) and Sentinel-2 fSCA retrievals. The scheme is applied to Arctic sites near Ny-Ålesund (79° N, Svalbard, Norway) where field measurements of fSCA and SWE distributions are available. The method is able to successfully recover accurate estimates of peak SSD on most of the occasions considered. Through the ES-MDA assimilation, the root-mean-square error (RMSE) for the fSCA, peak mean SWE and peak subgrid coefficient of variation is improved by around 75, 60 and 20 %, respectively, when compared to the prior, yielding RMSEs of 0.01, 0.09 m water equivalent (w.e.) and 0.13, respectively. The ES-MDA either outperforms or at least nearly matches the performance of other ensemble-based batch smoother schemes with regards to various evaluation metrics. Given the modularity of the method, it could prove valuable for a range of satellite-era hydrometeorological reanalyses.

Description: The CryoGrid community model is a flexible toolbox for simulating the ground thermal regime and the ice-water balance for permafrost and glaciers, extending a well-established suite of permafrost models (CryoGrid 1, 2, and 3). The CryoGrid community model can accommodate a wide variety of application scenarios, which is achieved by fully modular structures through object-oriented programming. Different model components, characterized by their process representations and parameterizations, are realized as classes (i.e., objects) in CryoGrid. Standardized communication protocols between these classes ensure that they can be stacked vertically. For example, the CryoGrid community model features several classes with different complexity for the seasonal snow cover, which can be flexibly combined with a range of classes representing subsurface materials, each with their own set of process representations (e.g., soil with and without water balance, glacier ice). We present the CryoGrid architecture as well as the model physics and defining equations for the different model classes, focusing on one-dimensional model configurations which can also interact with external heat and water reservoirs. We illustrate the wide variety of simulation capabilities for a site on Svalbard, with point-scale permafrost simulations using, e.g., different soil freezing characteristics, drainage regimes, and snow representations, as well as simulations for glacier mass balance and a shallow water body. The CryoGrid community model is not intended as a static model framework but aims to provide developers with a flexible platform for efficient model development. In this study, we document both basic and advanced model functionalities to provide a baseline for the future development of novel cryosphere models.