Temporal evolution of emerging contaminants, priority substances and heavy metals most frequently detected in reclaimed water and grounwater in a volcanic aquifer

[EN] The presence of emerging contaminants has been previously described in reclaimed water and groundwater of Gran Canaria (Spain). Despite of the environmental risk associated to irrigation with reclaimed water (R), this practice is necessary considering sustainability of the hydrological cycle in semiarid zones, especially regarding agricultural activity. The aim of this study was: i) to analyse the evolution during two years of contaminants of emerging concern, priority substances (2008/105/EC) and heavy metals in reclaimed water (R) and in a volcanic aquifer in the NE of Gran Canaria where a golf course has been irrigated with R since 1976 and ii) to relate this presence with physicochemical water properties and hydrogeological media. Reclaimed water and groundwater (GW) were monitoring quarterly from July 2009 to September 2011. Sorption and degradation processes in soil account for more compounds being detected in R. Diazinon and chlorfenvinphos were detected always in R and terbuthylazine, terbutryn and diuron at 90% of frequency. Considering all the samples, the most frequent compounds were chlorpyrifos ethyl, fluorene, phenanthrene and pyrene. Although their concentrations were frequently below 50 ngL-1, some contaminants, were occasionally detected at higher concentrations. Chlorpyrifos ethyl and diuron are priority substances detected frequently and at high concentrations so they must be included in monitoring studies. Geology and location seem to be related to the emerging compounds presence due to occasional contamination events (not related to R irrigation) and therefore not to an existence of a dangerous diffuse contamination level. Thus, it is preferable to select wells with less stable chemical water quality, in order to monitor the risk of emerging compounds presence. Considering the relationship between contaminant presence, chemical water quality, seasonal variation, hydrogeological characteristics and wells location we can conclude that chlorpyrifos ethyl and diuron were the most dangerous priority substances in terms of GW quality so they must be included in all of the monitoring studies, at least in Canary Islands.

Reconstruction of the aquifer history of Gran Canaria

[EN]A study on the recent history and current state of the aquifer in the Island of Gran Canaria (Canary Is., 28oN, 15oW) is performed. Though rainfall is scarce on the island, traditional agricultural practices and small population were able to keep the aquifer in a constant state for centuries. Nevertheless, at the beginning of the 20th Century, culture of several water-consuming species was introduced on a commercial basis due to the relative proximity of the Canaries to continental Europe and to the possibility of more than one yearly harvest. This led to generalised well digging (more than 300m deep in many cases) and to the appearance of a chronic hydraulic deficit, as well as to spoiling vastcoastal areas of the aquifer through intrusion of brackish water. In the mid 1960’s, coincident with the apex of agricultural exploitation, massive tourism appeared in the scene. This new activity soon became a susbstitute for Agriculture, but it attracted more new labour force to the island, and a fast growth of population was the main result. Moreover, new water use practices entered the scene. As a consequence, the main causes for the aquifer decline are population growth and extensive Agriculture practices in use during the last half of the 20th Century. Some remarks on sustainability issues in order to cope with Climate Change are also offered.

Shallow and deep deformation in northern Apennines region using seismological data

For its particular position and the complex geological history, the Northern Apennines has been considered as a natural laboratory to apply several kinds of investigations. By the way, it is complicated to joint all the knowledge about the Northern Apennines in a unique picture that explains the structural and geological emplacement that produced it. The main goal of this thesis is to put together all information on the deformation - in the crust and at depth - of this region and to describe a geodynamical model that takes account of it. To do so, we have analyzed the pattern of deformation in the crust and in the mantle. In both cases the deformation has been studied using always information recovered from earthquakes, although using different techniques. In particular the shallower deformation has been studied using seismic moment tensors information. For our purpose we used the methods described in Arvidsson and Ekstrom (1998) that allowing the use in the inversion of surface waves [and not only of the body waves as the Centroid Moment Tensor (Dziewonski et al., 1981) one] allow to determine seismic source parameters for earthquakes with magnitude as small as 4.0. We applied this tool in the Northern Apennines and through this activity we have built up the Italian CMT dataset (Pondrelli et al., 2006) and the pattern of seismic deformation using the Kostrov (1974) method on a regular grid of 0.25 degree cells. We obtained a map of lateral variations of the pattern of seismic deformation on different layers of depth, taking into account the fact that shallow earthquakes (within 15 km of depth) in the region occur everywhere while most of events with a deeper hypocenter (15-40 km) occur only in the outer part of the belt, on the Adriatic side. For the analysis of the deep deformation, i.e. that occurred in the mantle, we used the anisotropy information characterizing the structure below the Northern Apennines. The anisotropy is an earth properties that in the crust is due to the presence of aligned fluid filled cracks or alternating isotropic layers with different elastic properties while in the mantle the most important cause of seismic anisotropy is the lattice preferred orientation (LPO) of the mantle minerals as the olivine. This last is a highly anisotropic mineral and tends to align its fast crystallographic axes (a-axis) parallel to the astenospheric flow as a response to finite strain induced by geodynamic processes. The seismic anisotropy pattern of a region is measured utilizing the shear wave splitting phenomenon (that is the seismological analogue to optical birefringence). Here, to do so, we apply on teleseismic earthquakes recorded on stations located in the study region, the Sileny and Plomerova (1996) approach. The results are analyzed on the basis of their lateral and vertical variations to better define the earth structure beneath Northern Apennines. We find different anisotropic domains, a Tuscany and an Adria one, with a pattern of seismic anisotropy which laterally varies in a similar way respect to the seismic deformation. Moreover, beneath the Adriatic region the distribution of the splitting parameters is so complex to request an appropriate analysis. Therefore we applied on our data the code of Menke and Levin (2003) which allows to look for different models of structures with multilayer anisotropy. We obtained that the structure beneath the Po Plain is probably even more complicated than expected. On the basis of the results obtained for this thesis, added with those from previous works, we suggest that slab roll-back, which created the Apennines and opened the Tyrrhenian Sea, evolved in the north boundary of Northern Apennines in a different way from its southern part. In particular, the trench retreat developed primarily south of our study region, with an eastward roll-back. In the northern portion of the orogen, after a first stage during which the retreat was perpendicular to the trench, it became oblique with respect to the structure.

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.

The Influence of Fractures on Vulnerability to Pollution of a Carbonate Aquifer: Monte Conero (Italy)

The carbonate outcrops of the anticline of Monte Conero (Italy) were studied in order to characterize the geometry of the fractures and to establish their influence on the petrophysical properties (hydraulic conductivity) and on the vulnerability to pollution. The outcrops form an analog for a fractured aquifer and belong to the Maiolica Fm. and the Scaglia Rossa Fm. The geometrical properties of fractures such as orientation, length, spacing and aperture were collected and statistically analyzed. Five types of mechanical fractures were observed: veins, joints, stylolites, breccias and faults. The types of fractures are arranged in different sets and geometric assemblages which form fracture networks. In addition, the fractures were analyzed at the microscale using thin sections. The fracture age-relationships resulted similar to those observed at the outcrop scale, indicating that at least three geological episodes have occurred in Monte Conero. A conceptual model for fault development was based on the observations of veins and stylolites. The fracture sets were modelled by the code FracSim3D to generate fracture network models. The permeability of a breccia zone was estimated at microscale by and point counting and binary image methods, whereas at the outcrop scale with Oda’s method. Microstructure analysis revealed that only faults and breccias are potential pathways for fluid flow since all veins observed are filled with calcite. According this, three scenarios were designed to asses the vulnerability to pollution of the analogue aquifer: the first scenario considers the Monte Conero without fractures, second scenario with all observed systematic fractures and the third scenario with open veins, joints and faults/breccias. The fractures influence the carbonate aquifer by increasing its porosity and hydraulic conductivity. The vulnerability to pollution depends also on the presence of karst zones, detric zones and the material of the vadose zone.

Land use management in mountainous areas: combining ground-based and EO (Earth Observation) data to investigate the shallow landsliding susceptibility in the Duron valley (Trento, Italy).

Throughout the alpine domain, shallow landslides represent a serious geologic hazard, often causing severe damages to infrastructures, private properties, natural resources and in the most catastrophic events, threatening human lives. Landslides are a major factor of landscape evolution in mountainous and hilly regions and represent a critical issue for mountainous land management, since they cause loss of pastoral lands. In several alpine contexts, shallow landsliding distribution is strictly connected to the presence and condition of vegetation on the slopes. With the aid of high-resolution satellite images, it's possible to divide automatically the mountainous territory in land cover classes, which contribute with different magnitude to the stability of the slopes. The aim of this research is to combine EO (Earth Observation) land cover maps with ground-based measurements of the land cover properties. In order to achieve this goal, a new procedure has been developed to automatically detect grass mantle degradation patterns from satellite images. Moreover, innovative surveying techniques and instruments are tested to measure in situ the shear strength of grass mantle and the geomechanical and geotechnical properties of these alpine soils. Shallow landsliding distribution is assessed with the aid of physically based models, which use the EO-based map to distribute the resistance parameters across the landscape.

Hydrogeochemical characterization of Seehausen territory - Bremen

The aim of this thesis is to provide a geochemical characterization of the Seehausen territory (a neighborhood) of Bremen, Germany. In this territory it is hosted a landfill of dredged sediments coming both from Bremerhaven (North See) and Bremen harbor (directly on the river Weser). For this reason this work has been focused also on possible impacts of the landfill on the groundwaters (shallow and deep aquifer). The Seehausen landfill uses the dewatering technique to manage the dredged sediments: incoming sediments are put into dewatering fields until they are completely dried (it takes almost a year). Then they are randomly sampled and analyzed: if the pollutants content is acceptable, sediments are treated with other materials and used instead of raw material for embankment, bricks, etc., otherwise they are disposed in the landfill. During this work it has been made a study of the natural geology and hydrogeology of the whole area of interest, especially because it is characterized by ancient natural salt deposits. Then, together with the Geological Survey of Bremen and the Harbor Authority of Bremen there have been identified all useful piezometers for a monitoring net around the landfill. During the sampling campaign there have been collected data of the principal anions and cations, physical parameters and stable water isotopes. Data analysis has been focused particularly on Cl, Na, SO4 and EC because these parameters might be helpful to attribute geochemical trends to the landfill or to a natural background. Furthermore dataloggers have been installed for a month in some piezometers and EC, pressure, dissolved oxygen and temperature data have been collected. Finally there has been made a deep comparison between current and historical data (1996 – 2011) and between old interpolation maps and current ones in order to see time trends of the aquifer geochemistry.

Diversity of microbial communities at shallow-hydrothermal vents off Dominica, Lesser Antilles

Hydrothermal vents are often compared to desert oases, because of the presence of highly diverse and abundant biotic communities inhabiting these extreme environments. Nevertheless, the microbial communities associated with shallow-hydrothermal systems have been poorly studied. Hydrothermal activity at Dominica Island is quite well known under the geological and geochemical aspects, but no previous information existed about the microbial communities associated to this area. This thesis is therefore targeting the microbiology of hydrothermal sediments combining geochemical and molecular biological investigations, focusing on differences between hydrothermal vents and background (i.e. control) areas, and between hydrothermal sites. It was also intended to assess relationship between geochemical parameters and microbial diversity at the two hydrothermally impacted sites. Two shallow-sea hydrothermal vents located south-west off Dominica Island (Lesser Antilles) have been investigated in this study: Champagne Hot Springs and Soufriere Bay offshore vent. During this study, sediments for geochemical and molecular analyses were collected every 2 cm from the two impacted areas and from two control sites not associated with hydrothermal activity; in situ temperatures measurements were also taken every 5 cm deep in the sediment for all the sites. A geochemical characterization of the sediment porewater was performed through the analysis of several elements’ concentrations (i.e. H2S, Cl-, Br-, SO42-, Fe2+, Na+, K+, B+, Si+). Microbial communities at the different sites were studied by Automated Intergenic Spacer Analysis (ARISA). Inspection of the operational taxonomic units (OTUs) distribution was performed, as well as statistical analyses for communities’ structure and composition differences, and for changes of β-diversity along with sediment geochemistry. Data suggested that mixing between hydrothermal fluids and seawater results in distinct different environmental gradients and potential ecological niches between the two investigated hydrothermal vents, reflecting a difference in microbial community structures between them.

IMEX finite volume methods for the shallow water equations

Die Flachwassergleichungen (SWE) sind ein hyperbolisches System von Bilanzgleichungen, die adäquate Approximationen an groß-skalige Strömungen der Ozeane, Flüsse und der Atmosphäre liefern. Dabei werden Masse und Impuls erhalten. Wir unterscheiden zwei charakteristische Geschwindigkeiten: die Advektionsgeschwindigkeit, d.h. die Geschwindigkeit des Massentransports, und die Geschwindigkeit von Schwerewellen, d.h. die Geschwindigkeit der Oberflächenwellen, die Energie und Impuls tragen. Die Froude-Zahl ist eine Kennzahl und ist durch das Verhältnis der Referenzadvektionsgeschwindigkeit zu der Referenzgeschwindigkeit der Schwerewellen gegeben. Für die oben genannten Anwendungen ist sie typischerweise sehr klein, z.B. 0.01. Zeit-explizite Finite-Volume-Verfahren werden am öftersten zur numerischen Berechnung hyperbolischer Bilanzgleichungen benutzt. Daher muss die CFL-Stabilitätsbedingung eingehalten werden und das Zeitinkrement ist ungefähr proportional zu der Froude-Zahl. Deswegen entsteht bei kleinen Froude-Zahlen, etwa kleiner als 0.2, ein hoher Rechenaufwand. Ferner sind die numerischen Lösungen dissipativ. Es ist allgemein bekannt, dass die Lösungen der SWE gegen die Lösungen der Seegleichungen/ Froude-Zahl Null SWE für Froude-Zahl gegen Null konvergieren, falls adäquate Bedingungen erfüllt sind. In diesem Grenzwertprozess ändern die Gleichungen ihren Typ von hyperbolisch zu hyperbolisch.-elliptisch. Ferner kann bei kleinen Froude-Zahlen die Konvergenzordnung sinken oder das numerische Verfahren zusammenbrechen. Insbesondere wurde bei zeit-expliziten Verfahren falsches asymptotisches Verhalten (bzgl. der Froude-Zahl) beobachtet, das diese Effekte verursachen könnte.Ozeanographische und atmosphärische Strömungen sind typischerweise kleine Störungen eines unterliegenden Equilibriumzustandes. Wir möchten, dass numerische Verfahren für Bilanzgleichungen gewisse Equilibriumzustände exakt erhalten, sonst können künstliche Strömungen vom Verfahren erzeugt werden. Daher ist die Quelltermapproximation essentiell. Numerische Verfahren die Equilibriumzustände erhalten heißen ausbalanciert.rnrnIn der vorliegenden Arbeit spalten wir die SWE in einen steifen, linearen und einen nicht-steifen Teil, um die starke Einschränkung der Zeitschritte durch die CFL-Bedingung zu umgehen. Der steife Teil wird implizit und der nicht-steife explizit approximiert. Dazu verwenden wir IMEX (implicit-explicit) Runge-Kutta und IMEX Mehrschritt-Zeitdiskretisierungen. Die Raumdiskretisierung erfolgt mittels der Finite-Volumen-Methode. Der steife Teil wird mit Hilfe von finiter Differenzen oder au eine acht mehrdimensional Art und Weise approximniert. Zur mehrdimensionalen Approximation verwenden wir approximative Evolutionsoperatoren, die alle unendlich viele Informationsausbreitungsrichtungen berücksichtigen. Die expliziten Terme werden mit gewöhnlichen numerischen Flüssen approximiert. Daher erhalten wir eine Stabilitätsbedingung analog zu einer rein advektiven Strömung, d.h. das Zeitinkrement vergrößert um den Faktor Kehrwert der Froude-Zahl. Die in dieser Arbeit hergeleiteten Verfahren sind asymptotisch erhaltend und ausbalanciert. Die asymptotischer Erhaltung stellt sicher, dass numerische Lösung das "korrekte" asymptotische Verhalten bezüglich kleiner Froude-Zahlen besitzt. Wir präsentieren Verfahren erster und zweiter Ordnung. Numerische Resultate bestätigen die Konvergenzordnung, so wie Stabilität, Ausbalanciertheit und die asymptotische Erhaltung. Insbesondere beobachten wir bei machen Verfahren, dass die Konvergenzordnung fast unabhängig von der Froude-Zahl ist.

The effect of temperature on pressure tolerance of the shallow-water shrimp Palaemon serratus

Tolerance to low temperature and high pressure may allow shallow-water species to extend bathymetric range in response to changing climate, but adaptation to contrasting shallow-water environments may affect tolerance to these factors. The brackish shallow-water shrimp Palaemon varians demonstrates remarkable tolerance to elevated hydrostatic pressure and low temperature, but inhabits a highly variable environment: environmental adaptation may therefore make P. varians tolerances unrepresentative of other shallow-water species. Critical thermal maximum (CTmax), critical hydrostatic pressure maximum (CPmax), and acute respiratory response to hydrostatic pressure were assessed in the shallow-water shrimp Palaemon serratus, which inhabits a more stable intertidal habitat. P. serratus’ CTmax was 22.3°C when acclimated at 10°C, and CPmax was 5.9, 10.1, and 14.1 MPa when acclimated at 5, 10, and 15°C respectively: these critical tolerances were consistently lower than P. varians. Respiratory responses to acute hyperbaric exposures similarly indicated lower tolerance to hydrostatic pressure in P. serratus than in P. varians. Contrasting tolerances likely reflect physiological adaptation to differing environments and reveal that the capacity for depth-range extension may vary among species from different habitats.

The contribution of the Finite Volume Point Dilution Method to the determination of solute mass flux in an alluvial aquifer.

Groundwater represents one of the most important resources of the world and it is essential to prevent its pollution and to consider remediation intervention in case of contamination. According to the scientific community the characterization and the management of the contaminated sites have to be performed in terms of contaminant fluxes and considering their spatial and temporal evolution. One of the most suitable approach to determine the spatial distribution of pollutant and to quantify contaminant fluxes in groundwater is using control panels. The determination of contaminant mass flux, requires measurement of contaminant concentration in the moving phase (water) and velocity/flux of the groundwater. In this Master Thesis a new solute flux mass measurement approach, based on an integrated control panel type methodology combined with the Finite Volume Point Dilution Method (FVPDM), for the monitoring of transient groundwater fluxes, is proposed. Moreover a new adsorption passive sampler, which allow to capture the variation of solute concentration with time, is designed. The present work contributes to the development of this approach on three key points. First, the ability of the FVPDM to monitor transient groundwater fluxes was verified during a step drawdown test at the experimental site of Hermalle Sous Argentau (Belgium). The results showed that this method can be used, with optimal results, to follow transient groundwater fluxes. Moreover, it resulted that performing FVPDM, in several piezometers, during a pumping test allows to determine the different flow rates and flow regimes that can occurs in the various parts of an aquifer. The second field test aiming to determine the representativity of a control panel for measuring mass flus in groundwater underlined that wrong evaluations of Darcy fluxes and discharge surfaces can determine an incorrect estimation of mass fluxes and that this technique has to be used with precaution. Thus, a detailed geological and hydrogeological characterization must be conducted, before applying this technique. Finally, the third outcome of this work concerned laboratory experiments. The test conducted on several type of adsorption material (Oasis HLB cartridge, TDS-ORGANOSORB 10 and TDS-ORGANOSORB 10-AA), in order to determine the optimum medium to dimension the passive sampler, highlighted the necessity to find a material with a reversible adsorption tendency to completely satisfy the request of the new passive sampling technique.