Search

This paper aims to solve the limitations of traditional gravity physical property inversion methods such as insufficient depth resolution and difficulties in parameter selection, by proposing an improved 3D gravity inversion method based on deep learning. The deep learning network model is established using the fully convolutional U-net network. To enhance the generalization ability of the sample set, the large-scale training set and test set are generated by the random walk, based on the forward theory. Founded on the traditional loss function's definition, this paper introduces an improvement incorporating a physical constraint to measure the degree of data fitting between the predicted and the real gravity data. This improvement significantly boosted the accuracy of the deep learning inversion method, as verified through both a single model and an intricate combination model. Finally, we applied this improved inversion method to the gravity data from the Gamburtsev Subglacial Mountains in the interior of East Antarctica, obtaining a comprehensive 3D crustal density structure. The results provide new evidence for the presence of a dense crustal root situated beneath the central Gamburtsev Province near the Gamburtsev Suture. [ABSTRACT FROM AUTHOR]

The thickness of the Antarctic ice sheet is a crucial parameter for inferring glacier mass and its evolution process. In the literature, the gravity method has been proven to be one of the effective means for estimating ice sheet thickness. And it is a preferred approach when direct measurements are not available. However, few gravity inversion methods are valid in rugged terrain areas with undulating observation surfaces (UOSs). To solve this problem, this paper proposes an improved high-precision 3D density interface inversion method considering terrain and UOSs simultaneously. The proposed method utilizes airborne gravity data at their flight altitudes, instead of the continued data yield from the unstable downward continuation procedure. In addition, based on the undulating right rectangular prism model, the large reliefs of the terrain are included in the iterative inversion. The proposed method is verified on two synthetic examples and is successfully applied to real data in East Antarctica. [ABSTRACT FROM AUTHOR]

Antarctica is of great importance in terms of global warming, the sustainability of resources, and the conservation of biodiversity. However, due to 99.66% of the continent being covered in ice and snow, geological research and geoscientific study in Antarctica face huge challenges. Geophysical surveys play a crucial role in enhancing comprehension of the fundamental structure of Antarctica. This study used bibliometric analysis to analyze citation data retrieved from the Web of Science for the period from 1982 to 2022 with geophysical research on Antarctica as the topic. According to the analysis results, the amount of Antarctic geophysical research has been steadily growing over the past four decades as related research countries/regions have become increasingly invested in issues pertaining to global warming and sustainability, and international cooperation is in sight. Moreover, based on keyword clustering and an analysis of highly cited papers, six popular research topics have been identified: Antarctic ice sheet instability and sea level change, Southern Ocean and Sea Ice, tectonic activity of the West Antarctic rift system, the paleocontinental rift and reorganization, magmatism and volcanism, and subglacial lakes and subglacial hydrology. This paper provides a detailed overview of these popular research topics and discusses the applications and advantages of the geophysical methods used in each field. Finally, based on keywords regarding abrupt changes, we identify and examine the thematic evolution of the nexus over three consecutive sub-periods (i.e., 1990–1995, 1996–2005, and 2006–2022). The relevance of using geophysics to support numerous and diverse scientific activities in Antarctica becomes very clear after analyzing this set of scientific publications, as is the importance of using multiple geophysical methods (satellite, airborne, surface, and borehole technology) to revolutionize the acquisition of new data in greater detail from inaccessible or hard-to-reach areas. Many of the advances that they have enabled be seen in the Antarctic terrestrial areas (detailed mapping of the geological structures of West and East Antarctica), ice, and snow (tracking glaciers and sea ice, along with the depth and features of ice sheets). These valuable results help identify potential future research opportunities in the field of Antarctic geophysical research and aid academic professionals in keeping up with recent advances. [ABSTRACT FROM AUTHOR]

Remote sensing is a very powerful tool that has been used to identify, map and monitor Antarctic features and processes for nearly one century. Satellite remote sensing plays the main role for about the last five decades, as it is the only way to provide multitemporal views at continental scale. But the emergence of small consumer-grade unoccupied aerial vehicles (UAVs) over the past two decades has paved the way for data in unprecedented detail. This has been also verified by an increasing noticeable interest in Antarctica by the incorporation of UAVs in the field activities in diversified research topics. This paper presents a comprehensive review about the use of UAVs in scientific activities in Antarctica. It is based on the analysis of 190 scientific publications published in peer-reviewed journals and proceedings of conferences which are organised into six main application topics: Terrestrial, Ice and Snow, Fauna, Technology, Atmosphere and Others. The analysis encompasses a detailed overview of the activities, identifying advantages and difficulties, also evaluating future possibilities and challenges for expanding the use of UAV in the field activities. The relevance of using UAVs to support numerous and diverse scientific activities in Antarctica becomes very clear after analysing this set of scientific publications, as it is revolutionising the remote acquisition of new data with much higher detail, from inaccessible or difficult to access regions, in faster and cheaper ways. Many of the advances can be seen in the terrestrial areas (detailed 3D mapping; vegetation mapping, discrimination and health assessment; periglacial forms characterisation), ice and snow (more detailed topography, depth and features of ice-sheets, glaciers and sea-ice), fauna (counting penguins, seals and flying birds and detailed morphometrics) and in atmosphere studies (more detailed meteorological measurements and air-surface couplings). This review has also shown that despite the low environmental impact of UAV-based surveys, the increasing number of applications and use, may lead to impacts in the most sensitive Antarctic ecosystems. Hence, we call for an internationally coordinated effort to for planning and sharing UAV data in Antarctica, which would reduce environmental impacts, while extending research outcomes. [ABSTRACT FROM AUTHOR]

This is one of the first papers to study the ionospheric effects of two solar eclipses that occurred in South America and Antarctica under geomagnetic activity in different seasons (summer and autumn) and their impact on the equatorial ionization anomaly (EIA). The changes in total electron content (TEC) during the 15 February 2018 and 30 April 2022 partial solar eclipses will be analyzed. The study is based on more than 390 GPS stations, Swarm-A, and DMSP F18 satellite measurements, such as TEC, electron density, and electron temperature. The ionospheric behaviors over the two-fifth days on both sides of each eclipse were used as a reference for estimating TEC changes. Regional TEC maps were created for the analysis. Background TEC levels were significantly higher during the 2022 eclipse than during the 2018 eclipse because ionospheric levels depend on solar index parameters. On the days of the 2018 and 2022 eclipses, the ionospheric enhancement was noticeable due to levels of geomagnetic activity. Although geomagnetic forcing impacted the ionosphere, both eclipses had evident depletions under the penumbra, wherein differential vertical TEC (DVTEC) reached values <−40%. The duration of the ionospheric effects persisted after 24 UT. Also, while a noticeable TEC depletion (DVTEC ∼−50%) of the southern EIA crest was observed during the 2018 eclipse (hemisphere summer), an evident TEC enhancement (DVTEC > 30%) at the same crest was seen during the eclipse of 2022 (hemisphere autumn). Swarm-A and DMSP F18 satellite measurements and analysis of other solar eclipses in the sector under quiet conditions supported the ionospheric behavior. [ABSTRACT FROM AUTHOR]

Antarctic mapping with satellite images is an important basic task for polar environmental monitoring. Since the first Chinese marine satellite was launched in 2002, China has formed three series of more than 10 marine satellites in orbit. As global operational monitoring satellites of ocean color series, HY-1C and HY-1D have good coverage characteristics and imaging performance in polar regions, and they provide an effective tool for Antarctic monitoring and mapping. In this paper, Antarctic images acquired by the HY-1 satellite Coastal Zone Imager (CZI) sensor were used to study color matching in the mosaic process. According to the CZI characteristics for Antarctic imaging, experiments were carried out on the illuminance nonuniformity of a single image and color registration of multiple images. A gray-level segmentation color-matching method is proposed to solve the problem of image overstretching in the Antarctic image color-matching process. The results and statistical analysis show that the proposed method can effectively eliminate the color deviation between HY-1 Antarctic images, and the mosaic results have a good effect. [ABSTRACT FROM AUTHOR]

The paper presents the first retrievals of clean snow properties using spaceborne hyperspectral observations via the Environmental Mapping and Analysis Program (EnMAP). The location close to the Concordia station at the Dome C Plateau (Antarctica) was selected. At this location, the atmospheric effects (except molecular light scattering and absorption) are weak, and the simplified atmospheric correction scheme could be applied. The ice grain size, snow specific surface area, and snow spectral and broadband albedos were retrieved using single-view EnMAP measurements. In addition, we propose a technique to retrieve trace gas concentrations (e.g., water vapor and ozone) from EnMAP observations over the snow surfaces. A close correspondence of satellite and ground-measured parameters was found. [ABSTRACT FROM AUTHOR]

As the largest species of birds in Antarctica, penguins are called "biological indicators". Changes in the environment will cause population fluctuations. Therefore, developing a penguin census regularly will not only help carry out conservation activities but also provides a basis for studying climate change. Traditionally, scholars often use indirect methods, e.g., identifying penguin guano and establishing regression relationships to estimate the size of penguin colonies. In this paper, we explore the feasibility of automatic object detection algorithms based on aerial images, which locate each penguin directly. We build a dataset consisting of images taken at 400 m altitude over the island populated by Adelie penguins, which are cropped with a resolution of 640 × 640. To address the challenges of detecting minuscule penguins (often 10 pixels extent) amidst complex backgrounds in our dataset, we propose a new object detection network, named YoloPd (Yolo for penguin detection). Specifically, a multiple frequency features fusion module and a Bottleneck aggregation layer are proposed to strengthen feature representations for smaller penguins. Furthermore, the Transformer aggregation layer and efficient attention module are designed to capture global features with the aim of filtering out background interference. With respect to the latency/accuracy trade-off, YoloPd surpasses the classical detector Faster R-CNN by 8.5% in mean precision (mAP). It also beats the latest detector Yolov7 by 2.3% in F1 score with fewer parameters. Under YoloPd, the average counting accuracy reaches 94.6%, which is quite promising. The results demonstrate the potential of automatic detectors and provide a new direction for penguin counting. [ABSTRACT FROM AUTHOR]

The snow temperature above the ice sheet is one of the basic characteristic parameters of the ice sheet, which plays an important role in the study of the global climate. Because infrared and microwaves with different frequencies have different penetration depths in snow, it is possible to retrieve the snow temperature profiles by combining microwave and infrared brightness temperatures. This paper proposes a conjoint inversion algorithm to retrieve snow temperature profiles by combining multi-frequency microwave brightness temperature (BT) with infrared BT, in which different weight functions of microwave BT at different frequencies are adopted, and the atmosphere influence has also been corrected. The snow temperature profile data are retrieved based on AMSR2 microwave BT data and MODIS infrared BT data in 2017 and 2018, which are evaluated by comparing with the measured snow temperature at Dome-C station. The results confirm that the inverted snow temperature profiles are consistent with the field observation data from the Dome-C station. Multi-frequency microwave brightness temperature can be used to invert the snow temperature profiles; however, the inverted snow surface temperature is more accurate by combining the infrared BT with the microwave BT in the conjoint inversion algorithm. [ABSTRACT FROM AUTHOR]

The analytical solutions of the radiative transfer equation are needed for the solution of various applied atmospheric and snow optics problems. In this paper, we propose a simple analytical equation for the top-of-atmosphere (TOA) spectral reflectance. To simplify the problem under study we consider the case of Antarctica, where both snow and atmosphere are almost free of pollutants. This work is focused on the simulation of the moderate spectral resolution TOA measurements (1 nm or so) and the spectral range 400–1000 nm. The values of the coefficient of variance (CV) between the measured by the Ocean and Land Colour Instrument (OLCI) on board Sentinel-3A and modelled spectra are smaller than 10% for most cases in Antarctica. There are regions in Eastern Antarctica, where the values of CV are smaller than 5%. The areas with larger deviations between measured and retrieved spectra could be due to the presence of clouds or structures on the snow surface not captured by the proposed model. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Antarctic basal water storage variation (BWSV) refers to mass changes of basal water beneath the Antarctic ice sheet (AIS). Identifying these variations is critical for understanding Antarctic basal hydrology variations and basal heat conduction, yet they are rarely accessible due to a lack of direct observation. This paper proposes a layered gravity density forward/inversion iteration method to investigate Antarctic BWSV based on multi-source satellite observations and relevant models. During 2003–2009, BWSV increased at an average rate of 43 ± 23 Gt/yr, which accounts for 29% of the previously documented total mass loss rate (−76 ± 20 Gt/yr) of AIS. Major uncertainty arises from satellite gravimetry, satellite altimetry, the glacial isostatic adjustment (GIA) model, and the modelled basal melting rate. We find that increases in basal water mainly occurred in regions with widespread active subglacial lakes, such as the Rockefeller Plateau, Siple Coast, Institute Ice Stream regions, and marginal regions of East Antarctic Ice Sheet (EAIS), which indicates the increased water storage in these active subglacial lakes, despite the frequent water drainage events. The Amundsen Sea coast experienced a significant loss during the same period, which is attributed to the basal meltwater discharging into the Amundsen Sea through basal channels. [ABSTRACT FROM AUTHOR]

Recent studies have demonstrated that wideband microwave radiometers provide significant potential for profiling important subsurface polar firn characteristics necessary to understand the dynamics of the cryosphere and predict future changes in ice and snow coverage. Different frequencies within the wide spectra of radiometers result in different electromagnetic propagation losses and thus reveal characteristics at different depths in ice and snow. This paper, expanding on those investigations, explores the utilization of the Global Precipitation Measurement (GPM) constellation as a single wideband (6.93 GHz–91.655 GHz) spaceborne radiometer, covering the entire microwave spectrum from C-band to W-band, to profile subsurface properties of the Antarctic firn. Results of the initial analyses over Concordia and Vostok Stations in Antarctica indicate that GPM brightness temperature measurements provide critical information regarding the subsurface temperatures and physical properties of the firn from the surface down to several meters of depth. Considering the high spatiotemporal coverage of polar-orbiting spaceborne radiometers, these results are promising for future continent-level thermal and physical characterization of the Antarctic firn. [ABSTRACT FROM AUTHOR]

La Salle and Ebro Observatory have been involved in remote sensing projects in Antarctica for the last 11 years (approximately one solar cycle). The Ebro Observatory has been monitoring and analyzing the geomagnetic and the ionospheric activity in the Antarctic Spanish station Juan Carlos I (ASJI) (62.7°S, 299.6°E) for more than eighteen and ten years, respectively. La Salle has two main goals in the project. The first one is the data transmission and reception from Antarctica to Spain to obtain a historical series of measurements of channel sounding of this 12,760-km ionospheric HF (high frequency) radio link. The second one is the establishment of a stable data low power communication system between the ASJI and Cambrils, Spain (41.0°N, 1.0°E), to transmit the data from the remote sensors located on the island. In this paper, both narrowband and wideband soundings have been carried out to figure out the channel availability performed using a frequency range from 2 to 30 MHz with 0.5 MHz step during the 24 h of the day, encompassing wider channel measurements than previously done, in terms of hours and frequency. This paper presents the results obtained for the austral summer in 2014, using a monopole antenna at the transmitter and an inverted V on the receiver side. These results led us to the final physical layer design for the long haul link, dividing the day into two parts: daytime, with low data throughput design, and nighttime, reaching high data throughput. [ABSTRACT FROM AUTHOR]

Antarctica is covered by a thick ice sheet, and the application of geophysical methods is necessary to image the subglacial structures for studying the hydrologic systems and tectonic deformations in the Antarctic continent. The magnetotelluric (MT) method is one of the best approaches to obtain the subglacial electrical resistivities. However, only a very small volume of data has been collected so far in Antarctica using this method. In this paper, we report on a broadband MT profile collected at 10 sites in the south of the Chinese Zhongshan Station at Prydz Bay, and a 3D resistivity model was constructed by inversion of these data. This 3D model shows two low resistivity zones at a depth shallower than 6 km. They are consistent with the low-velocity zones in the previous shear-wave model and can be interpreted as the result of interconnected fluids associated with a downward migration of subglacial water. In addition, a distinct eastward dipping low resistivity zoneis present in the crust, which extends from the top surface to the lower crust. Since its location coincides with the highly positive magnetization named Amery Lineament, it is proposed that this eastward dipping low resistivity zoneextending in the entire crust probably reveals the geometry structure of the Amery Lineament at depth. Besides, it can be inferred from this new 3D resistivity model that the Amery Lineament is at least a crustal-scale structure, which probably outcrops on the land surface but was covered by the ice sheets in the study area. [ABSTRACT FROM AUTHOR]

Measurements of Very-Low-Frequency (VLF) transmitter signals have been widely used to investigate the effects of various space weather events on the D-region ionosphere, including nowcasting solar flares. Previous studies have established a method to nowcast solar flares using VLF measurements, but only using measurements from dayside propagation paths, and there remains limited focus on day–night mixed paths, which are important for method applicability. Between March and May of 2022, the Sun erupted a total of 56 M-class and 6 X-class solar flares, all of which were well captured by our VLF receiver in Antarctica. Using these VLF measurements, we reexamine the capability of the VLF technique to nowcast solar flares by including day–night mixed propagation paths and expanding the path coverage in longitude compared to that in previous studies. The amplitude and phase maximum changes are generally positively correlated with X-ray fluxes, whereas the time delay is negatively correlated. The curve-fitting parameters that we obtain for the X-ray fluxes and VLF signal maximum changes are consistent with those in previous studies for dayside paths, even though different instruments are used, supporting the flare-nowcasting method. Moreover, the present results show that, for day–night mixed paths, the amplitude and phase maximum changes also scale linearly with the logarithm of the flare X-ray fluxes, but the level of change is notably different from that for dayside paths. The coefficients used in the flare-nowcasting method need to be updated for mixed propagation paths. [ABSTRACT FROM AUTHOR]

Byers Peninsula is considered one of the largest ice-free areas in maritime Antarctica. Since 2006, the Spanish Polar Program has taken part in a large number of environmental studies involving the effects of climate change on biological life cycles, limnology, and microbiology. Soils from maritime Antarctica are generally weakly developed and have chemical, physical, and morphological characteristics that are strongly influenced by the parent material. However, biological activity during the short Antarctic summer promotes intense transference of nutrients and organic matter in areas occupied by different species of birds and marine mammals. Mapping and monitoring those areas that are highly occupied by various species could be very useful to create models prepared from satellite images of the edaphic properties. In this approach, deep learning and linear regression models of the soil properties and spectral indexes, which were considered as explicative variables, were used. We trained the models on soil properties closely related to biological activity such as dissolved organic carbon (DOC) and the iron fraction associated with the organic matter (Fe). We tested the best approach to model the spatial distribution of DOC, Fe, and pH by training the linear regression and deep learning models on Sentinel-2 and WorldView-2 images. The most robust models, the pH model built with the deep learning approach on Sentinel images (MAE of 0.51, RMSE of 0.70, and R2 with a residual of −0.49), the DOC model built with linear regression on Sentinel images (MAE of 189.39, RMSE of 342.23, and R2 with a residual of 0.0), and the organic Fe model built with deep learning (MAE of 116.20, RMSE of 209.93, and R2 of −0.05), were used to track possible areas with ornithogenic soils, as well as areas of Byers Peninsula that could be supporting the highest biological development. [ABSTRACT FROM AUTHOR]

Temperature analysis is of special interest in polar areas because temperature is an essential variable in the energy exchange between the Earth's surface and atmosphere. Although land surface temperature (LST) obtained using satellites and air temperature (Ta) have different physical meanings and are measured with different techniques, LST has often been successfully employed to estimate Ta. For this reason, in this work, we estimated Ta from LST MODIS collection 6 (C6) and used other predictor variables. Daily mean Ta was calculated from Spanish State Meteorological Agency (AEMET) stations data on the Livingston and Deception Islands, and from the PERMASNOW project stations on Livingston Island; both islands being part of the South Shetland Islands (SSI) archipelago. In relation to our previous work carried out in the study area with collection 5 (C5) data, we obtained higher R2 values (R2CV = 0.8, in the unique model with Terra daytime data) and lower errors (RMSECV = 2.2 °C, MAECV = 1.6 °C). We corroborated significant improvements in MODIS C6 LST data. We analyzed emissivity as a possible factor of discrepancies between C5 and C6, but we did not find conclusive results, therefore we could not affirm that emissivity is the factor that causes differences between one collection and another. The results obtained with the applied filters indicated that MODIS data can be used to study Ta in the area, as these filters contribute to the reduction of uncertainties in the modeling of Ta from satellites. [ABSTRACT FROM AUTHOR]

In this work, we develop a software suite for studies of atmosphere–underlying SNOW-spaceborne optical receiver light TRANsmission calculations (SNOWTRAN) with applications for the solution of forward and inverse radiative transfer problems in polar regions. Assuming that the aerosol load is extremely low, the proposed theory does not require the numerical procedures for the solution of the radiative transfer equation and is based on analytical equations for the spectral nadir reflectance and simple approximations for the local optical properties of atmosphere and snow. The developed model is validated using EnMAP and PRISMA spaceborne imaging spectroscopy data close to the Concordia research station in Antarctica. A new, fast technique for the determination of the snow grain size and assessment of the snowpack vertical inhomogeneity is then proposed and further demonstrated on EnMAP imagery over the Aviator Glacier and in the vicinity of the Concordia research station in Antarctica. The results revealed a large increase in precipitable water vapor at the Concordia research station in February 2023 that was linked to a warming event and a four times larger grain size at Aviator Glacier compared with Dome C. [ABSTRACT FROM AUTHOR]

Polynyas are an important factor in the Antarctic and Arctic climate, and their changes are related to the ecosystems in the polar regions. The phenomenon of polynyas is influenced by the combination of inherent persistence and dynamic factors. The dynamics of polynyas are greatly affected by temporal dynamical factors, and it is difficult to objectively reflect the internal characteristics of their formation. Separating the two factors effectively is necessary in order to explore their essence. The Special Sensor Microwave/Imager (SSM/I) passive microwave sensor has been making observations of Antarctica for more than 20 years, but it is difficult for existing current sea ice concentration (SIC) products to objectively reflect how the inherent persistence factors affect the formation of polynyas. In this paper, we proposed a long-term multiple spatial smoothing method to remove the influence of dynamic factors and obtain stable annual SIC products. A halo located on the border of areas of low and high ice concentration around the Antarctic coast, which has a strong similarity with the local seabed in outline, was found using the spatially smoothed SIC products and seabed. The relationship of the polynya location to the wind and topography is a long-understood relationship; here, we quantify that where there is an abrupt slope and wind transitions, new polynyas are best generated. A combination of image expansion and threshold segmentation was used to extract the extent of sea ice and coastal polynyas. The adjusted record of changes in the extent of coastal polynyas and sea ice in the Southern Ocean indicate that there is a negative correlation between them. [ABSTRACT FROM AUTHOR]

We conducted a comparison of sea ice concentration (SIC) in the Amundsen Sea using three satellite datasets: Hadley Centre's sea ice and sea surface temperature (HadISST1), Operational Sea Surface Temperature and Ice Analysis (OSTIA), and Advanced Microwave Scanning Radiometer 2 (AMSR2). HadISST1 has the longest time period, while AMSR2 has the shortest. In terms of grid resolution, HadISST1 has the coarsest resolution, while AMSR2 has the finest. The sea ice areas (SIAs) observed in HadISST1, OSTIA, and AMSR2 are similar. We studied the decadal variations in SICs by dividing the study period into four temporal segments. We investigated the differences between HadISST1 and OSTIA for each temporal segment. HadISST1 exhibited a more pronounced positive trend compared to OSTIA between 2005 and 2010. Additionally, we compared the interannual and seasonal variations in SICs between HadISST1 and OSTIA. Lastly, it should be noted that the Amundsen Sea polynya area varies across all three datasets. [ABSTRACT FROM AUTHOR]

Mapping Antarctic Specially Protected Areas (ASPAs) remains a critical yet challenging task, especially in extreme environments like Antarctica. Traditional methods are often cumbersome, expensive, and risky, with limited satellite data further hindering accuracy. This study addresses these challenges by developing a workflow that enables precise mapping and monitoring of vegetation in ASPAs. The processing pipeline of this workflow integrates small unmanned aerial vehicles (UAVs)—or drones—to collect hyperspectral and multispectral imagery (HSI and MSI), global navigation satellite system (GNSS) enhanced with real-time kinematics (RTK) to collect ground control points (GCPs), and supervised machine learning classifiers. This workflow was validated in the field by acquiring ground and aerial data at ASPA 135, Windmill Islands, East Antarctica. The data preparation phase involves a data fusion technique to integrate HSI and MSI data, achieving the collection of georeferenced HSI scans with a resolution of up to 0.3 cm/pixel. From these high-resolution HSI scans, a series of novel spectral indices were proposed to enhance the classification accuracy of the model. Model training was achieved using extreme gradient boosting (XGBoost), with four different combinations tested to identify the best fit for the data. The research results indicate the successful detection and mapping of moss and lichens, with an average accuracy of 95%. Optimised XGBoost models, particularly Model 3 and Model 4, demonstrate the applicability of the custom spectral indices to achieve high accuracy with reduced computing power requirements. The integration of these technologies results in significantly more accurate mapping compared to conventional methods. This workflow serves as a foundational step towards more extensive remote sensing applications in Antarctic and ASPA vegetation mapping, as well as in monitoring the impact of climate change on the Antarctic ecosystem. [ABSTRACT FROM AUTHOR]

The characteristics of sea ice leads (SILs) in the Weddell Sea are an important basis for understanding the mechanism of the atmosphere–ocean system in the Southern Ocean. In this study, we derived the sea ice surface temperature (IST) of the Weddell Sea from MODIS thermal images and then generated a daily SIL map for 2015 and 2022 by utilizing the iterative threshold method on the optimised MOD35 cloud-masked IST. The results showed that SIL variations in the Weddell Sea presented remarkable seasonal characteristics. The trend of the SIL area exhibited an initial rise followed by a decline from January to December, characterised by lower values in spring and summer and higher values in fall and winter. SILs in the Weddell Sea were predominantly concentrated between 70~78°S and 60~30°W. The coastal spatial distribution density of the SILs exceeded that of offshore regions, peaking near the Antarctic Peninsula and then near Queen Maud Land. The SIL variation was mainly influenced by dynamical factors, and there were strong positive correlations between the wind field, ocean currents, and sea-ice motion. [ABSTRACT FROM AUTHOR]

Polar amplification appears in response to greenhouse gas forcing, which has become a focus of climate change research. However, polar amplification has not been systematically investigated over the Earth's three poles (the Arctic, Antarctica, and the Third Pole). An index of polar amplification is employed, and the annual and seasonal variations of land surface temperature over the Earth's three poles are examined using MODIS (Moderate Resolution Imaging Spectroradiometer) observations for the period 2001–2018. As expected, the warming of the Arctic is most conspicuous, followed by the Third Pole, and is weakest in Antarctica. Compared to the temperature changes for the global land region, positive polar amplification appears in the Arctic and the Third Pole on an annual scale, whereas Antarctic amplification disappears, with a negative amplification index of −0.72. The polar amplification for the Earth's three poles shows seasonal differences. Strong Arctic amplification appears in boreal spring and winter, with a surface warming rate of more than 3.40 times the global mean for land regions. In contrast, the amplification of the Third Pole is most conspicuous in boreal summer. The two poles located in the Northern Hemisphere have the weakest amplification in boreal autumn. Differently from the positive amplification for the Arctic and the Third Pole in all seasons, the faster variations in Antarctic temperature compared to the globe only appear in austral autumn and winter, and the amplification signal is negative in these seasons, with an amplification index of −1.68 and −2.73, respectively. In the austral winter, the strong negative amplification concentrates on West Antarctica and the coast of East Antarctica, with an absolute value of amplification index higher than 5 in general. Generally, the polar amplification is strongest in the Arctic except from June to August, and Antarctic amplification is the weakest among the Earth's three poles. The Earth's three poles are experiencing drastic changes, and the potential influence of climate change should receive attention. [ABSTRACT FROM AUTHOR]

Ice motion is an essential element for accurately evaluating glacier mass balance. Interferometric synthetic aperture radar (InSAR) has been widely applied for monitoring ice motion with high precision and wide coverage in the Antarctic. However, the ionospheric effects can significantly impact InSAR-based ice-motion measurements. At low radar frequencies in particular, the ionospheric effects have been regarded as a serious source of noise in L-band SAR data. The split-spectrum method (SSM) is commonly used for correcting the ionospheric effects of the InSAR technique. However, it requires spatial filtering with the relatively large factors used to scale the sub-bands' interferograms, which often results in an unwrapped phase error. In this paper, a reformulation of the split-spectrum method (RSSM) is introduced to correct the ionospheric effects in the Grove Mountains of East Antarctica, which have slow ice flow and frequent ionosphere changes. The results show that RSSM can effectively correct the ionospheric effects of InSAR-based ice-motion measurements. To evaluate the ability of ionospheric correction using RSSM, the result of ionospheric correction derived from SSM is compared with the results of RSSM. In addition, ionosphere-corrected ice motion is also compared with GPS and MEaSUREs. The results show that the ionosphere-corrected ice velocities are in good agreement with GPS observations and MEaSUREs. The average ice velocity from the InSAR time series is compared to that from MEaSUREs, and the average ionosphere-corrected ice velocity error reduces 43.9% in SSM and 51.1% in RSSM, respectively. The ionosphere-corrected ice velocity error is the most significant, reducing 86.9% in SSM and 90.4% in RSSM from 1 November 2007 to 19 December 2007. The results show that the ability of RSSM to correct ionospheric effects is slightly better than that of SSM. Therefore, we deduce that the RSSM offers a feasible way to correct ionospheric effects in InSAR-based ice-motion measurements in Antarctica. [ABSTRACT FROM AUTHOR]

Due to rapid global warming, the relationship between the mass loss of the Antarctic ice sheet and rising sea levels are attracting widespread attention. The Lambert–Amery glacial system is the largest drainage system in East Antarctica, and its mass balance has an important influence on the stability of the Antarctic ice sheet. In this paper, the recent ice flux in the Lambert Glacier of the Lambert–Amery system was systematically analyzed based on recently updated remote sensing data. According to Landsat-8 ice velocity data from 2018 to April 2019 and the updated Bedmachine v2 ice thickness dataset in 2021, the contribution of ice flux approximately 140 km downstream from Dome A in the Lambert Glacier area to downstream from the glacier is 8.5 ± 1.9   Gt · a − 1 , and the ice flux in the middle of the convergence region is 18.9 ± 2.9   Gt · a − 1 . The ice mass input into the Amery ice shelf through the grounding line of the whole glacier is 19.9 ± 1.3   Gt · a − 1 . The ice flux output from the mainstream area of the grounding line is 19.3 ± 1.0   Gt · a − 1 . Using the annual SMB data of the regional atmospheric climate model (RACMO v2.3) as the quality input, the mass balance of the upper, middle, and lower reaches of the Lambert Glacier was analyzed. The results show that recent positive accumulation appears in the middle region of the glacier (about 74–78°S, 67–85°E) and the net accumulation of the whole glacier is 2.4 ± 3.5   Gt · a − 1 . Although the mass balance of the Lambert Glacier continues to show a positive accumulation, and the positive value in the region is decreasing compared with values obtained in early 2000. [ABSTRACT FROM AUTHOR]

Convolutional neural networks (CNNs) are becoming an increasingly popular approach for classification mapping of large complex regions where manual data collection is too time consuming. Stream boundaries in hyper-arid polar regions such as the McMurdo Dry Valleys (MDVs) in Antarctica are difficult to locate because they have little hydraulic flow throughout the short summer months. This paper utilizes a U-Net CNN to map stream boundaries from lidar derived rasters in Taylor Valley located within the MDVs, covering ∼770 km2. The training dataset consists of 217 (300 × 300 m2) well-distributed tiles of manually classified stream boundaries with diverse geometries (straight, sinuous, meandering, and braided) throughout the valley. The U-Net CNN is trained on elevation, slope, lidar intensity returns, and flow accumulation rasters. These features were used for detection of stream boundaries by providing potential topographic cues such as inflection points at stream boundaries and reflective properties of streams such as linear patterns of wetted soil, water, or ice. Various combinations of these features were analyzed based on performance. The test set performance revealed that elevation and slope had the highest performance of the feature combinations. The test set performance analysis revealed that the CNN model trained with elevation independently received a precision, recall, and F1 score of 0.94 ± 0.05 , 0.95 ± 0.04 , and 0.94 ± 0.04 respectively, while slope received 0.96 ± 0.03 , 0.93 ± 0.04 , and 0.94 ± 0.04 , respectively. The performance of the test set revealed higher stream boundary prediction accuracies along the coast, while inland performance varied. Meandering streams had the highest stream boundary prediction performance on the test set compared to the other stream geometries tested here because meandering streams are further evolved and have more distinguishable breaks in slope, indicating stream boundaries. These methods provide a novel approach for mapping stream boundaries semi-automatically in complex regions such as hyper-arid environments over larger scales than is possible for current methods. [ABSTRACT FROM AUTHOR]

The ice shelf is an important component of the Antarctic system, and the interaction between the ice sheet and the ocean often proceeds through mass variations of the ice shelf. The digital elevation model (DEM) of the ice shelf is particularly important for ice shelf elevation change and mass balance estimation. With the development of satellite altimetry technology, it became an important data source for DEM research of Antarctica. The National Aeronautics and Space Administration (NASA) Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) launched in 2018 is a significant improvement in along-track sampling rate and measurement accuracy compared with previous altimetry satellites. This study uses ordinary kriging interpolation to present new DEMs (ICESat-2 DEM hereinafter) for the three ice shelves (Ross, Filchner–Ronne and Amery) in Antarctica with ICESat-2 altimetry data. Two variogram models (linear and spherical) of ordinary kriging interpolation are compared in this paper. The result shows that the spherical model generally shows better performance and lower standard deviation (STD) than the linear models. The precision of the ultimate DEM was evaluated by NASA Operation IceBridge (OIB) data and compared with five previously published Antarctic DEM products (REMA, TanDEM-X PolarDEM, Slater DEM, Helm DEM, and Bamber DEM). The comparison reveals that the mean difference between ICESat-2 DEM of the Ross ice shelf and OIB is −0.016 m with a STD of 0.918 m, and the mean difference between ICESat-2 DEM of the Filchner–Ronne ice shelf and OIB is −0.533 m with a STD of 0.718 m. The three ICESat-2 DEMs show higher spatial resolution and elevation accuracy than five previously published Antarctic DEMs. [ABSTRACT FROM AUTHOR]

With global warming accelerating, polar amplification is one of the hot issues in climate research. However, most studies focus on Arctic amplification, and little attention has been paid to Antarctic amplification (AnA), and there is no relevant research based on MODIS (Moderate Resolution Imaging Spectroradiometer) land surface temperature observations. Compared with 128 stations' observations, MODIS can capture the variations in temperature over Antarctica. In addition, the temperature changes in Antarctica, East Antarctica, West Antarctica and the Antarctic Peninsula during the period 2001–2018 reflected by the MODIS and ERA5 are basically consistent, and the temperature changes in Antarctica are negatively correlated with the Southern Annular Mode. AnA occurs under all annual and seasonal scales, with an AnA index greater than 1.27 (1.31) from the MODIS (ERA5), and is strongest in the austral winter and weakest in summer. AnA displays regional differences, and the signal from the MODIS is similar to that from ERA5. The strongest amplification occurs in East Antarctica, with an AnA index greater than 1.45 (1.48) from the MODIS (ERA5), followed by West Antarctica, whereas the amplified signal is absent at the Antarctic Peninsula. In addition, seasonal differences can be observed in the sub regions of Antarctica. For West Antarctica, the greatest amplification appears in austral winter, and in austral spring for East Antarctica. The AnA signal also can be captured in daytime and nighttime observations, and the AnA in nighttime observations is stronger than that in daytime. Generally, the MODIS illustrates the appearance of AnA for the period 2001–2018, and the Antarctic climate undergoes drastic changes, and the potential impact should arouse attention. [ABSTRACT FROM AUTHOR]

Although extensive research of Moderate Resolution Imaging Spectroradiometer (MODIS) albedo data is available on the Greenland Ice Sheet, there is a lack of studies evaluating MODIS albedo products over Antarctica. In this paper, MOD10A1, MYD10A1, and MCD43 (C6) daily albedo products were compared with the in situ albedo data on Livingston Island, South Shetland Islands (SSI), Antarctica, from 2006 to 2015, for both all-sky and clear-sky conditions, and for the entire study period and only the southern summer months. This is the first evaluation in which MYD10A1 and MCD43 are also included, which can be used to improve the accuracy of the snow BRDF/albedo modeling. The best correlation was obtained with MOD10A1 in clear-sky conditions (r = 0.7 and RMSE = 0.042). With MCD43, only data from the backup algorithm could be used, so the correlations obtained were lower (r = 0.6). However, it was found that there was no significant difference between the values obtained for all-sky and for clear-sky data. In addition, the MODIS products were found to describe the in situ data trend, with increasing albedo values in the range between 0.04 decade−1 and 0.16 decade−1. We conclude that MODIS daily albedo products can be applied to study the albedo in the study area. [ABSTRACT FROM AUTHOR]

Over the last two decades, spaceborne polarimetric synthetic aperture radar (PolSAR) has been widely used to penetrate sea ice surfaces to achieve fully polarimetric high-resolution imaging at all times of day and in a range of weather conditions. Model-based polarimetric decomposition is a powerful tool used to extract useful physical and geometric information about sea ice from the matrix datasets acquired by PolSAR. The volume scattering of sea ice is usually modeled as the incoherent average of scatterings of a large volume of oriented ellipsoid particles that are uniformly distributed in 3D space. This uniform spatial distribution is often approximated as a uniform orientation distribution (UOD), i.e., the particles are uniformly oriented in all directions. This is achieved in the existing literature by ensuring the canting angle φ and tilt angle τ of particles uniformly distributed in their respective ranges and introducing a factor cos ⁡ τ in the ensemble average. However, we find this implementation of UOD is not always effective, while a real UOD can be realized by distributing the solid angles of particles uniformly in 3D space. By deriving the total solid angle of the canting-tilt cell spanned by particles and combining the differential relationship between solid angle and Euler angles φ and τ , a complete expression of the joint probability density function p φ , τ that can always ensure the uniform orientation of particles of sea ice is realized. By ensemble integrating the coherency matrix of φ , τ -oriented particle with p φ , τ , a generalized modeling of the volume coherency matrix of 3D uniformly oriented spheroid particles is obtained, which covers factors such as radar observation geometry, particle shape, canting geometry, tilt geometry and transmission effect in a multiplicative way. The existing volume scattering models of sea ice constitute special cases. The performance of the model in the characterization of the volume behaviors was investigated via simulations on a volume of oblate and prolate particles with the differential reflectivity Z D R , polarimetric entropy H and scattering α angle as descriptors. Based on the model, several interesting orientation geometries were also studied, including the aligned orientation, complement tilt geometry and reflection symmetry, among which the complement tilt geometry is specifically highlighted. It involves three volume models that correspond to the horizontal tilt, vertical tilt and random tilt of particles within sea ice, respectively. To match the models to PolSAR data for adaptive decomposition, two selection strategies are provided. One is based on Z D R , and the other is based on the maximum power fitting. The scattering power that reduces the rank of coherency matrix by exactly one without violating the physical realizability condition is obtained to make full use of the polarimetric scattering information. Both the models and decomposition were finally validated on the Gaofen-3 PolSAR data of a young ice area in Prydz Bay, Antarctica. The adaptive decomposition result demonstrates not only the dominant vertical tilt preference of brine inclusions within sea ice, but also the subordinate random tilt preference and non-negligible horizontal tilt preference, which are consistent with the geometric selection mechanism that the c-axes of polycrystallines within sea ice would gradually align with depth. The experiment also indicates that, compared to the strategy based on Z D R , the maximum power fitting is preferable because it is entirely driven by the model and data and is independent of any empirical thresholds. Such soft thresholding enables this strategy to adaptively estimate the negative Z D R offset introduced by the transmission effect, which provides a novel inversion of the refractive index of sea ice based on polarimetric model-based decomposition. [ABSTRACT FROM AUTHOR]

This study focuses on the development of a time-variable regional geo-potential model for Antarctica using the spherical cap harmonic analysis (SCHA) basis functions. The model is derived from line-of-sight gravity difference (LGD) measurements obtained from the GRACE-Follow-On (GFO) mission. The solution of a Laplace equation for the boundary values over a spherical cap is used to expand the geo-potential coefficients in terms of Legendre functions with a real degree and integer order suitable for regional modelling, which is used to constrain the geo-potential coefficients using LGD measurements. To validate the performance of the SCHA, it is first utilized with LGD data derived from a L2 JPL (Level 2 product of the Jet Propulsion Laboratory). The obtained LGD data are used to compute the local geo-potential model up to Kmax = 20, corresponding to the SH degree and order up to 60. The comparison of the radial gravity on the Earth's surface map across Antarctica with the corresponding radial gravity components of the L2 JPL is carried out using local geo-potential coefficients. The results of this comparison provide evidence that these basis functions for Kmax = 20 are valid across the entirety of Antarctica. Subsequently, the analysis proceeds using LGD data obtained from the Level 1B product of GFO by transforming these LGD data into the SCHA coordinate system and applying them to constrain the SCHA harmonic coefficients up to Kmax = 20. In this case, several independent LGD profiles along the trajectories of the satellites are devised to verify the accuracy of the local model. These LGD profiles are not employed in the inverse problem of determining harmonic coefficients. The results indicate that using regional harmonic basis functions, specifically spherical cap harmonic analysis (SCHA) functions, leads to a close estimation of LGD compared to the L2 JPL. The regional harmonic basis function exhibits a root mean square error (RMSE) of 3.71 × 10−4 mGal. This represents a substantial improvement over the RMSE of the L2 JPL, which is 6.36 × 10−4 mGal. Thus, it can be concluded that the use of local geo-potential coefficients obtained from SCHA is a reliable method for extracting nearly the full gravitational signal within a spherical cap region, after validation of this method. The SCHA model provides significant realistic information as it addresses the mass gain and loss across various regions in Antarctica. [ABSTRACT FROM AUTHOR]

Geological models are essential components in various applications. To generate reliable realizations, the geostatistical method focuses on reproducing spatial structures from training images (TIs). Moreover, uncertainty plays an important role in Earth systems. It is beneficial for creating an ensemble of stochastic realizations with high diversity. In this work, we applied a pattern classification distribution (PCD) method to quantitatively evaluate geostatistical modeling. First, we proposed a correlation-driven template method to capture geological patterns. According to the spatial dependency of the TI, region growing and elbow-point detection were launched to create an adaptive template. Second, a combination of clustering and classification was suggested to characterize geological realizations. Aiming at simplifying parameter specification, the program employed hierarchical clustering and decision tree to categorize geological structures. Third, we designed a stacking framework to develop the multi-grid analysis. The contribution of each grid was calculated based on the morphological characteristics of TI. Our program was extensively examined by a channel model, a 2D nonstationary flume system, 2D subglacial bed topographic models in Antarctica, and 3D sandstone models. We activated various geostatistical programs to produce realizations. The experimental results indicated that PCD is capable of addressing multiple geological categories, continuous variables, and high-dimensional structures. [ABSTRACT FROM AUTHOR]

A polynya is an area of open water or reduced concentration of sea ice surrounded by either concentrated sea ice or land ice. They are often seen as sites of intense ocean–atmosphere heat exchange and as ice production factories. Given their importance, it is crucial to quantify the accuracy of satellite-derived polynya information. Polynyas in their early evolution phase are generally narrow and occur at scales likely too fine to be detected by widely used passive microwave (PMW) radiometric sensors. We derived 40 m scale polynya information over the western Ross Sea from high-resolution Synthetic Aperture Radar (SAR) Sentinel-1 C-band data and examined discrepancies with larger-scale estimates. We utilized two automated algorithms, supervised (a rule-based approach) and unsupervised (a combination of texture analysis with k-means clustering), to accurately identify the polynya areas. We generated data for validation using Sentinel-1 data at instances where polynyas can be visually delineated. Results from PMW sensors (NSIDC and AMSR2) and SAR-based algorithms (rule-based and texture-based) are compared with manually delineated polynya areas obtained through Sentinel-1. Analysis using PMW sensors revealed that NSIDC overestimates larger polynyas and underestimates smaller polynyas compared to AMSR2. We were more accurately able to identify polynya presence and area using Sentinel-1 SAR observations, especially in clear cases and cases when PMW data miscalculates the polynya's presence. Of our SAR-based algorithms, the rule-based approach was more accurate than the texture-based approach at identifying clear polynyas when validated against manually delineated regions. Altogether, we emphasize the need for finer spatio-temporal resolution data for polynya studies. [ABSTRACT FROM AUTHOR]

The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2), launched in September 2018, has been widely used in forestry and surveying. A high-accuracy digital elevation model (DEM)/digital surface model (DSM) for terrain matching can effectively evaluate the ICESat-2 design requirements and provide essential data support for further study. The conventional terrain-matching methods regard the laser ground track as a whole, ignoring the individual differences caused by the interaction of photons during flight. Therefore, a novel terrain-matching method using a two-dimensional affine transformation model was proposed to describe the deformation of laser tracks. The least-square optimizes the model parameters with the high-accuracy terrain data to obtain the best matching result. The results in McMurdo Dry Valley (MDV), Antarctica, and Zhengzhou (ZZ), China, demonstrate that the proposed method can verify geolocation accuracy and indicate that the average horizontal accuracy of ICESat-2 V5 data is about 3.86 m in MDV and 4.67 m in ZZ. It shows that ICESat-2 has good positioning accuracy, even in mountainous areas with complex terrain. Additionally, the random forest (RF) model was calculated to analyze the influence of four factors on geographic location accuracy. The slope and signal-to-noise ratio (SNR) are considered the crucial factors affecting the accuracy of ICESat-2 data. [ABSTRACT FROM AUTHOR]

David Glacier and Drygalski Ice Tongue are massive glaciers in Victoria Land, Antarctica. The ice from the East Antarctic Ice Sheet is drained through the former, and then discharged into the western Ross Sea through the latter. David Drygalski is the largest outlet glacier in Northern Victoria Land, floating kilometers out to sea. The floating and grounded part of the David Glacier are the main focus of this article. During the XXI Italian Antarctic Expedition (2005–2006), within the framework of the National Antarctic Research Programme (PNRA), two GNSS stations were installed at different points: the first close to the grounding line of David Glacier, and the second approximately 40 km downstream of the first one. Simultaneous data logging was performed by both GNSS stations for 24 days. In the latest data processing, the kinematic PPP technique was adopted to evaluate the dominant diurnal components and the very small semi-diurnal variations in ice motion induced by the ocean tide and the mean ice flow rates of both GNSS stations. Comparison of the GNSS time series with predicted ocean tide calculated from harmonic coefficients of the nearest tide gauge stations, installed at Cape Roberts and Mario Zucchelli Station, highlight different local response of the glacier to ocean tide, with a minor amplitude of vertical motion at a point partially anchored at the bedrock close to the grounding line. During low tide, the velocity of the ice flow reaches its daily maximum, in accordance with the direction of seawater outflow from the fjord into the ocean, while the greatest daily tidal excursion generates an increase in the horizontal ice flow velocity. With the aim to extend the analysis in spatial terms, five COSMO-SkyMED Stripmap scenes were processed. The comparison of the co-registered offset tracking rates, obtained from SAR images, with the GNSS estimation shows good agreement. [ABSTRACT FROM AUTHOR]

(1) The Environmental Trace Gases Monitoring Instrument-2(EMI-2) is a high-quality spaceborne imaging spectrometer that launched in September 2021. To evaluate its radiometric calibration performance in-flight, the UV2 and VIS1 bands of EMI-2 were cross-calibrated by the corresponding bands (band3 and band4) of TROPOMI over the pseudo-invariant calibration site Dome C. (2) After angle limitation and cloud filtering of the Earth radiance data measured by EMI-2 and TROPOMI over Dome C, the top of atmosphere (TOA) reflectance time series were calculated. The spectral adjustment factors (SAF) were derived from the solar spectrum measured by the sensor to minimize the uncertainties caused by the different spectral response functions (SRF) of sensors. In addition, a correction method based on the radiative transfer model (RTM) SCIATRAN was used to suppress unaccounted angular dependence of atmospheric scattering. The radiation performance of EMI-2 is evaluated using the TOA reflectance ratio of EMI-2 and TROPOMI, combining the SAF correction and RTM-based correction methods. (3) It was shown that the time series trending of the TOA reflectance ratio between EMI-2 measurements and TROPOMI demonstrate flat characteristics and strong correlation. The mean reflectance ratios range from 0.998 to 1.09. The standard deviation of the reflection ratio is less than 3%. For 328 nm, 335 nm, 340 nm, 460 nm, and 490 nm, the mean values are close to one, and the relative radiometric bias estimated through EMI-2 and TROPOMI intercalibration is less than 3%, and for other wavelengths, the biases are less than 6%, except for 416 nm, which behaves higher than 7%. The cross-calibration results show that the radiometric calibration of EMI-2 is within the relative accuracy requirement. [ABSTRACT FROM AUTHOR]

Remote sensing, and unmanned aerial vehicles (UAVs) in particular, can be a valid tool for assessing the dynamics of cryotic features as frost blisters and to monitor the surface changes and the sublimation rates on perennially frozen lakes that host important ecosystems. In this paper, through the use of these remote sensing techniques, we aim to understand the type of groundwater supply of an Antarctic perennial frozen lake that encompasses two frost blisters (M1 and M2) through the temporal analysis of the features' elevation changes (frost blisters and lake ice level). The frozen lake is located at Boulder Clay (northern Victoria Land, Antarctica). We relied on several photogrammetric models, past satellite images and ground pictures to conduct differencing of digital elevation models, areal variations and pixel counting. In addition, in situ measurements of the ice sublimation or snow accumulation were carried out. The two frost blisters showed different elevation trends with M1 higher in the past (1996–2004) than recently (2014–2019), while M2 showed an opposite trend, similarly to the ice level. Indeed, the linear regression between M2 elevation changes and the ice level variation was statistically significant, as well as with the annual thawing degree days, while M1 did not show significant results. From these results we can infer that the groundwater supply of M1 can be related to a sublake open talik (hydraulic system) as confirmed also by pressurized brines found below M1, during a drilling in summer 2019. For M2 the groundwater flow is still not completely clear although the hydrostatic system seems the easiest explanation as well as for the uplift of the lake ice. [ABSTRACT FROM AUTHOR]

Supraglacial meltwater accumulation on ice sheets can be a main driver for accelerated ice discharge, mass loss, and global sea-level-rise. With further increasing surface air temperatures, meltwater-induced hydrofracturing, basal sliding, or surface thinning will cumulate and most likely trigger unprecedented ice mass loss on the Greenland and Antarctic ice sheets. While the Greenland surface hydrological network as well as its impacts on ice dynamics and mass balance has been studied in much detail, Antarctic supraglacial lakes remain understudied with a circum-Antarctic record of their spatio-temporal development entirely lacking. This study provides the first automated supraglacial lake extent mapping method using Sentinel-1 synthetic aperture radar (SAR) imagery over Antarctica and complements the developed optical Sentinel-2 supraglacial lake detection algorithm presented in our companion paper. In detail, we propose the use of a modified U-Net for semantic segmentation of supraglacial lakes in single-polarized Sentinel-1 imagery. The convolutional neural network (CNN) is implemented with residual connections for optimized performance as well as an Atrous Spatial Pyramid Pooling (ASPP) module for multiscale feature extraction. The algorithm is trained on 21,200 Sentinel-1 image patches and evaluated in ten spatially or temporally independent test acquisitions. In addition, George VI Ice Shelf is analyzed for intra-annual lake dynamics throughout austral summer 2019/2020 and a decision-level fused Sentinel-1 and Sentinel-2 maximum lake extent mapping product is presented for January 2020 revealing a more complete supraglacial lake coverage (~770 km2) than the individual single-sensor products. Classification results confirm the reliability of the proposed workflow with an average Kappa coefficient of 0.925 and a F 1 -score of 93.0% for the supraglacial water class across all test regions. Furthermore, the algorithm is applied in an additional test region covering supraglacial lakes on the Greenland ice sheet which further highlights the potential for spatio-temporal transferability. Future work involves the integration of more training data as well as intra-annual analyses of supraglacial lake occurrence across the whole continent and with focus on supraglacial lake development throughout a summer melt season and into Antarctic winter. [ABSTRACT FROM AUTHOR]

In Antarctica, the severe climatic conditions and the thick ice sheet that covers the largest and most internal part of the continent make it particularly difficult to systematically carry out geophysical and geodetic observations on a continental scale. It prevents the comprehensive understanding of both the onshore and offshore geology as well as the relationship between the inner part of East Antarctica (EA) and the coastal sector of Victoria Land (VL). With the aim to reduce this gap, in this paper multiple geophysical dataset collected since the 1980s in Antarctica by Programma Nazionale di Ricerche in Antartide (PNRA) were integrated with geodetic observations. In particular, the analyzed data includes: (i) Geodetic time series from Trans Antarctic Mountains DEFormation (TAMDEF), and Victoria Land Network for DEFormation control (VLNDEF) GNSS stations installed in Victoria Land; (ii) the integration of on-shore (ground points data and airborne) gravity measurements in Victoria Land and marine gravity surveys performed in the Ross Sea and the narrow strip of Southern Ocean facing the coasts of northern Victoria Land. Gravity data modelling has improved the knowledge of the Moho depth of VL and surrounding the offshore areas. By the integration of geodetic and gravitational (or gravity) potential results it was possible to better constrain/identify four geodynamic blocks characterized by homogeneous geophysical signature: the Southern Ocean to the N, the Ross Sea to the E, the Wilkes Basin to the W, and VL in between. The last block is characterized by a small but significant clockwise rotation relative to East Antarctica. The presence of a N-S to NNW-SSE 1-km step in the Moho in correspondence of the Rennick Geodynamic Belt confirms the existence of this crustal scale discontinuity, possibly representing the tectonic boundary between East Antarctica and the northern part of VL block, as previously proposed by some geological studies. [ABSTRACT FROM AUTHOR]

This paper presents an ice classification algorithm based on combined active and passive microwave radiometer data in Lützow-Holm Bay (LHB), East Antarctica. The ice classification algorithm is developed based on the threshold values of an advanced scatterometer (ASCAT) and advanced microwave scanning radiometer 2 (here, AMSR2). These values are calculated via the features of various ice types, including open ice, first-year (FY) ice, multi-year (MY) ice, MY ice including icebergs (MY IB), ice shelves, coastal ice sheets, and inland ice sheets. To verify the validity of the ice classification algorithm, the algorithm results are compared with visual observation data and satellite imagery. Except for the flaw polynya and area with surface melting, the FY ice, MY ice, and the ice shelf areas estimated here using the proposed ice classification algorithm match those discernible from the verification data. Inter-annual changes in the areal extents of FY ice, MY ice, and the ice shelves are investigated here using the proposed ice classification algorithm. Investigation of MY ice and ice shelf areas revealed that the breakup of MY ice induced a breakup of an ice shelf. A comparison of the FY ice and MY ice areas showed the replacement of these ice types. The proposed ice classification algorithm can detect ice breakup events as quantitative changes in the distribution and ice type. In future work, we plan to classify sea ice in other sea ice areas, applying the proposed algorithm throughout the Antarctic region. [ABSTRACT FROM AUTHOR]

Mapping spatiotemporal changes in the distribution of blue ice regions (BIRs) in Antarctica requires repeated, precise, and high-resolution baseline maps of the blue ice extent. This study demonstrated the design and application of a newly-developed semi-automatic method to map BIRs in the Antarctic environment using very high-resolution (VHR) WorldView-2 (WV-2) satellite images. We discussed the potential of VHR satellite data for the mapping of BIRs in the Antarctic environment using a customized normalized-difference blue-ice index (NDBI) method devised using yellow, green, and near-infrared spectral bands of WV-2 data. We compared the viability of the newly developed customized NDBI approach against state-of-the-art target detection (TD), spectral processing (SP) and pixel-wise supervised (PSC) feature extraction (FE) approaches. Four semi-automatic FE approaches (three existing plus one newly developed) consisting of 16 standalone FE methods (12 existing + four customized) were evaluated using an extensive quantitative and comparative assessment for mapping BIRs in the vicinity of Schirmacher Oasis, on the continental Antarctic coastline. The results suggested that the customized NDBI approach gave a superior performance and the highest statistical stability when compared with existing FE techniques. The customized NDBI generally rendered the lowest level of misclassification (average RMSE = 654.48 ± 58.26 m2), followed by TD (average RMSE = 987.81 ± 55.05 m2), SP (average RMSE = 1327.09 ± 127.83 m2) and PSC (average RMSE = 2259.43 ± 115.36 m2) for mapping BIRs. Our results indicated that the use of the customized NDBI approach can greatly improve the semi-automatic mapping of BIRs in the Antarctic environment. This study presents the first refined map of distribution of BIRs around the Schirmacher Oasis. The total area of blue ice in the study area was estimated to be 106.875 km2, approximately 61% of the study area. The WV-2 derived BIR map area presented in this study locally refined the existing BIR map derived using Landsat Enhanced Thematic Mapper Plus (ETM+) and the Moderate Resolution Imaging Spectroradiometer (MODIS)-based mosaic of Antarctica (MOA) dataset by ~31% (~33.40 km2). Finally, we discussed the practical challenges and future directions in mapping BIRs across Antarctica. [ABSTRACT FROM AUTHOR]

The Southern Patagonian Icefield (SPI) is the largest continuous ice mass in the Southern Hemisphere outside Antarctica. It has been shrinking since the Little Ice Age (LIA) period, with increasing rates in recent years. An uplift of crustal deformation in response to this deglaciation process has been expected. The goal of this investigation is to analyze the crustal deformation caused by ice retreat using time-series data from continuous GPS stations (2015–2020) in the northern area of the SPI. For this purpose, we installed two continuous GPS stations on rocky nunataks of the SPI (the GRCS near Greve glacier and the GBCS close by Cerro Gorra Blanca). In addition, ice elevation changes (2000–2019) were analyzed by the co-registration of the SRTM digital elevation model and ICESat elevation data points. The results of the vertical components are positive (36.55 ± 2.58 mm a−1), with a maximum at GBCS, indicating the highest rate of crustal uplift ever continuously recorded in Patagonia; in addition, the mean horizontal velocities reached 11.7 mm a−1 with an azimuth of 43°. The negative ice elevation changes detected in the region have also accelerated in the recent two decades, with a median Δ h (elevation change) of −3.36 ± 0.01 m a−1 in the ablation zone. The seasonality of the GPS signals was contrasted with the water levels of the main Patagonian lakes around the SPI, detecting a complex interplay between them. Hence, the study sheds light on the knowledge of the crustal uplift as evidence of the wastage experienced by the SPI glaciers. [ABSTRACT FROM AUTHOR]

Subglacial water bodies are critical components in analyzing the instability of the Antarctic ice sheet. Their detection and identification normally rely on geophysical and remote sensing methods such as airborne radar echo sounding (RES), ground seismic, and satellite/airborne altimetry and gravity surveys. In particular, RES surveys are able to detect basal terrain with a relatively high accuracy that can assist with the mapping of subglacial hydrology systems. Traditional RES processing methods for the identification of subglacial water bodies mostly rely on their brightness in radargrams and hydraulic flatness. In this study, we propose an automatic method with the main objective to differentiate the basal materials by quantitatively evaluating the shape of the A-scope waveform near the basal interface in RES radargrams, which has been seldom investigated. We develop an automatic algorithm mainly based on the traditional short-time Fourier transform (STFT) to quantify the shape of the A-scope waveform in radargrams. Specifically, with an appropriate window width applied on the main peak of each A-scope waveform in the RES radargram, STFT shows distinct and contrasting frequency responses at the ice-water interface and ice-rock interface, which is largely dependent upon their different reflection characteristics at the basal interface. We apply this method on 882 RES radargrams collected in the Antarctic's Gamburtsev Province (AGAP) in East Antarctica. There are 8822 identified A-scopes with the calculated detection value larger than the set threshold, out of the overall 1,515,065 valid A-scopes in these 882 RES radargrams. Although these identified A-scopes only takes 0.58% of the overall A-scope population, they show exceptionally continuous distribution to represent the subglacial water bodies. Through a comprehensive comparison with existing inventories of subglacial lakes, we successfully verify the validity and advantages of our method in identifying subglacial water bodies using the detection probability for other basal materials of theoretically the highest along-track resolution. The frequency signature obtained by the proposed joint time–frequency analysis provides a new corridor of investigation that can be further expanded to multivariable deep learning approaches for subglacial and englacial material characterization, as well as subglacial hydrology mapping. [ABSTRACT FROM AUTHOR]

An automatic geomagnetic station for monitoring the Earth's magnetic field variations was installed in December 2020 at Talos Dome, a remote site on the Antarctic Plateau, about 300 km away from the permanent geomagnetic observatory at Mario Zucchelli Station (MZS). Designed and assembled at the laboratory of electronics of the Istituto Nazionale di Geofisica e Vulcanologia (INGV) in Rome, this autonomous station is formed by a vector magnetometer specifically manufactured by Lviv Institute (Ukraine) for very low temperatures and a low-power system supplied by batteries charged by a wind generator and solar panel. Data, sampled at 1 Hz, are locally stored and can be downloaded once a year during the Antarctic summer expeditions. The goal was to integrate observatory data for better monitoring the geomagnetic field from an uncovered Antarctic area. In fact, it is well known that the distribution of geomagnetic observatories strongly favors the northern hemisphere, and each new instrumental installation in Antarctica should be considered as a useful attempt to balance the geomagnetic monitoring in the two hemispheres. The achieved goal was to obtain a long data series, keeping the station working even during the austral winter when the temperature can reach −60 °C; we recorded almost 11 months of data in one year and the station is still operating. Data from the new station, jointly with data from permanent observatories, improve the analysis of the magnetospheric dynamics and the ionosphere–magnetosphere coupling. Talos Dome, together with the Italian geomagnetic observatory at Mario Zucchelli Station and New Zealand geomagnetic observatory at Scott Base, constitutes a network along the 80°S geomagnetic parallel, which is interesting for studying the longitudinal propagation of geomagnetic signals of external origin. In this work we present the characteristics of the station and of the data it provides, with the aim of them for analysis in the framework of space weather. [ABSTRACT FROM AUTHOR]

The Gamburtsev Subglacial Mountains (GSMs) in the interior East Antarctic Craton are entirely buried under the massive East Antarctic ice sheet, with a ~50–60 km thick crust and ~200 km thick lithosphere, but little is known of the crustal structure and uplift mechanism. Here, we use airborne gravity and aeromagnetic anomalies for characteristic analysis and inverse calculations. The gravity and magnetic images show three distinct geophysical domains. Based on the gravity anomalies, a dense lower crustal root is modelled to underlie the GSMs, which may have formed by underplating during the continental collision of Antarctica and India. The high frequency linear magnetic characteristics parallel to the suture zone suggest that the upper crustal architecture is dominated by thrusts, consisting of a large transpressional fault system with a trailing contractional imbricate fan. A 2D model along the seismic profile is created to investigate the crustal architecture of the GSMs with the aid of depth to magnetic source estimates. Combined with the calculated crustal geometry and physical properties and the geological background of East Antarctica, a new evolutionary model is proposed, suggesting that the GSMs are underlain by part of a Pan-African age advancing accretionary orogen superimposed on Precambrian basement. [ABSTRACT FROM AUTHOR]

The Antarctic ice sheet and the continent both play critical roles in global sea level rise and climate change but they remain poorly understood because data collection is greatly limited by the remote location and hostile conditions there. Airborne platforms have been extensively used in Antarctica due to their capabilities and flexibility and have contributed a great deal of knowledge to both the ice sheet and the continent. The Snow Eagle 601 fixed-wing airborne platform has been deployed by China for Antarctic expeditions since 2015. Scientific instruments on the airplane include an ice-penetrating radar, a gravimeter, a magnetometer, a laser altimeter, a camera and a Global Navigation Satellite System (GNSS). In the past five austral seasons, the airborne platform has been used to survey Princess Elizabeth Land, the largest data gap in Antarctica, as well as other critical areas. This paper reviews the scientific operations of Snow Eagle 601 including airborne and ground-based scientific instrumentation, aviation logistics, field data acquisition and processing and data quality control. We summarize the progress of airborne surveys to date, focusing on scientific motivations, data coverage and national and international collaborations. Finally, we discuss potential regions for applications of the airborne platform in Antarctica and developments of the airborne scientific system for future work. [ABSTRACT FROM AUTHOR]

The ionosphere provides a channel that is able to propagate electromagnetic waves for long-haul communications, allowing Non-Line-Of-Sight (NLOS) radio communications in the high-frequency band (HF). Nonetheless, its performance as a communications channel mainly depends on solar activity, as Earth's diurnal and seasonal cycles modify the ionospheric ionization by the sun and consequently the channel performance. La Salle and the Observatori de l'Ebre have conducted oblique soundings of a 12,760 km ionospheric channel from Livingston Island (South Shetland Archipelago—Antarctica − 62.7 ∘ S, 299.6 ∘ E) to Cambrils (Spain − 41.0 ∘ N, 1.0 ∘ E) to evaluate this evidence and analyze the characteristics of this particular channel. The final goal of the project is to establish a stable communications link to be used as a backup for low throughput data transmission from the various geomagnetic and meteorological sensors located in the Spanish Antarctic Station. The aim of this paper is to confirm the relationship between the channel availability of narrowband and wideband communications and the resulting channel parameters with the solar and ionospheric activity for four consecutive sounding campaigns. [ABSTRACT FROM AUTHOR]

Protected areas (PAs) have been established worldwide for achieving long-term goals in the conservation of nature with the associated ecosystem services and cultural values. Globally, 15% of the world's terrestrial lands and inland waters, excluding Antarctica, are designated as PAs. About 4.12% of the global ocean and 10.2% of coastal and marine areas under national jurisdiction are set as marine protected areas (MPAs). Protected lands and waters serve as the fundamental building blocks of virtually all national and international conservation strategies, supported by governments and international institutions. Some of the PAs are the only places that contain undisturbed landscape, seascape and ecosystems on the planet Earth. With intensified impacts from climate and environmental change, PAs have become more important to serve as indicators of ecosystem status and functions. Earth's remaining wilderness areas are becoming increasingly important buffers against changing conditions. The development of remote sensing platforms and sensors and the improvement in science and technology provide crucial support for the monitoring and management of PAs across the world. In this editorial paper, we reviewed research developments using state-of-the-art remote sensing technologies, discussed the challenges of remote sensing applications in the inventory, monitoring, management and governance of PAs and summarized the highlights of the articles published in this Special Issue. [ABSTRACT FROM AUTHOR]

The HY-1C satellite, as part of China's optical satellite constellation for global ocean monitoring, monitors the ocean and coastal environment by the three broad visible bands of the Coastal Zone Imager (CZI) instrument. However, as a result of the sensor instrument noise, the atmospheric environment during imaging, and the shooting angle, the satellite images often show uneven illumination and inconsistent color between neighboring images. In this paper, according to the characteristics of the HY-1C CZI instrument, we propose a color consistency processing framework for coastal zone images of Antarctica. First of all, the high-frequency and low-frequency information of the image is separated by a statistical filter with simple clustering. The uneven lighting is then replaced by artificial lighting, which is globally uniform. Finally, the color difference between images is corrected by a color transfer method. In order to evaluate the color consistency results quantitatively, a new quantitative evaluation method is proposed. The experimental results for the coastal zone images of Antarctica show that the new processing framework can effectively eliminate the unevenness in the lighting and color. The mosaic results show a good performance in consistent lighting and tones, and the lack of visible mosaic lines proves the effectiveness of the proposed method. The quantitative evaluation analysis confirms the superiority of the proposed method over the Wallis method. [ABSTRACT FROM AUTHOR]