The mass balance of glaciers requires more detailed and continuous observations to understand their seasonal change in relation to climate. Here, we designed and installed an automated real-time monitoring platform at 4645 m a.s.l. on the Baishui River Glacier No.1 to collect continuous high-resolution observational data, and analyzed the seasonal dynamic from glacier movement and surface mass balance from glacier melting and snow accumulation. The results showed that the platform moved northeastward ~12.9 m at a rate of 0.06 ± 0.02 m d−1 between September 2021 and April 2022. The surface mass balance showed a varied temporal period. July and August were the main ablation periods, while ablation decreased and ceased in September. The glacier neither melted nor accumulated much between October and December, but began to have rapid snow accumulation in January. The glacier surface temperature varied with the air temperature and showed significant inter-seasonal differences among monsoon, post-monsoon and winter seasons. The surface mass balance also exhibited a strong response to the air temperature changes, with an average decrease of 1°C the point mass balance increased by 0.11 m w.e. from monsoon to post-monsoon and 0.22 m w.e. from post-monsoon to winter. Moreover, we found snowfall caused a decrease in the glacier surface temperature by increasing the surface albedo.

Passive seismology is becoming increasingly popular for glacier/ice-sheet structure investigations in Polar regions. Single-station passive seismic methods including P-wave receiver functions (PRFs), horizontal-to-vertical spectral ratio (HVSR) and a recently proposed autocorrelation method have been used to retrieve glacier/ice-sheet structures. Despite their successful applications, analysis regarding their detection abilities in different glaciological environments has not been reported. In this study, we compare ice thicknesses and vp/vs ratios obtained from the three methods using data collected at GAMSEIS and POLENET/ANET seismic arrays in Antarctica. Ice thickness estimates derived from the three methods are found to be consistent. Comparisons conducted under various model setups, including those involving tiled layers and sedimentary layers, show that the effectiveness of the autocorrelation method is not superior to the PRF method for retrieving ice-sheet structures. The autocorrelation method however can complement other methods as it only requires a single component seismic record.

In January 2016, static GPS measurements were carried out in a 30 × 30 km2 area centered around Kunlun station at Dome Argus (Dome A), East Antarctica, to acquire high-precision 3-D geodetic coordinates at 49 sites. By comparing the coordinates with previous GPS measurements in 2008 and 2013 at the same sites, we constructed a detailed and long-term record of the ice-surface velocity field, 2008–2016, around Dome A. During this time span, the estimated ice-surface velocity ranges from 0.8 ± 0.3 to 28.7 ± 1.6 cm a−1, with a mean of 10.4 ± 0.3 cm a−1. From 2013 to 2016, the surface elevation of most Dome A areas exhibits a rising trend, and the maximum increase of snow surface elevation is 84.8 cm. The mean snow surface elevation change rate at Dome A is estimated to be 6.6 ± 0.7 cm a−1. The difference of 1.0 cm a−1 between the snow surface change rate derived from GPS and pole-height change rate from surface mass balance is suspected to be a result of a combination of firn densification and basal melt under Dome A.

We conducted a 9-d seismic experiment in October 2015 at Laohugou Glacier No. 12. We identified microseismic signals using the short-term/long-term average trigger algorithm at four stations and classified them as long and short-duration events based on waveform, frequency, duration and magnitude characteristics. Both categories show systematical diurnal trends. The long-duration events are low-frequency tremor-like events that mainly occurred during the daytime with only several events per day. These events lasted tens of seconds to tens of minutes and are likely related to resonance of daytime meltwater. The dominant short-duration events mostly occurred during the night time with a peak occurrence frequency of ~360 h−1. Their short-duration (<0.2 s), high frequency (20–100 Hz) and dominance of Rayleigh waves are typical of events for near-surface crack opening. A strong negative correlation between the hourly event number and temperature change rate suggests that the occurrence of night-time events is controlled by the rate of night-time cooling. We estimated the near-surface tensile stress due to thermal contraction at night to be tens of kilopascals, which is enough to induce opening of surface cracks with pre-existing local stress concentrations, although we cannot exclude the effect of refreezing of meltwater produced during the day.

Using repeat GPS measurements during 2005–16, we calculated and updated two-dimensional high-resolution decadal ice surface velocity estimates along the traverse route from Zhongshan Station to and around Dome Argus, East Antarctica. Along the 71 sites of the transect, the magnitudes of ice velocity increased from near 0 in Dome Argus to 1, 10 and ~100 m a−1 at the sites DT416, DT333 and LT980, respectively. The comparison between GPS and interferometric synthetic aperture radar (InSAR) derived results agree well when the magnitude of the ice surface velocities is faster than 5 m a−1, and disagree for slower flow velocities. A scale value 1.15 and 0.12 can be applied to InSAR derived results over this region with ice surface velocity larger and <5 m a−1, respectively. We attributed the cause of the discrepancy to the insensitivity of InSAR to the magnitude of low ice surface velocities, thus confirming the importance of GPS fieldwork-based ground truth high-resolution ice velocity estimates to constrain ice-sheet dynamics.

Knowledge of the surface topography, velocity field and strain field at an ice-core site is critical to the accurate interpretation of ice-core records. At Dome Argus (Dome A), where a Chinese deep ice-core drilling project is being carried out, we have produced an accurate surface digital elevation model from GPS measurements in January 2013 at 47 sites. We identify two peaks at Dome A, with the northern peak ~7 cm higher than the southern peak. Repeat GPS measurements at 12 sites in 2008 and 2013 provide a surface velocity field around the dome. The surface velocity ranges from 3.1±2.6 to 29.4±1.2 cm a–1, with a mean of 11.1 ~2.4 cm a–1. The surface flow directions are near perpendicular to the surface elevation contours. Velocities from GPS are lower than derived from satellite radar interferometry (InSAR). From GPS velocities, the accuracy of velocity from the existing InSAR velocity field is determined, resulting in a standard deviation of 0.570 m a–1 in speed and 117.5º in direction. This result is consistent with the reported accuracy of InSAR, showing the value of in situ GPS measurements for assessing and correcting remote-sensing results. A surface strain field for the drilling site over Dome A is calculated from 24 strain triangles, showing north–south extension, east– west compression and vertical layer thinning.

This study explores an iterative method for simultaneously estimating the present-day glacial isostatic adjustment (GIA), ice mass change and elastic vertical crustal deformation of the Antarctic ice sheet (AIS) for the period October 2003–October 2009. The estimations are derived by combining mass measurements of the GRACE mission and surface height observations of the ICESat mission under the constraint of GPS vertical crustal deformation rates in the spatial domain. The influence of active subglacial lakes on GIA estimates are mitigated for the first time through additional processing of ICESat data. The inferred GIA shows that the strongest uplift is found in the Amundsen Sea Embayment (ASE) sector and subsidence mostly occurs in Adelie Terre and the East Antarctica inland. The total GIA-related mass change estimates for the entire AIS, West Antarctica Ice Sheet (WAIS), East Antarctica Ice Sheet (EAIS), and Antarctic Peninsula Ice Sheet (APIS) are 43 ± 38, 53 ± 24, −23 ± 29 and 13 ± 6 Gt a−1, respectively. The overall ice mass change of the AIS is −46 ± 43 Gt a−1 (WAIS: −104 ± 25, EAIS: 77 ± 35, APIS: −20 ± 6). The most significant ice mass loss and most significant elastic vertical crustal deformations are concentrated in the ASE and northern Antarctic Peninsula.