Evolução cenozóica da área de Longroiva-Vilariça (NE Portugal)

In the Longroiva-Vilariça area, the identification of Cenozoic lithostratigraphic units, the sedimentology and the characterization of its geometric relations with tectonic structures allowed the interpretation of the palaeogeographic main stages: 1) the greenwhitish Vilariça Arkoses (Middle Eocene to Oligocene ?) represent proximal sediments of a very low gradient drainage towards the eastern Spanish Tertiary Duero Basin; 2)Quintãs Formation (late Miocene ?) are brown-reddish coloured piedmont alluvial deposits, correlative of important vertical displacement (western tectonic block relative uplift) along the NNE-SSW indent-linked strike-slip Bragança-Vilariça-Longroiva fault zone, interpreted as a reactivated deep hercynian fracture, with left-lateral movement; 3) the red Sampaio Formation (Gelasian-early Pleistocene ?)was interpreted as downhill conglomeratic deposits related with important overtrusting along this fault zone (the definition of the present-day narrow graben configuration) and correlative of the atlantic hydrographic incision stage beginning; 4) conglomeratic terraces (middle and late Pleistocene ?); 5) alluvial plains and colluvial deposits (Holocene).

Litostratigrafia do Terciario da região Miranda do Corvo-Viseu (Bacia do Mondego, Portugal)

The lithostratigraphic units (2 groups, 5 formations and 4 members) of the Tertiary of the Miranda do Corvo-Viseu region (Central Portugal) are here described. For each unit the characterization and description (boundaries, diagnostic properties and atributes) were included. The stratigraphic, lithological, sequential and tectonic data allow correlations with other units of the same Tertiary basin located more to west and support the chronostratigraphic attribution. Sedimentologic characteristics of the deposits lead to the interpretation of the influence of tectonism, climate and eustasy during the sedimentary evolution of this Basin.

Budgeting rockfall and modeling sedimentary delivery in torrent systems

This thesis is a compilation of projects to study sediment processes recharging debris flow channels. These works, conducted during my stay at the University of Lausanne, focus in the geological and morphological implications of torrent catchments to characterize debris supply, a fundamental element to predict debris flows. Other aspects of sediment dynamics are considered, e.g. the coupling headwaters - torrent, as well as the development of a modeling software that simulates sediment transfer in torrent systems. The sediment activity at Manival, an active torrent system of the northern French Alps, was investigated using terrestrial laser scanning and supplemented with geostructural investigations and a survey of sediment transferred in the main torrent. A full year of sediment flux could be observed, which coincided with two debris flows and several bedload transport events. This study revealed that both debris flows generated in the torrent and were preceded in time by recharge of material from the headwaters. Debris production occurred mostly during winter - early spring time and was caused by large slope failures. Sediment transfers were more puzzling, occurring almost exclusively in early spring subordinated to runoffconditions and in autumn during long rainfall. Intense rainstorms in summer did not affect debris storage that seems to rely on the stability of debris deposits. The morpho-geological implication in debris supply was evaluated using DEM and field surveys. A slope angle-based classification of topography could characterize the mode of debris production and transfer. A slope stability analysis derived from the structures in rock mass could assess susceptibility to failure. The modeled rockfall source areas included more than 97% of the recorded events and the sediment budgets appeared to be correlated to the density of potential slope failure. This work showed that the analysis of process-related terrain morphology and of susceptibility to slope failure document the sediment dynamics to quantitatively assess erosion zones leading to debris flow activity. The development of erosional landforms was evaluated by analyzing their geometry with the orientations of potential rock slope failure and with the direction of the maximum joint frequency. Structure in rock mass, but in particular wedge failure and the dominant discontinuities, appear as a first-order control of erosional mechanisms affecting bedrock- dominated catchment. They represent some weaknesses that are exploited primarily by mass wasting processes and erosion, promoting not only the initiation of rock couloirs and gullies, but also their propagation. Incorporating the geological control in geomorphic processes contributes to better understand the landscape evolution of active catchments. A sediment flux algorithm was implemented in a sediment cascade model that discretizes the torrent catchment in channel reaches and individual process-response systems. Each conceptual element includes in simple manner geomorphological and sediment flux information derived from GIS complemented with field mapping. This tool enables to simulate sediment transfers in channels considering evolving debris supply and conveyance, and helps reducing the uncertainty inherent to sediment budget prediction in torrent systems. Cette thèse est un recueil de projets d'études des processus de recharges sédimentaires des chenaux torrentiels. Ces travaux, réalisés lorsque j'étais employé à l'Université de Lausanne, se concentrent sur les implications géologiques et morphologiques des bassins dans l'apport de sédiments, élément fondamental dans la prédiction de laves torrentielles. D'autres aspects de dynamique sédimentaire ont été abordés, p. ex. le couplage torrent - bassin, ainsi qu'un modèle de simulation du transfert sédimentaire en milieu torrentiel. L'activité sédimentaire du Manival, un système torrentiel actif des Alpes françaises, a été étudiée par relevés au laser scanner terrestre et complétée par une étude géostructurale ainsi qu'un suivi du transfert en sédiments du torrent. Une année de flux sédimentaire a pu être observée, coïncidant avec deux laves torrentielles et plusieurs phénomènes de charriages. Cette étude a révélé que les laves s'étaient générées dans le torrent et étaient précédées par une recharge de débris depuis les versants. La production de débris s'est passée principalement en l'hiver - début du printemps, causée par de grandes ruptures de pentes. Le transfert était plus étrange, se produisant presque exclusivement au début du printemps subordonné aux conditions d'écoulement et en automne lors de longues pluies. Les orages d'été n'affectèrent guère les dépôts, qui semblent dépendre de leur stabilité. Les implications morpho-géologiques dans l'apport sédimentaire ont été évaluées à l'aide de MNT et études de terrain. Une classification de la topographie basée sur la pente a permis de charactériser le mode de production et transfert. Une analyse de stabilité de pente à partir des structures de roches a permis d'estimer la susceptibilité à la rupture. Les zones sources modélisées comprennent plus de 97% des chutes de blocs observées et les bilans sédimentaires sont corrélés à la densité de ruptures potentielles. Ce travail d'analyses des morphologies du terrain et de susceptibilité à la rupture documente la dynamique sédimentaire pour l'estimation quantitative des zones érosives induisant l'activité torrentielle. Le développement des formes d'érosion a été évalué par l'analyse de leur géométrie avec celle des ruptures potentielles et avec la direction de la fréquence maximale des joints. Les structures de roches, mais en particulier les dièdres et les discontinuités dominantes, semblent être très influents dans les mécanismes d'érosion affectant les bassins rocheux. Ils représentent des zones de faiblesse exploitées en priorité par les processus de démantèlement et d'érosion, encourageant l'initiation de ravines et couloirs, mais aussi leur propagation. L'incorporation du control géologique dans les processus de surface contribue à une meilleure compréhension de l'évolution topographique de bassins actifs. Un algorithme de flux sédimentaire a été implémenté dans un modèle en cascade, lequel divise le bassin en biefs et en systèmes individuels répondant aux processus. Chaque unité inclut de façon simple les informations géomorpologiques et celles du flux sédimentaire dérivées à partir de SIG et de cartographie de terrain. Cet outil permet la simulation des transferts de masse dans les chenaux, considérants la variabilité de l'apport et son transport, et aide à réduire l'incertitude liée à la prédiction de bilans sédimentaires torrentiels. Ce travail vise très humblement d'éclairer quelques aspects de la dynamique sédimentaire en milieu torrentiel.

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Sedimentary Organic Matter in Cores of the Cananeia-Iguape Lagoonal-Estuarine System, Sao Paulo State, Brazil

Sedimentary organic matter is a good tool for environmental evaluation where the sediments are deposited. We determined the elemental and C- and N-isotopic compositions of 211 sub-surface sediment samples from 13 cores (ranging from 18 to 46cm), collected in the Cananeia-Iguape estuarine-lagoonal system. The aim of this research is to evaluate the environmental variations of this tropical coastal micro-tidal system over the last decades, through SOM distribution. The studied parameters show differences between the cores located in the northern (sandy-silt sediments) and southern (sand and silty-sand) portions. The whole area presents a mixed organic matter origin signature (local mangrove plants: < -25.6 parts per thousand PDB/ phytoplancton delta(13)C values: -19.4 parts per thousand PDB). The northern cores, which submitted higher sedimentation deposition (1.46cm year(-1)), are more homogenous, presenting lower delta(13)C (< -25.2 parts per thousand PDB) and higher C/N values (in general >14), directly related to the terrestrial input from Ribeira de Iguape River (24,000 km(2) basin). The southern portion presents lower sedimentation rates (0.38cm year(-1)) and is associated to a small river basin (1,340 km(2)), presenting values Of delta(13)C: -25.0 to 23.0 parts per thousand PDB and of C/N ratio: 11 to 15. In general, the elemental contents in the 15 cores may be considered from low to medium (< 2.0% C - < 0.1% N), compared to similar environments. Although a greater marine influence is observed in the southern system portion, the majority of the cores present an elevated increase of continental deposition, most likely related to the strong silting process that the area has been subjected to since the 1850s, when an artificial channel was built linking, directly, the Ribeira River to the estuarine-lagoonal system.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Sedimentary cyclicity in CRP drillcore, Victoria Land Basin, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-1.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

Global distribution and genetic discontinuities of mangroves - emerging patterns in the evolution of Rhizophora

Mangroves are often described as a group of plants with common features and common origins based mostly on their broad distributional patterns, together with an erroneous view of comparable abilities in long-distance dispersal. However, whilst mangroves have common needs to adapt to rigorous environmental constraints associated with regular seawater inundation, individual taxa have developed different strategies and characteristics. Since mangroves are a genetically diverse group of mostly flowering plants, they may also have evolved at quite different geological periods, dispersed at different rates from different locations and developed different adaptive strategies. Current distributions of individual taxa show numerous instances of unusual extant distribution which demonstrate finite dispersal limitations, especially across open water. Our preliminary assessment of broad distribution and discontinuities reveals important patterns. Discontinuities, in the absence of current dispersal barriers, may be explained by persistent past barriers. As we learn more about discontinuities, we are beginning to appreciate their immense implications and what they might tell us about past geological conditions and how these might have influenced the distribution and evolution of mangroves. In this article, we describe emerging patterns in genetic relationships and distributions based on both current knowledge and preliminary results of our studies of molecular and morphometric characteristics of Rhizophora species in the Indo West Pacific region.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Stability of planar flames as gasdynamic discontinuities

The stability of a steadily propagating planar premixed flame has been the subject of numerous studies since Darrieus and Landau showed that in their model flames are unstable to perturbations of any wavelength. Moreover, the instability was shown to persist even for very small wavelengths, i.e. there was no high-wavenumber cutoff of the instability. In addition to the Darrieus-Landau instability, which results from thermal expansion, analysis of the diffusional thermal model indicates that premixed flames may exhibit cellular and pulsating instabilities as a consequence of preferential diffusion. However, no previous theory captured all the instabilities including a high-wavenumber cutoff for each. In Class, Matkowsky & Klimenko (2003) a unified theory is proposed which, in appropriate limits and under appropriate assumptions, recovers all the relevant previous theories. It also includes additional new terms, not present in previous theories. In the present paper we consider the stability of a uniformly propagating planar flame as a solution of the unified model. The results are then compared to those based on the models of Darrieus-Landau, Sivashinsky and Matalon-Matkowsky. In particular, it is shown that the unified model is the only model to capture the Darrieus-Landau, cellular and pulsating instabilities including a high-wavenumber cutoff for each.

A new sedimentary benchmark for the Deccan Traps volcanism?

The origin of the Cretaceous-Paleogene boundary (KPB) mass extinction is still the center of acrimonious debates by opposing partisans of the bolide impact theory to those who favored a terrestrial origin linked to the Deccan Traps volcanism. Here we apply an original and high-resolution environmental magnetic study of the reference Bidart section, France. Our results show that the KPB is identified by an abrupt positive shift of the magnetic susceptibility (MS), also observed by others at the KPB elsewhere. In addition, an anomalous interval of very low MS, carried by an unknown Cl-bearing iron oxide similar to specular hematite, is depicted just below the KPB. Grain-size and morphology of the Cl-iron oxide are typically in the range of hematitic dust currently transported by winds from Sahara to Europe. This discovery is confirmed in the referenced Gubbio section (Italy) suggesting a global scale phenomenon. As a conjecture we suggest an origin by heterogeneous reaction between HCl-rich volcanic gas and liquid-solid aerosols within buoyant atmospheric plumes formed above the newly emitted Deccan flood basalts. Based on this hypothesis, our discovery provides a new benchmark for the Deccan volcanism and witnesses the nature and importance of the related atmospheric change.

Variation of magnetic properties in sedimentary rocks hosting the Foum Zguid dyke (southern Morocco): Combined effects of re-crystallization and Fe-metasomatism

The effects of dyke intrusion on the magnetic properties of host sedimentary rocks are still poorly understood. Therefore, we have evaluated bulk magnetic parameters of standard palaeomagnetic samples collected along several sections across the sediments hosting the Foum Zguid dyke in southern Morocco. The study has been completed with the evaluation of the magnetic fabric after laboratory application of sequential heating experiments. The present study shows that: (1) close to Fourn Zguid dykes, the variations of the bulk magnetic parameters and of the magnetic fabric is strongly related with re-crystallization and Fe-metasomatism intensity. (2) The thermal experiments on AMS of samples collected farther from the dyke and, thus, less affected by heating during dyke emplacement, indicate that 300-400 degrees C is the minimum experimental temperature necessary to trigger appreciable transformations of the pre-existing magnetic fabrics. For temperatures higher than ca. 580 degrees C, the magnetic fabric transformations are fully realized, with complete transposition of the initial fabric to a fabric similar to that of samples collected close to the dyke. Therefore, measured variations of the magnetic fabric can be used to evaluate re-crystallization temperatures experienced by the host sedimentary rock during dyke emplacement. The distinct magnetic behaviour observed along the cross-sections strongly suggests that samples collected farther from the dyke margins did not experience thermal episodes with temperatures higher than 300 degrees C after dyke emplacement. (3) AMS data shows a gradual variation of the magnetic fabric with distance from the dyke margin, from sub-horizontal K-3 away from the dyke to vertical K3 close to the dyke. Experimental heating shows that heat alone can be responsible for this strong variation. Therefore, such orientation changes should not be unequivocally interpreted as the result of a stress field (resulting from the emplacement of the dyke, for instance). (4) Magnetic studies prove to be a very sensitive tool to assess rock magnetic transformations, thermally and chemically induced by dyke intrusion in hosting sediments.

Coastal sedimentary environments and sea-level changes

A detailed knowledge of the 3-D arrangement and lateral facies relationships of the stacking patterns in coastal deposits is essential to approach many geological problems such as precise tracing of sea level changes, particularly during small scale fluctuations. These are useful data regarding the geodynamic evolution of basin margins and yield profit in oil exploration. Sediment supply, wave-and tidal processes, coastal morphology, and accommodation space generated by eustasy and tectonics govern the highly variable architecture of sedimentary bodies deposited in coastal settings. But these parameters change with time, and erosional surfaces may play a prominent role in areas located towards land. Besides, lateral shift of erosional or even depositional loci very often results in destruction of large parts of the sediment record. Several case studies illustrate some commonly found arrangements of facies and their distinguishing features. The final aim is to get the best results from the sedimentological analysis of coastal units.