Microtectonic and geochemical characterization of thrusting in a foreland basin: Example of the South-Pyrenean orogenic wedge (Spain)

In orogenic systems, thrust faults play a major role in stacking different tectonic units and may act as conduits for the expulsion of large amounts of fluid of different origins (metamorphic, diagenetic, meteoric). This study focuses on the Monte Perdido thrust unit emplaced in the Paleogene Jaca thrust-sheet-top basin, in the SW-central Pyrenees. We aim to decipher the mechanisms and P-T conditions of deformation in fault zones and characterize the related fluid involvement, through combined microstructural, geochemical and microthermometry analyses. Two thrust faults cutting platform limestones, marls and siliciclastic turbidites of the lower part of the basin-fill (Paleocene–lower Eocene) have been studied. The fault zones are characterized by metre-thick shear zones with highly deformed, foliated clay-rich sediments. Foliation is underlined by preferentially oriented phyllosilicates. Several generations of shear and extension calcite, quartz and chlorite-bearing veins attest to fluid-rock interactions during a multi-stage deformation. Microstructural observations and stable isotope analyses on calcite from veins and host sediments suggest that deformation was aseismic and dominated by diffusive mass transfer from pressure solution sites along cleavage and stylolites to the precipitation sites in veins, with mineralizing fluids in equilibrium with the host sediments. Our results suggest an essentially closed hydrologic system, and imply the absence of significant fluid flow along the studied fault zones. Microthermometric study on fluid inclusions present in calcite and quartz veins, and calcite-quartz oxygen isotopic fractionation determined for the first generation shear veins, allow a geothermal gradient of 34 °C/km to be estimated. Analytical results demonstrate an evolution of the fault zones in three stages. The first stage was related to the emplacement of the Monte Perdido thrust unit during the middle Eocene at a temperature of ～208 °C and a burial depth of ～5.7 km. The second stage corresponds to a fault reactivation at a temperature of ～240 °C and a burial depth of ～6.5 km. The latter deformation may have been related to folding of the Monte Perdido thrust unit during the emplacement of the underlying Gavarnie thrust unit during the late Eocene–early Oligocene, with deeper burial resulting from aggradation of the thrust-sheet-top basin-fill. The last event corresponds to the formation of a dilatant vein system likely related to the exhumation of the massif.     

Quantifying fracture geometry with X-ray tomography: Technique of Iterative Local Thresholding (TILT) for 3D image segmentation

This paper presents a new method—the Technique of Iterative Local Thresholding (TILT)—for processing 3D X-ray computed tomography (xCT) images for visualization and quantification of rock fractures. The TILT method includes the following advancements. First, custom masks are generated by a fracture-dilation procedure, which significantly amplifies the fracture signal on the intensity histogram used for local thresholding. Second, TILT is particularly well suited for fracture characterization in granular rocks because the multi-scale Hessian fracture (MHF) filter has been incorporated to distinguish fractures from pores in the rock matrix. Third, TILT wraps the thresholding and fracture isolation steps in an optimized iterative routine for binary segmentation, minimizing human intervention and enabling automated processing of large 3D datasets. As an illustrative example, we applied TILT to 3D xCT images of reacted and unreacted fractured limestone cores. Other segmentation methods were also applied to provide insights regarding variability in image processing. The results show that TILT significantly enhanced separability of grayscale intensities, outperformed the other methods in automation, and was successful in isolating fractures from the porous rock matrix. Because the other methods are more likely to misclassify fracture edges as void and/or have limited capacity in distinguishing fractures from pores, those methods estimated larger fracture volumes (up to 80 %), surface areas (up to 60 %), and roughness (up to a factor of 2). These differences in fracture geometry would lead to significant disparities in hydraulic permeability predictions, as determined by 2D flow simulations.     

Variability of organic δC and C/N in the Mersey Estuary, U.K. and its implications for sea-level reconstruction studies

Microfossil analysis (e.g. diatoms, foraminifera and pollen) represents the cornerstone of Holocene relative sea-level (RSL) reconstruction because their distribution in the contemporary inter-tidal zone is principally controlled by ground elevation within the tidal frame. A combination of poor microfossil preservation and a limited range in the sediment record may severely restrict the accuracy of resulting RSL reconstructions. Organic δ¹³C and C/N analysis of inter-tidal sediments have shown some potential as coastal palaeoenvironmental proxies. Here we assess their viability for reconstructing RSL change by examining patterns of organic δ¹³C and C/N values in a modern estuarine environment. δ¹³C and C/N analysis of bulk organic inter-tidal sediments and vegetation, as well as suspended and bedload organic sediments of the Mersey Estuary, U.K., demonstrate that the two main sources of organic carbon to surface saltmarsh sediments (terrestrial vegetation and tidal-derived particulate organic matter) have distinctive δ¹³C and C/N signatures. The resulting relationship between ground elevation within the tidal frame and surface sediment δ¹³C and C/N is unaffected by decompositional changes. The potential of this technique for RSL reconstruction is demonstrated by the analysis of part of an early Holocene sediment core from the Mersey Estuary. Organic δ¹³C and C/N analysis is less time consuming than microfossil analysis and is likely to provide continuous records of RSL change.     

The Meso-Cenozoic thermo-tectonic evolution of the Eastern Pyrenees: an <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar fission track and (U–Th)/He thermochronological study of the Canigou and Mont-Louis massifs

Vertical displacements on the SW–NE Têt fault (Eastern Pyrenees Axial Zone, France), which separates the Variscan Canigou-Carança and Mont-Louis massifs, were constrained using a thermochronologic multi-method approach. ⁴⁰Ar/³⁹Ar data from the granitic Mont-Louis massif record its Variscan cooling history and reveal no ages younger than Early Cretaceous, while the Canigou-Carança gneiss massif records systematically younger ⁴⁰Ar/³⁹Ar ages. These younger ⁴⁰Ar/³⁹Ar ages in the Canigou-Carança gneiss massif are the result of partial to total rejuvenation of argon isotopic systems related to a thermal flow coeval with the Cretaceous HT-BP metamorphism in the North Pyrenean Zone. Only the deepest rocks from the Canigou-Carança suffered this extensive Mid-Cretaceous thermal overprint probably due to differential burial around 4 km at that time. The post Mid-Cretaceous vertical displacements along the Têt fault are recorded by “low” temperature thermochronology using K-feldspar ⁴⁰Ar/³⁹Ar, zircon and apatite fission track and (U–Th)/He datings. The Mont-Louis granite samples experienced a long period of protracted cooling reflecting a lack of thermo-tectonic activity in this area from Late Palaeozoic to Early Cenozoic, followed by cooling from 55–60 Ma to Late Eocene at a mean rate of 15–20°C/Ma in the final stage. This cooling stage corresponds to Têt fault reactivation with a reversed component, promoting exhumation of the Mont-Louis roof zone contemporaneously with the south-vergent Pyrenean thrusting. In the Canigou-Carança massif, the main cooling event occurred from 32 to 18 Ma at a maximum rate of 30°C/Ma during Early Oligocene followed by a more moderate rate of 3°C/Ma from Late Oligocene to Early Burdigalian, coeval with the normal reactivation of the Têt fault in brittle conditions that accommodated the final exhumation of the massif during the opening of the Gulf of Lion.     

Suitability analysis for sand and gravel extraction site location in the context of a sustainable development in the surroundings of Zaragoza (Spain)

Zaragoza city is located in the central Ebro Basin, in the Iberian Peninsula. The fluvial terraces formed by the Ebro River present a valuable resource of sand and gravel deposits. However, taking advantage of these available resources implies conflicts with other land use interests like urban and industrial development as well as agricultural use, which has also traditionally occupied the alluvial terraces. These deposits represent a substantial groundwater resource that should be preserved for future generations. The development of spatial decision support systems (SDSS) has greatly assisted efforts for solving land-use conflicts. These systems combine the benefits of geographic information systems (GIS) and decision support methodologies and are therefore suitable to manage sustainable development of urban areas. In this contribution, an extraction suitability map taking into consideration a variety of environmental criteria is created with the help of a SDSS. The method used is the analytical hierarchy process which is integrated in ArcGIS. Areas most suitable to sand and gravel extraction are located in the high terraces, and in those terraces covered by pediments where the thickness of resource is relatively high. These areas are far from valuable natural areas, outside areas most vulnerable to groundwater contamination, and beneath soils with poor irrigation characteristics.     

High-frequency monitoring of hydrological and biogeochemical fluxes in forested catchments of southern Chile

The variability of rainfall-dependent streamflow at catchment scale modulates many ecosystem processes in wet temperate forests. Runoff in small mountain catchments is characterized by a quick response to rainfall pulses which affects biogeochemical fluxes to all downstream systems. In wet-temperate climates, water erosion is the most important natural factor driving downstream soil and nutrient losses from upland ecosystems. Most hydrochemical studies have focused on water flux measurements at hourly scales, along with weekly or monthly samples for water chemistry. Here, we assessed how water and element flows from broad-leaved, evergreen forested catchments in southwestern South America, are influenced by different successional stages, quantifying runoff, sediment transport and nutrient fluxes during hourly rainfall events of different intensities. Hydrograph comparisons among different successional stages indicated that forested catchments differed in their responses to high intensity rainfall, with greater runoff in areas covered by secondary forests (SF), compared to old-growth forest cover (OG) and dense scrub vegetation (CH). Further, throughfall water was greatly nutrient enriched for all forest types. Suspended sediment loads varied between successional stages. SF catchments exported 455 kg of sediments per ha, followed by OG with 91 kg/ha and CH with 14 kg/ha, corresponding to 11 rainfall events measured from December 2013 to April 2014. Total nitrogen (TN) and phosphorus (TP) concentrations in stream water also varied with rainfall intensity. In seven rainfall events sampled during the study period, CH catchments exported less nutrients (46 kg/ha TN and 7 kg/ha TP) than SF catchments (718 kg/ha TN and 107 kg/ha TP), while OG catchments exported intermediate sediment loads (201 kg/ha TN and 23 kg/ha TP). Further, we found significant effects of successional stage attributes (vegetation structure and soil physical properties) and catchment morphometry on runoff and sediment concentrations, and greater nutrients retention in OG and CH catchments. We conclude that in these southern hemisphere, broad-leaved evergreen temperate forests, hydrological processes are driven by multiple interacting phenomena, including climate, vegetation, soils, topography, and disturbance history.     

Helium isotopes in ferromanganese crusts from the central Pacific Ocean

Helium isotopes have been measured in samples of two ferromanganese crusts (VA13/2 and CD29-2) from the central Pacific Ocean. With the exception of the deepest part of crust CD29-2 the data can be explained by a mixture of implanted solar- and galactic cosmic ray-produced (GCR) He, in extraterrestrial grains, and radiogenic He in wind-borne continental dust grains. ⁴He concentrations are invariant and require retention of less than 12% of the in situ He produced since crust formation. Loss has occurred by recoil and diffusion. High ⁴He in CD29-2 samples older than 42 Ma are correlated with phosphatization and can be explained by retention of up to 12% of the in situ-produced ⁴He. ³He/⁴He of VA13/2 samples varies from 18.5 to 1852 R_a due almost entirely to variation in the extraterrestrial He contribution. The highest ³He/⁴He is comparable to the highest values measured in interplanetary dust particles (IDPs) and micrometeorites (MMs). Helium concentrations are orders of magnitude lower than in oceanic sediments reflecting the low trapping efficiency for in-falling terrestrial and extraterrestrial grains of Fe-Mn crusts. The extraterrestrial ³He concentration of the crusts rules out whole, undegassed 4-40 μm diameter IDPs as the host. Instead it requires that the extraterrestrial He inventory is carried by numerous particles with significantly lower He concentrations, and occasional high concentration GCR-He-bearing particles.     

Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Numerical modeling of three‐phase dissolution of underground cavities using a diffuse interface model

Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.