Volume XLII-3
Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci., XLII-3, 1657-1660, 2018
https://doi.org/10.5194/isprs-archives-XLII-3-1657-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.
Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci., XLII-3, 1657-1660, 2018
https://doi.org/10.5194/isprs-archives-XLII-3-1657-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

  30 Apr 2018

30 Apr 2018

SURFACE SNOW DENSITY OF EAST ANTARCTICA DERIVED FROM IN-SITU OBSERVATIONS

Y. Tian1,2, S. Zhang1,2, W. Du1,2, J. Chen1,2, H. Xie1,2, X. Tong1,2, and R. Li1,2 Y. Tian et al.
  • 1College of Surveying and Geo-Informatics, Tongji University, 1239 Siping Road, Shanghai, China
  • 2Center for Spatial Information Science and Sustainable Development, Tongji University, 1239 Siping Road, Shanghai, China

Keywords: Surface snow density, Matrix completion, Empirical Orthogonal Function, Kriging

Abstract. Models based on physical principles or semi-empirical parameterizations have used to compute the firn density, which is essential for the study of surface processes in the Antarctic ice sheet. However, parameterization of surface snow density is often challenged by the description of detailed local characterization. In this study we propose to generate a surface density map for East Antarctica from all the filed observations that are available. Considering that the observations are non-uniformly distributed around East Antarctica, obtained by different methods, and temporally inhomogeneous, the field observations are used to establish an initial density map with a grid size of 30 × 30 km2 in which the observations are averaged at a temporal scale of five years. We then construct an observation matrix with its columns as the map grids and rows as the temporal scale. If a site has an unknown density value for a period, we will set it to 0 in the matrix. In order to construct the main spatial and temple information of surface snow density matrix we adopt Empirical Orthogonal Function (EOF) method to decompose the observation matrix and only take first several lower-order modes, because these modes already contain most information of the observation matrix. However, there are a lot of zeros in the matrix and we solve it by using matrix completion algorithm, and then we derive the time series of surface snow density at each observation site. Finally, we can obtain the surface snow density by multiplying the modes interpolated by kriging with the corresponding amplitude of the modes. Comparative analysis have done between our surface snow density map and model results. The above details will be introduced in the paper.