East Antarctic Research Articles

Glacial-Holocene climate-driven shifts in lacustrine and terrestrial environments: Rock magnetic and geochemical evidence from East Antarctic Mochou Lake

Geomorphic reconstructions of the East Antarctic Ice Sheet history across Antarctica suggest diachronous and varying retreat patterns with changing climate. However, little is known of how terrestrial environments responded to these climate changes because the continental ice sheet covered most Antarctic coastal oases during the Last Glacial Maximum (LGM), eliminating terrestrial records. In this study, the environmental history spanning the last 25,400 years was reconstructed using geochronology, lithology, rock magnetism and geochemistry to understand the effect of glacial–deglacial climate variations on surface processes and Mochou Lake ecology. During the last glacial period (25.4 to 18.8 cal. kyr BP), intensified winds and the freeze-thaw action of ice actively weathered the high elevation western catchment that was a dominant sediment source to the lake. The glacial environmental conditions induced a persistent perennial lake ice-cover resulting in anoxia, low sedimentation, primary production and weak brackish conditions. Progressively, as deglacial conditions augmented around 18.8 cal. kyr BP, increased catchment meltwater flow enhanced sediment transport and also transformed the lake into a freshwater basin, although perennial ice-cover still prevailed. Further, the biotic development of the lake began later in the Mid-Holocene as compared to other lakes in the region, wherein the seasonal lake ice cycle and/or ice moat developed. Weak chemical weathering and pedogenesis also commenced in the catchment in response to the warming climate trend. This delayed response to the Holocene warming was possibly due to the local cooling effect of the Dälk Glacier's discharge that also led to sea-ice persistence in the Prydz Bay. This study also provides further evidence that parts of the Larsemann Hills remained free of erosive grounded ice during the LGM and subsequent deglaciation. The Mirror Peninsula, specifically, escaped ice overriding despite being enclosed by the Dälk Glacier, Lake Nella ice lobe and Amery Ice Shelf.

Surface Mass Balance Models Vs. Stake Observations: A Comparison in the Lake Vostok Region, Central East Antarctica

The surface mass balance (SMB) is very low over the vast East Antarctic Plateau, for example in the Vostok region, where the mean SMB is on the order of 20–35 kg m-2 a-1. The observation and modeling of spatio-temporal SMB variations are equally challenging in this environment. Stake measurements carried out in the Vostok region provide SMB observations over half a century (1970–2019). This unique data set is compared with SMB estimations of the regional climate models RACMO2.3p2 (RACMO) and MAR3.11 (MAR). We focus on the SMB variations over time scales from months to decades. The comparison requires a rigorous assessment of the uncertainty in the stake observations and the spatial scale dependence of the temporal SMB variations. Our results show that RACMO estimates of annual and multi-year SMB agree well with the observations. The regression slope between modelled and observed temporal variations is close to 1.0 for this model. SMB simulations by MAR are affected by a positive bias which amounts to 6 kg m-2 a-1 at Vostok station and 2 kg m-2 a-1 along two stake profiles between Lake Vostok and Ridge B. None of the models is capable to reproduce the seasonal distributions of SMB and precipitation. Model SMB estimates are used in assessing the ice-mass balance and sea-level contribution of the Antarctic Ice Sheet by the input-output method. Our results provide insights into the uncertainty contribution of the SMB models to such assessments.

Sizes and Shapes of Sea Ice Floes Broken by Waves–A Case Study From the East Antarctic Coast

The floe size distribution (FSD) is an important characteristics of sea ice, influencing several physical processes that take place in the oceanic and atmospheric boundary layers under/over sea ice, as well as within sea ice itself. Through complex feedback loops involving those processes, FSD might modify the short-term and seasonal evolution of the sea ice cover, and therefore significant effort is undertaken by the scientific community to better understand FSD-related effects and to include them in sea ice models. An important part of that effort is analyzing the FSD properties and variability in different ice and forcing conditions, based on airborne and satellite imagery. In this work we analyze a very high resolution (pixel size: 0.3 m) satellite image of sea ice from a location off the East Antarctic coast (65.6°S, 101.9°E), acquired on February 16, 2019. Contrary to most previous studies, the ice floes in the image have angular, polygonal shapes and a narrow size distribution. We show that the observed FSD can be represented as a weighted sum of two probability distributions, a Gaussian and a tapered power law, with the Gaussian part clearly dominating in the size range of floes that contribute over 90% to the total sea ice surface area. Based on an analysis of the weather, wave and ice conditions in the period preceding the day in question, we discuss the most probable scenarios that led to the breakup of landfast ice into floes visible in the image. Finally, theoretical arguments backed up by a series of numerical simulations of wave propagation in sea ice performed with a scattering model based on the Matched Eigenfunction Expansion Method are used to show that the observed dominating floe size in the three different regions of the image (18, 13 and 51 m, respectively) agree with those expected as a result of wave-induced breaking of landfast ice.

Linking Sub‐Tropical Evaporation and Extreme Precipitation Over East Antarctica: An Atmospheric River Case Study

AbstractWe investigate an intense snowfall event between 15 and 18 February 2011 over the East Antarctic coastal region which contributed to roughly 24% of the annual snow accumulation. The event was previously associated with an atmospheric river, and here we use both Eulerian and Lagrangian analysis to gain an understanding of the processes contributing to the atmospheric river signature. The planetary‐scale configuration during the event consisted of a persistent blocking situation resulting in a sustained meridional flow from the sub‐tropics to the Antarctic ice sheet between 20 and 50°E. Within this configuration, synoptic‐scale cyclogenesis contributed to slantwise ascent of moisture loaded air parcels toward Antarctica. Landfall of this cyclone’s warm sector coincided with the onset of Antarctic precipitation. Subsequently, a secondary cyclone developed along a pre‐existing baroclinic zone. The rapid intensification and propagation speed of this mesoscale cyclone alongside the warm, moist air mass resulted in strong moisture flux convergence ahead of the cyclone, providing additional poleward moisture transport. The poleward progression of warm moist air and a corresponding decrease of sea‐surface temperatures implied downward surface sensible and latent heat fluxes throughout the region of intense poleward moisture, roughly between 40 and 60°S. Hence, moisture uptake via surface evaporation was suppressed between the sub‐tropics and the polar continent, favoring long‐range transport. Identification of the surface moisture uptake region by tracing changes in moisture in air parcels confirmed the limited uptake of moisture during the poleward transport in this case study, with the primary moisture source for Antarctic precipitation located in the sub‐tropics.

Shallow-Layers-Detection Ice Sounding Radar for Mapping of Polar Ice Sheets

The accumulation rate is a key parameter in computing the mass balance of glaciers and ice sheets to estimate sea level rise. A shallow-layers-detection ice sounding radar (SLDISR) is developed to measure the accumulation rate and shape of near-surface internal layers with high resolution. With a transmitting frequency from 500 to 2000 MHz, this frequency-modulated continuous wave (FMCW) radar provides a range resolution of about 16 cm in free space by using a Hanning window and a penetrating depth about 150 m under polar ice. The spectral analysis and coherent integration techniques are used to obtain a high processing gain and to improve the signal-to-noise ratio of the system. A phase-locked loop with wideband yttrium iron garnet (YIG) oscillator is applied to generate a sweeping chirp signal as an input source for the transmitter. A stable, low-frequency reference chirp signal is generated with a direct digital synthesizer (DDS) integrated in field-programmable gate array (FPGA). To reduce the high-speed requirement to the analog-to-digital converter (ADC), dechirp technology is adopted at the RF section of the receiver. The implementation of the digital unit is based on an FPGA chip. The designed radar has been successfully deployed in Antarctica during the 31st Chinese Antarctic Research Expedition (CHINARE 31) and CHINARE 33, mainly over the East Antarctic Ice Sheet (EAIS). The echograms indicate the effectiveness of the radar system on detecting clear internal reflecting horizons (IRHs) over ice sheets.

Joint Inversion for Surface Accumulation Rate and Geothermal Heat Flow From Ice‐Penetrating Radar Observations at Dome A, East Antarctica. Part II: Ice Sheet State and Geophysical Analysis

AbstractDome A is the summit of the East Antarctic Ice Sheet, underlain by the rugged Gamburtsev Subglacial Mountains. The rugged basal topography produces a complex hydrological system featuring basal melt, water transport and storage, and freeze‐on. In a companion study, we used an inverse model to infer the spatial distributions of geothermal heat flow (GHF) and accumulation rate that best fit a variety of observational constraints. Here, we present and analyze the best‐fit state of the ice sheet in detail. Our modeled result agrees well with the observed water bodies and freeze‐on structures, while also predicting a significant amount of unobserved water and suggesting a change in stratigraphic interpretation that reduces the volume of the freeze‐on units. Our modeled stratigraphy agrees well with observations, and we predict that there will be two distinct patches of ice up to 1.5 Ma suitable for ice coring underneath the divide. Past divide migration could have interrupted stratigraphic continuity at the old ice patches, but various indirect lines of evidence suggest that the divide has been stable for about the last one and a half glacial cycles, which is an encouraging but not definitive sign for stability in the longer term. Finally, our GHF estimate is higher than previous estimates for this region, but consistent with possible heterogeneity in crustal heat production.

Multistage Evolution of Proterozoic Crust of East Antarctica by the Example of the Filla Terrane (Rauer Islands): New Geological and Isotope Data

Abstract —The paper presents new data on the Rauer Islands, one of the unique objects of the East Antarctic Shield. The interest in this area is triggered by its complex geologic structure, including both Archean and Proterozoic fragments of the Earth’s crust, and by its multiphase formation. A detailed scheme of the geologic structure of the area is proposed, new petrologic complexes are revealed, and the stages of tectonomagmatic activity at ~1400–1320 Ma and 1150 Ma are reliably dated. This serves as a factual basis for comparison the study area with other regions of East Antarctica. Based on the geological and isotope data obtained, the Meso–Neoproterozoic Filla Terrane in the area of the Rauer Islands is recognized. It is composed of metamorphic and primarily intrusive rocks, whose protoliths formed in the time interval 1400–950 Ma. Three periods of tectonothermal activity have been established in the Filla Terrane: Mid-Mesoproterozoic (1400–1320 Ma), Meso–Neoproterozoic (1150–886 Ma), and early Cambrian (536–504 Ma). The first period is the formation time of Mesoproterozoic crust, and it is time-correlated with the tectonogenesis phase in the adjacent Rayner province. The second period corresponds to the later phase of tectonothermal activity in the Rayner province. In the Filla Terrane, this period can be divided into two intervals, 1150–1100 Ma and 1010–886 Ma. The former interval is treated as intense crustal growth in the course of granitoid and mantle magmatism. The latter interval is a period of tectonothermal processes accompanied by intense deformations, high-temperature metamorphism, and intrusion of porphyritic granitoids. Apparently, the gap between the first and the second intervals is the time of deposition of the sedimentary protolith of paragneisses, which, together with the surrounding rocks, underwent high-temperature metamorphism and deformations at 950–914 Ma. The synchronous evolution of the Archean block and the Filla Terrane began at least within 1100–1000 Ma. The youngest, early Cambrian period of tectonic activity coincides with the development of local low-temperature mylonite zones and the intrusion of synkinematic pegmatite veins. Thus, the tectonothermal evolution of the Filla Terrane includes almost the same main phases of crustal growth and transformation as the Rayner province. This indicates that the Filla Terrane is a fragment of the Rayner province, which accreted to the Archean terrane at least in the late Mesoproterozoic.

The influence of Antarctic ice loss on polar motion: an assessment based on GRACE and multi-mission satellite altimetry

Increasing ice loss of the Antarctic Ice Sheet (AIS) due to global climate change affects the orientation of the Earth’s spin axis with respect to an Earth-fixed reference system (polar motion). Here the contribution of the decreasing AIS to the excitation of polar motion is quantified from precise time variable gravity field observations of the Gravity Recovery and Climate Experiment (GRACE) and from measurements of the changing ice sheet elevation from altimeter satellites. While the GRACE gravity field models need to be reduced by noise and leakage effects from neighboring subsystems, the ice volume changes observed by satellite altimetry have to be converted into ice mass changes. In this study we investigate how much individual gravimetry and altimetry solutions differ from each other. We show that due to combination of individual solutions systematic and random errors of the data processing can be reduced and the robustness of the geodetic derived AIS polar motion excitations can be increased. We investigate the interannual variability of the Antarctic polar motion excitation functions by means of piecewise linear trends. We find that the long-term behavior of the three ice sheet subregions: EAIS (East Antarctic Ice Sheet), WAIS (West Antarctic Ice Sheet) and APIS (Antarctic Peninsula Ice Sheet) is quite different. While APIS polar motion excitations show no significant interannual variations during the study period 2003-2015, the trend of the WAIS and EAIS polar motion excitations increased in 2006 and again in 2009 while it started slightly to decline in 2013. AIS mass changes explain about 45% of the observed magnitude of the polar motion vector (excluding glacial isosatic adjustment). They cause the pole position vector to drift along 59^{circ } East longitude with an amplitude of 2.7 mas/yr. Thus the contribution of the AIS has to be considered to close the budget of the geophysical excitation functions of polar motion.

Present and Future of Rainfall in Antarctica

AbstractWhile most precipitation in Antarctica falls as snow, little is known about liquid precipitation, although it can have ecological and climatic impacts. This study combines meteorological reports at 10 stations with the ERA5 reanalysis to provide a climatological characterization of rainfall occurrence over Antarctica. Along the East Antarctic coast, liquid precipitation occurs 22 days per year at most and coincides with maritime intrusions and blocking anticyclones. Over the north‐western Antarctic Peninsula, rainfall occurs more than 50 days per year on average and the recent summer cooling was accompanied by a decrease of −35 annual rainy days per decade between 1998 and 2015 at Faraday‐Vernadsky. Projections from seven latest‐generation climate models reveal that Antarctic coasts will experience a warming and more frequent and intense rainfall by the end of the century. Rainfall is expected to impact new regions of the continent, increasing their vulnerability to melting by the preconditioning of surface snow.

The Long‐Term Cooling Trend in East Antarctic Plateau Over the Past 2000 Years Is Only Robust Between 550 and 1550 CE

AbstractThe uncertainties in Antarctic climate reconstructions due to scarcity of proxy records have restricted the understanding of mechanisms of climate change, and further hindered the improvement of climate models. Here, we provide a new climate record derived from water stable isotopes in a Dome A, East Antarctica ice core. Together with six other ice core records, the Dome A record is used to investigate temperature changes in East Antarctic Plateau (EAP) during period 1–1900 CE. Our results show that, a previously reported long‐term cooling trend in EAP during the recent (pre‐1900 CE) 1900 years is only robust between 550 and 1550 CE. A combination of solar and volcanic forcing may have induced the EAP centennial‐scale cold events, and further caused the long‐term cooling trend from 550 to 1550 CE with a small contribution from orbital forcing.

Aerogeophysical characterization of Titan Dome, East Antarctica, and potential as an ice core target

Abstract. Based on sparse data, Titan Dome has been identified as having a higher probability of containing ice that would capture the middle Pleistocene transition (1.25 to 0.7 Ma). New aerogeophysical observations (radar and laser altimetry) collected over Titan Dome, located about 200 km from the South Pole within the East Antarctic Ice Sheet, were used to characterize the region (e.g., geometry, internal structure, bed reflectivity, and flow history) and assess its suitability as a paleoclimate ice core site. The radar coupled with an available ice core chronology enabled the tracing of dated internal reflecting horizons throughout the region, which also served as constraints on basal ice age modeling. The results of the survey revealed new basal topographic detail and better constrain the ice topographical location of Titan Dome, which differs between community datasets. Titan Dome is not expected to be relevant to the study of the middle Pleistocene transition due to a combination of past fast flow dynamics, the basal ice likely being too young, and the temporal resolution likely being too coarse if 1 Ma ice were to exist.

Deriving Melt Rates at a Complex Ice Shelf Base Using In Situ Radar: Application to Totten Ice Shelf

AbstractA phase‐sensitive radar (ApRES) was deployed on Totten Ice Shelf to provide the first in situ basal melt estimate at this dynamic East Antarctic ice shelf. Observations of internal ice dynamics at tidal time scales showed that early arrivals from off‐nadir reflectors obscure the true depth of the ice shelf base. Using the observed tidal deformation, the true base was found to lie at 1,910–1,950‐m depth, at 350–400 m greater range than the first reflection from an ice‐ocean interface. The robustness of the basal melt rate estimate was increased by using multiple basal reflections over the radar footprint, yielding a melt rate of 22 ± 2.1 m a−1. The ApRES estimate is over 40% lower than the three existing satellite estimates covering Totten Ice Shelf. This difference in basal melt is dynamically significant and highlights the need for independent melt rate estimates using complementary instrumentation and techniques that rely on different sets of assumptions.

Uncertainty in East Antarctic Firn Thickness Constrained Using a Model Ensemble Approach

AbstractMass balance assessments of the East Antarctic ice sheet (EAIS) are highly sensitive to changes in firn thickness, causing substantial disagreement in estimates of its contribution to sea‐level. To better constrain the uncertainty in recent firn thickness changes, we develop an ensemble of 54 model scenarios of firn evolution between 1992 and 2017. Using statistical emulation of firn‐densification models, we quantify the impact of firn compaction formulation, differing climatic forcing, and surface snow density on firn thickness evolution. At basin scales, the ensemble uncertainty in firn thickness change ranges between 0.2 and 1.0 cm yr−1 (15%–300% relative uncertainty), with the choice of climate forcing having the largest influence on the spread. Our results show the regions of the ice sheet where unexplained discrepancies exist between observed elevation changes and an extensive set of modeled firn thickness changes estimates, marking an important step toward more accurately constraining ice sheet mass balance.

Determination of Cd, Pb, and Cu in the Atmospheric Aerosol of Central East Antarctica at Dome C (Concordia Station).

Trace heavy metals Cd, Pb, and Cu were determined (by square wave anodic stripping voltammetry) in aerosol samples collected at Dome C (the Italo-French Station Concordia), a remote site of the Central East Antarctic plateau, for which no data are available until now. During the Austral Summer 2005–2006, three PM10 high-volume impactors were installed in two locations nearby of Concordia station: the first one very close and downwind of the station (about 50 m north), the other two (very close to each other) in a ‘distant’ site, upwind of the station and close to the astrophysics tent (not used in that expedition) at ~800 m south of Station Concordia. For each sample, the availability of the mass of the aerosol collected (obtained by differential weighing carried out on site), in addition to the volume of the filtered air, allowed us to express results both in terms of metal mass fractions in the aerosol and in the usual way of metal atmospheric concentrations. Metal contents increased in the order Cd < Pb < Cu with the following ranges of values: Cd 1.0–8.4 µg g−1 (0.09–3.1 pg m−3), Pb 96–470 µg g−1 (12–62 pg m−3), and Cu 0.17–20 mg g−1 (0.027–2.4 ng m−3). From the metal temporal profiles obtained we estimated the following background values for the area of Dome C, expressed both in mass fractions and in atmospheric concentrations: Cd 1.2 ± 0.2 µg g−1 (0.24 ± 0.13 pg m−3), Pb (here fixed as upper limit) 113 ± 13 µg g−1 (21 ± 8 pg m−3), and Cu 0.91 ± 0.48 mg g−1 (0.12 ± 0.07 ng m−3). The highest values were observed in the first part of the season, and particularly for the site close to the station, possibly related to sample contamination linked to intense activity at the Concordia station connected with the beginning of the expedition, including aircraft arrivals/departures. Increments of up to 10 times (and even 20 times for Cu) were recorded with respect to the background values. The metal excesses of the contaminated over background samples were found approximately, except for Cu, in the same proportion of the metal contents of the special Antarctic blend (SAB) diesel fuel, which is used almost exclusively at Concordia Station. The effect of the wind direction was also observed. Thus in the intermediate period of the campaign, when the wind direction reversed for several days with respect to the prevailing one, Cd and Pb metal contents decreased at the sampling point installed close to the station, now upwind of Concordia station, and increased at the ‘clean’ site astrophysics tent, turned downwind at the main station. No simple and easily interpretable effect of the wind direction was observed for Cu, which suggests that some other extemporaneous and not clearly identified factor may have intervened in this case. These results suggest that the human impact at Dome C influences mainly the zone very close to the station, but also the area in the neighborhood, including the supposed clean site of the astrophysics tent (about 800 m far from the station), when the wind direction reverses with respect to the prevailing one, leaving the site downwind of the station Concordia. Since no other data are reported for the Dome C area, our results are compared with literature data referred to the South Pole Station (the only other plateau site for which data are available) and several other coastal Antarctic sites, observing that our results (excluding Cu) are the lowest ever observed for Antarctic aerosol.

Atmospheric carbon dioxide variations across the middle Miocene climate transition

Abstract. The middle Miocene climate transition ∼ 14 Ma marks a fundamental step towards the current “ice-house” climate, with a ∼ 1 ‰ δ18O increase and a ∼ 1 ‰ transient δ13C rise in the deep ocean, indicating rapid expansion of the East Antarctic Ice Sheet associated with a change in the operation of the global carbon cycle. The variation of atmospheric CO2 across the carbon-cycle perturbation has been intensely debated as proxy records of pCO2 for this time interval are sparse and partly contradictory. Using boron isotopes (δ11B) in planktonic foraminifers from Ocean Drilling Program (ODP) Site 1092 in the South Atlantic, we show that long-term pCO2 varied at 402 kyr periodicity between ∼ 14.3 and 13.2 Ma and follows the global δ13C variation remarkably well. This suggests a close link to precessional insolation forcing modulated by eccentricity, which governs the monsoon and hence weathering intensity, with enhanced weathering and decreasing pCO2 at high eccentricity and vice versa. The ∼ 50 kyr lag of δ13C and pCO2 behind eccentricity in our records may be related to the slow response of weathering to orbital forcing. A pCO2 drop of ∼ 200 µatm before 13.9 Ma may have facilitated the inception of ice-sheet expansion on Antarctica, which accentuated monsoon-driven carbon cycle changes through a major sea-level fall, invigorated deep-water ventilation, and shelf-to-basin shift of carbonate burial. The temporary rise in pCO2 following Antarctic glaciation would have acted as a negative feedback on the progressing glaciation and helped to stabilize the climate system on its way to the late Cenozoic ice-house world.

Spatial Distribution of Crusts in Antarctic and Greenland Snowpacks and Implications for Snow and Firn Studies

The occurrence of snowpack features has been used in the past to classify environmental regimes on the polar ice sheets. Among these features are thin crusts with high density, which contribute to firn stratigraphy and can have significant impact on firn ventilation as well as on remotely inferred properties like accumulation rate or surface mass balance. The importance of crusts in polar snowpack has been acknowledged, but nonetheless little is known about their large-scale distribution. From snow profiles measured by means of microfocus X-ray computer tomography we created a unique dataset showing the spatial distribution of crusts in snow on the East Antarctic Plateau as well as in northern Greenland including a measure for their local variability. With this method, we are able to find also weak and oblique crusts, to count their frequency of occurrence and to measure the high-resolution density. Crusts are local features with a small spatial extent in the range of tens of meters. From several profiles per sampling site we are able to show a decreasing number of crusts in surface snow along a traverse on the East Antarctic Plateau. Combining samples from Antarctica and Greenland with a wide range of annual accumulation rate, we find a positive correlation (R2 = 0.89) between the logarithmic accumulation rate and crusts per annual layer in surface snow. By counting crusts in two Antarctic firn cores, we can show the preservation of crusts with depth and discuss their temporal variability as well as the sensitivity to accumulation rate. In local applications we test the robustness of crusts as a seasonal proxy in comparison to chemical records like impurities or stable water isotopes. While in regions with high accumulation rates the occurrence of crusts shows signs of seasonality, in low accumulation areas dating of the snowpack should be done using a combination of volumetric and stratigraphic elements. Our data can bring new insights for the study of firn permeability, improving of remote sensing signals or the development of new proxies in snow and firn core research.

A 43-ka mineral magnetic record of environmental variations from lacustrine sediments of Schirmacher Oasis, East Antarctica

Glacial-interglacial climatic and environmental variations recorded in lake-bed sediments of Schirmacher Oasis (SO), East Antarctica, are rare; limited data exist on the past behavior of the East Antarctic Ice Sheet (EAIS). In this work, we present the mineral magnetic properties of a radiocarbon-dated sediment core, spanning the past 43,000 years, from a land-locked lake L-49 in SO. We also performed Scanning Electron Microscopic-Energy Dispersive X-ray Spectroscopic (SEM-EDS) studies on magnetic extracts to confirm the magnetic mineralogy. The values of magnetic susceptibility (χlf) and saturation isothermal remanent magnetization (SIRM) are relatively low during the Holocene but high during the last glacial period. Sediments of the last glacial period are composed mainly of coarse-grained (multi-domain - MD) magnetite/titanomagnetite (high values of S20 – an MD proxy, low susceptibility of anhysteretic remanent magnetism (χARM)/SIRM ratio values and high S-ratio values) indicating a stronger intensity of physical weathering in the Lake L-49 catchment. The mineralogy is confirmed by SEM-EDS data for magnetic extracts. Principal Component Analysis (PCA) of the data suggests three processes affecting magnetic properties of the lake sediments, i.e., ferrimagnetic input, biogenic processes and aeolian input. During deglaciation, a shift in the magnetic properties is documented: a decrease in χlf values from 18.6 cal ka BP onwards suggests melting of the snow and ice in the catchment thus leading to a significant production of fine-grained iron oxide minerals. At ~12.77 cal ka BP, high values of χlf and S20 ratio indicates cooling which is related to the Antarctic Cold Reversal (ACR). The Early Holocene Climate Optimum is distinctly imprinted in our record as a sudden decrease in χlf values and an increase in magnetic grain size ratio values (indicating finer iron oxide grains), suggesting a rapid thawing of the ice-sheet in response to a warm climate. The Mid-Holocene Hypsithermal (from 5.8 to 3.5 cal ka BP) and Neoglaciation (3.6 cal ka BP onwards) that are evident in this study correspond with reconstructions from lake and marine sediment, and ice-core records from Antarctica.

Geochemical zones and environmental gradients for soils from the central Transantarctic Mountains, Antarctica

Abstract. Previous studies have established links between biodiversity and soil geochemistry in the McMurdo Dry Valleys, Antarctica, where environmental gradients are important determinants of soil biodiversity. However, these gradients are not well established in the central Transantarctic Mountains, which are thought to represent some of the least hospitable Antarctic soils. We analyzed 220 samples from 11 ice-free areas along the Shackleton Glacier (∼ 85∘ S), a major outlet glacier of the East Antarctic Ice Sheet. We established three zones of distinct geochemical gradients near the head of the glacier (upper), its central part (middle), and at the mouth (lower). The upper zone had the highest water-soluble salt concentrations with total salt concentrations exceeding 80 000 µg g−1, while the lower zone had the lowest water-soluble N:P ratios, suggesting that, in addition to other parameters (such as proximity to water and/or ice), the lower zone likely represents the most favorable ecological habitats. Given the strong dependence of geochemistry on geographic parameters, we developed multiple linear regression and random forest models to predict soil geochemical trends given latitude, longitude, elevation, distance from the coast, distance from the glacier, and soil moisture (variables which can be inferred from remote measurements). Confidence in our random forest model predictions was moderately high with R2 values for total water-soluble salts, water-soluble N:P, ClO4-, and ClO3- of 0.81, 0.88, 0.78, and 0.74, respectively. These modeling results can be used to predict geochemical gradients and estimate salt concentrations for other Transantarctic Mountain soils, information that can ultimately be used to better predict distributions of soil biota in this remote region.

Carbonaceous Aerosol in Polar Areas: First Results and Improvements of the Sampling Strategies

While more and more studies are being conducted on carbonaceous fractions—organic carbon (OC) and elemental carbon (EC)—in urban areas, there are still too few studies about these species and their effects in polar areas due to their very low concentrations; further, studies in the literature report only data from intensive campaigns, limited in time. We present here for the first time EC–OC concentration long-time data records from the sea-level sampling site of Ny-Ålesund, in the High Arctic (5 years), and from Dome C, in the East Antarctic Plateau (1 year). Regarding the Arctic, the median (and the interquartile range (IQR)) mass concentrations for the years 2011–2015 are 352 (IQR: 283–475) ng/m3 for OC and 4.8 (IQR: 4.6–17.4) ng/m3 for EC, which is responsible for only 3% of total carbon (TC). From both the concentration data sets and the variation of the average monthly concentrations, the influence of the Arctic haze on EC and OC concentrations is evident. Summer may be interested by high concentration episodes mainly due to long-range transport (e.g., from wide wildfires in the Northern Hemisphere, as happened in 2015). The average ratio of EC/OC for the summer period is 0.05, ranging from 0.02 to 0.10, and indicates a clean environment with prevailing biogenic (or biomass burning) sources, as well as aged, highly oxidized aerosol from long-range transport. Contribution from ship emission is not evident, but this result may be due to the sampling time resolution. In Antarctica, a 1 year-around data set from December 2016 to February 2018 is shown, which does not present a clear seasonal trend. The OC median (and IQR) value is 78 (64–106) ng/m3; for EC, it is 0.9 (0.6–2.4) ng/m3, weighing for 3% on TC values. The EC/OC ratio mean value is 0.20, with a range of 0.06–0.35. Due to the low EC and OC concentrations in polar areas, correction for the blank is far more important than in campaigns carried out in other regions, largely affecting uncertainties in measured concentrations. Through the years, we have thus developed a new sampling strategy that is presented here for the first time: samplers were modified in order to collect a larger amount of particulates on a small surface, enhancing the capability of the analytical method since the thermo-optical analyzer is sensitive to carbonaceous aerosol areal density. Further, we have recently coupled such modified samplers with a sampling strategy that makes a more reliable blank correction of every single sample possible.

Beryllium isotopes in sediments from Lake Maruwan Oike and Lake Skallen, East Antarctica, reveal substantial glacial discharge during the late Holocene

Constraining East Antarctic Ice Sheet (EAIS) evolution during the Holocene is important for exploring the forcing mechanisms behind ice sheet retreat and to constrain numerical ice sheet models that aid predictions of future sea-level rise. Beryllium (Be) isotope analysis of bedrock and marine sediments have offered unparalleled insight into Antarctic ice sheet history since the Pliocene, but much of EAIS remains poorly studied. Here, we report the reactive (authigenic) 10Be abundance, 9Be abundance and 10Be/9Be ratios of Antarctic lake sediments, for the first time, from Lake Maruwan Oike and Lake Skallen along Soya Coast of Lützow-Holm Bay, East Antarctica. Beryllium isotope records reveal melting of local glaciers associated with higher subglacial erosion between ∼4.1 and ∼3.6 ka BP. Comparison to marine records from the Antarctic continental shelf suggests this was part of a circum-Antarctic phenomena that led to widespread glacial discharge from other sectors of the East and West Antarctic Ice Sheet. We suggest the incursion of relatively warm Circumpolar Deep Water (CDW) into Lützow-Holm Bay during the Late Holocene led to frontal and basal melting of ice sheets along the Soya Coast, supporting the notion of Antarctic ice sheet instabilities as a contributor to global sea-level rise since the Mid Holocene.