The internal layers of ice sheets from ice-penetrating radar (IPR) investigation preserve critical information about the ice-flow field and englacial conditions. This paper presents a new detailed analysis of spatial distribution characteristics of internal layers and subglacial topography of the East Antarctic ice sheet (EAIS) from Zhongshan Station to Dome A. The radar data of 1244 km along a traverse between Zhongshan Station and Dome A of EAIS were collected during the 29th Chinese National Antarctic Research Expedition (CHINARE 29, 2012/2013). In this study, the Internal Layering Continuity Index (ILCI) and basal roughness were taken as indicators to provide an opportunity to evaluate the past internal environment and dynamics of the ice sheet. Except for the upstream of Lambert Glacier, the fold patterns of internal layers are basically similar to that of the bed topography. The relatively flat basal topography and the decrease of ILCI with increasing depth provide evidence for identifying previous rapid ice flow areas that are unavailable to satellites, especially in the upstream of Lambert Glacier. Continuous internal layers of Dome A, recording the spatial change of past ice accumulation and ice-flow history over 160 ka, almost extend to the bed, with high ILCI and high basal roughness of the corresponding bed topography. There are three kinds of basal roughness patterns along the traverse, that is, “low ξ_t low η”, “low ξ_t high η”, and “high ξ_t high η”, where ξ_t represents the amplitude of the undulations, and quantifies the vertical variation of the bedrock, and η measures the frequency variation of fluctuations and the horizontal irregularity of the profile. The characteristics of internal layers and basal topography of the traverse between Zhongshan Station and Dome A provide new information for understanding the ancient ice-flow activity and the historical evolution of EAIS.

The airborne ice-penetrating radar (IPR) is an effective method used for ice sheet exploration and is widely applied for detecting the internal structures of ice sheets and for understanding the mechanism of ice flow and the characteristics of the bottom of ice sheets. However, because of the ambient influence and the limitations of the instruments, IPR data are frequently overlaid with noise and interference, which further impedes the extraction of layer features and the interpretation of the physical characteristics of the ice sheet. In this paper, we first applied conventional filtering methods to remove the feature noise and interference in IPR data. Furthermore, machine learning methods were introduced in IPR data processing for noise removal and feature extraction. Inspired by a comparison of the filtering methods and machine learning methods, we propose a fusion method combining both filtering methods and machine-learning-based methods to optimize the feature extraction in IPR data. Field data tests indicated that, under different conditions of IPR data, the application of different methods and strategies can improve the layer feature extraction.

Understanding englacial and subglacial structures is a fundamental method of inferring ice sheets’ historical evolution and surface mass balance. The internal layer continuity index and the basal roughness are key parameters and indicators for the speculation of the relationship between the ice sheet’s internal structure or bottom and ice flow. Several methods have been proposed in the past two decades to quantitatively calculate the continuity index of ice layer geometry and the roughness of the ice–bedrock interface based on radar echo signals. These methods are mainly based on the average of the absolute value of the vertical gradient of the echo signal amplitude and the standard deviation of the horizontal fluctuation of the bedrock interface. However, these methods are limited by the amount and quality of unprocessed radar datasets and have not been widely used, which also hinders further research, such as the analysis of the englacial reflectivity, the subglacial conditions, and the history of the ice sheets. In this paper, based on geophysical processing methods for radar image denoising and deep learning for ice layer and bedrock interface extraction, we propose a new method for calculating the layer continuity index and basal roughness. Using this method, we demonstrate the ice-penetrating radar data processing and compare the imaging and calculation of the radar profiles from Dome A to Zhongshan Station, East Antarctica. We removed the noise from the processed radar data, extracted ice layer continuity features, and used other techniques to verify the calculation. The potential application of this method in the future is illustrated by several examples. We believe that this method can become an effective approach for future Antarctic geophysical and glaciological research and for obtaining more information about the history and dynamics of ice sheets from their radar-extracted internal structure.

In the Arctic and Antarctic scientific expeditions, ice-penetrating radar is an effective method for studying the bedrock under the ice sheet and ice information within the ice sheet. Because of the low conductivity of ice and the relatively uniform composition of ice sheets in the polar regions, ice-penetrating radar is able to obtain deeper and more abundant data than other geophysical methods. However, it is still necessary to suppress the noise in radar data to obtain more accurate and plentiful effective information. In this paper, the entirely non-recursive Variational Mode Decomposition (VMD) is applied to the data noise reduction of ice-penetrating radar. VMD is a decomposition method of adaptive and quasi-orthogonal signals, which decomposes airborne radar data into multiple frequency-limited quasi-orthogonal Intrinsic Mode Functions (IMFs). The IMFs containing noise are then removed according to the information distribution in the IMF’s components and the remaining IMFs are reconstructed. This paper employs this method to process the real ice-penetrating radar data, which effectively eliminates the interference noise in the data, improves the signal-to-noise ratio and obtains the clearer inner layer structure of ice. It is verified that the method can be applied to the noise reduction processing of polar ice-penetrating radar data very well, which provides a better basis for data interpretation. At last, we present fine ice structure within the ice sheet based on VMD denoised radar profile.

Basal melt in the Dome A region will influence the deep-ice-core drilling at Kunlun Station. The melting point (wet bedrock) has a higher reflectivity than the surrounding area, which can be assessed using radar echoes from the bedrock. This paper uses a linear absorption model to determine wet and dry ice–bedrock interfaces around the Kunlun drilling site. In the determination process, an artificial intelligence model was applied to extract the ice–bedrock interface for inferring the ice thickness. Additionally, the various depth-averaged attenuation rates were used to identify the maximal range of basal melting. We mapped the patterns of the wet points on the bottom of the ice sheet and the modeled basal temperature to verify the results of the wet bed conditions. According to these maps of wet bed conditions, the areas with basal melting around the drilling site primarily appear in valley walls with low basal temperatures and are linked with hydraulic potential and bedrock elevation. Identifying the ice–bedrock interface is challenging in the valley bottom area, where the melting points are less than at the valley walls. Additionally, the melting proportions of 11.8% and 3.62% were calculated from two ice-penetrating radar data in this research. The mapped melting points around the site of Kunlun ice core drilling suggest complex ice flow effects and can be used to better interpret archives of old ice for paleoclimate research.