Groundwater flow modelling in the central zone of Hanoi, Vietnam

On the basis of a review of the Quaternary sedimentary architecture in the area of Hanoi city (Vietnam), a numerical, deterministic and three-dimensional groundwater flow model has been built for a simulation between 1995 and 2004. The sedimentary architecture has been constructed on the basis of the data from 32 drillings covering the entire Quaternary sequence (but with little sedimentological detail), as well as hydrographical and hydrogeological data. Both steady- and transient-state conditions were tested. Results calculated by the model seem to indicate that the conceptual hypotheses adopted are reasonable for the modelled domain and period. The simulation allows for calculation of the regional groundwater flow trends. It is also used for assessing the relative importance of the various recharge sources of the shallow aquifer system in Hanoi, and for estimating the interactions between groundwater and the Red River.     

Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.     

LA-MC-ICPMS and SHRIMP U-Pb dating of complex zircons from Quaternary tephras from the French Massif Central: Magma residence time and geochemical implications

Analyses of zircon grains from the Queureuilh Quaternary tephras (pumice) provide new information about their pre-eruptive history. U-Pb dating was performed in situ using two methods: SHRIMP and LA-MC-ICPMS equipped with a multi-ion counting system. Both methods provided reliable ²⁰⁷Pb/²⁰⁶Pb and ²⁰⁶Pb/²³⁸U ratios as well as U and Th abundances required for U-Pb Concordia intercept age determination, after initial ²³⁰Th disequilibrium correction. The new LA-MC-ICPMS method was validated by dating a reference zircon (61.308B) and zircons from a phonolitic lava dated independently with the two techniques. A time resolution of about 20 kyr for 1 Ma zircon crystals was achieved for both methods.The clear euhedral zircon population from Queureuilh tephras is quite complex from several points of view: (1) some grains are reddish or yellowish while others are colorless; (2) the U and Th composition changes by more than an order of magnitude and Th/U is generally high (∼1-2); (3) there are three discrete ages recorded at 2.35 ± 0.04, 1.017 ± 0.008 and 0.640 ± 0.010 Ma.From the previously determined ⁴⁰Ar/³⁹Ar age at 0.571 ± 0.060 Ma [Duffell H. (1999) Contribution géochronologique à la stratigraphie volcanique du Massif des Monts Dore par la méthode ⁴⁰Ar/³⁹Ar. D.E.A. Univ. Clermont-Ferrand, 56 p.], the discontinuous zircon age populations, the color of the grains and their composition, we favor the following model as explanation: The oldest, less numerous group of reddish zircons represents xenocrystic grains resulting from assimilation of the local material during magma ascent. A primitive magma chamber, perhaps deep in crustal level, was formed at 1.0 Ma. The related magma, previously characterized by high Th/U ratio (2.2 ± 1.1), underwent rejuvenation during ascent to a new chamber at shallow depth and/or during injection of more mafic magmas. During this stage, at 0.64 Ma, the colorless zircon grains of lower Th/U ratio (1.3 ± 0.5) crystallized. This last stage defined the magma residence time of 70 kyr prior to eruption dated by the ⁴⁰Ar/³⁹Ar method. However, if the primitive magma is considered, the magma residence time as a whole from this first stage reached 446 kyr.In the light of the complex history of such magmas, which commonly involves recycling of zircon grains that precipitated tens to hundreds of kyr earlier than eruptions, the use of Zr concentration in geochemical modeling of whole rock compositional data can be problematic.     

Archaeological reconstruction of medieval lead production: Implications for ancient metal provenance studies and paleopollution tracing by Pb isotopes

The identification of metal provenance is often based on chemical and Pb isotope analyses of materials from the operating chain, mainly ores and metallic artefacts. Such analyses, however, have their limits. Some studies are unable to trace metallic artefacts or ingots to their ore sources, even in well-constrained archaeological contexts. Possible reasons for this difficulty are to be found among a variety of limiting factors: (i) problems of ore signatures, (ii) mixing of different ores (alloys), (iii) the use of additives during the metallurgical process, (iv) metal recycling and (v) possible Pb isotopic fractionation during metal production. This paper focuses on the issue of Pb isotope fractionation during smelting to address the issue of metal provenance. Through an experimental reconstruction of argentiferous Pb production in the medieval period, an attempt was made to better understand and interpret the Pb isotopic composition of ore smelting products. It is shown that the measured differences (outside the total external uncertainties of 0.005 (2*sd) for ²⁰⁶Pb/²⁰⁴Pb ratios) in Pb signatures measured between ores, slag and smoke are not due to Pb mass fractionation processes, but to (1) ore heterogeneity (Δ²⁰⁶Pb/²⁰⁴Pb_slag-ores = 0.066) and (2) the use of additives during the metallurgical process (Δ²⁰⁶Pb/²⁰⁴Pb_slag-ores = 0.083). Even if these differences are due to causes (1) and/or (2), smoke from the ore reduction appears to reflect the ore mining area without a significant disturbance of its Pb signature for all the isotopic ratios (Δ²⁰⁶Pb/²⁰⁴Pb_smokes-ores = 0.026). Thus, because the isotopic heterogeneity of the mining district and additives is averaged in slags, slag appears as the most relevant product to identify ancient metal provenance. Whereas aiming at identifying a given mine seems beyond the possibilities provided by the method, searching for the mining district through analysis of the smelting workshop materials should provide a more appropriate approach in cases where no archaeological evidence of ancient mining is available. Furthermore, smoke Pb isotopic composition does not seem to be significantly affected by the metallurgical process. Paleopollution recorded in peat deposits could help to detect ancient mining production and workshops. Integrated collaboration between mining archaeologists and geochemists appears crucial to achieve this goal.     

Study of the initiation and propagation of excavation damaged zones around openings in argillaceous rock

The study of the creation and evolution of the excavation disturbed zone (EDZ) in argillaceous rocks is a major issue for the safety of nuclear wastes underground repositories. In this context, the argillaceous Tournemire site has provided a unique opportunity to study the evolution of the EDZ with time thanks to the existence of three openings of different ages. A thorough characterization of the EDZ has been conducted by different means such as visual observation, analysis of samples extracted from drilled boreholes, EDZ permeability measurements, etc. On the basis of these measurements, a conceptual model of the EDZ initiation and propagation at the Tournemire site has been proposed. In order to validate this model, numerical simulations of increasing complexity have been carried out. In a first attempt, the response of the rock mass to the excavation phase, followed by seasonal cyclic variations of temperature and relative humidity inside the opening, has been simulated by means of a purely mechanical analysis, using a simple elastic material model. The EDZ has been estimated by post-processing the calculated stress states, using a Mohr–Coulomb failure criterion. The results obtained show that no EDZ could be predicted unless adopting a low cohesion value for the rock mass. Moreover, the deferred nature of the EDZ formation in Tournemire could not be reproduced. These limitations have then been suppressed by using a coupled viscoplastic-damaging mechanical model, the parameters of which have been identified from different laboratory experiments. With this model, a time evolution of the EDZ could be predicted, but the EDZ pattern could not match the one observed in situ. Finally, in view of the importance of the hydraulic couplings, unsaturated hydro-mechanical calculations have been carried out to investigate the effect of the numerous seasonal variations cycles and the resulting shrinkage.     

Weak oceanic heat transport as a cause of the instability of glacial climates

The stability of the thermohaline circulation of modern and glacial climates is compared with the help of a two dimensional ocean—atmosphere—sea ice coupled model. It turns out to be more unstable as less freshwater forcing is required to induce a polar halocline catastrophy in glacial climates. The large insulation of the ocean by the extensive sea ice cover changes the temperature boundary condition and the deepwater formation regions moves much further South. The nature of the instability is of oceanic origin, identical to that found in ocean models under mixed boundary conditions. With similar strengths of the oceanic circulation and rates of deep water formation for warm and cold climates, the loss of stability of the cold climate is due to the weak thermal stratification caused by the cooling of surface waters, the deep water temperatures being regulated by the temperature of freezing. Weaker stratification with similar overturning leads to a weakening of the meridional oceanic heat transport which is the major negative feedback stabilizing the oceanic circulation. Within the unstable regime periodic millennial oscillations occur spontaneously. The climate oscillates between a strong convective thermally driven oceanic state and a weak one driven by large salinity gradients. Both states are unstable. The atmosphere of low thermal inertia is carried along by the oceanic overturning while the variation of sea ice is out of phase with the oceanic heat content. During the abrupt warming events that punctuate the course of a millennial oscillation, sea ice variations are shown respectively to damp (amplify) the amplitude of the oceanic (atmospheric) response. This sensitivity of the oceanic circulation to a reduced concentration of greenhouse gases and to freshwater forcing adds support to the hypothesis that the millennial oscillations of the last glacial period, the so called Dansgaard—Oeschger events, may be internal instabilities of the climate system.     

On the Impact of Anthropogenic Heat Fluxes on the Urban Boundary Layer: A Two-Dimensional Numerical Study

The heat generated in buildings and the manner in which this heat is exchanged with the ambient environment can play an important role in urban climate. Recent studies have shown that anthropogenic heat from air-conditioning facilities can increase the exterior ambient temperature and should be taken into account for a more complete urban heat island (UHI) mitigation study. For this purpose, the first part of the present work is focused on the coupling of a new building energy model (BEM) and an urban canopy parameterisation (UCP). The new scheme is implemented in a finite volume mesoscale model (MM) and tested in a two-dimensional (2D) configuration of a city over flat terrain. A sensitivity study is performed with respect to different parameters in order to test the simulation system and enhance the understanding of the possible impacts of the BEM on the exterior microclimate.     

Holocene salt‐marsh sedimentary infilling and relative sea‐level changes in West Brittany (France) using foraminifera‐based transfer functions

In order to reconstruct former sea‐levels and to better characterize the history of Holocene salt‐marsh sedimentary infillings in West Brittany (western France), local foraminifera‐based transfer functions were developed using weighted‐average‐partial‐least‐squares (WAPLS) regression, based on a modern data set of 26 and 51 surface samples obtained from salt‐marshes in the bay of Tresseny and the bay of Brest, respectively. Fifty cores were retrieved from Tresseny, Porzguen, Troaon and Arun salt‐marshes, which were litho‐ and biostratigraphically analysed in order to reconstruct palaeoenvironmental changes. A total of 26 AMS ¹⁴C age determinations were performed within the sediment successions. The Holocene evolution of salt‐marsh environments can be subdivided into four stages: (i) a development of brackish to freshwater marshes (from c. 6400 to 4500 cal. a BP); (ii) salt‐marsh formation behind gravel barriers in the bay of Brest (from 4500 to 2900 cal. a BP); (iii) salt‐marsh erosion and rapid changes of infilling dynamics due to the destruction of coastal barriers by storm events (c. 2900?2700 cal. a BP); (iv) renewed salt‐marsh deposition and small environmental changes (from 2700 cal. a BP to present). From the application of transfer functions to fossil assemblages, 14 new sea‐level index points were obtained, indicating a mean relative sea‐level rise of around 0.90±0.12 mm a^?1 since 6300 cal. a BP.