The El Masnou Infralittoral sedimentary environment (Barcelona province, NW Mediterranean Sea): morphology and Holocene seismic stratigraphy

A detailed analysis of the morphology and the Holocene seismic and sequence stratigraphy and architecture of the infralittoral sedimentary environment of the El Masnou coast (Catalonia, NW Mediterranean Sea) was carried out using multibeam bathymetry and GeoPulse seismic data. This environment extends down to 26-30 m water depth, and is defined morphologically by two depositional wedges whose seafloor is affected by erosive furrows, slides, fields of large- and small-scale wavy bedforms, and dredging trenches and pits. Erosive terraces are also identified in the transition domain toward the inner continental shelf. The Holocene stratigraphy of the infralittoral environment is defined by two major seismic sequences (lower and upper), each one formed by internal seismic units. The sequences and units are characterised by downlapping surfaces made up of deposits formed by progradation of coastal lithosomes. The stratigraphy and stratal architecture, displaying a retrogradational arrangement with progradational patterns of minor order, were controlled by different sea-level positions. The stratigraphic division represents the coastal response to the last fourth-order transgressive and highstand conditions, modulated by small-scale sea-level oscillations (≈1-2 m) of fith to sixth order. This study also highlights the advantage of an integrated analysis using acoustic/seismic methods for practical assessment of the anthropogenic effects on infralittoral domains based on the association of marine geological observations.

Tectonic-palaeoenvironmental forcing of clay-mineral assemblages in nonmarine settings: the Oligocene-Miocene As Pontes Basin (Spain)

Two small, alluvial-lacustrine subbasins developed during the early restraining overstep stages of the Oligocene-Miocene As Pontes strike-slip Basin (NW Spain). Later, the basin evolved into a restraining bend stage and an alluvial-swamp-dominated depositional framework developed. The palaeobiological record demonstrates that the Oligocene-Miocene palaeoclimate in NW Spain was subtropical, warm and humid to subhumid.

An interactive web-GIS tool for risk analysis: a case study in the Fella River basin, Italy

This paper presents a prototype of an interactive web-GIS tool for risk analysis of natural hazards, in particular for floods and landslides, based on open-source geospatial software and technologies. The aim of the presented tool is to assist the experts (risk managers) in analysing the impacts and consequences of a certain hazard event in a considered region, providing an essential input to the decision-making process in the selection of risk management strategies by responsible authorities and decision makers. This tool is based on the Boundless (OpenGeo Suite) framework and its client-side environment for prototype development, and it is one of the main modules of a web-based collaborative decision support platform in risk management. Within this platform, the users can import necessary maps and information to analyse areas at risk. Based on provided information and parameters, loss scenarios (amount of damages and number of fatalities) of a hazard event are generated on the fly and visualized interactively within the web-GIS interface of the platform. The annualized risk is calculated based on the combination of resultant loss scenarios with different return periods of the hazard event. The application of this developed prototype is demonstrated using a regional data set from one of the case study sites, Fella River of northeastern Italy, of the Marie Curie ITN CHANGES project.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

Deep incision in an Aptian carbonate succession indicates major sea-level fall in the Cretaceous

Long-term relative sea-level cycles (0 5 to 6 Myr) have yet to be fully understood for the Cretaceous. During the Aptian, in the northern Maestrat Basin (Eastern Iberian Peninsula), fault-controlled subsidence created depositional space, but eustasy governed changes in depositional trends. Relative sea-level history was reconstructed by sequence stratigraphic analysis. Two forced regressive stages of relative sea-level were recognized within three depositional sequences. The first stage is late Early Aptian age (intra Dufrenoyia furcata Zone) and is characterized by foreshore to upper shoreface sedimentary wedges, which occur detached from a highstand carbonate platform, and were deposited above basin marls. The amplitude of relative sea-level drop was in the order of tens of metres, with a duration of <1 Myr. The second stage of relative sea-level fall occurred within the Late Aptian and is recorded by an incised valley that, when restored to its pre-contractional attitude, was >2 km wide and cut 115 m down into the underlying Aptian succession. With the subsequent transgression, the incision was back-filled with peritidal to shallow subtidal deposits. The changes in depositional trends, lithofacies evolution and geometric relation of the stratigraphic units characterized are similar to those observed in coeval rocks within the Maestrat Basin, as well as in other correlative basins elsewhere. The pace and magnitude of the two relative sea-level drops identified fall within the glacio-eustatic domain. In the Maestrat Basin, terrestrial palynological studies provide evidence that the late Early and Late Aptian climate was cooler than the earliest part of the Early Aptian and the Albian Stage, which were characterized by warmer environmental conditions. The outcrops documented here are significant because they preserve the results of Aptian long-term sea-level trends that are often only recognizable on larger scales (i.e. seismic) such as for the Arabian Plate.

Deep incision in an Aptian carbonate succession indicates major sea-level fall in the Cretaceous

Long-term relative sea-level cycles (0 5 to 6 Myr) have yet to be fully understood for the Cretaceous. During the Aptian, in the northern Maestrat Basin (Eastern Iberian Peninsula), fault-controlled subsidence created depositional space, but eustasy governed changes in depositional trends. Relative sea-level history was reconstructed by sequence stratigraphic analysis. Two forced regressive stages of relative sea-level were recognized within three depositional sequences. The first stage is late Early Aptian age (intra Dufrenoyia furcata Zone) and is characterized by foreshore to upper shoreface sedimentary wedges, which occur detached from a highstand carbonate platform, and were deposited above basin marls. The amplitude of relative sea-level drop was in the order of tens of metres, with a duration of <1 Myr. The second stage of relative sea-level fall occurred within the Late Aptian and is recorded by an incised valley that, when restored to its pre-contractional attitude, was >2 km wide and cut 115 m down into the underlying Aptian succession. With the subsequent transgression, the incision was back-filled with peritidal to shallow subtidal deposits. The changes in depositional trends, lithofacies evolution and geometric relation of the stratigraphic units characterized are similar to those observed in coeval rocks within the Maestrat Basin, as well as in other correlative basins elsewhere. The pace and magnitude of the two relative sea-level drops identified fall within the glacio-eustatic domain. In the Maestrat Basin, terrestrial palynological studies provide evidence that the late Early and Late Aptian climate was cooler than the earliest part of the Early Aptian and the Albian Stage, which were characterized by warmer environmental conditions. The outcrops documented here are significant because they preserve the results of Aptian long-term sea-level trends that are often only recognizable on larger scales (i.e. seismic) such as for the Arabian Plate.

Factors affecting Hg (II) adsorption in soils from the Rio Negro basin (Amazon)

Mercury (II) adsorption studies in top soils (top 10 cm) from the Rio Negro basin show this process depends strongly on some selected parameters of the aqueous phase in contact with the soils. Maximum adsorption occurred in the pH range 3.0-5.0 (>90%). Dissolved organic matter shows an inhibitory effect on the availability of Hg (II) to be adsorbed by the soils, whereas a higher chloride content of the solution resulted in a lower adsorption of Hg (II) at pH 5.0. Soils with higher organic matter content were less affected by changes in the salinity. An increase in the initial Hg (II) concentration increased the amount of Hg (II) adsorbed by the soil and decreased the time needed to reach equilibrium. A Freundlich isotherm provided a good model for Hg (II) adsorption in the two types of soil studied. The kinetics of Hg (II) adsorption on Amazonian soils showed to be very fast and followed pseudo-second order kinetics. An environmental implication of these results is discussed under the real scenario present in the Negro River basin, where acidic waters are in contact with a soil naturally rich in mercury.

Modelling syntectonic sedimentation: combining a discrete element model of tectonic. Deformation and process-based sedimentary Model in 3D.

This paper presents a new numerical program able to model syntectonic sedimentation. The new model combines a discrete element model of the tectonic deformation of a sedimentary cover and a process-based model of sedimentation in a single framework. The integration of these two methods allows us to include the simulation of both sedimentation and deformation processes in a single and more effective model. The paper describes briefly the antecedents of the program, Simsafadim-Clastic and a discrete element model, in order to introduce the methodology used to merge both programs to create the new code. To illustrate the operation and application of the program, analysis of the evolution of syntectonic geometries in an extensional environment and also associated with thrust fault propagation is undertaken. Using the new code, much more complex and realistic depositional structures can be simulated together with a more complex analysis of the evolution of the deformation within the sedimentary cover, which is seen to be affected by the presence of the new syntectonic sediments.

Organochlorine pesticides in Piracicaba river basin (São Paulo/Brazil): a survey of sediment, bivalve and fish

The main objective of this paper was to evaluate the level of occurrence of the organochlorine compounds in samples of sediments, bivalves and two fish species collected in the Piracicaba River basin (São Paulo, Brazil). The isomers alpha and gamma of HCH and Heptachlor were most frequently detected in samples of sediments and specimens of bivalve and fish. Therefore, although the levels of these compounds found were not critically high, they are still found in the environment. This fact suggests that they are still being used, despite the fact that the use of these compounds was outlawed more than twenty years ago.

Impacts of climate change on water resources in the Mediterranean basin : a case study in Catalonia, Spain

Most climate change projections show important decreases in water availability in the Mediterranean region by the end of this century. We assess those main climate change impacts on water resources in three medium-sized catchments with varying climatic conditions in north-eastern Spain. A combination of hydrological modelling and climate projections with B1 and A2 IPCC emission scenarios is performed to infer future stream flows. The largest reduction (22-48% for 2076-2100) of stream flows is expected in the headwaters of the two wettest catchments, while lower decreases (22-32% for 2076-2100) are expected in the drier one. In all three catchments, autumn and summer are the seasons with the most notable projected decreases in stream flow, 50% and 34%, respectively (2076-2100). Thus, ecological flows might be noticeably impacted by climate change in the catchments, especially in the headwaters of those wet catchments.

Fontes e distribuição de matéria orgânica sedimentar no complexo estuarino-lagunar de Mundaú-Manguaba/AL, utilizando esterois e alcoóis como indicadores

The origin of the sedimentary organic matter in the Mundaú-Manguaba estuarine system was evaluated through the distribution and composition of sterols, alcohols and ancillary data. The muddy sediments of the lagoons are enriched in organic matter, derived from a mixture of autochthonous and allochthonous inputs. Mundaú exhibited moderate contamination by sewage, with coprostanol concentration as high as 4.4 µg g-1. The channels are characterized by sand and organic matter-poor sediments. The sediments from some rivers in the drainage basin accumulated organic matter derived only from terrestrial vegetation. Stanol/sterol ratios evidenced the preferential diagenesis of the autochthonous fraction of organic matter.

Terpanos pentacíclicos como indicadores de heterogeneidades composicionais em reservatório de petróleo biodegradado

Reservoir geochemistry has been effectively used to evaluate compositional heterogeneity in petroleum reservoirs by employing oil fingerprinting technique to assess reservoir continuity, primarily in paraffin oil accumulations. In-reservoir biodegradation has resulted in vast deposits of heavy oils globally; therefore, the application of reservoir geochemistry to such accumulations becomes necessary. We recommend the use of pentacyclic terpanes instead of lower molecular compounds, which are less resistant to biodegradation. Using oil fingerprinting technique in this novel way enabled the differentiation of oil composition among sedimentary facies of tar sands (Pirambóia Formation, Paraná Basin).

Reactive dissolution of sedimentary rocks in flue gas desulfurization : modeling and experimental investigation

Svavel förekommer i kol och olja och oxideras vid förbränning till svaveldioxid (SO2). Årligen utsläpps stora mängder svaveldioxid som åstadkommer sura regn, minskning av stratosfäriskt ozon och sjukdomar. Av dessa orsaker är det nödvändigt att minska utsläppen av svaveldioxid. Den teknologi som rör planering av reaktorer och processer för rökgasavsvavling (FGD) har utvecklats kraftigt och idag använder man olika typer av lösningar. De vanligaste alternativen för FGD är våtskrubber- och semitorra skrubbersystem (Spray Dry Scrubbers, SDS) och injektionsprocesser för absorbenter. SDS-processer har studerats bland annat av Ahlbeck [4] och Klingspor [5] och injektionsprocesser av Alvfors [6]. Kalksten, som i huvusak består av kalciumkarbonat, används i rökgasavsvavling på grund av sin förmåga att binda svavel i form av sulfatsalter. Den vanligaste rökgasavsvavlingsmetoden är våtskrubbning, där det sedimentära stenmaterialets upplösningshastighet är en av de faktorer som påverkar resultatet mest. Utvärdering av kalkstensreaktivitet är därför speciellt viktig vid planering och drift av anläggningar för rökgasavsvavling. Målsättningen med detta arbete var att modellera upplösningen av olika typers kalksten för att få en kvantitativ utvärdering av kvaliteten på de analyserade proverna. Därtill testades även karbonatbiprodukter från stålindustrin för att utvärdera möjligheter att använda andra råmaterial. Det transienta förloppet har analyserats, varvid upplösningshastigheten modellerades bl.a. i avseende på tid och pH. Under arbetets gång har antalet empiriska korrelationer minskats till fördel för fysikaliska modeller av diffusiva och konvektiva masstransportfenomen. En målsättning var att skapa en effektiv och snabb metod för att testa olika absorbenter för rökgasavsvavling under transienta förlopp. I arbetet användes PSD-analys, gjordes pH-mätningar och andra utvärderingar av de fysikaliska parametrar som ingår i beräkningarna. On-line mätningar för de icke-stationära variablerna tid och pH ger möjlighet att eliminera osäkerheter. Vissa modeller kan vara komplicerade. En modell för upplösningshastigheten med mer detaljerad utvärdering av parametrar och färre approximationer är därför nödvändig då man vill utvärdera reaktionshastigheten för fasta partiklar i sur miljö. Arbetet utfördes under fyra år och fem peer review-artiklar ingår i avhandlingen.

METHODOLOGY FOR AUTOMATICALLY DELIMITING PERMANENT PRESERVATION AREAS ALONG WATER COURSES - THE USE OF GIS IN THE HYDROLOGICAL BASIN OF THE SERGIPE RIVER, BRAZIL

ABSTRACT Permanent Preservation Areas (PPAs) along watercourses have been the focus of numerous studies, not only because of the fragility and ecological relevance of riverine vegetation, but also because of the inefficiency demonstrated in conforming to the legislation protecting it. One of the major difficulties encountered in terms of guaranteeing the effective conservation of these riverside areas is the absence of methodologies that can be used to define them rapidly and accurately without manually determining the widths of the rivers or assigning only uniform linear values for the entire watercourse. The present work sought to develop a spatial analysis methodology capable of automatically defining permanent preservation areas along watercourses using geographic information system (GIS) software. The present study was undertaken in the Sergipe River basin, "considering the river itself and its principal affluents. We used the database of the Digital Atlas of Hydrological Resources (SEMARH/SE), and the delimitations of the PPAs were performed using ArcGIS 10.1 and the XToolPro 9.0 extension. A total of 5,003.82 hectares of Permanent Preservation Areas were delimited along the margins of the rivers analyzed, with a margin of error of <1% in delimiting the widths of the rivers within the entire area considered. The methodology described here can be used to define PPAs efficiently, relatively rapidly, and with very small margins of error, thus representing a technological advance in terms of using GIS for land management.