Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

Geometry and kinematics of recent deformation in the Mondy–Tunka area (south-westernmost Baikal rift zone, Mongolia–Siberia)

We used satellite imagery and field data to investigate the south‐westernmost Baikal rift zone. We focus our study in the Mondy and Ikhe Ukhgun valleys, site of an Mw = 6.9 seismic event in 1950. Surface deformations are observed along the E–W‐trending Mondy strike‐slip fault and along the Ikhe Ukhgun thrust. The Mondy fault system is 80 km long and is composed of four segments 10–15 km long. These segments are characterized by subvertical planes with left‐lateral movements. The Ikhe Ukhgun thrust is 20 km long, dips 40° to the south and shows reverse movement with a left‐lateral component. These observations are consistent with the present‐day regional NNE–SSW compression and with the focal mechanism of the 1950 Mondy earthquake that was recently re‐evaluated. These features, like those observed in the Tunka basin, demonstrate a recent change of regional strain regime from transtension to transpression that we place before the Late Pleistocene.     

Sensitivity of a Greenland ice sheet model to ice flow and ablation parameters: consequences for the evolution through the last climatic cycle

Sensitivity experiments are conducted to test the influence of poorly known model parameters on the simulation of the Greenland ice sheet by means of a three dimensional numerical model including the mechanical and thermal processes within the ice. Two types of experiments are performed: steady-state climatic conditions and simulations over the last climatic cycle with a climatic forcing derived from the GRIP record. The experiments show that the maximum altitude of the ice sheet depends on the ice flow parameters (deformation and sliding law coefficients, geothermal flux) and that it is low when the ice flow is fast. On the other hand, the maximum altitude is not sensitive to the ablation strength and consequently during the climatic cycle it is driven by changes in accumulation rate. The ice sheet extension shows the opposite sensitivity: it is barely affected by ice flow velocity and the ice covered area is smaller for large ablation coefficients. For colder climates, when there is no ablation, the ice sheet extension depends on the sea level. An interesting result is that the variations with time of the altitude at the ice divide (Summit) do not depend on the parameters we tested. The present modelled ice sheets resulting from the climatic cycle experiments are compared with the present measured ice sheet in order to find the set of parameters that gives the best fit between modelled and measured geometry. It seems that, compared to the parameter set most commonly used, higher ablation rate coefficents must be used. Received: 19 September 1995 / Accepted: 30 May 1996     

Luminescence dating of a gigantic palaeolandslide in the Gobi-Altay mountains, Mongolia

In an attempt to date a palaeolandslide that took place along the Baga Bogd Massif, in Mongolia, the infrared stimulated luminescence (IRSL) method has been applied to lacustrine silty sediments directly overlying the landslide mass. The IRSL age estimates obtained on alkali feldspar grains (>40 μm) and polymineral fine grains (4–11 μm) provide a minimum age for the landslide event. The IRSL ages on alkali feldspars corrected for long-term fading using the protocol of Mejdahl (1988, 1989) suggest that the palaeolandslide occurred at the beginning of the Last Interglacial. These are in good agreement with the ¹⁰Be cosmogenic dates obtained on faulted and abandoned alluvial fans in the Gobi-Altay mountains. This study demonstrates for the first time that the IRSL dating method can successfully be applied for establishing landslide chronologies.     

Origin and mechanism of Subantarctic Mode Water formation and transformation in the Southern Indian Ocean

The sources and pathways of mode waters and lower thermocline waters entering the subtropical gyre of the Indian Ocean are examined. A Lagrangian analysis is performed on an eddy-admitting simulation of the Global Ocean performed by the DRAKKAR Group (NEMO/OPA), which captures the main observed features. We trace the subducted mode water’s pathways, identify their formation regions and trace whether their source waters come from the Atlantic, Pacific or Indian sectors of the Southern Ocean. Three main sites for mode waters ventilation in the Indian sector are identified with different circulation pathways and source water masses: (a) just north of Kerguelen, where 4.2 Sv of lighter Subantarctic Mode Waters (SAMW); σ ₀ ∼ 26.5) are exported—originating in the Atlantic and Agulhas Retroflection regions; (b) SW of Australia, where 6.5 Sv of medium SAMW (σ ₀ ∼ 26.6) are ventilated—originating in the southern and denser Agulhas Retroflection region; (c) SW of Tasmania and along the South Australian coast, where 3 Sv of denser SAMW (σ ₀ ∼ 26.75) are ventilated—originating from three sources: Leeuwin Current waters, Tasman Sea (Pacific) waters and Antarctic Surface Waters. In all cases, modelled mode waters were last ventilated in the Indian Ocean just north of the deepest winter-mixed layers. For the waters subducted SW of Australia, the last ventilation site extends even further north. Waters ventilated in the deepest mixed layers north of the Subantarctic Front are then re-ventilated 5 years later southwest of Australia. The model results raise new hypotheses that revisit the classical picture of the SAMW formation and transformation, where a large homogeneous mixed layer is subducted and ‘slides’ equatorward, essentially maintaining the T/S characteristics acquired at the surface. Firstly, the last ventilation of the modelled mode waters is not in the region of the deepest mixed layers, as previously thought, but further north in regions of moderate meso-scale eddy activity. Secondly, the model shows for the first time a significant source region for Indian Ocean mode waters coming from deep winter-mixed layers along the south Australian coast. Finally, this analysis shows how the mode water characteristics are modified after subduction, due to internal eddy mixing. The simulation shows that resolved eddies have a strong impact on the mixed layer properties and that isopycnal eddy mixing also contributes to the generation of more homogeneous mode water characteristics in the interior.     

A Spatio‐Temporal Graph Model for Marine Dune Dynamics Analysis and Representation

Defining a model for the representation and the analysis of spatio‐temporal dynamics remains an open domain in geographical information sciences. In this article we investigate a spatio‐temporal graph‐based model dedicated to managing and extracting sets of geographical entities related in space and time. The approach is based on spatial and temporal local relations between neighboring entities during consecutive times. The model allows us to extract sets of connected entities distant in time and space over long periods and large spaces. From GIS concepts and qualitative reasoning on space and time, we combine the graph model with a dedicated spatial database. It includes information on geometry and geomorphometric parameters, and on spatial and temporal relations. This allows us to extend classical measurements of spatial parameters, with comparisons of entities linked by complex relations in space and time. As a case study, we show how the model suggests an efficient representation of dunes dynamics on a nautical chart for safe navigation.     

Physical processes contributing to the water mass transformation of the Indonesian Throughflow

The properties of the waters that move from the Pacific to the Indian Ocean via passages in the Indonesian archipelago are observed to vary with along-flow-path distance. We study an ocean model of the Indonesian Seas with reference to the observed water property distributions and diagnose the mechanisms and magnitude of the water mass transformations using a thermodynamical methodology. This model includes a key parameterization of mixing due to baroclinic tidal dissipation and simulates realistic water property distributions in all of the seas within the archipelago. A combination of air–sea forcing and mixing is found to significantly change the character of the Indonesian Throughflow (ITF). Around 6 Sv (approximately 1/3 the model net ITF transport) of the flow leaves the Indonesian Seas with reduced density. Mixing transforms both the intermediate depth waters (transforming 4.3 Sv to lighter density) and the surface waters (made denser despite the buoyancy input by air–sea exchange, net transformation?=?2 Sv). The intermediate transformation to lighter waters suggests that the Indonesian transformation contributes significantly to the upwelling of cold water in the global conveyor belt. The mixing induced by the wind is not driving the transformation. In contrast, the baroclinic tides have a major role in this transformation. In particular, they are the only source of energy acting on the thermocline and are responsible for creating the homostad thermocline water, a characteristic of the Indonesian outflow water. Furthermore, they cool the sea surface temperature by between 0.6 and 1.5°C, and thus allow the ocean to absorb more heat from the atmosphere. The additional heat imprints its characteristics into the thermocline. The Indonesian Seas cannot only be seen as a region of water mass transformation (in the sense of only transforming water masses in its interior) but also as a region of water mass formation (as it modifies the heat flux and induced more buoyancy flux). This analysis is complemented with a series of companion numerical experiments using different representations of the mixing and advection schemes. All the different schemes diagnose a lack of significant lateral mixing in the transformation.     

Origin of the absarokite–banakite association of the Damavand volcano (Iran): trace elements and Sr,Nd, Pb isotope constraints

The activity of the Damavand volcano (Central Alborz, northern Iran) began 1.8 Ma ago and continued up to 7 ka BP. Although the volcanic suite is clearly of shoshonitic affinity, only two petrographic types can be distinguished in the studied lavas: (1) weakly differentiated absarokites (49 < %SiO₂ < 51), scattered around the volcano but with a regional extension, (2) highly differentiated banakites (59 < %SiO₂ < 63), which form the bulk of the 4,000 m thick volcanic pile. All lavas are alkalic (3.7 < %K₂O < 5), REE and LILE-rich (e.g., 85 < La < 148 ppm; 9 < Th < 32 ppm) and show highly fractionated REE patterns (69 < La/Yb < 115) and pronounced Nb–Ta negative anomalies. The absarokites are characterised by Sr (0.7045–0.7046) and Nd (0.51266–0.51269) isotope compositions close to the Bulk Earth values, and distinct from those of the banakites (0.7047 < ⁸⁷Sr/⁸⁶Sr < 0.7049, 0.51258 < ¹⁴³Nd/¹⁴⁴Nd < 0.51262). The Pb isotope ratios are also slightly lower in the absarokites than in the banakites (18.71 < ²⁰⁶Pb/²⁰⁴Pb < 18.77, 15.62 < ²⁰⁷Pb/²⁰⁴Pb < 15.63, 38.85 < ²⁰⁸Pb/²⁰⁴Pb < 38.91, and 18.77 < ²⁰⁶Pb/²⁰⁴Pb < 18.84, 15.62 < ²⁰⁷Pb/²⁰⁴Pb < 15.64, 38.94 < ²⁰⁸Pb/²⁰⁴Pb < 39.06, respectively). Overall, there is a clear tendency towards higher Sr, Pb and lower Nd isotope ratios with increasing degree of differentiation. This study suggests that the absarokites result from a low degree of partial melting (∼5%) of a highly metasomatized mantle source, which inherited its characteristics from an old subduction setting. The initiation of volcanic activity 1.8 Ma ago results from variations in the lithospheric thermal regime, probably related to lithospheric delamination as proposed for Anatolia (Pearce et al. 1990). The banakites are mainly generated by extensive fractional crystallisation (∼70%) of the absarokitic magma, with a limited amount (a few percents) of assimilation of an old crustal component, in the form of bulk assimilation or AFC processes, which both can explain the Sr, Nd and Pb isotope data.     

Mapping Coastal Aquifers by Joint Inversion of DC and TEM Soundings-Three Case Histories

Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.