Late Quaternary record of sea-level changes in the Antarctic

The Late Quaternary sediment sequence of the continental margin in the eastern Weddell Sea is well suited for palaeoenvironmental reconstructions. Two cores from the upper slope, which contain the sedimentary record of the last 300 ky, have been sedimentologically investigated. Age models are based on lithostratigraphy and are correlated with the stable isotope record. As a result of a detailed analysis of the clay mineral composition, grain size distributions and structures, this sedimentary record provides the first marine evidence that the Antarctic ice sheet extended to the shelf edge during the last glacial.The variations in volume and size of the ice sheet were also simulated in numerical models. Changes in accumulation rate and ice temperature are of some importance, but the model revealed that fluctuations are primarily driven by changes in eustatic sea-level and that the ice edge extended to the shelf edge during the last glacial maximum. This causal relationship implies that the maximum ice extension strongly depends on the magnitude and duration of the sea-level depression during a glacial period. The results of the sedimentological investigations and of the numerical models show that the Antarctic ice sheet follows glacial events in the northern hemisphere by teleconnections of sea level. Correspondence to: H. Grobe     

A 3-D model for the Antarctic ice sheet: a sensitivity study on the glacial-interglacial contrast

On the longer climatic time scales, changes in the elevation and extent of the Antarctic ice sheet have an important role in modulating global atmospheric and oceanographic processes, and contribute significantly to world-wide sea levels. In this paper, a 3-D time-dependent thermomechanical model for the entire ice sheet is presented, that is subsequently used to examine the effects of glacial-interglacial shifts in environmental boundary conditions on its geometry. The model takes into account a coupled ice shelf, grounding-line dynamics, basal sliding and isostatic bed adjustment and considers the fully coupled velocity and temperature fields. Ice flow is calculated on a fine mesh (40 km horizontal grid size and 10 layers in the vertical) for grounded and floating ice and a stress transition zone in between at the grounding line, where all stress components contribute in the effective stress in the flow law. There is free interaction between ice sheet and ice shelf, so that the entire geometry is internally generated. A simulation of the present ice sheet reveals that the model is able to yield realistic results. A series of sensitivity experiments are then performed, in which lower temperatures, reduced accumulation rates and lower global sea level stands are imposed, either singly or in combination. By comparing results of pairs of experiments, the effects of each of these environmental changes can be determined. In agreement with glacial-geological evidence, we found that the most pronounced changes show up in the West Antarctic ice sheet configuration. They appear to be essentially controlled by variations in eustatic sea level, whereas typical glacial-interglacial changes in temperature and ice deposition rates tend to balance one another. These findings support the hypothesis that the Antarctic ice sheet basically follows glacial episodes in the northern hemisphere by means of sea-level teleconnections. Grounding occurs more readily in the Weddell sea than in the Ross sea and long time scales appear to be involved: it may take up to 30–40000 years for these continental shelf areas to become completely grounded after an initial stepwise perturbation in boundary conditions. According to these reconstructions, a steady state Antarctic ice sheet may contribute some 16 m to global sea level lowering at maximum glaciation.     

Calibrating a glaciological model of the Greenland ice sheet from the Last Glacial Maximum to present-day using field observations of relative sea level and ice extent

We constrain a three-dimensional thermomechanical model of Greenland ice sheet (GrIS) evolution from the Last Glacial Maximum (LGM, 21 ka BP) to the present-day using, primarily, observations of relative sea level (RSL) as well as field data on past ice extent. Our new model (Huy2) fits a majority of the observations and is characterised by a number of key features: (i) the ice sheet had an excess volume (relative to present) of 4.1 m ice-equivalent sea level at the LGM, which increased to reach a maximum value of 4.6 m at 16.5 ka BP; (ii) retreat from the continental shelf was not continuous around the entire margin, as there was a Younger Dryas readvance in some areas. The final episode of marine retreat was rapid and relatively late (c. 12 ka BP), leaving the ice sheet land based by 10 ka BP; (iii) in response to the Holocene Thermal Maximum (HTM) the ice margin retreated behind its present-day position by up to 80 km in the southwest, 20 km in the south and 80 km in a small area of the northeast. As a result of this retreat the modelled ice sheet reaches a minimum extent between 5 and 4 ka BP, which corresponds to a deficit volume (relative to present) of 0.17 m ice-equivalent sea level. Our results suggest that remaining discrepancies between the model and the observations are likely associated with non-Greenland ice load, differences between modelled and observed present-day ice elevation around the margin, lateral variations in Earth structure and/or the pattern of ice margin retreat.     

The Greenland ice sheet and greenhouse warming

Increased melting on glaciers and ice sheets and rising sea level are often mentioned as important aspects of the anticipated greenhouse warming of the earth's atmosphere. This paper deals with the sensitivity of Greenland's ice mass budget and presents a tentative projection of the Greenland component of future sea level rise for the next few hundred years. To do this, the ‘Villach II temperature scenario’ is prescribed,output from a comprehensive mass balance model is used to drive a high-resolution 3-D thermomechanic model of the ice sheet.The mass balance model consists of two parts: the accumulation part is based on presently observed values and is forced by changes in mean anr tempeerature. The ablation model is based on the degree-day method and accounts for daily and annual temperature cycle, a different degree-day factor for ice and snow melting and superimposed ice formation. Under present-day climatic conditions, the following total mass balance results (in ice equivalent per years): 599.3 × 10⁹ m³ of accumulation, 281.7 × 10⁹m³ of runoff assuming a balanced budget, 317.6 × 10⁹m³ of iceberg calving. A 1K uniform warming is then calculated to increase the runoff by 119.5 × 10⁹ m³. Since accumulation also increases by 32 × 10⁹ m³, this leads to reduction of the total mass balance by 887.5 × 10⁹ m³ of ice, corresponding to a sea level rise of 0.22 mm/yr. For temperature increase larger than 2.7 K, runoff, exceeds accumulation, and if ice sheet dynamics were to remain unchanged, this would add an extra amount of 0.8 mmyr to the worl's oceans.Imposing the Villach II scenario (warming up to 4.23 K) and accumulating mass balance changes forward in time (static response) would then result in a global sea level rise of 7.1 cm by 2100 AD, but this figure may go up to as much as 40 cm per century in case the warming is doubled. In a subsequent dynamic model involving the ice flow, the ice sheet is found to produce a counteracting effect by dynamically producing steeper slopes at the margin, thereby reducing the area over which runoff can take place. This effect is particularly apparent in the northeastern part of the ice sheet, and is also more pronounced for the smaller temperature perturbations. Nevertheless, all these experiments certainly highlight the vulnerability of the Greenland ice sheet with respect to a climatic warming.     

Predicted present-day evolution patterns of ice thickness and bedrock elevation over Greenland and Antarctica

This paper discusses predicted evolution patterns of present-day changes of ice thickness, surface elevation, and bedrock elevation over the Greenland and Antarctic continents. These were obtained from calculations with dynamic 3-D ice sheet models which were coupled to a visco-elastic solid Earth model. The experiments were initialized over the last two glacial cycles and subsequently averaged over the last 200 years to obtain the current evolution. The calculations indicate that the Antarctic Ice Sheet is still adjusting to the last glacial-interglacial transition yielding a decreasing ice volume and a rising bedrock elevation of the order of several centimetres per year. The Greenland Ice Sheet was found to be close to a stationary state with a mean thickness change of only a few millimetres per year, but the calculations revealed large spatial differences. Predicted patterns over Greenland are characterized by a small thickening over the ice sheet interior and a general thinning of the ablation area. In Antarctica, almost all of the predicted changes are concentrated in the West Antarctic Ice Sheet, which is still retreating at both the Weddell and Ross Sea margins. Over most of both ice sheets, the model indicates that the surface elevation trend is dominated by ice thickness changes rather than by bedrock elevation changes.     

Contrasting isotopic signatures between anthropogenic and geogenic Zn and evidence for post-depositional fractionation processes in smelter-impacted soils from Northern France

Zinc isotopes have been studied along two smelter-impacted soil profiles sampled near one of the largest Pb and Zn processing plants in Europe located in northern France, about 50 km south of Lille. Mean δ⁶⁶Zn values along these two soil profiles range from +0.22 ± 0.17‰ (2σ) to +0.34 ± 0.17‰ (2σ) at the lowest horizons and from +0.38 ± 0.45‰ (2σ) to +0.76 ± 0.14‰ (2σ) near the surface. The δ⁶⁶Zn values in the lowest horizons of the soils are interpreted as being representative of the local geochemical background (mean value +0.31 ± 0.38‰), whereas heavier δ⁶⁶Zn values near the surface of the two soils are related to anthropogenic Zn. This anthropogenic Zn occurs in the form of franklinite (ZnFe₂O₄)-bearing slag grains originating from processing wastes at the smelter site and exhibiting δ⁶⁶Zn values of +0.81 ± 0.20‰ (2σ). The presence of franklinite is indicated by EXAFS analysis of the topsoil samples from both soil profiles as well as by micro-XANES analysis of the surface horizon of a third smelter-impacted soil from a distant site. These results indicate that naturally occurring Zn and smelter-derived Zn exhibit significantly different δ⁶⁶Zn values, which suggests that zinc isotopes can be used to distinguish between geogenic and anthropogenic sources of Zn in smelter-impacted soils. In addition to a possible influence of additional past sources of light Zn (likely Zn-sulfides and Zn-sulfates directly emitted by the smelter), the light δ⁶⁶Zn values in the surface horizons compared to smelter-derived slag materials are interpreted as resulting mainly from fractionation processes associated with biotic and/or abiotic pedological processes (Zn-bearing mineral precipitation, Zn complexation by organic matter, and plant uptake of Zn). This conclusion emphasizes the need for additional Zn isotopic studies before being able to use Zn isotopes to trace sources and pathways of this element in surface environments.     

Two-dimensional modeling of fine sediment transport with mixed sediment and consolidation: Application to the Gironde Estuary,France

In order to optimize ship navigation in the macrotidal Gironde Estuary, a recent project funded by the port of Bordeaux aims at better understand and forecast hydrodynamic and fine sediment transport within the estuary. In the framework of this project, a two-dimensional hydro-sedimentary model is built. The model includes hydrodynamic forcings, mixed-sediment transport, and consolidation processes. The harmonic analysis of the astronomical tides reveals a strong distortion of the tidal wave inducing the growth of overtide constituents and the non-significant effect of tide-surge interactions in annual-scale prediction. Depending on hydrological conditions, river discharge can considerably alter the model accuracy due to the migration of the turbidity maximum zone modifying the bottom roughness. Comparison with measurements shows the ability of the model to reproduce suspended-sediment concentrations in the central Estuary. Sensitivity of the model to sediment features has also been discussed in regard of suspended-sediment concentrations and fluid mud deposits. The model will be further coupled with ship squat predictions and a morphodynamic model.     

The Greenland ice sheet through the last glacial-interglacial cycle

The evolution of the Greenland ice sheet during the last 150,000 years, in response to a climate history derived from a Greenland ice-margin oxygen-18 record, is simulated by means of a three-dimensional, time-dependent ice-sheet model. The calculations indicate that the ice sheet displayed considerable thinning and ice-margin retreat during the last interglacial (the Eemian) and during a warm interstadial c. 100,000 yr B.P., resulting in a splitting up of the ice sheet into a central-northern and a southern part. However, the ice sheet in Central Greenland survived the warm stages with almost unchanged surface elevations as compared with the present.     

Effect of self-weight consolidation on a hydro-sedimentological model for the Río de la Plata estuary

The effect of the consolidation process on the morphodynamics and fine sediment dynamics of the Río de la Plata estuary is explored through a circulation-wave-sediment transport model. The consolidation model is calibrated based on settling column experimental data. Different simulations are done in order to initialize the mud layer distribution and to investigate the impact of different erosion parameter assumptions on the modeled sediment dynamics. Finally a two-year simulation is done with and without the consolidation process and realistic hydrodynamic forcings. Considering the consolidation process, the model correctly reproduces measured vertical density profiles in the Montevideo Bay access channel. The simulated suspended sediment dynamics behavior in Montevideo Bay with the consolidation process provides a more realistic deposition pattern in regard to the dredging activities.