Hyper-resolution ensemble-based snow reanalysis in mountain regions using clustering

Abstract

Spatial variability in high-relief landscapes is immense, and grid-based models cannot be run at spatial resolutions to explicitly represent important physical processes. This hampers the assessment of the current and future evolution of important issues such as water availability or mass movement hazards. Here, we present a new processing chain that couples an efficient sub-grid method with a downscaling tool and a data assimilation method with the purpose of improving numerical simulation of surface processes at multiple spatial and temporal scales in ungauged basins. The novelty of the approach is that while we add 1-2 orders of magnitude of computational cost due to ensemble simulations, we save 4-5 orders of magnitude over explicitly simulating a high-resolution grid. This approach makes data assimilation at large spatio-temporal scales feasible. In addition, this approach utilizes only freely available global datasets and is therefore able to run globally. We demonstrate marked improvements in estimating snow height and snow water equivalent at various scales using this approach that assimilates retrievals from a MODIS snow cover product. We propose that this as a suitable method for a wide variety of operational and research applications where surface models need to be run at large scales with sparse to non-existent ground observations and with the flexibility to assimilate diverse variables retrieved by Earth observation missions.