Implementing spatially and temporally varying snow densities into the GlobSnow snow water equivalent retrieval

Snow water equivalent (SWE) is a valuable characteristic of snow cover, and it can be estimated using passive spaceborne radiometer measurements. The radiometer-based GlobSnow SWE retrieval methodology, which assimilates weather station snow depth observations into the retrieval, has improved the reliability and accuracy of SWE retrieval when compared to stand-alone radiometer SWE retrievals. To further improve the GlobSnow SWE retrieval methodology, we investigate implementing spatially and temporally varying snow densities into the retrieval procedure. Thus far, the GlobSnow SWE retrieval has used a constant snow density throughout the retrieval despite differing locations, snow depth, or time of winter. This constant snow density is a known source of inaccuracy in the retrieval. Four different versions of spatially and temporally varying snow densities are tested over a 10-year period (2000-2009). These versions use two different spatial interpolation techniques: ordinary Kriging interpolation and inverse distance weighted regression (IDWR). All versions were found to improve the SWE retrieval compared to the baseline GlobSnow v3.0 product, although differences between versions are small. Overall, the best results were obtained by implementing IDWR-interpolated densities into the algorithm, which reduced RMSE (root mean square error) and MAE (mean absolute error) by about 4 mm (8 % improvement) and 5 mm (16 % improvement) when compared to the baseline GlobSnow product, respectively. Furthermore, implementing varying snow densities into the SWE retrieval improves the magnitude and seasonal evolution of the Northern Hemisphere snow mass estimate compared to the baseline product and a product post-processed with varying snow densities.