Radiocarbon dating of modern groundwater: the role of the unsaturated zone

Biological and physical processes occurring in soils may lead to significant isotopic changes between the isotopic compositions of atmospheric CO2 and of soil CO2. Also, during water and gas transport from the soil surface to the water table, isotopic changes likely occur due to numerous physical processes such as gas production and diffusion, water advection, and gas-water-rock interactions. In most cases, these changes are not included in the correction models developed for groundwater dating, whereas they can significantly impact the calculation of the 14C age. We explore the role of these processes using: i) experimental data from two aquifer sites (Fontainebleau sands and Astian sands, France), ii) a distributed model to simulate the 14C activities of soil CO2, and iii) numerical simulations in order to highlight the role of the physical processes.¦The 13C content in soil CO2 showed seasonal variations and highlighted the competition between CO2 production and CO2 diffusion. Their respective contributions played a significant role in defining the isotopic composition of CO2 at the water table. On both study sites, variations of the 14C activity in soil CO2 reflect the competition between the fluxes of root derived-CO2 and organic matter derived-CO2. Since the nuclear weapon tests in the fifties and sixties, soil CO2 became significantly depleted in 14C compared to modern atmospheric CO2. Models that take into account this 14C depletion in soil CO2 for dating modern groundwater would lead to apparent younger 14C ages than models that only consider the 14C activity in atmospheric CO2. Moreover, since 2000-2005, the inverse effect is observed as soil CO2 is enriched in 14C compared to atmospheric CO2.¦Therefore, we conclude that the isotopic composition of CO2 at the water table have to be taken into account for the dating of modern groundwater. This requires a systematic sampling of soil CO2 and the measurement of its 13C and 14C contents. We used this information in a numerical simulation to calculate the evolution of isotopic composition of CO2 from the soil surface to the water table. This simulation integrated physical processes in the unsaturated zone (e.g. CO2 production and diffusion, water advection, etc.) and gas-water-rock interactions.

Spatio-temporal variations in reproductive patterns and population structure of Pasiphaea multidentata (Decapoda: Caridea) in the Blanes canyon and adjacent margin, North-western Mediterranean Sea

The bathyal faunal communities of the NW Mediterranean slopes have been studied consistently in the last two decades, with a special focus on population structure, trophic dynamics and benthopelagic coupling of commercial deep-sea decapod crustaceans and fishes (reviewed in Sardà et al. 2004) and associated species (Cartes and Sardà, 1993; Company and Sardà, 1997, 2000; Cartes et al., 2001; Company et al., 2001, 2003, 2004). One of the major topographic features in the North-western Mediterranean slope is the presence of submarine canyons. Canyons play a major role in funnelling energy and organic matter from the shelf to bathyal and abyssal depths (Puig et al., 2000), but the implications of this enhanced organic supply in the deep-sea benthic communities is still mostly unknown. Trophic supply can follow two major pathways – vertical deposition in the water column (Billett et al., 1983; Baldwin et al., 1998; Lampitt et al., 2001) or down-slope advection on the margins (Puig et al., 2001; Bethoux et al., 2002; Canals et al., 2006) – and can be a limiting factor in the deep-sea, being especially important in the oligotrophic Mediterranean Sea (Sardà et al., 2004). Differences in the quantity, quality and timing of organic matter input to the deep seafloor have been used to explain patterns of biomass and abundance in benthic communities (Levin et al., 1994; Gooday & Turley, 1990; Billett et al., 2001; Galéron et al., 2001; Puig et al., 2001; Gage, 2003) as well as other biological process and in particular the existence of seasonal reproduction (Tyler et al., 1994; Company et al., 2004 (MEPS). Reproduction is a highly energetic process tightly linked to food availability and quality.

Low-pressure, water-assisted anatexis of basic dykes in a contact metamorphic aureole, Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite-gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative delta O-18 whole-rock values (+0.2 to -3.4 parts per thousand), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest delta O-18 values (-2.8, -4.2, - 4.4 and -7.6 parts per thousand, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic delta O-18 values up to the anatectic zone, where the values are lower and more heterogeneous. Low delta O-18 values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.

A block solver for the exponentally fitted IIPG-0 method

We consider an exponentially fitted discontinuous Galerkin method for advection dominated problems and propose a block solver for the resulting linear systems. In the case of strong advection the solver is robust with respect to the advection direction and the number of unknowns.

Application of low temperature thermochronology to understanding cooling and exhumation in the Central Swiss Alps

THESIS ABSTRACT : Low-temperature thermochronology relies on application of radioisotopic systems whose closure temperatures are below temperatures at which the dated phases are formed. In that sense, the results are interpreted as "cooling ages" in contrast to "formation ages". Owing to the low closure-temperatures, it is possible to reconstruct exhumation and cooling paths of rocks during their residence at shallow levels of the crust, i.e. within first ~10 km of depth. Processes occurring at these shallow depths such as final exhumation, faulting and relief formation are fundamental for evolution of the mountain belts. This thesis aims at reconstructing the tectono-thermal history of the Aar massif in the Central Swiss Alps by means of zircon (U-Th)/He, apatite (U-Th)/He and apatite fission track thermochronology. The strategy involved acquisition of a large number of samples from a wide range of elevations in the deeply incised Lötschen valley and a nearby NEAT tunnel. This unique location allowed to precisely constrain timing, amount and mechanisms of exhumation of the main orographic feature of the Central Alps, evaluate the role of topography on the thermochronological record and test the impact of hydrothermal activity. Samples were collected from altitudes ranging between 650 and 3930 m and were grouped into five vertical profiles on the surface and one horizontal in the tunnel. Where possible, all three radiometric systems were applied to each sample. Zircon (U-Th)/He ages range from 5.1 to 9.4 Ma and are generally positively correlated with altitude. Age-elevation plots reveal a distinct break in slope, which translates into exhumation rate increasing from ~0.4 to ~3 km/Ma at 6 Ma. This acceleration is independently confirmed by increased cooling rates on the order of 100°C/Ma constrained on the basis of age differences between the zircon (U-Th)/He and the remaining systems. Apatite fission track data also plot on a steep age-elevation curve indicating rapid exhumation until the end of the Miocene. The 6 Ma event is interpreted as reflecting tectonically driven uplift of the Aar massif. The late Miocene timing implies that the increase of precipitation in the Pliocene did not trigger rapid exhumation in the Aar massif. The Messinian salinity crisis in the Mediterranean could not directly intensify erosion of the Aar but associated erosional output from the entire Alps may have tapered the orogenic wedge and caused reactivation of thrusting in the Aar massif. The high exhumation rates in the Messinian were followed by a decrease to ~1.3 km/Ma as evidenced by ~8 km of exhumation during last 6 Ma. The slowing of exhumation is also apparent from apatite (U-Th)1He age-elevation data in the northern part of the Lötschen valley where they plot on a ~0.5km/Ma line and range from 2.4 to 6.4 Ma However, from the apatite (U-Th)/He and fission track data from the NEAT tunnel, there is an indication of a perturbation of the record. The apatite ages are youngest under the axis of the valley, in contrast to an expected pattern where they would be youngest in the deepest sections of the tunnel due to heat advection into ridges. The valley however, developed in relatively soft schists while the ridges are built of solid granitoids. In line with hydrological observations from the tunnel, we suggest that the relatively permeable rocks under the valley floor, served as conduits of geothermal fluids that caused reheating leading to partial Helium loss and fission track annealing in apatites. In consequence, apatite ages from the lowermost samples are too young and the calculated exhumation rates may underestimate true values. This study demonstrated that high-density sampling is indispensable to provide meaningful thermochronological data in the Alpine setting. The multi-system approach allows verifying plausibility of the data and highlighting sources of perturbation. RÉSUMÉ DE THÈSE : La thermochronologie de basse température dépend de l'utilisation de systèmes radiométriques dont la température de fermeture est nettement inférieure à la température de cristallisation du minéral. Les résultats obtenus sont par conséquent interprétés comme des âges de refroidissement qui diffèrent des âges de formation obtenus par le biais d'autres systèmes de datation. Grâce aux températures de refroidissement basses, il est aisé de reconstruire les chemins de refroidissement et d'exhumation des roches lors de leur résidence dans la croute superficielle (jusqu'à 10 km). Les processus qui entrent en jeu à ces faibles profondeurs tels que l'exhumation finale, la fracturation et le faillage ainsi que la formation du relief sont fondamentaux dans l'évolution des chaînes de montagne. Ces dernières années, il est devenu clair que l'enregistrement thermochronologique dans les orogènes peut être influencé par le relief et réinitialisé par l'advection de la chaleur liée à la circulation de fluides géothermaux après le refroidissement initial. L'objectif de cette thèse est de reconstruire l'histoire tectono-thermique du massif de l'Aar dans les Alpes suisses Centrales à l'aide de trois thermochronomètres; (U-Th)/He sur zircon, (U-Th)/He sur apatite et les traces de fission sur apatite. Afin d'atteindre cet objectif, nous avons récolté un grand nombre d'échantillons provenant de différentes altitudes dans la vallée fortement incisée de Lötschental ainsi que du tunnel de NEAT. Cette stratégie d'échantillonnage nous a permis de contraindre de manière précise la chronologie, les quantités et les mécanismes d'exhumation de cette zone des Alpes Centrales, d'évaluer le rôle de la topographie sur l'enregistrement thermochronologique et de tester l'impact de l'hydrothermalisme sur les géochronomètres. Les échantillons ont été prélevés à des altitudes comprises entre 650 et 3930m selon 5 profils verticaux en surface et un dans le tunnel. Quand cela à été possible, les trois systèmes radiométriques ont été appliqués aux échantillons. Les âges (U-Th)\He obtenus sur zircons sont compris entre 5.l et 9.4 Ma et sont corrélés de manière positive avec l'altitude. Les graphiques représentant l'âge et l'élévation montrent une nette rupture de la pente qui traduisent un accroissement de la vitesse d'exhumation de 0.4 à 3 km\Ma il y a 6 Ma. Cette accélération de l'exhumation est confirmée par les vitesses de refroidissement de l'ordre de 100°C\Ma obtenus à partir des différents âges sur zircons et à partir des autres systèmes géochronologiques. Les données obtenues par traces de fission sur apatite nous indiquent également une exhumation rapide jusqu'à la fin du Miocène. Nous interprétons cet évènement à 6 Ma comme étant lié à l'uplift tectonique du massif de l'Aar. Le fait que cet évènement soit tardi-miocène implique qu'une augmentation des précipitations au Pliocène n'a pas engendré cette exhumation rapide du massif de l'Aar. La crise Messinienne de la mer méditerranée n'a pas pu avoir une incidence directe sur l'érosion du massif de l'Aar mais l'érosion associée à ce phénomène à pu réduire le coin orogénique alpin et causer la réactivation des chevauchements du massif de l'Aar. L'exhumation rapide Miocène a été suivie pas une diminution des taux d'exhumation lors des derniers 6 Ma (jusqu'à 1.3 km\Ma). Cependant, les âges (U-Th)\He sur apatite ainsi que les traces de fission sur apatite des échantillons du tunnel enregistrent une perturbation de l'enregistrement décrit ci-dessus. Les âges obtenus sur les apatites sont sensiblement plus jeunes sous l'axe de la vallée en comparaison du profil d'âges attendus. En effet, on attendrait des âges plus jeunes sous les parties les plus profondes du tunnel à cause de l'advection de la chaleur dans les flancs de la vallée. La vallée est creusée dans des schistes alors que les flancs de celle-ci sont constitués de granitoïdes plus durs. En accord avec les observations hydrologiques du tunnel, nous suggérons que la perméabilité élevée des roches sous l'axe de la vallée à permi l'infiltration de fluides géothermaux qui a généré un réchauffement des roches. Ce réchauffement aurait donc induit une perte d'Hélium et un recuit des traces de fission dans les apatites. Ceci résulterait en un rajeunissement des âges apatite et en une sous-estimation des vitesses d'exhumation sous l'axe de la vallée. Cette étude à servi à démontrer la nécessité d'un échantillonnage fin et précis afin d'apporter des données thermochronologiques de qualité dans le contexte alpin. Cette approche multi-système nous a permi de contrôler la pertinence des données acquises ainsi que d'identifier les sources possibles d'erreurs lors d'études thermochronologiques. RÉSUMÉ LARGE PUBLIC Lors d'une orogenèse, les roches subissent un cycle comprenant une subduction, de la déformation, du métamorphisme et, finalement, un retour à la surface (ou exhumation). L'exhumation résulte de la déformation au sein de la zone de collision, menant à un raccourcissement et un apaissessement de l'édifice rocheux, qui se traduit par une remontée des roches, création d'une topographie et érosion. Puisque l'érosion agit comme un racloir sur la partie supérieure de l'édifice, des tentatives de corrélation entre les épisodes d'exhumation rapide et les périodes d'érosion intensive, dues aux changements climatiques, ont été effectuées. La connaissance de la chronologie et du lieu précis est d'une importance capitale pour une quelconque reconstruction de l'évolution d'une chaîne de montagne. Ces critères sont donnés par un retraçage des changements de la température de la roche en fonction du temps, nous donnant le taux de refroidissement. L'instant auquel les roches ont refroidit, passant une certaine température, est contraint par l'application de techniques de datation par radiométrie. Ces méthodes reposent sur la désintégration des isotopes radiogéniques, tels que l'uranium et le potassium, tous deux abondants dans les roches de la croûte terrestre. Les produits de cette désintégration ne sont pas retenus dans les minéraux hôtes jusqu'au moment du refroidissement de la roche sous une température appelée 'de fermeture' , spécifique à chaque système de datation. Par exemple, la désintégration radioactive des atomes d'uranium et de thorium produit des atomes d'hélium qui s'échappent d'un cristal de zircon à des températures supérieures à 200°C. En mesurant la teneur en uranium-parent, l'hélium accumulé et en connaissant le taux de désintégration, il est possible de calculer à quel moment la roche échantillonnée est passée sous la température de 200°C. Si le gradient géothermal est connu, les températures de fermeture peuvent être converties en profondeurs actuelles (p. ex. 200°C ≈ 7km), et le taux de refroidissement en taux d'exhumation. De plus, en datant par système radiométrique des échantillons espacés verticalement, il est possible de contraindre directement le taux d'exhumation de la section échantillonnée en observant les différences d'âges entre des échantillons voisins. Dans les Alpes suisses, le massif de l'Aar forme une structure orographique majeure. Avec des altitudes supérieures à 4000m et un relief spectaculaire de plus de 2000m, le massif domine la partie centrale de la chaîne de montagne. Les roches aujourd'hui exposées à la surface ont été enfouies à plus de 10 km de profond il y a 20 Ma, mais la topographie actuelle du massif de l'Aar semble surtout s'être développée par un soulèvement actif depuis quelques millions d'années, c'est-à-dire depuis le Néogène supérieur. Cette période comprend un changement climatique soudain ayant touché l'Europe il y a environ 5 Ma et qui a occasionné de fortes précipitations, entraînant certainement une augmentation de l'érosion et accélérant l'exhumation des Alpes. Dans cette étude, nous avons employé le système de datation (U-TH)/He sur zircon, dont la température de fermeture de 200°C est suffisamment basse pour caractériser l'exhumation du Néogène sup. /Pliocène. Les échantillons proviennent du Lötschental et du tunnel ferroviaire le plus profond du monde (NEAT) situé dans la partie ouest du massif de l'Aar. Considérés dans l'ensemble, ces échantillons se répartissent sur un dénivelé de 3000m et des âges de 5.1 à 9.4 Ma. Les échantillons d'altitude supérieure (et donc plus vieux) documentent un taux d'exhumation de 0.4 km/Ma jusqu'à il y a 6 Ma, alors que les échantillons situés les plus bas ont des âges similaires allant de 6 à 5.4 Ma, donnant un taux jusqu'à 3km /Ma. Ces données montrent une accélération dramatique de l'exhumation du massif de l'Aar il y a 6 Ma. L'exhumation miocène sup. du massif prédate donc le changement climatique Pliocène. Cependant, lors de la crise de salinité d'il y a 6-5.3 Ma (Messinien), le niveau de la mer Méditerranée est descendu de 3km. Un tel abaissement de la surface d'érosion peut avoir accéléré l'exhumation des Alpes, mais le bassin sud alpin était trop loin du massif de l'Aar pour influencer son érosion. Nous arrivons à la conclusion que la datation (U-Th)/He permet de contraindre précisément la chronologie et l'exhumation du massif de l'Aar. Concernant la dualité tectonique-érosion, nous suggérons que, dans le cas du massif de l'Aar, la tectonique prédomine.

A methodology for the adaptation of a PMP at the determination of a PMF

This paper presents a new type of very fine grid hydrological model based on the spatiotemporal repartition of a PMP (Probable Maximum Precipitation) and on the topography. The goal is to estimate the influence of this rain on a PMF (Probable Maximum Flood) on a catchment area in Switzerland. The spatiotemporal distribution of the PMP was realized using six clouds modeled by the advection-diffusion equation. The equation shows the movement of the clouds over the terrain and also gives the evolution of the rain intensity in time. This hydrological modeling is followed by a hydraulic modeling of the surface and subterranean flow, done considering the factors that contribute to the hydrological cycle, such as the infiltration, the resurgence and the snowmelt. These added factors make the developed model closer to reality and also offer flexibility in the initial condition that is added to the factors concerning the PMP, such as the duration of the rain, the speed and direction of the wind. All these initial conditions taken together offer a complete image of the PMF.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.

Partial melting of basic dykes in the contact aureoles of two shallow-level, gabbro-pyroxenite intrusions, Fuerteventura

THESIS ABSTRACT Low-pressure anatexis of basic dykes gave rise to unusual, zebra-like migmatites, in the contact metamorphic aureoles of two layered gabbro-pyroxenite intrusions, PXl and PX2, in the root zone of an ocean island, Fuerteventura Basal Complex (Canary Islands). This thesis focuses on the understanding of processes attributing to the partial melting and formation of these migmatites, characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations with perfectly preserved igneous textures. The presence of fluids are required to decrease the solidus of basic igneous lithologies, to allow partial melting in such aloes-pressure (1-2 kb) environment. An oxygen isotope study was thus carried out on dykes inside and beyond the PX2 aureole, in order to decipher the nature and origin of such fluids. Low or negative δ18O values were obtained for whole rocks and mineral-separates, decreasing towards the contact, with the intrusion itself retaining fairly high values. This trend has been attributed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. A comparison of whole rock and mineral major- and trace- element data allowed the redistribution of elements .between different mineral phases and generations, during contact metamorphism and partial melting to be assessed. Certain trace-elements, e.g. Zr, Hf, Y, and REEs, were internally redistributed during contact metamorphic recrystallisation, causing- the enrichment of neocrystallised diopsides compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post- solidus reequilibration on cooling of the migmatite terrain. The mineralogical, textural and geochemical evolution of dykes in a contact metamorphic aureole, is recorded around an apophysis of the PX1 intrusion, where there is evidence of incipient partial melting. Hydrothermal mineral pseudomorphs in the outer parts of the aureole are progressively replaced by dry mineral assemblages, with increasingly recrystallised diopside and evidence of partial melting -the extent of which varies from one lithology to another. The appearance of more mafic lithologies towards the intrusion, with lower whole rock SiO2 and mobile element abundances, e.g. Rb, Cs, K, has been explained by the migration and accumulation of feldspathic material into leucosomes outside the samples. A micro-structural study of leucosomes and leucocratic pods, with the aid of high-resolution X-ray computed micro-tomography (HRXµCT), allowing the visualization and quantification of shapes and orientations, was carried out in order to better understand the processes of melt segregation in the PX1 aureole. Leucocratic pods, representing former amygdales, are considered as natural strain ellipsoids. Their short axes are oriented perpendicular to leucosome planes, which sub-parallel the intrusive contact. Leucosomes thus effectively represent foliation planes. This implies that the direction of maximum shortening, during migmatisation, was perpendicular to the orientation of leucosomes, contradicting earlier models that suggest leucosomes represent tension veins. RESUME DE LA THESE Un phénomène rare de fusion partielle de filons basiques à basse pression a été étudié dans les auréoles de contact de deux intrusions litées de gabbro-pyroxénite, PX1 et PX2, localisées dans le soubassement de l'île volcanique de Fuerteventura aux Canaries. Cette anatexie a engendré des migmatites finement zébrées d'aspect très inhabituel, dont les processus de formation ont été étudiés dans le présent travail. Ces roches sont caractérisées par un réseau dense de veinules leucocrates d'épaisseur millimétrique, dont les textures ignées sont parfaitement préservées. La fusion partielle de roches basiques à basse pression (1-2 kbar) requiert la présence d'eau afin d'abaisser le solidus du système à des températures géologiquement réalistes. Une étude comparative des isotopes de l'oxygène a ainsi été menée sur des filons respectivement affectés et non affectés par le métamorphisme de contact, afin de confirmer la présence de ces fluides, de déterminer l'importance de leur interaction avec les roches et leur origine. Des valeurs de δ180 basses ou négatives ont été mesurées sur roche totale et minéraux séparés, décroissantes en direction du contact, alors que l'intrusion elle-même a conservé des valeurs élevées. Ce gradient a été attribué à l'advection d'eau météorique durant la mise en place du magma, les températures les plus élevées favorisant d'autant plus la circulation des fluides et les échanges isotopiques durant la recristallisation des roches. Cette recristallisation engendré une redistribution chimique complète des éléments entre les différentes générations de minéraux résultant du métamorphisme de contact et de l'anatexie, mise en évidence par microanalyse. Certains éléments traces comme Zr, Hf, Y et les REE ont été concentrés dans le diopside néoformé consécutivement à la déstabilisation de minéraux primaires riches en ces éléments comme la kaersutite ou le sphène. Les compositions en éléments majeurs et traces des minéraux des mélanosomes et leucosomes des migmatites sont pratiquement identiques, indiquant une rééquilibration syn- à postsolidus lors du refroidissement de l'auréole de contact. La transformation progressive des filons basiques au niveau de leur minéralogie, textures et composition chimique a pu être observée en détail à l'approche du contact d'une apophyse de l'intrusion PX1. La paragenèse magmatique initiale n'est jamais préservée, les faciès les plus distants du contact étant constitués d'un assemblage pseudomorphique hydrothermal. Ce dernier est progressivement remplacé par des assemblages anhydres incluant du diopside néoformé, puis apparaissent les premiers signes de fusion partielle, dont l'importance varie fortement d'une lithologie à l'autre. L'apparition de faciès plus basiques en direction du contact, avec des teneurs réduites en SiO2 et en éléments incompatibles tels Rb, Cs, K, a été attribuée à l'échappement de leucosomes feldspathiques hors du système. Une étude microstructurale de la distribution spatiale du matériel leucocrate au sein des migmatites par microtomographie X de haute résolution (HRXµCT) a été menée pour mieux comprendre les processus de ségrégation des liquides dans l'auréole de PX1. De petites entités ovoïdes, représentant d'anciennes structures amygdalaires au sein des filons, ont été considérées comme des ellipsoïdes marqueurs de la déformation finie. Leur petit axe est orienté perpendiculairement aux plans définis par les leucosomes, eux-mêmes subparallèles au contact intrusif. Les leucosomes matérialisent donc des plans de clivage. Ainsi, la direction de raccourcissement maximum durant la fusion partielle était perpendiculaire à l'orientation des leucosomes, contrairement à ce qui a été dit dans de précédentes publications, qui suggéraient que les leucosomes représentaient des veines de tension. RESUME DE LA THESE (POUR LE GRAND PUBLIC) L'observation directe du soubassement d'une île volcanique est une occasion rare, accessible dans le «complexe de base » de l'île canarienne de Fuerteventura. Ce dernier a enregistré divers phénomènes magmatiques, métamorphiques et de fusion partielle induits par l'intrusion répétée de magmas alimentant des appareils volcaniques sus jacents, sous forme de petits plutons, essaims de filons et complexes annulaires de gabbros alcalins, pyroxénites, syénites et carbonatites. Dans ce contexte de flux de chaleur élevé, des filons basiques ont subi une fusion partielle au contact de deux intrusions de gabbro-pyroxénite, un phénomène extrêmement rare à une profondeur aussi réduite, estimée à quelque 3-6 km. Les produits de cette fusion partielle sont des liquides très riches en feldspath, concentrés en un réseau dense de veinules blanches (leucosomes) au sein du matériau résiduel sombre non fondu (mélanosome) pour former ce qu'on appelle des migmatites. Outre les aspects pétrologiques liés à la formation de ces migmatites, l'intérêt majeur du phénomène réside dans le fait qu'il puisse représenter la source des magmas évolués parfois observés sur les îles océaniques. A des pressions aussi faibles que dans le soubassement de Fuerteventura (1-2 Kbar), la présence de fluides abondants est nécessaire pour abaisser la température de début de fusion des roches (solidus) à des valeurs géologiquement réalistes. Des expériences ont montré que même en présence de plusieurs %-poids d'eau, une température de 1000°C était encore nécessaire pour obtenir une proportion de liquide équivalente à celle observée sur le terrain, soit 25%. Or les magmas alcalins des îles océaniques, bien qu'hydratés, n'en contiennent de loin pas autant, ce qui implique une source d'eau externe. Une étude isotopique de l'oxygène a été entreprise afin de tester cette hypothèse. Les valeurs obtenues en 5180 sont basses ou négatives et indiquent l'influence d'eau d'origine météorique. Cette eau de pluie se serait infiltrée le long des filons depuis la surface du volcan et les aurait complètement hydrothermalisés en profondeur (situation encore visible à l'extérieur de l'auréole de contact), leur permettant ainsi de stocker l'eau nécessaire à leur fusion partielle ultérieure. L'interaction entre eau de pluie et filons a été d'autant plus importante que ces derniers étaient proches du contact avec l'intrusion, ce qui suggère que la circulation de ces eaux et leur interaction avec les roches a été favorisée par la chaleur fournie par l'intrusion elle-même. Un autre aspect de ce travail s'est focalisé sur la redistribution des éléments traces au sein des minéraux des filons basiques durant le métamorphisme de contact et la fusion partielle. Ainsi, le pyroxène de seconde génération est-il sensiblement enrichi en traces telles Zr, Hf, Y et les terres rares, par rapport au pyroxène magmatique originel, en relation avec la déstabilisation de minéraux primaires riches en ces éléments tels le sphène et la kaersutite. Cependant, les compositions en éléments majeurs et traces des minéraux recristallisés des migmatites sont pratiquement les mêmes dans les leucosomes et les mélanosomes, suggérant une rééquilibration chimique complète durant le refroidissement de ces lithologies. Si certaines migmatites se sont comportées en système chimiquement fermé (hormis l'eau météorique), d'autres filons ont manifestement perdu une partie de leurs leucosomes, ainsi qu'en témoigne leur composition progressivement appauvrie en silice et autres éléments incompatibles mobiles, tels K et Rb à l'approche du contact de l'intrusion. Parallèlement à cette évolution chimique, les paragenèses hydrothermales distantes du contact sont progressivement remplacées par des paragenèses anhydres, puis par l'apparition des premiers leucosomes, tandis que les textures magmatiques initiales sont complètement effacées au profit d'une combinaison de textures magmatique dans les leucosomes et en mosaïque dans les mélanosomes. Enfin, la distribution spatiale des liquides de fusion partielle a été étudiée par microtomographie X de haute résolution, dans des filons contenant des entités ovoïdes leucocrates, sans doute d'anciennes amygdales à zéolites. Ces dernières ont été considérées comme des ellipsoïdes de la déformation finie. L'orientation de leur petit axe, perpendiculaire au plan défini par les veinules de leucosomes, indique que ces derniers représentent des plans de clivage perpendiculaires à la direction de raccourcissement maximum. Ainsi, la ségrégation des liquides de fusion partielle se serait faite dans les plans de compression et non dans des plans de dilatation, contrairement à ce que laisserait penser le sens commun.

An observational study of the 7 September 2005 Barcelona tornado outbreak

This paper presents an observational study of the tornado outbreak that took place on the 7 September 2005 in the Llobregat delta river, affecting a densely populated and urbanised area and the Barcelona International airport (NE Spain). The site survey confirmed at least five short-lived tornadoes. Four of them were weak (F0, F1) and the other one was significant (F2 on the Fujita scale). They started mostly as waterspouts and moved later inland causing extensive damage estimated in 9 million Euros, three injured people but fortunately no fatalities. Large scale forcing was provided by upper level diffluence and low level warm air advection. Satellite and weather radar images revealed the development of the cells that spawned the waterspouts along a mesoscale convergence line in a highly sheared and relatively low buoyant environment. Further analysis indicated characteristics that could be attributed indistinctively to non-supercell or to mini-supercell thunderstorms.

Improving QPF by blending techniques at the Meteorological Service of Catalonia

The current operational very short-term and short-term quantitative precipitation forecast (QPF) at the Meteorological Service of Catalonia (SMC) is made by three different methodologies: Advection of the radar reflectivity field (ADV), Identification, tracking and forecasting of convective structures (CST) and numerical weather prediction (NWP) models using observational data assimilation (radar, satellite, etc.). These precipitation forecasts have different characteristics, lead time and spatial resolutions. The objective of this study is to combine these methods in order to obtain a single and optimized QPF at each lead time. This combination (blending) of the radar forecast (ADV and CST) and precipitation forecast from NWP model is carried out by means of different methodologies according to the prediction horizon. Firstly, in order to take advantage of the rainfall location and intensity from radar observations, a phase correction technique is applied to the NWP output to derive an additional corrected forecast (MCO). To select the best precipitation estimation in the first and second hour (t+1 h and t+2 h), the information from radar advection (ADV) and the corrected outputs from the model (MCO) are mixed by using different weights, which vary dynamically, according to indexes that quantify the quality of these predictions. This procedure has the ability to integrate the skill of rainfall location and patterns that are given by the advection of radar reflectivity field with the capacity of generating new precipitation areas from the NWP models. From the third hour (t+3 h), as radar-based forecasting has generally low skills, only the quantitative precipitation forecast from model is used. This blending of different sources of prediction is verified for different types of episodes (convective, moderately convective and stratiform) to obtain a robust methodology for implementing it in an operational and dynamic way.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Reconstructing snow avalanches in the Southeastern Pyrenees

A regional study of snow avalanche processes was undertaken in the SE Pyrenees. Dendrogeomorphology was used to date and reconstruct large-scale snow avalanche events that occurred in the last four decades. Dendrochronological analyses yielded the dates of nine winters when avalanches occurred in the recent past in six studied avalanche paths. Some of these avalanches were already known, but others had not been documented. In one case, the existing avalanche path map was improved with the dendrogeomorphological information of a larger past event. As a result of the dendrogeomorphological analyses, evidence for three regional-scale major avalanche years was identified in the SE Pyrenees from 1971 to 2004: 1971¿1972, 1995¿1996 and 2002¿2003. The specific synoptic atmospheric situations and the most likely nivometeorological and snowpack conditions that released these major avalanches were determined using weather data for the seasons of major avalanche releases. In 1971¿1972 the snow avalanche episode was characterized by a deep trough crossing the Pyrenees. In 1995¿1996 a variety of meteorological situations produced several episodes of major avalanches. In 2002¿2003 the more significant of two episodes was attributed to a north advection pumping an arctic air mass over the Pyrenees. The 1995¿1996 avalanche season proved to be the most notable in the four past decades in the Pyrenees.

Long-term and long range migration of radioisotopes in terrestrial systems : mechanisms of radioisotopes (Cs, Sr, Pu, Am) mobilization in soils

Summary Artificial radionuclides were released in the environment during the atmospheric nuclear weapon tests and after accidental events involving nuclear industries. As a primary receptor of the deposition, the soil is a very sensitive compartment and understanding the interaction and migration of radionuclides within soils allows the development of scenario for the contamination risk of the population and of the environment. Most available field studies on radionuclides in soils only concern one or two isotopes, mostly 137Cs, and few physico-chemical soil parameters. The purpose of this study was a broader understanding of the radioecology of an Alpine valley. In a first part, we aimed to describe the depth distribution of 137Cs, 90Sr, 239+240Pu, and 241Am within different alpine soils and to identify some stable elements as indicators for accumulating layers. In the central part of the study, the goal was to investigate the repartition of ^Sr and 239Pu between the truly dissolved fraction and the colloidal fraction of the soil solutions and to identify the nature of colloids involved in the adsorption of ^Sr and 239Pu. These results were integrated in an "advection- sorption" transport model seeking to explain the migration of 239Pu and 90Sr within the soils and to assess the importance of colloidal transport for these two isotopes. A further aspect studied was the role of the competition between the radioisotopes (137Cs and 90Sr) and their stable chemical analogues (K and Ca) with respect to plant uptake by different plant species. The results on the depth distribution within the soils showed that 137Cs was mostly retained in the topsoil, to the exception of an organic-rich soil (Histosol 2) receiving important surface runoff, where migration down to a depth of 30 cm was observed. 137Cs depth distribution within the soils was similar to unsupported 210Pb depth distribution. The plant uptake of 137Cs clearly depended on the concentration of exchangeable potassium in the soils. Moreover, we showed that the 137Cs uptake by certain species of the taxonomic orders Poales and Rosales was more sensitive to the increase in exchangeable Κ compared to other orders. Strontium-90 was much more mobile in the soils than 137Cs and depth migration and accumulation in specific AI- and Fe-rich layers were found down to 30 cm. Copper and Ni showed accumulations in these same layers, indicating their potential to be used as indicators for the migration of ^Sr within the soils. In addition, we observed a 90Sr activity peak in the topsoil that can be attributable to recycling of 90Sr by plant uptake. We demonstrated for the first time that a part of 90Sr (at least 40%) was associated with the colloids in organic-rich soil solutions. Therefore, we predict a significant effect of the colloidal migration of ^Sr in organic-rich soil solutions. The plant uptake results for 90Sr indicated a phylogenetic effect between Non-Eudicot and Eudicots: the order Poales concentrating much less 90Sr than Eudicots do. Moreover, we were able to demonstrate that the sensitivity of the 90Sr uptake by 5 different Alpine plant species to the amount of exchangeable Ca was species-independent. Plutonium and 241Am accumulated in the second layer of all soils and only a slight migration deeper than 20 cm was observed. Plutonium and 241Am showed a similar depth distribution in the soils. The model results suggested that the present day migration of 239Pu was very slow and that the uptake by plants was negligible. 239Pu activities between 0.01 to 0.08 mBq/L were measured in the bulk soil solutions. Migration of 239Pu with the soil solution is dominated by colloidal transport. We reported strong evidences that humic substances were responsible of the sorption of 239Pu to the colloidal fraction of the soil solutions. This was reflected by the strong correlation between 239Pu concentrations and the content of (colloidal) organic matter in the soil solution. Résumé Certains radioéléments artificiels ont été disséminés dans l'environnement suite aux essais atmosphériques de bombes nucléaires et suite à des accidents impliquant les industries nucléaires. En tant que récepteur primaire de la déposition, le sol est un compartiment sensible et des connaissances sur les interactions et la migration des radioéléments dans le sol permettent de développer des modèles pour estimer la contamination de la population et de l'environnement. Actuellement, la plupart des études de terrain sur ce sujet concernent uniquement un ou deux radioéléments, surtout le 137Cs et peu d'études intègrent les paramètres du sol pour expliquer la migration des radioéléments. Le but général de cette étude était une compréhension étendue de la radio-écologie d'une vallée alpine. Notre premier objectif était de décrire la distribution en profondeur de 137Cs, ^Sr, 239+240pu et 241Am dans différents sols alpins en relation avec des éléments stables du sol, dans le but d'identifier des éléments stables qui pourraient servir d'indicateurs pour des horizons accumulateurs. L'objectif de la deuxième partie, qui était la partie centrale de l'étude, était d'estimer le pourcentage d'activité sous forme colloïdale du 239Pu et du 90Sr dans les solutions des sols. De plus nous avons déterminé la nature des colloïdes impliqués dans la fixation du ^Sr et 239Pu. Nous avons ensuite intégré ces résultats dans un modèle de transport développé dans le but de décrire la migration du 239Pu et 90Sr dans le sol. Finalement, nous avons étudié l'absorption de 137Cs et 90Sr par les plantes en fonction de l'espèce et de la compétition avec leur élément analogue stable (K et Ca). Les résultats sur la migration en profondeur du 137Cs ont montré que ce radioélément était généralement retenu en surface, à l'exception d'un sol riche en matière organique dans lequel nous avons observé une nette migration en profondeur. Dans tous les sols, la distribution en profondeur du 137Cs était corrélée avec la distribution du 210Pb. L'absorption du 137Cs par les plantes, était dépendante de la concentration en Κ échangeable dans le sol, le potassium étant un compétiteur. De plus, nous avons observé que les espèces ne réagissaient pas de la même manière aux variations de la concentration de Κ échangeable. En effet, les espèces appartenant aux ordres des Poales et des Rosales étaient plus sensibles aux variations de potassium échangeable dans le sol. Dans tous les sols Le 90Sr était beaucoup plus mobile que le 137Cs. En effet, nous avons observé des accumulations de 90Sr dans des horizons riches en Fe et Al jusqu'à 30 cm de profondeur. De plus, le Cu et le Ni montraient des accumulations dans les mêmes horizons que le 90Sr, indiquant qu'il pourrait être possible d'utiliser ces deux éléments comme analogues pour la migration du 90Sr. D'après le modèle développé, le pic de 90Sr dans les premiers centimètres du sol peut être attribué à du recyclage par les plantes. Le 90Sr en solution était principalement sous forme dissoute dans des solutions de sols peu organique (entre 60 et 100% de 90Sr dissous). Par contre, dans des solutions organiques, un important pourcentage de 90Sr (plus de 40%) était associé aux colloïdes. La migration colloïdale du 90Sr peut donc être significative dans des solutions organiques. Comme pour le 137Cs, l'absorption du 90Sr par les plantes dépendait de la concentration de son analogue chimique dans la fraction échangeable du sol. Par contre, les espèces de plantes étudiées avaient la même sensibilité aux variations de la concentration du calcium échangeable. Le plutonium et l'américium étaient accumulés dans le deuxième horizon du sol et nous avons observé seulement une faible migration plus profondément que 20 cm. Selon le modèle, la migration actuelle du plutonium est très lente et l'absorption par les plantes semble négligeable. Nous avons mesuré entre 0.01 et 0.08 mBq/L de 239Pu dans les solutions de sol brutes. La migration du plutonium par la solution du sol est due principalement aux colloïdes, probablement de nature humique. Résumé grand public Dans les années 1950 à 1960, l'environnement a été contaminé par des éléments radioactifs (radioéléments) artificiels provenant des essais des armes atomiques et de l'industrie nucléaire. En effet, durant ces années, les premiers essais de bombes atomiques se faisaient dans l'atmosphère, libérant de grandes quantités d'éléments radioactifs. De plus certains accidents impliquant l'industrie nucléaire civile ont contribué à la dissémination d'éléments radioactifs dans l'environnement. Ce fut par exemple le cas de l'accident de la centrale atomique de Tchernobyl en 1986 qui a causé une importante contamination d'une grande partie de l'Europe par le 137Cs. Lorsqu'ils sont libérés dans l'atmosphère, les radioéléments sont dispersés et transportés par les courants atmosphériques, puis peuvent être déposés dans l'environnement, principalement par les précipitations. Une fois déposés sur le sol, les radioéléments vont interagir avec les composants du sol et migrer plus ou moins vite. La connaissance des interactions des éléments radioactifs avec le sol est donc importante pour prédire les risques de contamination de l'environnement et de l'homme. Le but général de ce travail était d'évaluer la migration de différents éléments radioactifs (césium-137, strontium-90, plutonium et américium-241) à travers le sol. Nous avons choisi un site d'étude en milieu alpin (Val Piora, Tessin, Suisse), contaminé en radioéléments principalement par les retombées de l'accident de Tchernobyl et des essais atmosphériques de bombes atomiques. Dans un premier temps, nous avons caractérisé la distribution en profondeur des éléments radioactifs dans le sol et l'avons comparée à divers éléments stables. Cette comparaison nous a permit de remarquer que le cuivre et le nickel s'accumulaient dans les mêmes horizons du sol que le strontium-90 et pourraient donc être utilisés comme analogue pour la migration du strontium-90 dans les sols. Dans la plupart des sols étudiés, la migration du césium-137, du plutonium et de l'américium-241 était lente et ces radioéléments étaient donc accumulés dans les premiers centimètres du sol. Par contre, le strontium-90 a migré beaucoup plus rapidement que les autres radioéléments si bien qu'on observe des accumulations de strontium-90 à plus de 30 cm de profondeur. Les radioéléments migrent dans la solution du sol soit sous forme dissoute, soit sous forme colloïdale, c'est-à-dire associés à des particules de diamètre < Ιμηι. Cette association avec des colloïdes permet à des radioéléments peu solubles, comme le plutonium, de migrer plus rapidement qu'attendu. Nous avons voulu savoir quelle était la part de strontium-90 et plutonium associés à des colloïdes dans la solution du sol. Les résultats ont montré que le plutonium en solution était principalement associé à des colloïdes de type organique. Quant au strontium-90, ce dernier était en partie associé à des colloïdes dans des solutions de sol riches en matière organique, par contre, il était principalement sous forme dissoute dans les solutions de sol peu organiques. L'absorption de radioéléments par les plantes représente une voie importante pour le transfert vers la chaîne alimentaire, par conséquent pour la contamination de l'homme. Nous avons donc étudié le transfert du césium-137 et du strontium-90 de plusieurs sols vers différentes espèces de plantes. Les résultats ont montré que l'absorption des radioéléments par les plantes était liée à la concentration de leur analogue chimique (calcium pour le strontium-90 et potassium pour le césium- 137) dans la fraction échangeable du sol. De plus certaines espèces de plantes accumulent significativement moins de strontium-90.