Using data from three automatic weather stations (LGB69, Eagle and Dome A) from distinctly different climatological zones along the CHINARE (Chinese National Antarctic Research Expedition) traverse route from Zhongshan Station to Dome A, we investigated the characteristics of meteorological conditions and subsurface heat conduction. Spatial analysis indicated decreasing trends in air temperature, relative humidity and wind speed from the coastal katabatic wind zone to the inland plateau region, and air temperatures clearly showed a strong daily variability in winter, suggesting the effect from the fluctuation in the Antarctic atmospheric system. We also analyzed the optimal response time of the 1 and 3 m depth snow temperatures to the 0.1 m depth snow temperature for each site under clear/overcast and day/night situations. This showed an important enhancement to the heat transfer from shortwave radiation penetration. Using an iterative optimization method, we estimated the subsurface heat conduction variations along the transect. This was ~3–5 W m–2. Multiple maxima in daily mean subsurface fluxes were found in winter, with a typical value above 2 W m–2, while a single minimum value under –2 W m–2 was found in summer. On an annual scale, a larger mean loss of subsurface heat conduction was observed in the inland plateau compared to in the coastal katabatic area. Finally, we discussed the possible influences of turbulent and radiant transport on the vertical heat response and confirmed the wind enhancement on the growth of thermal conductivity. This preliminary study provides a brief perspective and an important reference for studying subsurface heat conduction in inland areas of Antarctica.

Stable isotope ratios (δ18O and δD) in Antarctic snow and ice are basic proxy indices of climate in ice core studies. The relation between the ratios has important indicative significance for moisture sources. In general, the fractionation characteristics of the two isotopes vary with different meteorological and topographical conditions. This paper presents the spatial and temporal distribution of meteoric water line (MWL) slopes along a traverse from the Zhongshan Station (ZSS) to Dome A in East Antarctica. It is found that the slopes decrease with the increasing distance inland from the coast and the lowest slope occurred at Dome A, where the long-range transported moisture dominates and clear sky snowing have an influence. The slopes in different layers of the snowpack showed a decreasing trend with depth and this is attributed to the fractionation during the interstitial sublimation and re-condensation processes of the water vapor. Frost flower development on the interior plateau surface can greatly alter the depth evolution of the MWL slope. The coastal snow pits also go through the post-depositional smoothing effect, but their influences are not so significant as the inland regions.

Using frequency-modulated continuous wave radar data from the 32nd Chinese Antarctic Research Expedition in 2015/16, subsurface profiles were obtained along an East Antarctic inland traverse from Zhongshan station to Dome A, and four distinct regions were selected to analyze the spatiotemporal variability in historical surface mass balance (SMB). Based on depth, density, and age data from ice cores along the traverse, the radar data were calibrated to yield average SMB data. The zone 49–195 km from the coast has the highest SMB (235 kg m−2 a−1). The 780–892 km zone was most affected by the Medieval Warm Period and the Little Ice Age, and the SMB during ad 1454–1836 (71 kg m−2 a−1) was only one-quarter of that in the 20th century. The SMB in the 1080–1157 km zone fluctuates the most, possibly due to erosion or irregular deposition of snow by katabatic winds in low SMB areas with surface elevation fluctuations. Dome A (1157–1236 km) has the lowest SMB (29 kg m−2 a−1) and did not decrease during Little Ice Age. Understanding the spatiotemporal variability of SMB in a larger space can help us understand the complex climate history of Antarctica.

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.

The stable oxygen isotope composition, major ions and isotopic compositions of strontium (Sr), neodymium (Nd) and lead (Pb) in insoluble dust from recent surface snow samples along the transect from the Zhongshan and Progress stations (located on the Amery Ice Shelf, East Antarctica) to Dome A (Summit, Antarctica) were analysed. No previous isotopic fingerprinting studies have been conducted for this transect. These data were used to document the dust provenances in Antarctica along the transect up to the highest site, Dome A, for the first time. The insoluble dust in snow samples along the coast displays an overall crust-line isotopic signature that is characterised by highly radiogenic 87Sr/86Sr values and less radiogenic 143Nd/144Nd values. These signatures are comparable with those of samples collected near the ice-free areas of the Zhongshan and Progress stations. Spatial differences are statistically significant along the transect, and the Sr, Nd and Pb isotope components in insoluble dust from two continuous snow samples at Dome A exhibit marked differences, indicating that additional dust reaches the East Antarctic Plateau. The isotopic characteristics of insoluble dust from this transect indicate that the long-distance natural dust and anthropogenic pollutants in these samples primarily originate from Australia.

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.

Stakes at 2 km intervals were installed in January 1997 and remeasured in February 1998, January 1999, January 2005 and during the 2007/08 austral summer along a 1248 km traverse route from Zhongshan station to Dome A, East Antarctica. Based on topographical parameters, meteorological features and the records of ∼650 stakes and six stake arrays, the route is divided into five zones. We find that the snow accumulation rate decreases with increasing altitude as one progresses inland, except in the zone 800–1128 km from the coast, where the average annual accumulation rate is higher than in the zone 524–800 km from the coast. The Dome A zone (1128–1248 km) has the lowest accumulation rate (35 kg m−2 a−1, 2005–08) due to having the highest elevation and being furthest from the coast. The surface mass balance in the region 202–1128 km from the coast exhibits no temporal change from 1999–2005 to 2005–08, but there is a change in the accumulation distribution. The zone from 202 to 524 km shows a decrease in surface mass balance from 84 kg m−2 a−1 in 1999–2005 to 67 kg m−2 a−1 in 2005–08, while the zone between 800 and 1128 km shows an increase from 67 kg m−2 a−1 in 1999–2005 to 75 kg m−2 a−1 in 2005–08.

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.

Knowledge of present-day ice temperature and velocity is important in order to determine how fast a glacier will respond to present and future climate change. We apply a two-dimensional higher-order thermomechanical flowband model to simulate present-day ice temperature and velocity along the main flowline of East Rongbuk Glacier, Qomolangma (Mount Everest), Himalaya. We use recent (2005–11) observational data to validate the numerical model. Modelled and observed ice surface velocities exhibit good agreement. Modelled ice temperatures agree well with observed values in two shallow boreholes that are ∼18 m deep. The model suggests that the ablation zone of East Rongbuk Glacier, km 4 to km 10 from the glacier head, is underlain by temperate ice.