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.     

External controls on the distribution,fabrics and mineralization of modern microbial mats in a coastal hypersaline lagoon,Cayo Coco (Cuba)

Active, carbonate‐mineralizing microbial mats flourish in a tropical, highly evaporative, marine‐fed lagoonal network to the south of Cayo Coco Island (Cuba). Hypersaline conditions support the development of a complex sedimentary microbial ecosystem with diverse morphologies, a variable intensity of mineralization and a potential for preservation. In this study, the role of intrinsic (i.e. microbial) and extrinsic (i.e. physicochemical) controls on microbial mat development, mineralization and preservation was investigated. The network consists of lagoons, forming in the interdune depressions of a Pleistocene aeolian substratum; they developed due to a progressive increase in sea‐level since the Holocene. The hydrological budget in the Cayo Coco lagoonal network changes from west to east, increasing the salinity. This change progressively excludes grazers and increases the saturation index of carbonate minerals, favouring the development and mineralization of microbial mats in the easternmost lagoons. Detailed mapping of the easternmost lagoon shows four zones with different flooding regimes. The microbial activity in the mats was recorded using light–dark shifts in conjunction with microelectrode O₂ and HS^? profiles. High rates of O₂ production and consumption, in addition to substantial amounts of exopolymeric substances, are indicative of a potentially strong intrinsic control on mineralization. Seasonal, climate‐driven water fluctuations are key for mat development, mineralization, morphology and distribution. Microbial mats show no mineralization in the permanently submersed zone, and moderate mineralization in zones with alternating immersion and exposure. It is suggested that mineralization is also driven by water‐level fluctuations and evaporation. Mineralized mats are laminated and consist of alternating trapping and binding of grains and microbially induced magnesium calcite and dolomite precipitation. The macrofabrics of the mats evolve from early colonizing Flat mats to complex Cerebroid or Terrace structures. The macrofabrics are influenced by the hydrodynamic regime: wind‐driven waves inducing relief terraces in windward areas and flat morphologies on the leeward side of the lagoon. Other external drivers include: (i) storm events that either promote (for example, by bioclasts covering) or prevent (for example, by causing erosion) microbial mat preservation; and (ii) subsurface degassing, through mangrove roots and desiccation cracks covered by Flat mats (i.e. forming Hemispheroids and Cerebroidal structures). These findings provide in‐depth insights into understanding fossil microbialite morphologies that formed in lagoonal settings.     

Monte-Carlo Simulations of Cosmic-ray and Internal Radiation Effects on ISOCAM on Board ISO

One of the main limitations to the sensitivity of the infrared camera ISOCAM on-board the Infrared Space Observatory (ISO) comes from responsivity variations and glitches caused by the impacts of charged particles in photo-detectors. Glitch rate measurements, glitch properties and removal methods have already been addressed during the first ISO detector workshop(Madrid, 1998) and published in a special issue of Experimental Astronomy. It appeared that glitch rate and most of glitch properties could be reproduced by Monte-Carlo simulations. This is very interesting in order to predict before launch the effect of charged particles in photo-detectors operated in space. This paper presents results of Monte-Carlo simulations of radiation effects on ISOCAM detectors. Glitch rates, spatial and energetic properties of glitches have been computed and are compared with measured values. This revised version was published online in July 2006 with corrections to the Cover Date.     

On different types of adjustment usable to calculate the parameters of the stream power law

Model parameterization through adjustment to field data is a crucial step in the modeling and the understanding of the drainage network response to tectonic or climatic perturbations. Using as a test case a data set of 18 knickpoints that materialize the migration of a 0.7-Ma-old erosion wave in the Ourthe catchment of northern Ardennes (western Europe), we explore the impact of various data fitting on the calibration of the stream power model of river incision, from which a simple knickpoint celerity equation is derived. Our results show that statistical least squares adjustments (or misfit functions) based either on the stream-wise distances between observed and modeled knickpoint positions at time t or on differences between observed and modeled time at the actual knickpoint locations yield significantly different values for the m and K parameters of the model. As there is no physical reason to prefer one of these approaches, an intermediate least-rectangles adjustment might at first glance appear as the best compromise. However, the statistics of the analysis of 200 sets of synthetic knickpoints generated in the Ourthe catchment indicate that the time-based adjustment is the most capable of getting close to the true parameter values. Moreover, this fitting method leads in all cases to an m value lower than that obtained from the classical distance adjustment (for example, 0.75 against 0.86 for the real case of the Ourthe catchment), corresponding to an increase in the non-linear character of the dependence of knickpoint celerity on discharge.     

Snow cover monitoring in the Northern Patagonia Icefield using MODIS satellite images (2000–2006)

The snow cover of the Northern Patagonia Icefield (NPI) was monitored after applying the Normalized Difference Snow Index (NDSI) and the Red/NIR band ratio to 134 Moderate Resolution Imaging Spectroradiometer (MODIS) images captured between 2000 and 2006. The final results show that the snow cover extent of the NPI fluctuates a lot in winter, in addition to its seasonal behaviour. The minimum snow cover extent of the period (3600 km²) was observed in March 2000 and the maximum (11,623 km²) in August 2001. We found that temperature accounts for approximately 76% of the variation of the snow cover extent over the entire icefield. We also show two different regimes of winter snow cover fluctuations corresponding to the eastern and the western sides of the icefield. The seasonality of the snow cover on the western side was determined by temperature rather than precipitation, while on the east side the seasonality of the snow cover was influenced by the seasonal behaviour of both temperature and precipitation. This difference can be explained by the two distinct climates: coastal and continental. The fluctuations in the winter snow cover extent were more pronounced and less controlled by temperature on the western side than on the eastern side of the icefield. Snow cover extent was correlated with temperature R² = 0.75 and R² = 0.74 for the western and eastern sides, respectively. Since limited meteorological data are available in this region, our investigation confirmed that the change in snow cover is an interesting climatic indicator over the NPI providing important insights in mass balance comprehension. Since snow and ice were distinguished snow cover fluctuations can be associated to fluctuations in the snow accumulation area of the NPI. In addition, days with minimum snow covers of summer season can be associated to the period in which Equilibrium Line Altitude (ELA) is the highest.     

New insight on the genus Nolaniceras Casey, 1961 (Ammonoidea,Cretaceous) and its consequences on the biostratigraphy of the Aptian Stage

Nolaniceras nolani ( Seunes, 1887) has been widely quoted in the Upper Aptian literature over the years. Re-examination of the holotype of the species shows that it has always been misidentified and that this taxon is represented by a single specimen, the holotype, of uncertain age. As a consequence, its use as an index species for the Upper Aptian should be reconsidered and abandoned.     

Leaching of lead metallurgical slags and pollutant mobility far from equilibrium conditions

Lead metallurgical slags are partially vitrified materials containing residual amounts of Zn, Pb, Cr, Cd and As. These hazardous materials are generally buried on heaps exposed to weathering. In this study, leaching behavior of lead blast furnace slags has been tested using pure water and open flow experiments. It appears that in such far from equilibrium and slightly acidic conditions, the main phase to be altered is the vitreous phase. As for lunar, basaltic and nuclear glasses, alkalis/proton exchanges prevail and lead to the formation of a non-protective altered layer enriched in Si, Fe and Al. The composition of the altered layer is quite constant except for Si whose concentration decreases towards the leachate interface. Owing to their sizes, micrometric Pb droplets are not always totally dissolved at the slag surface. Nevertheless, nanometric Pb droplets are instantaneously dissolved while a surrounding altered layer is formed. This leads to high Pb releases in open flow systems. Leachate chemistry and dissolution rates of the vitreous phase are closely comparable to previous leaching tests with basaltic and nuclear glasses in conditions far from equilibrium. Moreover, this study confirms that Fe is a stable element in such conditions.     

Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.     

Implementing landslide path dependency in landslide susceptibility modelling

Landslide susceptibility modelling—a crucial step towards the assessment of landslide hazard and risk—has hitherto not included the local, transient effects of previous landslides on susceptibility. In this contribution, we implement such transient effects, which we term “landslide path dependency”, for the first time. Two landslide path dependency variables are used to characterise transient effects: a variable reflecting how likely it is that an earlier landslide will have a follow-up landslide and a variable reflecting the decay of transient effects over time. These two landslide path dependency variables are considered in addition to a large set of conditioning attributes conventionally used in landslide susceptibility. Three logistic regression models were trained and tested fitted to landslide occurrence data from a multi-temporal landslide inventory: (1) a model with only conventional variables, (2) a model with conventional plus landslide path dependency variables, and (3) a model with only landslide path dependency variables. We compare the model performances, differences in the number, coefficient and significance of the selected variables, and the differences in the resulting susceptibility maps. Although the landslide path dependency variables are highly significant and have impacts on the importance of other variables, the performance of the models and the susceptibility maps do not substantially differ between conventional and conventional plus path dependent models. The path dependent landslide susceptibility model, with only two explanatory variables, has lower model performance, and differently patterned susceptibility map than the two other models. A simple landslide susceptibility model using only DEM-derived variables and landslide path dependency variables performs better than the path dependent landslide susceptibility model, and almost as well as the model with conventional plus landslide path dependency variables—while avoiding the need for hard-to-measure variables such as land use or lithology. Although the predictive power of landslide path dependency variables is lower than those of the most important conventional variables, our findings provide a clear incentive to further explore landslide path dependency effects and their potential role in landslide susceptibility modelling.     

Do landslides follow landslides? Insights in path dependency from a multi-temporal landslide inventory

Landslides are a major category of natural disasters, causing loss of lives, livelihoods and property. The critical roles played by triggering (such as extreme rainfall and earthquakes), and intrinsic factors (such as slope steepness, soil properties and lithology) have previously successfully been recognized and quantified using a variety of qualitative, quantitative and hybrid methods in a wide range of study sites. However, available data typically do not allow to investigate the effect that earlier landslides have on intrinsic factors and hence on follow-up landslides. Therefore, existing methods cannot account for the potentially complex susceptibility changes caused by landslide events. In this study, we used a substantially different alternative approach to shed light on the potential effect of earlier landslides using a multi-temporal dataset of landslide occurrence containing 17 time slices. Spatial overlap and the time interval between landslides play key roles in our work. We quantified the degree to which landslides preferentially occur in locations where landslides occurred previously, how long such an effect is noticeable, and how landslides are spatially associated over time. We also investigated whether overlap with previous landslides causes differences in landslide geometric properties. We found that overlap among landslides demonstrates a clear legacy effect (path dependency) that has influence on the landslide affected area. Landslides appear to cause greater susceptibility for follow-up landslides over a period of about 10  years. Follow-up landslides are on average larger and rounder than landslides that do not follow earlier slides. The effect of earlier landslides on follow-up landslides has implications for understanding of the landslides evolution and the assessment of landslide susceptibility.