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

• Published in: The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences
• Main Authors: Tian, Y., Zhang, S., Du, W., Chen, J., Xie, H., Tong, X., Li, R.
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2018
• Subjects: article; Verlagsveröffentlichung; Antarctic; East Antarctica; The Antarctic; Antarc*; Antarctica; Ice Sheet
• Online Access: https://doi.org/10.5194/isprs-archives-XLII-3-1657-2018; https://noa.gwlb.de/receive/cop_mods_00006505; https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00006462/isprs-archives-XLII-3-1657-2018.pdf; https://www.int-arch-photogramm-remote-sens-spatial-inf-sci.net/XLII-3/1657/2018/isprs-archives-XLII-3-1657-2018.pdf

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.