Biogeochemistry and dynamics of settling particle fluxes at the Antikythira Strait (Eastern Mediterranean)

For the first time, a 12-month trap experiment was conducted on both sides of the strait between Crete and Antikythira Island (Eastern Mediterranean Sea) from June 1994 to June 1995 as part of the PELAGOS experiment. Analyses of major chemical constituents, including carbohydrates and stable lead isotopes and Scanning Electron Microscope studies were performed on the trap samples. Total mass fluxes varied between 1 and 1273 mg m⁻² d⁻¹. The lowest fluxes observed were in summer and autumn 1994, when stratification of the water column was at its deepest. In general, mass fluxes exhibited very low values throughout this experiment confirming the strong oligotrophy of this area. The mean contents of the major constituents (carbonates, opal, lithogenic fraction) were quite similar during the survey and between traps, with the exception of organic carbon contents, which were highest (7–10%) in summer 1994, i.e. during the period of lowest mass fluxes. During the first 6-month deployment (summer–autumn 1994) there was an important mass flux peak, which was depleted in organic carbon, at the Ionian near-bottom trap. This event coincided with a violent wind episode, which may have caused the resuspension of particles, which were then transported down the steep continental slope on the Ionian side of the strait. A smaller peak in mass flux occurred at the Aegean near-bottom trap, coincident with rainfall. Both these events indicate that environmental factors can control flux variations in an oligotrophic environment. During the second 6-month deployment (winter–spring 1995) there was another important increase in mass fluxes, which occurred at all three traps, although in the Ionian traps mass flux peaks were delayed by one to two sampling intervals. The distance between the two mooring sites gives a rough estimate of a minimum horizontal advection speed of 2 cm s⁻¹ for this particulate transfer from the Aegean to the Ionian area. This estimate is in good agreement with the measured current velocities.     

Petrography,Geochemistry and Physical Properties of a Continuous Extrusive Section from the Sarami Massif,Semail Ophiolite

As opposed to ocean crustal drilling that often yields a discontinuous core, ophiolites provide a unique opportunity to study continuous sections of oceanic basalts. In order to provide high frequency constraints on the construction of the upper oceanic crust, a continuous 280 m thick volcanic transect was sampled in the Semail ophiolite in Oman. The analyzed section is located in the Sarami Massif, in the central part of the ophiolite, and exposed along Wadi Shaffan. A multidisciplinary study was carried out after sampling in the field. Core measurements including porosity, grain density, compressional velocity, magnetic susceptibility and electrical resistivity for over 100 samples from this transect have been made. Geochemical analyses including major, trace and Rare Earth Elements were also performed on 23 selected samples from this transect. The Wadi Shaffan transect appears as composed of two main petrological and geochemical members providing the main lithostratigraphic signal. The boundary between the two sequences is marked by the presence of a massive flow unit equating to the most primitive lava and across which significant changes in mineralogy and texture are observed. The physical properties equally record significant changes. In all, the section is characterized by chemical compositions coherent with that of V1-Geotime volcanism and appears as built through two main sequences of volcanic activity. The magnetic susceptibility profile correlated to geochemical variations present the most detailed evolution allowing to individualize a succession of at least 4 volcanic phases. These phases are interpreted as magmatic cycles characterized by differentiation processes, allowing one to place high frequency constraints on the volcanic behavior of fast spreading ridges. Some of the physical properties also reflect this series of differentiation trends, as noticed in the past from downhole measurements recorded in the context of deep ocean drilling into the crust.     

Structure of the Nemrut caldera (Eastern Anatolia, Turkey) and associated hydrothermal fluid circulation

Plio-Quaternary volcanism played an important role in the present physical state of Eastern Anatolia. Mount Nemrut, situated to the west of Lake Van is one of the main volcanic centers in the region, with a spectacular summit caldera 8.5 × 7 km in diameter. The most recent eruptions of the volcano were in 1441, 1597 and 1692. Nemrut Lake covers the western half of the caldera; it is a deep, half-bowl-shaped lake with a maximum depth of 176 m. Numerous eruption centers are exposed within the caldera as a consequence of magma–water interaction. Current activity of Nemrut caldera is revealed as hot springs, fumaroles and a small, hot lake.Self-potential and bathymetric surveys carried out in the caldera were used to characterize the structure of the caldera and the associated hydrothermal fluid circulation. In addition, analyses based on digital elevation models and satellite imagery were used to improve our knowledge about the structure of the caldera. According to SP results, the flanks of the volcano represent “the hydrogeologic zone”, whereas the intra-caldera region is an “active hydrothermal area” where the fluid circulation is controlled by structural discontinuities. There is also a northern fissure zone which exhibits hydrothermal signatures. Nemrut caldera collapsed piecemeal, with three main blocks. Stress controlling the collapse mechanism seems to be highly affected by the regional neotectonic regime. In addition to the historical activity, current hydrothermal and hydrogeologic conditions in the caldera, in which there is a large lake and shallow water table, increase the risk of the quiescent volcano.     

Efficient approximations for the simulation of density driven flow in porous media

Due to the non-linear coupling between flow and transport equations the simulation of real density driven flow problems requires a lot of computational time and/or heavy equipments. We suggest some approximations and numerical recipes to reduce the CPU costs for these strongly non-linear coupled equations without loss in accuracy.     

Innovative Contaminant Mass Flux Monitoring in an Aquifer Subject to Tidal Effects

Exposure from groundwater contamination to aquatic receptors residing in receiving surface water is dependent upon the rate of contaminated groundwater discharge. Characterization of groundwater fluxes is challenging, especially in coastal environments where tidal fluctuations result in transient groundwater flows towards these receptors. This can also be further complicated by the high spatial heterogeneity of subsurface deposits enhanced by anthropogenic influences such as the mixing of natural sediments and backfill materials, the presence of subsurface built structures such as sheet pile walls or even occurrence of other sources of contaminant discharge. In this study, the finite volume point dilution method (FVPDM) was successfully used to characterize highly transient groundwater flows and contaminant mass fluxes within a coastal groundwater flow system influenced by marked tides. FVPDM tests were undertaken continuously for more than 48 h at six groundwater monitoring wells, in order to evaluate groundwater flow dynamics during several tide cycles. Contaminant concentrations were measured simultaneously which allowed calculating contaminant mass fluxes. The study highlighted the importance of the aquifer heterogeneity, with groundwater fluxes ranging from 10⁻⁷ to 10⁻³ m/s. Groundwater flux monitoring enabled a significant refinement of the conceptual site model, including the fact that inversion of groundwater fluxes was not observed at high tide. Results indicated that contaminant mass fluxes were particularly higher at a specific monitoring well, by more than three orders of magnitude, than at other wells of the investigated aquifer. This study provided crucial information for optimizing further field investigations and risk mitigation measures.     

Random partitioning and adaptive filters for multiple-point stochastic simulation

Multiple-point geostatistical simulation is used to simulate the spatial structures of geological phenomena. In contrast to conventional two-point variogram based geostatistical methods, the multiple-point approach is capable of simulating complex spatial patterns, shapes, and structures normally observed in geological media. A commonly used pattern based multiple-point geostatistical simulation algorithms is called FILTERSIM. In the conventional FILTERSIM algorithm, the patterns identified in training images are transformed into filter score space using fixed filters that are neither dependent on the training images nor on the characteristics of the patterns extracted from them. In this paper, we introduce two new methods, one for geostatistical simulation and another for conditioning the results. At first, new filters are designed using principal component analysis in such a way to include most structural information specific to the governing training images resulting in the selection of closer patterns in the filter score space. We then propose to combine adaptive filters with an overlap strategy along a raster path and an efficient conditioning method to develop an algorithm for reservoir simulation with high accuracy and continuity. We also combine image quilting with this algorithm to improve connectivity a lot. The proposed method, which we call random partitioning with adaptive filters simulation method, can be used both for continuous and discrete variables. The results of the proposed method show a significant improvement in recovering the expected shapes and structural continuity in the final simulated realizations as compared to those of conventional FILTERSIM algorithm and the algorithm is more than ten times faster than FILTERSIM because of using raster path and using small overlap specially when we use image quilting.     

Methane Occurrences in Aquifers Overlying the Barnett Shale Play with a Focus on Parker County,Texas

Clusters of elevated methane concentrations in aquifers overlying the Barnett Shale play have been the focus of recent national attention as they relate to impacts of hydraulic fracturing. The objective of this study was to assess the spatial extent of high dissolved methane previously observed on the western edge of the play (Parker County) and to evaluate its most likely source. A total of 509 well water samples from 12 counties (14,500 km²) were analyzed for methane, major ions, and carbon isotopes. Most samples were collected from the regional Trinity Aquifer and show only low levels of dissolved methane (85% of 457 unique locations <0.1 mg/L). Methane, when present is primarily thermogenic (δ¹³C 10th and 90th percentiles of ?57.54 and ?39.00‰ and C1/C2+C3 ratio 10th, 50th, and 90th percentiles of 5, 15, and 42). High methane concentrations (>20 mg/L) are limited to a few spatial clusters. The Parker County cluster area includes historical vertical oil and gas wells producing from relatively shallow formations and recent horizontal wells producing from the Barnett Shale (depth of ～1500 m). Lack of correlation with distance to Barnett Shale horizontal wells, with distance to conventional wells, and with well density suggests a natural origin of the dissolved methane. Known commercial very shallow gas accumulations (<200 m in places) and historical instances of water wells reaching gas pockets point to the underlying Strawn Group of Paleozoic age as the main natural source of the dissolved gas.     

Late Weichselian relative sea-level changes and ice sheet history in southeast Greenland

Relative sea-level (RSL) observations from the margins of the Greenland Ice Sheet (GIS) provide information regarding the timing and rate of deglaciation and constraints on geophysical models of ice sheet evolution. In this paper we present the first RSL record for the southeast sector of the GIS based on field observations completed close to Ammassalik. The local marine limit is c. 69 m above sea-level (asl) and is dated to c. 11 k cal. yrs BP (thousand calibrated years before present) and is a minimum date for ice free conditions at the study site. RSL fell to c. 24 m asl by 9.5 k cal. yrs BP and continued to fall at a decreasing rate to reach close to present by 6.5 k cal. yrs BP. Our chronology agrees with radiocarbon dates from offshore cores that indicate ice free conditions on the adjacent mid-shelf by 15 k cal. yrs BP. We compare the new RSL data with predictions generated using two recently published glaciological models of the GIS that differ in the amount and timing of ice loading and unloading over our study area. These two GIS models are coupled to the same Earth viscosity model and background (global) ice model to aid in the data-model comparison. Neither model provides a close fit to the RSL observations. Based on a preliminary sensitivity study using a suite of Earth viscosity models, we conclude that the poor data-model fit is most likely due to an underestimate of the local ice unloading. An improved fit could be achieved by delaying the retreat of a thicker ice sheet across the continental shelf. A thick ice sheet extending well onto the continental shelf is in agreement with other recent observations elsewhere in east and south Greenland.