Seasonal variations of the near surfacelayer parameters over the Antarctic ice sheet in Princess Elizabeth Land,East Antarctica

Analysis of sensible heat flux （ Qh ）, latent heat flux （ Qe ）, Richardson number （Ri） ,bulk transport coefficient （Cd） and katabatic windsare presented by using the meteorological data in the near surface layer from an automatic weather station （AWS） in Princess Elizabeth Land, East Antarctica ice sheet and the data of corresponding period at Zhongshan station in 2002. It shows that annual mean air temperature at LGB69 is -25.6℃, which is 16.4℃ lower than that at Zhongshan, where the elevation is lower and located on the coast. The temperature lapse rate is about 1.0℃/110 m for the initial from coast to inland. The turbulence heat flux at LGB69 displays obvious seasonal variations with the average sensible heat flux -17.9 W/m^2 and latent heat flux -0.9 W/m^2. The intensity （Qh ＋ Qe ） of coolling source is - 18.8 W/m^2 meaning the snow surface layer obtains heat from atmosphere. The near surface atmosphere is near-neutral stratified with bulk transport coefficients （Cd） around 2.8 ×10^-3 ,and it is near constant when the wind speed higher than 8 m/s. The speed and the frequency of easterly Katabatic winds at LGB69 were higher than that at Zhongshan Station.

Some recent advances in remote sensing-based monitoring of changes in the Greenland Ice Sheet

The mass balance of the Greenland Ice Sheet(GrIS)plays a crucial role in global sea level change.Since the 1960s,remote sensing missions have been providing extensive and continuous observation data for change monitoring of the GrIS.In this paper,we present our recent research results from remote sensing-based GrIS change monitoring.First,historical satellite data are processed and used to fill data gaps and are combined with existing partial maps,completing an ice velocity map of the GrIS from the 1960s to 1980s.This map provides valuable data for estimating the historical mass balance of Greenland.Second,the monthly gravimetry-based mass balance of the GrIS from 2002 to 2020 is estimated by combining Gravity Recovery and Climate Experiment(GRACE)and GRACE Follow On(GRACE-FO)data.It is found that the GrIS has lost a total mass of approximately 4443±75 Gt during this period.Third,based on Global Land Ice Measurements from Space(GLIMS),an updated Greenland glacier inventory is achieved utilizing data collected between 2006 and 2020.This inventory provides more detailed and up-to-data glacier boundaries of Greenland.Overall,these advances provide essential data support for estimating the mass balance of the GrIS,contributing to the advancement of research on global sea level change.

Migration and role of zinc in biogeochemical cycles in the Antarctic Ice Sheet and the Southern Ocean

Zinc(Zn),a widespread metal in the Earth’s crust,serves as a crucial nutrient in the Southern Ocean’s primary production.Studies on Zn in Antarctic snow and ice offer insights into the origins of this metal and its transport routes,as well as its impact on the biogeochemical processes within the Antarctic atmosphere–land–ocean system.This review examines research on the spatial and temporal distribution of Zn in Antarctic snow and ice,as well as in Southern Ocean waters.It includes an overview of advanced methods for sampling and analyzing Zn,along with explanations for the observed variations.The review also discusses various sources of Zn as a nutrient to the Southern Ocean.Finally,it addresses prospective issues related to the use of Zn isotopes in identifying atmospheric sources and their biogeochemical effects on the development of the Southern Ocean ecosystem.

Oscillating Performance of A Rigid Buoy Floating Between Two Ice Sheets with Finite Lengths

Global warming has led to major melting of ice in the polar Arctic,making it possible to open Arctic shipping lanes.In this case,the large number of ice sheets are extremely dangerous for ship navigation,so in this paper,a body floating on water confined between two finite ice sheets is investigated.The linearized potential flow theory is adopted,and water is considered an incompressible ideal fluid with a finite depth of the fluid domain.The ice sheets are treated as elastic plates,and the problem is solved by matching eigenfunction expansion.The fluid domain is divided into subregions on the basis of the water surface conditions,and the velocity potential of the subdomains is expanded via the separated variable method.By utilizing the continuity of pressure and velocity at the interfaces of two neighboring regions,a system of linear equations is established to obtain the unknown coefficients in the expansion,which in turn leads to analytical solutions for different motion modes in different regions.The effects of different structural drafts,and different lengths of ice sheets on both sides,etc.,on the hydrodynamic characteristics of floats are analyzed.The amplitude of motion of the float is explored,as is the wave elevation between the ice sheets and the float.

Methods for ice velocity mapping in West Antarctica using historical optical satellite images

The Antarctic Ice Sheet harbors more than 90%of the Earth ice mass,with significant losses experienced through dynamic thinning,particularly in West Antarctica.The crucial aspect of investigating ice mass balance in historical periods preceding 1990 hinges on the utilization of ice velocities derived from optical satellite images.We employed declassified satellite images and Landsat images with normalized cross correlation based image matching,adopting an adaptive combination of skills and methods to overcome challenges encountered during the mapping of historical ice velocity in West Antarctica.A basin-wide synthesis velocity map encompassing the coastal regions of most large-scale glaciers and ice shelves in West Antarctica has already been successfully generated.Our results for historical ice velocities cover over 70%of the grounding line in most of the West Antarctic basins.Through adjustments,we uncovered overestimations in ice velocity measurements over an extended period,transforming our ice velocity map into a spatially deterministic,temporally average version.Among all velocity measurements,Thwaites Glacier exhibited a notable spatial variation in the fastest ice flowline and velocity distribution.Overestimation distributions on Thwaites Glacier displayed a clear consistency with the positions of subsequent front calving events,offering insights into the instabilities of ice shelves.

Ice thickness,internal layers,and surface and subglacial topography in the vicinity of Chinese Antarctic Taishan station in Princess Elizabeth Land,East Antarctica

We present the results of two ground-based radio-echo-sounding（RES） and GPS surveys performed in the vicinity of new Chinese Taishan station,Princess Elizabeth Land,East Antarctica,obtained in two austral summers during CHINARE 21（2004/2005） and CHINARE 29（2012/2013）.The radar surveys measured ice thickness and internal layers using 60- and 150-MHz radar systems,and GPS measurements showed smooth surface slopes around the station with altitudes of 2607-2636 m above sea level（a.s.l.）.Radar profiles indicate an average ice thickness of 1900 m,with a maximum of 1949 m and a minimum of 1856 m,within a square area measuring approximately 2 km × 2 km in the vicinity of the station.The ice thickness beneath the station site is 1870 m.The subglacial landscape beneath the station is quiet sharp and ranges from 662 to 770 m a.s.l.,revealing part of a mountainous topography.The ice volume in the grid is estimated to be 7.6 km^3.Along a 60-MHz radar profile with a length of 17.6 km at the region covering the station site,some disturbed internal layers are identified and traced;the geometry of internal layers within the englacial stratigraphy may imply a complex depositional process in the area.

Distribution of ice thickness and subglacial topography of the "Chinese Wall" around Kunlun Station,East Antarctica

As fundamental parameters of the Antarctic Ice Sheet,ice thickness and subglacial topography are critical factors for studying the basal conditions and mass balance in Antarctica.During CHINARE 24（the 24 th Chinese National Antarctic Research Expedition,2007/08）,the research team used a deep ice-penetrating radar system to measure the ice thickness and subglacial topography of the ＂Chinese Wall＂ around Kunlun Station,East Antarctica.Preliminary results show that the ice thickness varies mostly from 1600 m to 2800 m along the ＂Chinese Wall＂,with the thickest ice being 3444 m,and the thinnest ice 1255 m.The average bedrock elevation is 1722 m,while the minimum is just 604 m.Compared with the northern side of the ice divide,the ice thickness is a little greater and the subglacial topography lower on the southern side,which is also characterized by four deep valleys.We found no basal freeze-on ice in the Gamburtsev Subglacial Mountains area,subglacial lakes,or water bodies along the ＂Chinese Wall＂.Ice thickness and subglacial topography data extracted from the Bedmap 2 database along the ＂Chinese Wall＂ are consistent with our results,but their resolution and accuracy are very limited in areas where the bedrock fluctuates intensely.The distribution of ice thickness and subglacial topography detected by ice-penetrating radar clarifies the features of the ice sheet in this ＂inaccessible＂ region.These results will help to advance the study of ice sheet dynamics and the determination of future locations of the GSM＇s geological and deep ice core drilling sites in the Dome A region.

Greenland Ice Sheet Contribution to Future Global Sea Level Rise based on CMIP5 Models

Sea level rise （SLR） is one of the major socioeconomic risks associated with global warming. Mass losses from the Greenland ice sheet （GrIS） will be partially responsible for future SLR, although there are large uncertainties in modeled climate and ice sheet behavior. We used the ice sheet model SICOPOLIS （Simulation COde for POLythermal Ice Sheets） driven by climate projections from 20 models in the fifth phase of the Coupled Model Intercomparison Project （CMIP5） to estimate the GrlS contribution to global SLR. Based on the outputs of the 20 models, it is estimated that the GrIS will contribute 0-16 （0-27） cm to global SLR by 2100 under the Representative Concentration Pathways （RCP） 4.5 （RCP 8.5） scenarios. The projected SLR increases further to 7-22 （7-33） cm with 2~basal sliding included. In response to the results of the multimodel ensemble mean, the ice sheet model projects a global SLR of 3 cm and 7 cm （10 cm and 13 cm with 2~basal sliding） under the RCP 4.5 and RCP 8.5 scenarios, respectively. In addition, our results suggest that the uncertainty in future sea level projection caused by the large spread in climate projections could be reduced with model-evaluation and the selective use of model outputs.

The extraction of elevation information of ice-sheet surface on south area of the Larsemann Hills in East Antarctica

Based on the setting up of the function relation between the radiant brightness of ice surface and the elevation control point in Antarctica, the experiments for the extraction of elevation information by using thermal infrared image of Landsat TM band 6, and the ice surface topographical maps in area nearby Larsemann Hills have been performed.

Research progress in geophysical exploration of the Antarctic ice sheet

The Antarctic ice sheet is an important target of Antarctic research.Thickness and structure,including intraice and subice,are closely related to the mass balance of the ice sheet,and play an important role in the study of global sea level and climate change.Subglacial topography is an important basis for studying ice sheet dynamics and ice sheet evolution.This paper briefly reviews the geophysical detection methods and research status of the Antarctic ice sheet:(1)Conventional methods such as ice radar are the main methods for studying the ice sheet today,and passive source seismic methods such as the receiver function method,H/V method and P-wave coda autocorrelation method have good development prospects;(2)the high-resolution(1 km)ice thickness and subglacial topographic database BEDMAP2 established based on various data has greatly improved the ability to detect internal isochronous layers,anisotropic layers,and temperature changes within ice and has advanced research on ice sheet evolution;and(3)ice radar,numerical simulation and core drilling are the main methods to study subglacial lakes and sediments.More than 400 subglacial lakes have been confirmed,and more than 12000 simulation results have been obtained.Research on the Antarctic ice sheet faces enormous challenges and is of great urgency.Aiming at hot issues,such as Antarctic geological evolution,glacial retreat,ice sheet melting and their relationships with global climate change,it is the frontier and trend of future Antarctic ice sheet research to carry out multidisciplinary and multicountry comprehensive geophysical exploration based on the traditional ice radar method combined with passive seismic methods,especially new technologies such as short-period dense array technology,unmanned aerial vehicles and artificial intelligence.This is expected to further promote Antarctic research.

Accumulation over the Greenland Ice Sheet as Represented in Reanalysis Data

Annual precipitation,evaporation,and calculated accumulation from reanalysis model outputs have been investigated for the Greenland Ice Sheet （GrIS）,based on the common period of 1989-2001.The ERA-40 and ERA-interim reanalysis data showed better agreement with observations than do NCEP-1 and NCEP-2 reanalyses.Further,ERA-interim showed the closest spatial distribution of accumulation to the observation.Concerning temporal variations,ERA-interim showed the best correlation with precipitation observations at five synoptic stations,and the best correlation with in situ measurements of accumulation at nine ice core sites.The mean annual precipitation averaged over the whole GrIS from ERA-interim （363 mm yr 1） and mean annual accumulation （319 mm yr 1） are very close to the observations.The validation of accumulation calculated from reanalysis data against ice-core measurements suggests that further improvements to reanalysis models are needed.

Chinese radioglaciological studies on the Antarctic ice sheet: progress and prospects

Chinese radioglaciological studies on the Antarctic ice sheet(AIS) began in 2004/05 when the 21 st Chinese National Antarctic Research Expedition(CHINARE 21) team arrived at Dome A for the first time and radio echo sounding(RES) was conducted along the inland traverse and in the Dome A region. Subsequently, more field surveys were conducted along the traverse and in the Dome A region using different radar systems targeting different scientific purposes, such as revealing the landscape of the Gamburtsev Subglacial Mountains by detailed grid RES, or locating a deep ice core drilling site by mapping and studying internal structures, bedrock topography and subglacial conditions in the Dome A region. Furthermore, the evolution of the AIS was inferred from the typical mountain glaciation topography beneath Dome A, and the age of the deep ice core at Kunlun Station was estimated through numerical modeling. Recently, the Snow Eagle 601 airplane was acquired and an airborne geophysical system was constructed to survey the AIS in Princess Elizabeth Land during CHINARE 32(2015/16) and CHINARE 33(2016/17) in order to fill the large data gap there. In this paper, we review both the recent progress of Chinese radioglaciological science in Antarctica and future proposed work.

Ice sheet controls on fine-grained deposition at the southern Mendeleev Ridge since the penultimate interglacial

Clay minerals deposited at the southern Mendeleev Ridge in the Arctic Ocean have a unique provenance,which can be used to reconstruct changes in the local sedimentary environment.We show that sediments in core ARC7-E23 record high-frequency changes in clay minerals since the penultimate interglacial.The clay minerals,grain size,and ice-rafted debris indicate the extent of the East Siberia Ice Sheet(ESIS).During the glacial periods of Marine Isotope Stage 2(MIS2)and MIS4,the southern Mendeleev Ridge was likely covered by an ESIS-extended ice shelf,blocking almost all sediment input from the Canadian Arctic and Laptev Sea,but allowing transport of fine-grained sediments from the East Siberian and Chukchi Sea shelves.After ESIS retreat,the Beaufort Gyre and Transpolar Drift became the primary transport mechanism for the distally sourced sediments.Climate conditions in MIS3 enhanced both the oceanic circulation and sediment transport.

Greenland Ice Sheet Elevation Change in Winter and Influence of Atmospheric Teleconnections in the Northern Hemisphere

The relationship between the variability of the surface elevation of the Greenland Ice Sheet (GIS) in winter and sea level pressure is identified through analysis of data from satellite-borne radar altimeters, together with meteorological data fields during 1993 2005. We found that both the North Pacific Oscillation (NPO) and the North Atlantic Oscillation (NAO), the two major teleconnection patterns of the atmospheric surface pressure fields in the Northern Hemisphere, significantly influence the GIS winter elevation change. Further, it is suggested that the NPO may affect the GIS accumulation by influencing the NAO, particularly by changing the intensity and location of the Icelandic Low.

Evolution of the North Atlantic Current and Barents Ice Sheet as revealed by grain size populations in the northern Norwegian Sea during the last 60 ka

The grain size distribution of bulk sediment samples was decomposed in a core to reconstruct paleoceanographic evolution over the past 60 ka in the northern Norwegian Sea.The results show that sediments consisted of 3–4 grain populations derived from the North Atlantic Current(NAC)and Barents Ice Sheet(BIS).The grain size data suggest three palaeoceanographic evolution stages:(1)an environment affected by BIS and NAC and changed with the interstadial/stadial transition in phase with the Greenland ice-core record at 60–31 ka BP,during which discharge of icebergs and the content of the coarsest population containing ice-rafted debris(IRD)in the sediments increased significantly during stadial,while the fine silt population containing volcanic glasses increased with the enhancement of NAC during the interstadial;(2)an extreme environment controlled by BIS at31–13 ka BP.BIS reached to its maximum at about 31 ka BP and the turbid plumes that formed at the leading edge of BIS contributed to a significant increase in the clayey population in sediments.Icebergs drained into the northern Norwegian Sea with periodical calving of the BIS at 31–19 ka BP.Subsequently,the ablation of the BIS discharged massive floods with clayey sediments and icebergs into the Norwegian Sea at 19–13 ka BP,resulting in a constant increase in clay and IRD in sediments;and(3)a marine environment similar to the present one under the strong influence of NAC following the complete melting of the BIS after 13 ka BP,NAC is the dominant transport agent and no IRD occurred in sediments.The fine silt populations containing volcanic glasses transported by NAC significantly increased.

A frozen collision belt beneath ice:an overview of seismic studies around the Gamburtsev Subglacial Mountains,East Antarctica

The Gamburtsev Antarctic Mountains Seismic Experiment （GAMSEIS, 2007-2010） was jointly conducted by the United States, China, and Japan during and after the International Polar Year 2007-2008. Broadband seismic stations were deployed across the ice-covered Gamburtsev Subglacial Mountains （GSM） and other previously unexplored areas in the interior of East Antarctica. Using GAMSEIS data, published results not only have revealed the deep structure of Antarctica, but also improved our understanding of the tectonic evolution of Antarctica and the supercontinent Gondwana, and of the relationship between geothermal heat flux and glaciers. This contribution draws together the major findings from recent studies, and also offers further investigation into the relationship between tectonic history and the East Antarctic Ice Sheet. The elevation of the GSM is largely supported by thickened crust, with Moho depths of ~60 km near the crest of the range. The GSM are underlain by thick （〉200 km） and cold continental lithosphere that likely formed after collision of two ancient crustal blocks during the pan-African orogeny. Heat flux results obtained from seismic inversion support a model of ice sheet basal melting that depends more strongly on bedrock topography than on geothermal heat flux, while ice surface and ice thickness are inversely correlated with heat flux.

Persistently explore new understandings of Greenland Ice Sheet ablation

Greenland and Antarctic ice sheets are the largest potential contributors to global sea level rise(GSLR),amounting to more than 64 m of sea level equivalence(SLE).Between the two,Greenland Ice Sheet(GrIS)alone comprises about 7 mSLE,with a much faster speed of ablation than the Antarctic Ice Sheet.

Mass change of the Antarctic ice sheet inferred from ICESat and CryoSat-2

This study examined the mass change of the Antarctic ice sheet(AIS) based on ICESat and CryoSat-2 observations. We estimated the AIS exhibited mass losses of-101±15 Gt·aduring the ICESat period(Sept–Nov 2003 to Sept–Oct 2009) and-186±55 Gt·aduring the CryoSat-2 period(Jan 2011 to Dec 2015). Mass losses occurred mainly in the sectors of the Amundsen and Bellingshausen seas. Benefitting from the 30-d subcycle of CryoSat-2, we obtained monthly estimates of mass evolution. Considerable annual variations were observed in the mass evolution sequences and the climatological monthly mass evolution. Seasonal mass evolutions in the sectors of the Bellingshausen and Amundsen seas were found most representative of the annual variation. The geographical distribution characteristics of interannual AIS mass evolution were revealed by the annual average mass evolution sequences. During Jan 2011 to Dec 2015, the ice sheets in the sectors of the Bellingshausen and Amundsen seas, and the Totten Glacier, experienced increasingly rapid areal mass loss. An area of mass gain with a moderate rate of increase was found between Dronning Maud Land and Enderby Land. Rapid mass accumulation has occurred in a limited area of the Kamb Ice Stream.

Preliminary research on the transmission path of nssSO_4^(2-) and NO_3^- in Antarctic ice sheet

The main sources of nssSO 2- 4 and NO - 3 were summarized in this paper. By analyzing the spatial distribution features of major ions in Antarctic ice sheet and studying on the different time of the same volcanic event recorded by different ice cores from different regions in Antarctica, this paper intends to study the transmission path of nssSO 2- 4 and NO - 3. Results show that nssSO 2- 4 and NO - 3 are transmitted to the ice sheet through long distance and high altitude. The procedure of the transmission is that nssSO 2- 4 and NO - 3 are transmitted to the level between the top of troposphere and the bottom of stratosphere, then subsided to the ice sheet surface and spread to other regions.

An Improved Method for Modeling Spatial Distribution of δD in Surface Snow over Antarctic Ice Sheet

Using the recent compilation of the isotopic composition data of surface snow of Antarctic ice sheet, we proposed an improved interpolation method of δD, which utilizes geographical factors （i.e., latitude and altitude） as the primary predictors and incorporates inverse distance weighting （IDW） technique. The method was applied to a high-resolution digital elevation model （DEM） to produce a grid map of multi-year mean δD values with lkm spatial resolution for Antarctica. The mean absolute deviation between observed and estimated data in the map is about 5.4‰, and the standard deviation is 9‰. The resulting δD pattern resembles well known characteristics such as the depletion of the heavy isotopes with increasing latitude and distance from coast line, but also reveals the complex topographic effects.