Spatial Variability of Glaciochemistry along a Transect from Zhongshan Station to LGB69, Antarctica

: The spatial glaciochemical variability of snow samples collected along a transect from Zhongshan Station to Lambert Glacier Basin 69 (LGB69) in Antarctica was investigated. Sea-salt ion concentrations exponentially decreased with increasing distance from the coast and/or altitude. The observed high sea-salt ion concentrations within 20.6 km of the coast may be related to preferential wet or dry deposition of sea-salt aerosols. Methanesulfonic acid (MSA), non-sea-salt sulfate (nssSO 42 − ), and calcium (Ca 2+ ) concentrations decreased along the transect. The mean MSA/nssSO 42 − value of the surface snow samples (0.34 ± 0.08) indicates that coastal sea areas are their likely source regions. The non-sea-salt Ca 2+ (nssCa 2+ )/Ca 2+ percentages of the surface snow and LGB69 snow pit samples reveal that continental dust is the primary Ca 2+ source. The δ D and δ 18 O values decreased from the coast inland. The variation of deuterium excess (d-excess) along the transect was stable and d-excess values in the two snow pit samples were low and similar, which indicates that the moisture source region between Zhongshan Station and LGB69 is a coastal sea area. These results reveal the spatial distribution patterns and sources of ions and stable isotopes, as well as factors that inﬂuence the deposition of ions and the composition of stable isotopes, which provide important insight for further studies of ice cores drilled in Antarctic coastal regions.


Introduction
The ionic and isotopic compositions of Antarctic snow can serve as important indicators of environmental and climatic information [1]. These indicators include sea-salt ions, sulfur components, and stable isotopes. Sea-salt ions (mainly Cl − , Na + , and Mg 2+ ) come from sea ice and/or sea water [2,3], which are strongly influenced by sea-ice extent and meteorological conditions [4,5]. Marine sulfur components (non-sea-salt sulfate (nssSO 4 2− ) and methanesulfonic acid (MSA)) can be attributed to marine biogenic activity via dimethyl sulfide (DMS) oxidation [6]. DMS is the only source of MSA, and thus MSA is the best alternative indicator for DMS [7][8][9]. In contrast, nssSO 4 2− has various sources, such as volcanic eruptions, anthropogenic emissions, and crustal erosion [10]. The nssSO 4 2− from volcanic eruptions can be used to date ice cores [10,11] based on the eruption history [12][13][14]. Stable isotopes (δD and δ 18 O) reflect climate variations over time and space and are commonly used to quantify changes of air temperature in past environments. Deuterium excess (d-excess) is a second-order parameter used to identify moisture source regions [15][16][17].
High-resolution climate records are often obtained from ice cores drilled at sites with high annual accumulation rates. These data are used to interpret temporal variations Previous investigations of stable isotopes from surface snow samples along the transect from Zhongshan Station to Dome A showed that δD and δ 18 O values negatively correlated with distance from the coast and altitude and that the δD/δ 18 O ratio was lower than that of the Global Meteoric Water Line (GMWL) and the mean value for Antarctica (7.75) [24,26]. This implies that more intense kinetic fractionation processes occur during snow formation in the interior region of the transect [26,27]. Previous studies indicated a significant increase of Antarctic d-excess values at altitudes near or above 2000 m owing to different moisture origins [15,32]. The d-excess values along the transect from Zhongshan Station to Dome A exhibited a positive trend in the inland direction, with a step change as a function of either distance to the coast or latitude characterized by higher values for elevations abovẽ 2750 m. It has been speculated that a specific transport path brings moisture toward Dome A either from warm, distant moisture sources, or from sea-ice edge areas [17,26,27].
Ding et al. [33] showed that the coastal region (68-202 km) of the transect from Zhongshan Station to Dome A is characterized by a relatively high average annual accumulation rate (157 kg m −2 a −1 ), which indicates the ice cores drilled in the coastal region of the transect can provide high-resolution chemical ion and isotope records [34]. Moreover, the transect from Zhongshan Station to Dome A covers a distance of 1256 km and a wide range of environments. Different regions of the transect have diverse environments (e.g., the coast region has high accumulation rate, whereas very low accumulation rate is observed in the inland region [33]), which cause various distribution characteristics of chemical components [27]. Thus, it is necessary to collect more data and study each region individually. At present, some investigations have been carried out to reveal the spatial patterns of ions and stable isotopes along the entire transect from Zhongshan Station to Dome A [24,26,27]. However, due to limited snow chemistry data of coastal region in the past, there have only been a few systematic studies on the spatial distribution, sources, and factors that influence the deposition of the chemical components in the coastal region of the transect from Zhongshan Station to Dome A [24,26,27].
In this study, the spatial variability of the ionic and isotopic compositions of samples from surface snow and two snow pits obtained along a transect from Zhongshan Station to LGB69 was analyzed. In particular, surface snow and snow pit samples were collected in the steep coastal escarpment region within 20.6 km distance from Zhongshan Station, which were less documented in previous research [24,26,27]. The results provide a detailed illustration of the spatial distribution patterns and sources of the ions and stable isotopes, as well as factors that influence the deposition of ions and the composition of stable isotopes. This study not only fills in a gap of snow chemistry investigation within 20.6 km distance from the coast along the transect from Zhongshan Station to Dome A, but also contributes to the interpretation of long-term chemical records obtained from ice cores drilled in the coastal regions of Princess Elizabeth Land.

Sample Collection
Seventeen surface snow samples and two snow pit samples were collected along the transect from Zhongshan Station to LGB69 during the austral summer of the 18th CHINARE (2001-2002) ( Figure 1). The end site LGB69 is located at 70 • 50 06" S and 77 • 04 29" E in Princess Elizabeth Land of East Antarctica, approximately 172 km (route distance) from Zhongshan Station, at an altitude of 1860 m. The annual accumulation rate at LGB69 is as high as 255 kg m −2 a −1 based on snow pit dating and water equivalent depth. The transect from Zhongshan Station to LGB69 shows a clear difference in altitude changes between the coast and the inland region. Within 20.6 km (distance from the coast, unless otherwise mentioned), the altitude increases rapidly from 0 to 586 m, indicating a steep coastal escarpment region. Beyond 20.6 km, the altitude rises gradually from 586 to 1852 m ( Table 1). The slope was calculated as the altitude change/distance change between the two sites. Figure 2 shows the slope variation from the coast inland, as listed in Table 1. The steep coastal escarpment region within 20.6 km has an average slope of 28.4 m km −1 , whereas the interior region between 20.6 and 170 km has an average slope of 8.4 m km −1 .
Surface snow samples were collected at approximately 10-km intervals along the transect from Zhongshan Station to LGB69. One surface snow sample was collected at each site. In total, 17 surface snow samples were obtained to analyze the spatial variability of their ionic and isotopic compositions. Table 1 provides the geographical information of the surface snow samples.
Two sites were selected for the snow pit sampling: LHA006 near the coast and LGB69 located inland. The snow pit samples were collected at 2-cm intervals to determine their ionic and isotopic compositions. Samples from two snow pits were also collected to assess the spatial variability of their ionic and isotopic compositions. Table 2 lists the geographical and sampling information of the two snow pits.

Sample Analysis
The chemical samples were measured in the 100-class clean laboratory of the British Antarctic Survey. Anions and cations were analyzed simultaneously and a selection of chemical concentrations (F − , Cl − , NO 3 − , SO 4 2− , MSA, Na + , K + , Ca 2+ , and Mg 2+ ) was obtained using two Dionex ion chromatographs (DX 500, Autosampler, Chromeleon 6.6 chromatograph data system). Strict anti-contamination measures were adopted during the entire process of sample cutting, melting, and analysis. 18.2 MΩ·cm ELGA UHP water was used to clean the sample processing equipment and containers and to prepare blank samples and working standards. Analytical columns (AS18 and CS12A) and guard columns (AG18 and CG12A) were respectively used in the anion and cation measurements using an injection volume of 200 µL. The eluent conditions were chosen to result in similar run length of 12 min for both cations and anions, to allow simultaneous injection from the same autosampler vial. For the anions, isocratic elution used a 20 mM potassium hydroxide eluent at a flow rate of 1.0 mL per minute, while the cations were eluted isochratically with 20 mM methane sulphonic acid eluent at a flow rate of 1.0 mL per minute. Five multi-ion standards were prepared to bracket the range of concentrations measured in the traverse samples. The blank samples were regularly monitored during the analysis and the main chemical components of all blank samples were lower than the detection limit.
The stable isotope samples were measured on the liquid water isotope analyzer (DLT-100, Los Gatos Research, Inc., Mountain View, CA, USA) at the British Antarctic Survey. Hydrogen and oxygen isotopes were simultaneously analyzed. The DLT-100 was calibrated using working standards that bracketed the transect samples which had been prepared and calibrated against primary V-SMOW (Vienna Standard Mean Ocean Water) and V-SLAP (Vienna Standard Antarctic Precipitation) standards by the NERC Isotope Geosciences Laboratory (Keyworth, UK) using an Isoprime Isotope Ratio Mass Spectrometer. The accuracy of the δD or δ 18 O analysis (±0.8‰ and ±0.1‰, respectively) is comparable to that of traditional stable isotope mass spectrometry [35]. Some comparison samples were measured on the Isoprime mass spectrometers and the results are in good agreement with those obtained using the DLT-100.

Spatial Variability of Ionic Compositions
The concentrations of Cl − , Na + , and Mg 2+ were high in the surface snow samples, ranging from 952.48 to 1659.08, 970.78 to 1480.13, and 46.22 to 76.56 µg L −1 , respectively, between 7.4 and 20.6 km (285-586 m altitude). However, these concentrations exponentially decreased with increasing distance from the coast and altitude. The relatively low concentrations of Cl − , Na + , and Mg 2+ (ranged from 14.67 to 184.64, 7.93 to 150.00, and 1.37 to 22.05 µg L −1 , respectively) were observed between 31.5 and 170 km (737-1852 m altitude) ( Figure 3). The detailed Cl − , Na + , and Mg 2+ concentrations of the surface snow samples along the transect from Zhongshan Station to LGB69 are shown in Table 3. Moreover, because of significant correlation between distance from coast and altitude (R = 0.991, p <0.01), the variation trends of ions and stable isotopes with altitude are similar to those with distance from coast and not shown in the figures. High Cl − , Na + , and Mg 2+ concentrations occurred specifically in the steep coastal escarpment region within 20.6 km ( Figure 3). This may be related to the preferential wet or dry deposition of sea-salt aerosols. Sea-salt aerosols can be scavenged from the atmosphere by wet or dry deposition depending on the meteorological conditions. Most sea-salt aerosols can be deposited by snow precipitation in regions with high annual accumulation rates. High relative humidity conditions (60-70%) are conducive to the wet deposition of sea-salt aerosols [36]. In the present study, high accumulation rates and high relative humidity conditions (>60%) exist in the steep coastal escarpment region [33,37], leading to the preferential wet deposition of sea-salt aerosols. Alternatively, when local meteorological conditions are inadequate for wet deposition, the steep faces of the coastal escarpment can intercept the maritime air mass with abundant sea-salt aerosols. Katabatic winds are highly unidirectional and persistent and prevail in the steep coastal escarpment regions of Antarctica [37][38][39]. The intercepted maritime air masses interact strongly with the katabatic winds over the steep coastal escarpments, which leads to turbulence, thereby inducing the dry deposition of sea-salt aerosols [36]. Consequently, high Cl − , Na + , and Mg 2+ concentrations were found specifically in the steep coastal escarpment region within 20.6 km, which may be because of the preferential wet or dry deposition of sea-salt aerosols. Cl − , Na + , and Mg 2+ were also typically strongly positively correlated with one another (Table 4), which indicates they are from the same source. The sea-salt ions exhibited significant negative correlations with distance from the coast or altitude ( Table 4), suggesting that the spatial distributions of sea-salt ions from Zhongshan Station to LGB69 are mainly affected by distance from the coast and/or altitude. Cl − /Na + ratios in the surface snow samples have been used to determine the effects of sources or post-depositional processes on the deposition of sea-salt components [5,40]. In the present study, Cl − /Na + ratios were calculated from all of the surface snow samples, yielding an average ratio of 1.67 ± 0.44 from Zhongshan Station to LGB69, which was close to the ratio of sea water (1.81, W/W [41]). This suggests that sea salt is the main source of Cl − and Na + along the transect. Figure 4 shows the variation of Cl − /Na + ratios from the coast inland. However, a few low Cl − /Na + ratios occurred in the steep coastal escarpment region within 20.6 km. The lowest Cl − /Na + ratio (0.98) was recorded at site L156 at 20.6 km (586 m altitude), followed by 1.07 from site LH9601 at 7.4 km (285 m altitude), and the average of Cl − /Na + ratios was 1.06 ± 0.07 between 7.4 and 20.6 km. Beyond 20.6 km, the Cl − /Na + ratios ranged from 1.23 to 2.57 with an average of 1.80 ± 0.36, which was almost equal to that of the sea water ratio (Table 5). Previous studies at coastal Antarctic sites suggested that at relatively low accumulation rates (<40 kg m −2 a −1 ), the post-depositional loss of HCl becomes significant and reduces the Cl − concentrations in snow, resulting in a reduced Cl − /Na + ratio [5,42]. However, the present study region is characterized by a relatively high average snow accumulation rate (157 kg m −2 a −1 ) [33]. For example, the Cl − /Na + ratio of LHA006 at 17.9 km was only 1.12 (Table 5), whereas its annual accumulation rate reached as high as 200 kg m −2 a −1 ( Table 2). The post-depositional loss of HCl can therefore not account for the low Cl − /Na + ratios at these sites. Alternatively, Cl − depletion may also be related to the reaction of sea-salt aerosols (mainly NaCl) with atmospheric acids (mainly HNO 3 and H 2 SO 4 ). This reaction leads to the formation of gas-phase HCl, which lowers the initial Cl − /Na + ratios of sea-salt aerosols [5,42]. However, there was no low Cl − /Na + ratio occurred beyond 20.6 km. Thus, the reaction between sea-salt aerosols and atmospheric acids cannot be used to explain the low Cl − /Na + ratios that only occurred in the steep coastal escarpment region within 20.6 km. In the present surface snow samples, the observed low Cl − /Na + ratios within 20.6 km may be related to the steep slopes of the coastal escarpment. The variations of surface slope of the Antarctic ice sheet have a strong impact on the wind speed and direction [43]. The steep slopes of the coastal escarpment can lead to wind redistribution, inducing snow redistribution. In the process of snow redistribution by wind scouring, the Na + concentrations can be increased by the continuous input of sea-salt aerosols. However, unlike Na + , Cl − is depleted owing to the reaction between sea-salt aerosols and atmospheric acids during wind redistribution. The steep slopes of the coastal escarpment thus result in wind redistribution-induced snow redistribution, leading to an increase of Na + concentrations, thereby lowering the Cl − /Na + ratios [36]. Consequently, the low Cl − /Na + ratios within 20.6 km can be attributed to the steep slopes of the coastal escarpment. The spatial distributions of sea-salt ions in the two snow pits were also analyzed. The concentrations of Cl − , Na + , and Mg 2+ of the LHA006 snow pit ranged from 262.10 to 2660.70, 182.40 to 2205.96, and 27.66 to 223.38 µg L −1 , respectively, whereas those of the LGB69 snow pit ranged from 5.90 to 268.10, 2.02 to 195.51, and 0.53 to 28.75 µg L −1 , respectively ( Table 6). Given that LHA006 is located in the steep coastal escarpment region, the mean concentrations of these three sea-salt ions were considerably higher at LHA006 than at LGB69, which can be attributed to the preferential wet or dry deposition of seasalt aerosols. The snow pit results showed a clear difference in the mean Cl − /Na + ratios between LHA006 (1.18) and LGB69 (1.79). As with the Cl − /Na + ratio in the surface snow sample at LHA006, the mean Cl − /Na + ratio of the LHA006 snow pit samples was also low and lower than that of the LGB69 snow pit samples. Thus, the low mean Cl − /Na + ratio of the LHA006 snow pit samples may also be related to the increased Na + concentrations driven by the steep slopes of the coastal escarpment.
The MSA concentrations showed an overall decreasing trend with increasing distance from the coast and/or altitude. High MSA concentrations obtained between 7.4 and 68 km (285-1041 m altitude) ranged between 2.77 and 145.82 µg L −1 . However, the data from the interior region (78-170 km or 1157-1852 m altitude) showed relatively low MSA concentrations (ranged from 17.60 to 48.09 µg L −1 ) ( Figure 5). The detailed MSA concentrations of the surface snow samples along the transect are shown in Table 3. The MSA concentrations showed a significant relationship with distance from the coast or altitude (Table 4). This implies that distance from the coast and/or altitude are the main factors that affect the spatial distribution of MSA. The variation of nssSO 4 2− concentrations was similar to that of the MSA concentrations ( Figure 5). However, a few surface snow samples with negative and low nssSO 4 2− / SO 4 2− percentages were observed in the steep coastal escarpment region within 20.6 km (Table 5), as discussed later. The nssSO 4 2− concentrations showed high values (ranged from 118.45 to 335.38 µg L −1 ) between 31.5 and 68 km (737-1041 m altitude). The relatively low nssSO 4 2− concentrations (ranged from 58.48 to 162.24 µg L −1 ) were obtained from 78 to 170 km (1157-1852 m altitude) ( Figure 5). Table 3 shows the detailed nssSO 4 2− concentrations of the surface snow samples along the transect. When the surface snow samples with negative and low nssSO 4 2− /SO 4 2− percentages within 20.6 km were excluded, the nssSO 4 2− concentrations exhibited a significant negative correlation with distance from the coast (R = −0.624, p < 0.05) or altitude (R = −0.642, p < 0.05), implying that the variation of nssSO 4 2− is dependent both on distance from the coast and/or altitude. Moreover, the high concentrations of MSA and nssSO 4 2− occurred between 7.4 and 68 km, whereas those of sea-salt ions appeared only between 7.4 and 20.6 km. This is because super-micrometric sea-salt aerosol particles have relatively shorter atmospheric residence times and undergo more efficient atmospheric scavenging process than sub-micrometric MSA and nssSO 4 2− aerosol particles in the coastal regions [5,47,48]. The nssSO 4 2− concentrations exhibited a significant correlation with MSA concentration (Table 4), which suggests that both of them are from the atmospheric oxidation of DMS, which is produced by phytoplankton activity in the ocean.
In general, nssSO 4 2− contributes mostly to the sulfate budget [27,30]. However, in the present study, the spatial distribution pattern of nssSO 4 2− /SO 4 2− ratios was characterized by their gradual increase and stabilization from the coast inland, and a few negative and low ratios of nssSO 4 2− /SO 4 2− were observed in the steep coastal escarpment region within 20.6 km ( Table 5). Two main mechanisms have been proposed for the negative and low nssSO 4 2− /SO 4 2− percentages in coastal Antarctica. First, negative and low nssSO 4 2− /SO 4 2− percentages have been speculated to both originate from the underestimation of nssSO 4 2− . This implies that their sea-salt ions are likely from fractionated sea-salt aerosols produced by salty frost flowers on the sea ice surface at coastal Antarctica sites [2,19,[49][50][51]. The SO 4 2− /Na + ratio of salty frost flowers is lower than that of sea water owing to the brine fractionation processes on the sea ice surface, which occur when the ice temperature decreases. Mirabilite (Na 2 SO 4 ·10 H 2 O) begins to precipitate at −8 • C but NaCl precipitation only begins at −26 • C. Thus, the use of original SO 4 2− /Na + ratio from sea water in the nssSO 4 2− calculations may cause the negative and low concentrations of nssSO 4 2− in the snow samples. The low nssSO 4 2− /SO 4 2− percentages can also be attributed to the atmospheric scavenging of sea-salt aerosol particles by wet and dry removal processes. Super-micrometric ssSO 4 2− aerosol particles have relatively shorter atmospheric residence times and undergo more efficient atmospheric scavenging process than sub-micrometric biogenic nssSO 4 2− aerosol particles in the coastal regions [5,47,48]. In the present study, it is unlikely that the effects of salty frost flower fractionation can be used to explain the negative and low nssSO 4 2− /SO 4 2− percentages that occurred only in the steep coastal escarpment region within 20.6 km. Given that the nssSO 4 2− calculation was based on Na + , the increase of Na + concentrations driven by the steep slopes can affect the nssSO 4 2− calculations, resulting in negative or low nssSO 4 2− concentrations. Therefore, the negative and low nssSO 4 2− /SO 4 2− percentages in the steep coastal escarpment region within 20.6 km can be attributed to the increased Na + concentrations driven by the steep slopes of the coastal escarpment and the effects of atmospheric scavenging processes.
An inverse relationship has been reported between site temperature and MSA/nssSO 4 2− ratios in the marine atmospheric boundary layer of the southern hemisphere [52]. Higher MSA/nssSO 4 2− ratios (0.15-0.93, W/W) in the marine atmospheric boundary layer in summer were observed near coastal Antarctica, whereas lower MSA/nssSO 4 2− ratios (0.0024-0.12, W/W) were observed at mid and low latitudes [53]. Thus, the MSA/nssSO 4 2− ratios in snow samples from Antarctica can be used as a tracer of the latitudes of marine source regions of MSA and nssSO 4 2− at a given location, i.e., high MSA/nssSO 4 2− ratios in snow samples suggest that they are from a high-latitude marine source [54].
In the present study, a few abnormally negative and high MSA/nssSO 4 2− ratios were observed in the steep coastal escarpment region within 20.6 km. A negative MSA/nssSO 4 2− ratio (−0.01) was recorded at the LH9601 site at 7.4 km (285 m altitude). The highest MSA/nssSO 4 2− ratio (5.53) was recorded at the LHA006 site at 17.9 km (520 m altitude), followed by 1.09 at the L156 site at 20.6 km (586 m altitude) ( Figure 6 and Table 5). As mentioned previously, the increased Na + concentrations driven by the steep slopes of the coastal escarpment can result in negative and low nssSO 4 2− concentrations. Consequently, abnormally negative and high MSA/nssSO 4 2− ratios occurred in the steep coastal escarpment region within 20.6 km. When the abnormally negative and high MSA/nssSO 4 2− ratios within 20.6 km were removed, the mean MSA/nssSO 4 2− value of the surface snow samples from Zhongshan Station to LGB69 was 0.34 ± 0.08 (Figure 6), which was in agree-ment with the MSA/nssSO 4 2− ratios observed in the marine atmospheric boundary layer near coastal Antarctica, implying that the source regions of the MSA and nssSO 4 2− from Zhongshan Station to LGB69 are likely coastal sea areas. Ca 2+ concentrations also showed an overall decreasing trend with increasing distance from the coast and/or altitude. High Ca 2+ concentrations (ranged from 7.85 to 53.40 µg L −1 ) were measured between 7.4 and 55.2 km (285-955 m altitude). However, the data from the interior region between 68 and 170 km (1041-1852 m altitude) exhibited relatively low Ca 2+ concentrations (ranged from 1.44 to 14.76 µg L −1 ) (Figure 7). The detailed Ca 2+ concentrations of the surface snow samples along the transect are shown in Table 3. The Ca 2+ concentrations showed a significant correlation with distance from the coast or altitude ( Table 4), implying that the spatial distribution of Ca 2+ from Zhongshan Station to LGB69 is also mainly affected by distance from the coast and/or altitude. The main sources of Ca 2+ in Antarctic snow are continental dust and sea spray [55]. In the present study, the calculated nssCa 2+ /Ca 2+ percentages in the surface snow samples from Zhongshan Station to LGB69 (Table 5) suggest that continental dust is the primary source for Ca 2+ and the contribution from sea spray is relatively minor. However, similar to nssSO 4 2− , negative nssCa 2+ concentrations were also observed in the steep coastal escarpment region within 20.6 km (Table 3), which may also be because of the increased Na + concentrations driven by the steep slopes of the coastal escarpment.
The spatial distributions of non-sea-salt ions in the two snow pits were also simultaneously analyzed. The concentrations of nssSO 4 2− and MSA of the LHA006 snow pit ranged from −301.69 to 580.68 and 1.36 to 345.86 µg L −1 , respectively, and those of the LGB69 snow pit ranged from 2.24 to 96.27 and 0 to 18.08 µg L −1 , respectively (Table 6). However, the mean nssSO 4 2− /SO 4 2− percentage in the LHA006 snow pit (14%) was substantially lower than the LGB69 snow pit value (76%), which was similar to that at Hercules Névé (75%) [5] and lower than that at Talos Dome (88%) [4]. Moreover, the mean MSA/nssSO 4 2− ratio in the LHA006 snow pit (1.38) was considerably higher than that of the LGB69 snow pit (0.22). These large differences in the mean nssSO 4 2− /SO 4 2− percentage and mean MSA/nssSO 4 2− ratio between the LHA006 and LGB69 snow pits can be explained by the increased Na + concentrations driven by the steep slopes of the coastal escarpment and effects of atmospheric scavenging processes.
The concentrations of Ca 2+ of the LHA006 snow pit ranged from 8.87 to 110.00 µg L −1 , whereas those of the LGB69 snow pit ranged from 2.13 to 101.53 µg L −1 (Table 6). However, the calculated mean concentration of nssCa 2+ showed a negative value (−1.46 ± 17.20 µg L −1 ) in the LHA006 snow pit samples ( Table 6). As with the negative nssCa 2+ concentration in the surface snow sample at LHA006, the negative mean concentration of nssCa 2+ in the LHA006 snow pit samples can also be attributed to the increased Na + concentrations driven by the steep slopes of the coastal escarpment. The mean nssCa 2+ /Ca 2+ percentage of the LGB69 snow pit samples (88%) shows that the nssCa 2+ contributes mostly to the total Ca 2+ .

Spatial Variability of Isotopic Compositions
The spatial distributions of δD, δ 18 O, and d-excess from in the surface snow samples were analyzed. The second-order parameter d-excess was calculated following [56]: Both the δD and δ 18 O values in the surface snow samples decreased with increasing distance from the coast and/or altitude (Figure 8), exhibiting significant negative correlations ( Table 7). The detailed δD and δ 18 O values in the surface snow samples along the transect are shown in Table 3. The lapse rates of δD and δ 18 O with distance from the coast were respectively −54.17‰/100 km and −6.84‰/100 km ( Table 8). The topographic features were revealed by the contour lines in Figure 1 and showed that the slope from Zhongshan Station to LGB69 was steeper than that between LGB69 and Dome A, causing the rapid decreases of δD and δ 18 O values. The lapse rates of δD and δ 18 O with distance from the coast along the transect from Zhongshan Station to LGB69 were therefore considerably larger than those along the transect from Zhongshan Station to Dome A. The lapse rates of δD and δ 18 O with altitude from Zhongshan Station to LGB69 were respectively −5.99‰/100 m and −0.75‰/100 m, and lower than those from Zhongshan Station to Dome A (Table 8). This can still be attributed to the topographic difference between the transect from Zhongshan Station to LGB69 and that from LGB69 to Dome A.
The relationship between δD and δ 18 O values in the surface snow samples is given as δD = (7.90 ± 0.17) δ 18 O − (0.09 ± 3.33) (R 2 = 0.99, n = 17),   The slope between δD and δ 18 O was 7.90 ± 0.17 and lower than that of GMWL, but higher than the values reported by Li et al. (7.78 ± 0.04) [27] and Xiao et al. (7.50 ± 0.1) [26]. This indicates that more intense kinetic fractionation processes occur during snow formation in the inland between LGB69 and Dome A The d-excess values did not show a significant relationship with distance from the coast (R = 0.111, p = 0.67) or altitude (R = 0.064, p = 0.80) ( Table 7), indicating that the spatial distribution of d-excess from Zhongshan Station to LGB69 is mainly affected by the moisture source region. The gradients of d-excess with distance from the coast and altitude were 0.53‰/100 km and 0.03‰/100 m, respectively, which were substantially lower than values reported by Xiao et al. [26] ( Table 8). The d-excess varied from −4.47‰ to 4.97‰ (Table 8), a considerably smaller range than those reported by Li et al. [27] and Xiao et al. [26], from Zhongshan Station to Dome A. The variation of d-excess was stable, which suggests that the transect from Zhongshan Station to LGB69 has the same moisture source region.
The mean values of δD, δ 18 O, and d-excess in the snow pit samples were simultaneously analyzed. The δD and δ 18 Table 9). The d-excess values at LHA006 ranged between −1.25 and 10.88‰ with a mean value of 3.32 ± 2.53‰. The d-excess values at LGB69 ranged between −6.64 and 9.33‰ with a mean value of 2.56 ± 3.59‰ (Table 9). Previous studies noted that the d-excess values in Antarctica are low (~5‰) and stable below 2000 m altitude, which implies that the moisture source regions are coastal sea areas [32,57]. In our study, the altitudes of both LHA006 and LGB69 are below 2000 m. Furthermore, the mean d-excess values in the LHA006 and LGB69 snow pits were low and similar, indicating that the precipitation at both LHA006 and LGB69 originates from the same moisture source region (i.e., coastal sea area).

Conclusions
This study reports the spatial distribution patterns and sources of ions and stable isotopes as well as factors that influence the deposition of major ions and the composition of stable isotopes in surface snow and snow pit samples obtained along a transect from Zhongshan Station to LGB69, East Antarctica.
The spatial distributions of sea-salt ions (Cl − , Na + , and Mg 2+ ) along the transect showed that sea-salt ions exponentially decreased with increasing distance from the coast or altitude. The distance from the coast and/or altitude are the two major factors that influence the spatial distributions of sea-salt ions. High Cl − , Na + , and Mg 2+ concentrations that occurred specifically in the steep coastal escarpment region within 20.6 km can be attributed to the preferential wet or dry deposition of sea-salt aerosols. The average Cl − /Na + ratio of the surface snow samples from Zhongshan Station to LGB69 was 1.67 ± 0.44, which suggests that sea salt is the main source of Cl − and Na + along the transect. However, the steep slopes of the coastal escarpment can result in wind redistribution-induced snow redistribution, leading to increased Na + concentrations. Thus, the Cl − /Na + ratios, nssSO 4 2− , MSA/nssSO 4 2− ratios, and nssCa 2+ values calculated based on Na + showed a few abnormal values in the steep coastal escarpment region within 20.6 km. MSA, nssSO 4 2− , and Ca 2+ concentrations showed overall decreasing trends with increasing distance from the coast and/or altitude. The distance from the coast and/or altitude are the major influencing factors of the spatial distributions of MSA, nssSO 4 2− , and Ca 2+ . The similar spatial variation and significant correlation between MSA and nssSO 4 2− in the surface snow indicate their common source, i.e., biogenic activity. When the abnormally negative and high MSA/nssSO 4 2− ratios within 20.6 km were removed, the mean MSA/nssSO 4 2− value of surface snow samples was 0.34 ± 0.08, which suggests that their source regions are likely coastal sea areas. The nssCa 2+ /Ca 2+ percentages in the surface snow and LGB69 snow pit samples indicated that continental dust is the primary source for Ca 2+ along the transect from Zhongshan Station to LGB69. In addition, δD and δ 18 O values decreased from the coast inland. Both δD and δ 18 O in the surface snow exhibited significant negative correlations with distance from the coast and/or altitude. The lapse rates of δD and δ 18 O from Zhongshan Station to LGB69 were larger than those from Zhongshan Station to Dome A because of the steeper slope of the former transect. The d-excess values did not show a significant relationship with distance from the coast or altitude, indicating that the moisture source region is the major factor that affected the spatial distribution of d-excess from Zhongshan Station to LGB69. The variation of d-excess in the surface snow samples was stable and the mean d-excess values in the LHA006 and LGB69 snow pits were low and similar (3.32 ± 2.53‰ and 2.56 ± 3.59‰, respectively), which indicates that the moisture source region between Zhongshan Station and LGB69 is the coastal sea area.
The glaciochemistry results of the surface snow and snow pit samples along the transect from Zhongshan Station to LGB69 suggest that the ionic and moisture source region is the coastal sea area, which will be useful for interpreting the temporal variability of chemical ions and isotopes from ice cores drilled in the coastal regions of Princess Elizabeth Land.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author, upon reasonable request.