The East Anatolia-Lesser Caucasus ophiolite:An exceptional case of large-scale obduction,synthesis of data and numerical modelling

We present new,geological,metamorphic,geochemical and geochronological data on the East Anatolian-Lesser Caucasus ophiolites.These data are used in combination with a synthesis of previous data and numerical modelling to unravel the tectonic emplacement of ophiolites in this region.All these data allow the reconstruction of a large obducted ophiolite nappe,thrusted for>100 km and up to 250 km on the Anatolian-Armenian block.The ophiolite petrology shows three distinct magmatic series,highlighted by new isotopic and trace element data:(1)The main Early Jurassic Tholeiites(ophiolite s.s.)bear LILEenriched,subduction-modified,MORB chemical composition.Geology and petrology of the Tholeiite series substantiates a slow-spreading oceanic environment in a time spanning from the Late Triassic to the Middle-Late Jurassic.Serpentinites,gabbros and plagiogranites were exhumed by normal faults,and covered by radiolarites,while minor volumes of pillow-lava flows infilled the rift grabens.Tendency towards a subduction-modified geochemical signature suggests emplacement in a marginal basin above a subduction zone.(2)Late Early Cretaceous alkaline lavas conformably emplaced on top of the ophiolite.They have an OIB affinity.These lavas are featured by large pillow lavas interbedded a carbonate matrix.They show evidence for a large-scale OIB plume activity,which occurred prior to ophiolite obduction.(3)Early-Late Cretaceous calc-alkaline lavas and dykes.These magmatic rocks are found on top of the obducted nappe,above the post-obduction erosion level.This series shows similar Sr-Nd isotopic features as the Alkaline series,though having a clear supra-subduction affinity.They are thus interpreted to be the remelting product of a mantle previously contaminated by the OIB plume.Correlation of data from the Lesser Caucasus to western Anatolia shows a progression from back-arc to arc and fore-arc,which highlight a dissymmetry in the obducted oceanic lithosphere from East to West.The metamorphic P-T-t paths of the obduction sole lithologies define a southward propagation of the ophiolite:(1)P-T-t data from the northern Sevan-Akera suture zone(Armenia)highlight the presence and exhumation of eclogites(1.85±0.02 GPa and 590±5℃)and blueschists below the ophiolite,which are dated at ca.94 Ma by Ar-Ar on phengite.(2)Neighbouring Amasia(Armenia)garnet amphibolites indicate metamorphic peak conditions of 0.65±0.05 GPa and 600±20 C with a U-Pb on rutile age of 90.2±5.2 Ma and Ar-Ar on amphibole and phengite ages of 90.8±3.0 Ma and 90.8±1.2 Ma,respectively.These data are consistent with palaeontological dating of sediment deposits directly under(Cenomanian,i.e.>93.9 Ma)or sealing(Coniacian-Santonian,i.e.,≤89.8 Ma),the obduction.(3)At Hinis(NE Turkey)PT-t conditions on amphibolites(0.66±0.06 GPa and 660±20℃,with a U-Pb titanite age of80.0±3.2 Ma)agree with previous P-T-t data on granulites,and highlight a rapid exhumation below a top-to-the-North detachment sealed by the Early Maastrichtian unconformity(ca.70.6 Ma).Amphibolites are cross-cut by monzonites dated by U-Pb on titanite at 78.3±3.7 Ma.We propose that the HT-MP metamorphism was coeval with the monzonites,about 10 Ma after the obduction,and was triggered by the onset of subduction South of the Anatolides and by reactivation or acceleration of the subduction below the Pontides-Eurasian margin.Numerical modelling accounts for the obduction of an"old"~80 Myr oceanic lithosphere due to a significant heating of oceanic lithosphere through mantle upwelling,which increased the oceanic lithosphere buoyancy.The long-distance transport of a currently thin section of ophiolites(<1 km)onto the Anatolian continental margin is ascribed to a combination of northward mantle extensional thinning of the obducted oceanic lithosphere by the Hinis detachment at ca.80 Ma,and southward gravitational propagation of the ophiolite nappe onto its foreland basin.     

Physical characterization of core samples recovered from Gulf of Mexico

Seventy whole rounds from conventional cores obtained during drilling to 300 mbsf at Atwater Valley and Keathley Canyon in the Gulf of Mexico in April and May 2005 were tested to determine geophysical and geomechanical parameters (liquid and plastic limit, porosity, specific surface, pH, sediment electrical conductivity, P- and S-wave velocities and undrained shear strength). Available data from a pressure core are included as well. Results show that the sediments are high specific surface plastic clays, and exhibit pronounced time-dependent stiffness recovery. Strains during coring disturb specimens, yet, the water content retains the effective stress history and permits gaining stiffness and strength information from conventional cores. Remolding is exacerbated when gas expands upon decompression; the limited pressure core data available show the advantages of preserving the pore fluid pressure during core recovery and testing. Valuable parameters for sediment characterization and engineering analysis are extracted from the data using pre-existing soil models.     

Development of a grid‐based version of the SWAT landscape model

Integrated river basin models should provide a spatially distributed representation of basin hydrology and transport processes to allow for spatially implementing specific management and conservation measures. To accomplish this, the Soil and Water Assessment Tool (SWAT) was modified by integrating a landscape routing model to simulate water flow across discretized routing units. This paper presents a grid‐based version of the SWAT landscape model that has been developed to enhance the spatial representation of hydrology and transport processes. The modified model uses a new flow separation index that considers topographic features and soil properties to capture channel and landscape flow processes related to specific landscape positions. The resulting model is spatially fully distributed and includes surface, lateral and groundwater fluxes in each grid cell of the watershed. Furthermore, it more closely represents the spatially heterogeneous distributed flow and transport processes in a watershed. The model was calibrated and validated for the Little River Watershed (LRW) near Tifton, Georgia (USA). Water balance simulations as well as the spatial distribution of surface runoff, subsurface flow and evapotranspiration are examined. Model results indicate that groundwater flow is the dominant landscape process in the LRW. Results are promising, and satisfactory output was obtained with the presented grid‐based SWAT landscape model. Nash–Sutcliffe model efficiencies for daily stream flow were 0.59 and 0.63 for calibration and validation periods, and the model reasonably simulates the impact of the landscape position on surface runoff, subsurface flow and evapotranspiration. Additional revision of the model will likely be necessary to adequately represent temporal variations of transport and flow processes in a watershed. Copyright © 2014 John Wiley & Sons, Ltd.     

Geoelectrical Classification of Gypsum Rocks

Gypsum rocks are widely exploited in the world as industrial minerals. The purity of the gypsum rocks (percentage in gypsum mineral in the whole rock) is a critical factor to evaluate the potential exploitability of a gypsum deposit. It is considered than purities higher than 80% in gypsum are required to be economically profitable. Gypsum deposits have been studied with geoelectrical methods; a direct relationship between the electrical resistivity values of the gypsum rocks and its lithological composition has been established, with the presence of lutites being the main controlling factor in the geoelectrical response of the deposit. This phenomenon has been quantified in the present study, by means of a combination of theoretical calculations, laboratory measurements and field data acquisition. Direct modelling has been performed; the data have been inverted to obtain the mean electrical resistivity of the models. The laboratory measurements have been obtained from artificial gypsum-clay mixture pills, and the electrical resistivity has been measured using a simple electrical circuit with direct current power supply. Finally, electrical resistivity tomography data have been acquired in different evaporite Tertiary basins located in North East Spain; the selected gypsum deposits have different gypsum compositions. The geoelectrical response of gypsum rocks has been determined by comparing the resistivity values obtained from theoretical models, laboratory tests and field examples. A geoelectrical classification of gypsum rocks defining three types of gypsum rocks has been elaborated: (a) Pure Gypsum Rocks (>75% of gypsum content), (b) Transitional Gypsum Rocks (75–55%), and (c) Lutites and Gypsum-rich Lutites (<55%). From the economic point of view, the Pure Gypsum Rocks, displaying a resistivity value of >800 ohm.m, can be exploited as industrial rocks. The methodology used could be applied in other geoelectrical rock studies, given that this relationship between the resistive particles embedded within a conductive matrix depends on the connectivity of the matrix particles.     

The extreme initial kinetic energy allowed by a collapsing turbulent core

We present high-resolution hydrodynamical simulations aimed at following the gravitational collapse of a gas core, in which a turbulent spectrum of velocity is implemented only initially. We determine the maximal value of the ratio of kinetic energy to gravitational energy, denoted here by \((\frac{E_{\mathrm{kin}} }{E_{\mathrm{grav}}} )_{\max}\), so that the core (i) will collapse around one free-fall time of time evolution or (ii) will expand unboundedly, because it has a value of \(\frac{E_{\rm kin}}{E_{\mathrm{grav}}}\) larger than \(( \frac{E_{\mathrm{kin}}}{E_{\mathrm{grav}}} )_{\mathrm{max}}\). We consider core models with a uniform or centrally condensed density profile and with velocity spectra composed of a linear combination of one-half divergence-free turbulence type and the other half of a curl-free turbulence type. We show that the outcome of the core collapse are protostars forming either (i) a multiple system obtained from the fragmentation of filaments and (ii) a single primary system within a long filament. In addition, some properties of these protostars are also determined and compared with those obtained elsewhere.