Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area, Central Swiss Alps, based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic 10Be and in situ14C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ14C data are overall in good agreement with the 10Be ages and confirm continuous exposure throughout the Holocene. However, in situ14C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.

Impact of a Reduced Arctic Sea Ice Cover on Ocean and Atmospheric Properties

The Arctic sea ice cover declined over the last few decades and reached a record minimum in 2007, with a slight recovery thereafter. Inspired by this the authors investigate the response of atmospheric and oceanic properties to a 1-yr period of reduced sea ice cover. Two ensembles of equilibrium and transient simulations are produced with the Community Climate System Model. A sea ice change is induced through an albedo change of 1 yr. The sea ice area and thickness recover in both ensembles after 3 and 5 yr, respectively. The sea ice anomaly leads to changes in ocean temperature and salinity to a depth of about 200 m in the Arctic Basin. Further, the salinity and temperature changes in the surface layer trigger a “Great Salinity Anomaly” in the North Atlantic that takes roughly 8 yr to travel across the North Atlantic back to high latitudes. In the atmosphere the changes induced by the sea ice anomaly do not last as long as in the ocean. The response in the transient and equilibrium simulations, while similar overall, differs in specific regional and temporal details. The surface air temperature increases over the Arctic Basin and the anomaly extends through the whole atmospheric column, changing the geopotential height fields and thus the storm tracks. The patterns of warming and thus the position of the geopotential height changes vary in the two ensembles. While the equilibrium simulation shifts the storm tracks to the south over the eastern North Atlantic and Europe, the transient simulation shifts the storm tracks south over the western North Atlantic and North America. The authors propose that the overall reduction in sea ice cover is important for producing ocean anomalies; however, for atmospheric anomalies the regional location of the sea ice anomalies is more important. While observed trends in Arctic sea ice are large and exceed those simulated by comprehensive climate models, there is little evidence based on this particular model that the seasonal loss of sea ice (e.g., as occurred in 2007) would constitute a threshold after which the Arctic would exhibit nonlinear, irreversible, or strongly accelerated sea ice loss. Caution should be exerted when extrapolating short-term trends to future sea ice behavior.

Retrieving the paleoclimatic signal from the deeper part of the EPICA Dome C ice core

An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (δD-δ18Oice, δ18Oatm, total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the "basal clean ice facies". Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies. It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate "from within", and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice–bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future "oldest ice" drilling location in Antarctica.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

The role of seasonality of mineral dust concentration and size on glacial/interglacial dust changes in the EPICA Dronning Maud Land ice core

We present a record of particulate dust concentration and size distribution in subannual resolution measured on the European Project for Ice Coring in Antarctica (EPICA) Dronning Maud Land (EDML) ice core drilled in the Atlantic sector of the East Antarctic plateau. The record reaches from present day back to the penultimate glacial until 145,000 years B.P. with subannual resolution from 60,000 years B.P. to the present. Mean dust concentrations are a factor of 46 higher during the glacial (~850–4600 ng/mL) compared to the Holocene (~16–112 ng/mL) with slightly smaller dust particles during the glacial comparedto the Holocene and with an absolute minimum in the dust size at 16,000 years B.P. The changes in dust concentration are mainly attributed to changes in source conditions in southern South America. An increase in the modal value of the dust size suggests that at 16,000 years B.P. a major change in atmospheric circulation apparently allowed more direct transport of dust particles to the EDML drill site. We find a clear in-phase relation of the seasonal variation in dust mass concentration and dust size during the glacial (r(conc,size) = 0.8) but no clear phase relationship during the Holocene (0 ice sheet using the size as an indicator for transport intensity, we find that the effect of the changes in the seasonality of the source emission strength and the transport intensity on the dust decrease over Transition 1 can significantly contribute to the large decrease of dust concentration from the glacial to the Holocene.

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

Fire ecology north and south of the Alps since the last ice age

Wildfires are very rare in central Europe, which is probably why fire effects on vegetation have been neglected by most central European ecologists and palaeoecologists. Presently, reconstructions of fire history and fire ecology are almost absent. We analysed sediment cores from lakes on the Swiss Plateau (Lobsigensee and Soppensee) for pollen and charcoal to investigate the relationship between vegetation and fire. Microscopic charcoal evidence suggests increasing regional fire frequencies during the Neolithic (7350-4150 cal. BP, 5400-2200 BC) and the subsequent prehistoric epochs at Lobsigensee, whereas at Soppensee burnings remained rather rare until modern times. Neolithic peaks of charcoal at 6200 and 5500 cal. BP (4250 and 3550 BC) coincided with declines of pollen of fire-sensitive taxa at both sites (e.g., Ulmus, Tilia, Hedera, Fagus), suggesting synchronous vegetational responses to fire at regional scales. However, correlation analysis between charcoal and pollen for the period 6600-4400 cal. BP (4650-2650 BC) revealed no significant link between fire and vegetation at Soppensee, whereas at Lobsigensee increases of Corylus and decreases of Fagus were related to fire events. Fire impact on vegetation increased during the subsequent epochs at both sites. Correlation analyses of charcoal and pollen data for the period 4250-1150 cal. BP (2300 BC -AD 800) suggest that fires were intentionally set to disrupt forests and to provide open areas for arable and pastoral farming (e.g., significant positive correlations between charcoal and Cerealia, Plantago lanceolata, Asteroideae). These results are compared with southern European records (Lago di Origlio, Lago di Muzzano), which are situated in particularly fire-prone environments. After the Mesolithic period (I1 200-7350 cal. BP, 9250-5400 BC), charcoal influx was higher by an order of magnitude in the south, suggesting more frequent fires. Neolithic fires caused similar though more pronounced responses of vegetation in the south (e.g., expansions of Corylus). Post-Neolithic land-use practices involving (controlled) burning culminated in both regions at about 2550 cal. BP (c. 600 BC). However, fire-caused disappearances of entire forest communities were confined to the southern sites. Such differences in fire effects among the sites are explained by the dissimilar importance of fire as a result of different climatic conditions and cultural activities. Our results imply that the remaining (fire-sensitive) fragments of central European vegetation north of the Alps are especially endangered by increasing fire frequencies resulting from predicted climatic change.

An ice-rich flow origin for the banded terrain in the Hellas basin, Mars

The interior of Hellas Basin displays a complex landscape and a variety of geomorphological domains. One of these domains, the enigmatic banded terrain covers much of the northwestern part of the basin. We use high-resolution (CTX and HiRISE) Digital Terrain Models to show that most of the complex viscous flowing behavior exhibited by the banded terrain is controlled by topography and flow-like interactions between neighboring banded terrain. Furthermore, the interior of the basin hosts several landforms suggestive of the presence of near-surface ice, which include polygonal patterns with elongated pits, scalloped depressions, isolated mounds and collapse structures. We suggest that thermal contraction cracking and sublimation of near-surface ice are responsible for the formation and the development of most of the ice-related landforms documented in Hellas. The relatively pristine form, lack of superposed craters, and strong association with the banded terrain, suggest an Amazonian (<3 Ga) age of formation for these landforms. Finally, relatively high surface pressures (above the triple point of water) expected in Hellas and summer-time temperatures often exceeding the melting point of water ice suggest that the basin may have recorded relatively “temperate” climatic conditions compared to other places on Mars. Therefore, the potentially ice-rich banded terrain may have deformed with lower viscosity and stresses compared to other locations on Mars, which may account for its unique morphology.

Radiometric 81Kr dating identifies 120,000-year-old ice at Taylor Glacier, Antarctica

We present successful 81 Kr-Kr radiometric dating of ancient polarice. Krypton was extracted from the air bubbles in four∼350-kg polar ice samples from Taylor Glacier in the McMurdo Dry Valleys, Antarctica, and dated using Atom Trap Trace Analysis (ATTA). The 81 Kr radiometric ages agree with independent age estimates obtained from stratigraphic dating techniques with a mean abso-lute age offset of 6±2.5 ka. Our experimental methods and sampling strategy are validated by (i) 85 Kr and 39 Ar analyses that show the samples to be free of modern air contamination and (ii)air content measurements that show the ice did not experience gas loss. We estimate the error in the 81 Kr ages due to past geomagnetic variability to be below 3 ka. We show that ice from the previous interglacial period (Marine Isotope Stage 5e, 130–115 ka before present) can be found in abundance near the surface of Taylor Glacier. Our study paves the way for reliable radiometric dating of ancient ice in blue ice areas and margin sites where large samples are available, greatly enhancing their scientific value as archives of old ice and meteorites. At present, ATTA 81Kr analysis requires a 40–80-kg ice sample; as sample requirements continue to decrease, 81 Kr dating of ice cores is a future possibility.

Lake ice phenology from AVHRR data for European lakes: An automated two-step extraction method

Several lake ice phenology studies from satellite data have been undertaken. However, the availability of long-term lake freeze-thaw-cycles, required to understand this proxy for climate variability and change, is scarce for European lakes. Long time series from space observations are limited to few satellite sensors. Data of the Advanced Very High Resolution Radiometer (AVHRR) are used in account of their unique potential as they offer each day global coverage from the early 1980s expectedly until 2022. An automatic two-step extraction was developed, which makes use of near-infrared reflectance values and thermal infrared derived lake surface water temperatures to extract lake ice phenology dates. In contrast to other studies utilizing thermal infrared, the thresholds are derived from the data itself, making it unnecessary to define arbitrary or lake specific thresholds. Two lakes in the Baltic region and a steppe lake on the Austrian–Hungarian border were selected. The later one was used to test the applicability of the approach to another climatic region for the time period 1990 to 2012. A comparison of the extracted event dates with in situ data provided good agreements of about 10 d mean absolute error. The two-step extraction was found to be applicable for European lakes in different climate regions and could fill existing data gaps in future applications. The extension of the time series to the full AVHRR record length (early 1980 until today) with adequate length for trend estimations would be of interest to assess climate variability and change. Furthermore, the two-step extraction itself is not sensor-specific and could be applied to other sensors with equivalent near- and thermal infrared spectral bands.

Evidence for biological shaping of hair ice

An unusual ice type, called hair ice, grows on the surface of dead wood of broad-leaf trees at temperatures slightly below 0 °C. We describe this phenomenon and present physical, chemical, and biological investigations to gain insight in the properties and processes related to hair ice. Tests revealed that the biological activity of a winter-active fungus is required in the wood for enabling the growth of hair ice. We confirmed the fungus hypothesis originally suggested by Wegener (1918) by reproducing hair ice on wood samples. Treatment by heat and fungicide suppresses the formation of hair ice. Fruiting bodies of Asco- and Basidiomycota are identified on hair-ice-carrying wood. One species, Exidiopsis effusa (Ee), was present on all investigated samples. Both hair-ice-producing wood samples and those with killed fungus show essentially the same temperature variation, indicating that the heat produced by fungal metabolism is very small, that the freezing rate is not influenced by the fungus activity, and that ice segregation is the common mechanism of ice growth on the wood surface. The fungus plays the role of shaping the ice hairs and preventing them from recrystallisation. Melted hair ice indicates the presence of organic matter. Chemical analyses show a complex mixture of several thousand CHO(N,S) compounds similar to fulvic acids in dissolved organic matter (DOM). The evaluation reveals decomposed lignin as being the main constituent. Further work is needed to clarify its role in hair-ice growth and to identify the recrystallisation inhibitor.

Sublimation of water ice mixed with silicates and tholins: Evolution of surface texture and reflectance spectra, with implications for comets

The surfaces of many objects in the Solar System comprise substantial quantities of water ice sometimes mixed with minerals and/or organic molecules. The sublimation of the ice changes the structural and optical properties of these objects. We present laboratory data on the evolution of the structure and the visible and near-infrared spectral reflectance of icy surface analogues of cometary ices, made of water ice, complex organic matter (tholins) and silicates, as they undergo sublimation under low temperature (<-70°C) and pressure (10-⁵mbar) conditions inside the SCITEAS simulation chamber. As the water ice sublimated, we observed in situ the formation of a porous sublimation lag deposit, or sublimation mantle, at the top of the ice. This mantle is a network of filaments made of the non-volatile particles. Organics or phyllosilicates grains, able to interact via stronger inter-particulate forces than olivine grains, can form a foam-like structure having internal cohesiveness, holding olivine grains together. As this mantle builds-up, the band depths of the sub-surface water ice are attenuated until complete extinction under only few millimeters of mantle. Optically thick sublimation mantles are mainly featureless in the near infrared. The absorption bands of the minerals present in the mantle are weak, or even totally absent if minerals are mixed with organics which largely dominate the VIS–NIR reflectance spectrum. During sublimation, ejections of large fragments of mantle, triggered by the gas flow, expose ice particles to the surface. The contrast of brightness between mantled and ice-exposed areas depends on the wavelength range and the dust/ice ratio considered. We describe how the chemical nature of the non-volatiles, the size of their particles, the way they are mixed with the ice and the dust/ice mass ratio influence the texture, activity and spectro-photometric properties of the sublimation mantles. These data provide useful references for interpreting remote-sensing observations of comets and also icy satellites or trans-neptunian objects.