The Ponto-Caspian basin as a final trap for southeastern Scandinavian Ice-Sheet meltwater

This paper provides new data on the evolution of the Caspian Sea and Black Sea from the Last Glacial Maximum until ca. 12 cal kyr BP. We present new analyses (clay mineralogy, grain-size, Nd isotopes and pollen) applied to sediments from the river terraces in the lower Volga, from the middle Caspian Sea and from the western part of the Black Sea. The results show that during the last deglaciation, the Ponto-Caspian basin collected meltwater and fine-grained sediment from the southern margin of the Scandinavian Ice Sheet (SIS) via the Dniepr and Volga Rivers. It induced the deposition of characteristic red-brownish/chocolate-coloured illite-rich sediments (Red Layers in the Black Sea and Chocolate Clays in the Caspian Sea) that originated from the Baltic Shield area according to Nd data. This general evolution, common to both seas was nevertheless differentiated over time due to the specificities of their catchment areas and due to the movement of the southern margin of the SIS. Our results indicate that in the eastern part of the East European Plain, the meltwater from the SIS margin supplied the Caspian Sea during the deglaciation until ∼13.8 cal kyr BP, and possibly from the LGM. That led to the Early Khvalynian transgressive stage(s) and Chocolate Clays deposition in the now-emerged northern flat part of the Caspian Sea (river terraces in the modern lower Volga) and in its middle basin. In the western part of the East European Plain, our results confirm the release of meltwater from the SIS margin into the Black Sea that occurred between 17.2 and 15.7 cal kyr BP, as previously proposed. Indeed, recent findings concerning the evolution of the southern margin of the SIS and the Black Sea, show that during the last deglaciation, occurred a westward release of meltwater into the North Atlantic (between ca. 20 and 16.7 cal kyr BP), and a southward one into the Black Sea (between 17.2 and 15.7 cal kyr BP). After the Red Layers/Chocolate Clays deposition in both seas and until 12 cal kyr BP, smectite became the dominant clay mineral. The East European Plain is clearly identified as the source for smectite in the Caspian Sea sediments. In the Black Sea, smectite originated either from the East European Plain or from the Danube River catchment. Previous studies consider smectite as being only of Anatolian origin. However, our results highlight both, the European source for smectite and the impact of this source on the depositional environment of the Black Sea during considered period.

Nd-Sr isotopes, diatom and biogenic opal content of Pliocene sediments from IODP Site 318-U1361

Warm intervals within the Pliocene epoch (5.33-2.58 million years ago) were characterized by global temperatures comparable to those predicted for the end of this century (Haywood and Valdes, doi:10.1016/S0012-821X(03)00685-X) and atmospheric CO2 concentrations similar to today (Seki et al., 2010, doi:10.1016/j.epsl.2010.01.037; Bartoli et al., 2011, doi:10.1029/2010PA002055; Pagani et al., 2010, doi:10.1038/ngeo724). Estimates for global sea level highstands during these times (Miller et al., 2012, doi:10.1130/G32869.1) imply possible retreat of the East Antarctic ice sheet, but ice-proximal evidence from the Antarctic margin is scarce. Here we present new data from Pliocene marine sediments recovered offshore of Adélie Land, East Antarctica, that reveal dynamic behaviour of the East Antarctic ice sheet in the vicinity of the low-lying Wilkes Subglacial Basin during times of past climatic warmth. Sedimentary sequences deposited between 5.3 and 3.3 million years ago indicate increases in Southern Ocean surface water productivity, associated with elevated circum-Antarctic temperatures. The geochemical provenance of detrital material deposited during these warm intervals suggests active erosion of continental bedrock from within the Wilkes Subglacial Basin, an area today buried beneath the East Antarctic ice sheet. We interpret this erosion to be associated with retreat of the ice sheet margin several hundreds of kilometres inland and conclude that the East Antarctic ice sheet was sensitive to climatic warmth during the Pliocene.

Comprehensive characterization of snowfall in terms of microphysical features, quantitative precipitation estimation and scavenging properties by in situ and remote sensing observations at an Antarctic coastal site

Snow plays a crucial role in the Earth's hydrological cycle and energy budget, making its monitoring necessary. In this context, ground-based radars and in situ instruments are essential thanks to their spatial coverage, resolution, and temporal sampling. Deep understanding and reliable measurements of snow properties are crucial over Antarctica to assess potential future changes of the surface mass balance (SMB) and define the contribution of the Antarctic ice sheet on sea-level rise. However, despite its key role, Antarctic precipitation is poorly investigated due to the continent's inaccessibility and extreme environment. In this framework, this Thesis aims to contribute to filling this gap by in-depth characterization of Antarctic precipitation at the Mario Zucchelli station from different points of view: microphysical features, quantitative precipitation estimation (QPE), vertical structure of precipitation, and scavenging properties. For this purpose, a K-band vertically pointing radar collocated with a laser disdrometer and an optical particle counter (OPC) were used. The radar probed the lowest atmospheric layers with high vertical resolution, allowing the first trusted measurement at only 105 m height. Disdrometer and OPC provided information on the particle size distribution and aerosol concentrations. An innovative snow classification methodology was designed by comparing the radar reflectivity (Ze) and disdrometer-derived reflectivity by means of DDA simulations. Results of classification were exploited in QPE through appropriate Ze-snow rate relationships. The accuracy of the resulting QPE was benchmarked against a collocated weighing gauge. Vertical radar profiles were also investigated to highlight hydrometeors' sublimation and growth processes. Finally, OPC and disdrometer data allowed providing the first-ever estimates of scavenging properties of Antarctic snowfall. Results presented in this Thesis give rise to advances in knowledge of the characteristics of snowfall in Antarctica, contributing to a better assessment of the SMB of the Antarctic ice sheet, the major player in the global sea-level rise.

Sedimentary cyclicity in CRP drillcore, Victoria Land Basin, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-1.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

Twentieth-century global-mean sea-level rise: is the whole greater than the sum of the parts?

Confidence in projections of global-mean sea level rise (GMSLR) depends on an ability to account for GMSLR during the twentieth century. There are contributions from ocean thermal expansion, mass loss from glaciers and ice sheets, groundwater extraction, and reservoir impoundment. Progress has been made toward solving the “enigma” of twentieth-century GMSLR, which is that the observed GMSLR has previously been found to exceed the sum of estimated contributions, especially for the earlier decades. The authors propose the following: thermal expansion simulated by climate models may previously have been underestimated because of their not including volcanic forcing in their control state; the rate of glacier mass loss was larger than previously estimated and was not smaller in the first half than in the second half of the century; the Greenland ice sheet could have made a positive contribution throughout the century; and groundwater depletion and reservoir impoundment, which are of opposite sign, may have been approximately equal in magnitude. It is possible to reconstruct the time series of GMSLR from the quantified contributions, apart from a constant residual term, which is small enough to be explained as a long-term contribution from the Antarctic ice sheet. The reconstructions account for the observation that the rate of GMSLR was not much larger during the last 50 years than during the twentieth century as a whole, despite the increasing anthropogenic forcing. Semiempirical methods for projecting GMSLR depend on the existence of a relationship between global climate change and the rate of GMSLR, but the implication of the authors' closure of the budget is that such a relationship is weak or absent during the twentieth century.

Simulations of glacial/interglacial cycles with simple box-models. The ocean CO2 source during deglaciations

Máster en Oceanografía

Sea Level Variations from Decades to Millennia: A few Case Studies

The main goals of this Ph.D. study are to investigate the regional and global geophysical components related to present polar ice melting and to provide independent cross validation checks of GIA models using both geophysical data detected by satellite mission, and geological observations from far field sites, in order to determine a lower and upper bound of uncertainty of GIA effect. The subject of this Thesis is the sea level change from decades to millennia scale. Within ice2sea collaboration, we developed a Fortran numerical code to analyze the local short-term sea level change and vertical deformation resulting from the loss of ice mass. This method is used to investigate polar regions: Greenland and Antarctica. We have used mass balance based on ICESat data for Greenland ice sheet and a plausible mass balance for Antarctic ice sheet. We have determined the regional and global fingerprint of sea level variations, vertical deformations of the solid surface of the Earth and variations of shape of the geoid for each ice source mentioned above. The coastal areas are affected by the long wavelength component of GIA process. Hence understanding the response of the Earth to loading is crucial in various contexts. Based on the hypothesis that Earth mantle materials obey to a linear rheology, and that the physical parameters of this rheology can be only characterized by their depth dependence, we investigate the Glacial Isostatic Effect upon the far field sites of Mediterranean area using an improved SELEN program. We presented new and revised observations for archaeological fish tanks located along the Tyrrhenian and Adriatic coast of Italy and new RSL for the SE Tunisia. Spatial and temporal variations of the Holocene sea levels studied in central Italy and Tunisia, provided important constraints on the melting history of the major ice sheets.

Paleoglaciologys Grand Unsolved Problem

The paleoglaciological concept that during the Pleistocene glacial hemi-cycles a super-large, structurally complex ice sheet developed in the Arctic and behaved as a single dynamic system, as the Antarctic ice sheet does today, has not yet been subjected to concerted studies designed to test the predictions of this concept. Yet, it may hold the keys to solutions of major problems of paleoglaciology, to understanding climate and sea-level changes. The Russian Arctic is the least-known region exposed to paleoglaciation by a hypothetical Arctic ice sheet but now it is more open to testing the concept. Implementation of these tests is a challenging task, as the region is extensive and the available data are controversial. Well-planned and coordinated field projects are needed today, as well as broad discussion of the known evidence, existing interpretations and new field results. Here we present the known evidence for paleoglaciation of the Russian Arctic continental shelf and reconstruct possible marine ice sheets that could have produced that evidence.

Basal Melting Along the Floating Part of Byrd Glacier

A mass balance calculation was made for the floating part of Byrd Glacier, using 1978-79 ice elevation and velocity data, over the 45 km of Byrd Glacier from its grounding line to where it leaves its fjord and merges with the Ross Ice Shelf. Smoothed basal melting rates were relatively uniform over this distance and averaged 12 +/- 2 m yr(-1).

Monsoonal Circulation of the McMurdo Dry Valleys, Ross Sea Region, Antarctica: Signal from the Snow Chemistry

McMurdo Dry Valleys (MDV, Ross Sea region, Antarctica) precipitation exhibits extreme seasonality in ion concentration, 3-5 orders of magnitude between summer and winter precipitation. To identify aerosol sources and investigate causes for the observed amplitude in concentration variability, four snow pits were sampled along a coast-Polar Plateau transect across the MDV. The elevation of the sites ranges from 50 to 2400 m and the distance from the coast from 8 to 93 km. Average chemistry gradients along the transect indicate that most species have either a predominant marine or terrestrial source in the MDV. Empirical orthogonal function analysis on the snow-chemistry time series shows that at least 57% of aerosol deposition occurs concurrently. A conceptual climate model, based on meteorological observations, is used to explain the strong seasonality in the MDV. Our results suggest that radiative forcing of the ice-free valleys creates a surface low-pressure cell during summer which promotes air-mass flow from the Ross Sea. The associated precipitating air mass is relatively warm, humid and contains a high concentration of aerosols. During winter, the MDV are dominated by air masses draining off the East Antarctic ice sheet, that are characterized by cold, dry and low concentrations of aerosols. The strong differences between these two air-mass sources create in the MDV a polar version of the monsoonal flow, with humid, warm summers and dry, cold winters.

Authigenic uranium, mass accumulation rates, benthic d13C and SST during the last interglacial period of ODP Site 177-1094

Southern Ocean sediments reveal a spike in authigenic uranium 127,000 years ago, within the last interglacial, reflecting decreased oxygenation of deep water by Antarctic Bottom Water (AABW). Unlike ice age reductions in AABW, the interglacial stagnation event appears decoupled from open ocean conditions and may have resulted from coastal freshening due to mass loss from the Antarctic ice sheet. AABW reduction coincided with increased North Atlantic Deep Water (NADW) formation, and the subsequent reinvigoration in AABW coincided with reduced NADW formation. Thus, alternation of deep water formation between the Antarctic and the North Atlantic, believed to characterize ice ages, apparently also occurs in warm climates.

Sedimentology and Lead and Strontium isotope ratios of Miocene sediments from the Kerguelen Plateau

Characterization of sediment from Ocean Drilling Program Site 745, representing the East Kerguelen Ridge sediment drift, addresses important issues surrounding the timing of Miocene to present East Antarctic ice sheet stability and oceanic environmental change. Our results show three periods of greatly enhanced accumulation of Antarctic-derived sediment, at 6.4-5.9 Ma, 4.9-4.4 Ma and 1.1-0.8 Ma, potentially indicative of warmer, less stable ice sheets at these times. Conversely, the accumulation of Antarctic-derived material is comparatively less during the middle of the Pliocene warm epoch (4.8-3.2 Ma). The deep flow forming the Kerguelen drift was stronger during the latest Miocene and earliest Pliocene and has decreased in intensity continuously since then.

Foraminiferal study of ODP Hole 189-1171

Sequence boundary ages determined in shallow-water sediments obtained from ODP (Ocean Drilling Program) Leg 189 Site 1171 (South Tasman Rise) compare well with other stratigraphic records (New Jersey, United States, and northwestern Europe) and d18O increases from deep-sea records, indicating that significant (>10 m) eustatic changes occurred during the early to middle Eocene (51-42 Ma). Sequence boundaries were identified and dated using lithology, bio- and magnetostratigraphy, water-depth changes, CaCO3 content, and physical properties (e.g., photospectrometry). They are characterized by a sharp bioturbated surface, low CaCO3 content, and an abrupt increase in glauconite above the surface. Foraminiferal biofacies and planktonic/benthic foraminiferal ratios were used to estimate water-depth changes. Ages of six sequence boundaries (50.9, 49.2, 48.5-47.8, 47.1, 44.5, and 42.6 Ma) from Site 1171 correlate well to the timings of d18O increases and sequence boundaries identified from other Eocene studies. The synchronous nature of sequence boundary development from globally distal sites and d18O increases indicates a global control and that glacioeustasy was operating in this supposedly ice-free world. This is supported by previous modeling studies and atmospheric pCO2 estimates showing that the first time pCO2 levels decreased below a threshold that would support the development of an Antarctic ice sheet occurred at ca. 51 Ma. Estimates of sea-level amplitudes range from ~20 m for the early Eocene (51-49 Ma) and ~25 m to ~45 m for the middle Eocene (48-42 Ma) using constraints established for Oligocene d18O records.

Stable isotopes and calcium carbonate at DSDP Holes 85-574 and 85-574A

Detailed stable isotopic and calcium carbonate records (with a sampling resolution of 3000 yr.) from the middle Miocene section of hydraulic piston corer (HPC) Hole 574A provide a sequence that records the major shift in the oxygen isotopic composition of the world's oceans that occurred at about 14 Ma. The data suggest that this transition was rapid and spans about 30,000 yr. of sediment deposition. In intervals before and after the shift, the mean d18O values are characterized by a constant mean with a high degree of variability. The degree of variability in both the d18O and d13C records is comparable to that observed for the Pliocene and earliest Pleistocene and does not show a significant change before or after the major shift in the d18O record. Whereas the oxygen isotopic record is characterized by relatively stable mean values before and after the middle Miocene event, the d13C record shows a number of significant offsets in the mean value separated by intervals of high-frequency variations. Time and frequency domain analysis of all records from Hole 574A indicate that the frequency components shown to be related to orbital changes in the Pleistocene record are also present in the middle Miocene. The high variability observed in the Site 574 isotopic records places important constraints on models describing the role of formation of the Antarctic ice sheet during the middle Miocene climatic transitions. Thus, HPC Hole 574A provides a valuable sequence for detailed study of climatic variability during an important time in the Earth's history, although we cannot provide a definitive explanation of the major oxygen isotopic event of the middle Miocene.