A case study of land cover change (1950-2003) and runoff in a Mediterranean catchment

Mediterranean environments have been subject to major land cover change since the end of the second world war. Housing, agricultural activities, forests, green spaces and other land uses have shifted due to urbanisation and tourism. These changes influence runoff, and municipal authorities often cannot estimate the net impact of complex land cover transitions. During this period, elected representatives have become increasingly sensitive to the risks of flooding and have implemented a number of channel management strategies. The main objective of this case study was to analyse the impact of land cover change on total storm runoff between 1950 and 2003 in a Mediterranean catchment near St Tropez, France. A secondary objective was to compare these changes to the impacts of channel management on bankfull discharge. Aerial photographs were used to classify land cover in 3 urban categories, vineyards and bare soil, forests, and green spaces. Stream discharge was estimated using a distributed event based total runoff approach. After validating the model for a large winter event (114 mm) for 1982, runoff was calculated for the same event for 1950 and 2003. Land cover changes occurred mainly in the alluvial plain area. Total gauge catchment urban area increased from 30.1 ha to 393.8 between 1950 and 2003 at the expense mainly of agricultural land, but this was compensated in part by an increase in grassed area. Some of the loss in vineyards was replaced by clearing forested land on the first hills close to the plain. Bank stabilisation and channel maintenance since the 1980’s reduced surface roughness and increased channel area, thereby greatly increasing bankfull discharge. While the impact of urbanisation on runoff was small, channel management effects increased bankfull discharge substantially. Flood damage from extreme events was not studied here.     

Exploration of the Outer Solar System by Stellar Occultations

Stellar occultations are a powerful method for exploring the outer solar system, where faintness and small angular diameters prevent us from exploring in details objects like satellites, rings, or Kuiper Belt Objects. Unique kilometric spatial resolutions or better can be reached through that method. Occultations usually observe identified objects whose trajectory is known, though the occultation events might be difficult to predict. It is also possible to explore populations of small objects populations whose density in the sky plane is large enough to search for serendipitous occultations. Various instrumental methods exist for both predicted and serendipitous occultation, both needing fast photometric recordings of target stars.     

The fresnel interferometric imager

The Fresnel Interferometric Imager has been proposed to the European Space Agency (ESA) Cosmic Vision plan as a class L mission. This mission addresses several themes of the CV Plan: Exoplanet study, Matter in extreme conditions, and The Universe taking shape. This paper is an abridged version of the original ESA proposal. We have removed most of the technical and financial issues, to concentrate on the instrumental design and astrophysical missions. The instrument proposed is an ultra-lightweight telescope, featuring a novel optical concept based on diffraction focussing. It yields high dynamic range images, while releasing constraints on positioning and manufacturing of the main optical elements. This concept should open the way to very large apertures in space. In this two spacecraft formation-flying instrument, one spacecraft holds the focussing element: the Fresnel interferometric array; the other spacecraft holds the field optics, focal instrumentation, and detectors. The Fresnel array proposed here is a 3.6 ×3.6 m square opaque foil punched with 10⁵ to 10⁶ void “subapertures”. Focusing is achieved with no other optical element: the shape and positioning of the subapertures (holes in the foil) is responsible for beam combining by diffraction, and 5% to 10% of the total incident light ends up into a sharp focus. The consequence of this high number of subapertures is high dynamic range images. In addition, as it uses only a combination of vacuum and opaque material, this focussing method is potentially efficient over a very broad wavelength domain. The focal length of such diffractive focussing devices is wavelength dependent. However, this can be corrected. We have tested optically the efficiency of the chromatism correction on artificial sources (500 < λ < 750 nm): the images are diffraction limited, and the dynamic range measured on an artificial double source reaches 6.2 10^− 6. We have also validated numerical simulation algorithms for larger Fresnel interferometric arrays. These simulations yield a dynamic range (rejection factor) close to 10^− 8 for arrays such as the 3.6 m one we propose. A dynamic range of 10^− 8 allows detection of objects at contrasts as high as than 10^− 9 in most of the field. The astrophysical applications cover many objects in the IR, visible an UV domains. Examples are presented, taking advantage of the high angular resolution and dynamic range capabilities of this concept.     

“Galileo Galilei” (GG) a small satellite to test the equivalence principle of Galileo, Newton and Einstein

“Galileo Galilei” (GG) is a small satellite designed to fly in low Earth orbit with the goal of testing the Equivalence Principle—which is at the basis of the General Theory of Relativity—to 1 part in 10¹⁷. If successful, it would improve current laboratory results by 4 orders of magnitude. A confirmation would strongly constrain theories; proof of violation is believed to lead to a scientific revolution. The experiment design allows it to be carried out at ambient temperature inside a small 1-axis stabilized satellite (250 kg total mass). GG is under investigation at Phase A-2 level by ASI (Agenzia Spaziale Italiana) at Thales Alenia Space in Torino, while a laboratory prototype (known as GGG) is operational at INFN laboratories in Pisa, supported by INFN (Istituto Nazionale di fisica Nucleare) and ASI. A final study report will be published in 2009.     

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.     

Pliocene lahar deposits in the Coastal Cordillera of central Chile: Implications for uplift, avalanche deposits, and porphyry copper systems in the Main Andean Cordillera

Lahar deposits occur within a shallow marine sedimentary succession of the Pliocene La Cueva Formation in the Coastal Cordillera of central Chile (33°40′–34°15′S). Provenance studies of the abundant volcanic material in the lahar deposits suggest that they derive from denudation by mass wasting of Oligocene–Miocene volcanic rocks on the western slopes of the Main Andean Cordillera at the same latitude. Pliocene rock debris deposits preserved in the region of El Teniente (34°S) and scattered along the westernmost part of the Andes of central Chile indicate catastrophic erosive events related to the rapid uplift of the cordilleran block. This rock debris was deposited by avalanches and transformed further downslope into lahars by dilution with stream water. Lahars were channeled along the ancient drainage system that reached a shallow Pliocene sea at the site of the present Coastal Cordillera. The exceedingly rapid exhumation of active porphyry systems during the Early Pliocene in this part of the Andes may have played a role in affecting hydrothermal processes, brecciation, and diatreme formation at the porphyry systems of El Teniente and Río Blanco–Los Bronces.