Landscape Evolution Modelling of naturally dammed rivers

Natural damming of upland river systems, such as landslide or lava damming, occurs worldwide. Many dams fail shortly after their creation, while other dams are long‐lived and therefore have a long‐term impact on fluvial and landscape evolution. This long‐term impact is still poorly understood and landscape evolution modelling (LEM) can increase our understanding of different aspects of this response. Our objective was to simulate fluvial response to damming, by monitoring sediment redistribution and river profile evolution for a range of geomorphic settings. We used LEM LAPSUS, which calculates runoff erosion and deposition and can deal with non‐spurious sinks, such as dam‐impounded areas. Because fluvial dynamics under detachment‐limited and transport‐limited conditions are different, we mimicked these conditions using low and high erodibility settings, respectively. To compare the relative impact of different dam types, we evaluated five scenarios for each landscape condition: one scenario without a dam and four scenarios with dams of increasing erodibility. Results showed that dam‐related sediment storage persisted at least until 15 000 years for all dam scenarios. Incision and knickpoint retreat occurred faster in the detachment‐limited landscape than in the transport‐limited landscape. Furthermore, in the transport‐limited landscape, knickpoint persistence decreased with increasing dam erodibility. Stream capture occurred only in the transport‐limited landscape due to a persisting floodplain behind the dam and headward erosion of adjacent channels. Changes in sediment yield variation due to stream captures did occur but cannot be distinguished from other changes in variation of sediment yield. Comparison of the model results with field examples indicates that the model reproduces several key phenomena of damming response in both transport‐limited and detachment‐limited landscapes. We conclude that a damming event which occurred 15 000 years ago can influence present‐day sediment yield, profile evolution and stream patterns. Copyright © 2014 John Wiley & Sons, Ltd.     

Marble-hosted ruby deposits from Central and Southeast Asia: Towards a new genetic model

Marble-hosted ruby deposits represent the most important source of colored gemstones from Central and South East Asia. These deposits are located in the Himalayan mountain belt which developed during Tertiary collision of the Indian plate northward into the Eurasian plate. They are spatially related to granitoid intrusions and are contained in platform carbonates series that underwent high-grade metamorphism. All occurrences are located close to major tectonic features formed during Himalayan orogenesis, directly in suture zones in the Himalayas, or in shear zones that guided extrusion of the Indochina block after the collision in South East Asia. Ar–Ar dating of micas syngenetic with ruby and U–Pb dating of zircon included in ruby gives evidence that these deposits formed during Himalayan orogenesis, and the ages document the extensional tectonics that were active, from Afghanistan to Vietnam, between the Oligocene and the Pliocene.The petrography shows that ruby-bearing marbles formed in the amphibolite facies (T = 610 to 790 °C and P ~ 6 kbar). A fluid inclusion study defines the conditions of gem ruby formation during the retrograde metamorphic path (620 < T < 670 °C and 2.6 < P < 3.3 kbar) for the deposits of Jegdalek, Hunza and northern Vietnam.Whole rock analyses of non-ruby-bearing marbles indicate that they contain enough aluminum and chromiferous elements to produce all the ruby crystals that they contain. In addition, (C, O)-isotopic analyses of carbonates from the marbles lead to the conclusion that the marbles acted as a metamorphic closed fluid system that were not infiltrated by externally-derived fluids. The carbon isotopic composition of graphite in marbles reveals that it is of organic origin and that it exchanged C-isotopes with the carbonates during metamorphism. Moreover, the O-isotopic composition of ruby was buffered by metamorphic CO₂ released during devolatilisation of marble and the H-isotopic composition of mica is consistent with a metamorphic origin for water in equilibrium with the micas. The (C, O, H)-isotopic compositions of minerals associated with marble-hosted ruby are all in agreement with the hypothesis, drawn from the unusual chemistry of CO₂–H₂S–COS–S₈–AlO(OH)-bearing fluids contained in fluid inclusions, that gem ruby formed at P ~ 3 kbar and 620 < T < 670 °C, during thermal reduction of evaporite by organic matter, at high temperature-medium pressure metamorphism of platform carbonates during the Tertiary India–Asia collision. The carbonates were enriched in Al- and chromiferous-bearing detrital minerals, such as clay minerals that were deposited on the platform with the carbonates, and in organic matter. Ruby formed during the retrograde metamorphic path, mainly by destabilization of muscovite or spinel. The metamorphic fluid system was rich in CO₂ released from devolatilisation of carbonates, and in fluorine, chlorine and boron released by molten salts (NaCl, KCl, CaSO₄). Evaporites are key to explaining the formation of these deposits. Molten salts mobilized in situ Al and metal transition elements contained in marbles, leading to crystallization of ruby.     

Zooplankton and forage fish species off Peru: Large-scale bottom-up forcing and local-scale depletion

The Humboldt Current System, like all upwelling systems, has dramatic quantities of plankton-feeding fish, which suggested that their population dynamics may ‘drive’ or ‘control’ ecosystem dynamics. With this in mind we analysed the relationship between forage fish populations and their main prey, zooplankton populations. Our study combined a zooplankton sampling program (1961–2005) with simultaneous acoustic observations on fish from 40 pelagic surveys (1983–2005) conducted by the Peruvian Marine Research Institute (IMARPE) and landing statistics for anchoveta (Engraulis ringens) and sardine (Sardinops sagax) along the Peruvian coast from 1961 to 2005. The multi-year trend of anchoveta population abundance varied consistently with zooplankton biovolume trend, suggesting bottom-up control on anchovy at the population scale (since oceanographic conditions and phytoplankton production support the changes in zooplankton abundance). For a finer-scale analysis (km) we statistically modelled zooplankton biovolume as a function of geographical (latitude and distance from the 200-m isobath), environmental (sea surface temperature), temporal (year, month and time-of-day) and biological (acoustic anchovy and sardine biomass within 5 km of each zooplankton sample) covariates over all survey using both classification and regression trees (CART) and generalized additive models (GAM). CART showed local anchoveta density to have the strongest effect on zooplankton biovolume, with significantly reduced levels of biovolume for higher neighbourhood anchoveta biomass. Additionally, zooplankton biovolume was higher offshore than on the shelf. GAM results corroborated the CART findings, also showing a clear diel effect on zooplankton biovolume, probably due to diel migration or daytime net avoidance. Apparently, the observed multi-year population scale bottom-up control is not inconsistent with local depletion of zooplankton when anchoveta are locally abundant, since the latter effect was observed over a wide range of overall anchoveta abundance.     

Mesoscale eddies off Peru in altimeter records: Identification algorithms and eddy spatio-temporal patterns

Relatively little is known about coherent vortices in the eastern South-Pacific along the Peruvian coast, even with regard to basic facts about their frequency of occurrence, longevity and structure. This study addresses these issues with nearly 15 years of relatively high-resolution satellite altimetry measurements.We first compare two distinct automated methods for eddy identification. The objective validation protocol shows that the rarely-used geometrical or “winding-angle method”, based on the curvature of the streamline functions, is more accurate than the commonly-used “Okubo–Weiss algorithm”, which defines a vortex as a simple connected region with values of Okubo–Weiss parameter weaker than a given threshold.We then investigate vortices off Peru using more than 20,000 mesoscale eddies identified by the winding-angle method. Coherent eddies, characterized by a high ratio of vorticity to deformation rate, are typically formed along the coast and propagate westward at 3–6 cm s⁻¹. The vortices have a mean radius of 80 km, increasing northward, and are most frequently observed off of Chimbote (9°S) and south of San Juan (15°S). The mean eddy lifetime is about 1 month, but if eddies survive at least 2 months, the probability for surviving an additional week (or month) is constant at 90% (or 67%). Anticyclonic eddies tend to propagate northwestward whereas cyclonic vortices migrate southwestward. In general, cyclones and anticyclones are similar, except for eddies surviving at least 6 months. In this case, after a similar 3–4 months of radius and amplitude growth, amplitudes (or sizes) decay particularly rapidly for anticyclonic (or cyclonic) eddies. In terms of intensity, cyclonic eddies show a rapid decay during the first 3 months before arriving at a quasi-constant value, whereas anticyclones exhibit steady decline. Finally, eddy temporal variations were examined at seasonal and interannual scales in the “coastal” region favorable to the formation of energetic mesoscale structures. On seasonal scales, eddy activity is maximal in fall and minimum in spring. At interannual scales, the eddy activity index was maximal during the strong El Niño of 1997–1998 but another strong maximum of eddy activity also occurred late in 2004. These temporal variations are probably associated with the intensification of the upwelling thermal front and with the passage of coastal-trapped waves which generate baroclinic instabilities. Further investigation of the mechanisms involved on the eddy genesis is needed.     

Evidence for excess argon during high pressure metamorphism in the dora maira massif (Western Alps,Italy), using an ultra-violet laser ablation microprobe 40Ar-39Ar technique

 Ultra-high pressure eclogite/amphibolite grade metamorphism of the Dora Maira Massif in the western Alps is a well established and intensively studied event. However, the age of peak metamorphism and early cooling remains controversial. The ⁴⁰Ar-³⁹Ar step-heating and laser spot ages from high pressure phengites yield plateau ages as old as 110 Ma which have been interpreted as the time of early cooling after the high pressure event. Recent U/Pb and Sm/Nd results challenge this assertion, indicating a much younger age for the event, around 45 Ma, and hence a radically different timing for the tectonic evolution of the western Alps. In a new approach to the problem, samples from the undeformed Hercynian metagranite, Brossasco, were studied using an ultra-violet laser ablation microprobe technique for ⁴⁰Ar-³⁹Ar dating. The new technique allowed selective in situ analysis, at a spatial resolution of 50 μm, of quartz, phengite, biotite and K-feldspar. The results demonstrate the frequent occurrence of excess argon with high ⁴⁰Ar-³⁶Ar ratios (1000–10000) and a strong relationship between apparent ages and metamorphic textures. The highest excess argon ratios are always associated with high closure temperature minerals or large diffusion domains within single mineral phases. The best interpretation of this relationship seems to be that excess argon was incorporated in all phases during the high pressure event, then mixed with an atmospheric component during rapid cooling and retrogression, producing a wide range of argon concentrations and ⁴⁰Ar/³⁶Ar ratios. Step-heating analysis of minerals with this mixture would produce linear arrays on a ³⁶Ar/⁴⁰Ar versus ³⁹Ar/⁴⁰Ar correlation diagram, leading to geologically meaningless plateau ages, older than the true closure age. In the present case, some ages in the range 60–110 Ma could be explained by the presence of excess argon incorporated around 40–50 Ma ago. Similar results found in other high pressure terrains in the Alps may reconcile the argon geochronometer with other systems such as Rb/Sr, U/Pb or Sm/Nd. This study therefore calls for an increasing use of high resolution in situ sampling techniques to clarify the meaning of ⁴⁰Ar/³⁹Ar ages in many high pressure terrains. Received: 6 January 1994/Accepted: 4 April 1995     

The quantification of low-probability–high-consequences events: part I. A generic multi-risk approach

Dynamic risk processes, which involve interactions at the hazard and risk levels, have yet to be clearly understood and properly integrated into probabilistic risk assessment. While much attention has been given to this aspect lately, most studies remain limited to a small number of site-specific multi-risk scenarios. We present a generic probabilistic framework based on the sequential Monte Carlo Method to implement coinciding events and triggered chains of events (using a variant of a Markov chain), as well as time-variant vulnerability and exposure. We consider generic perils based on analogies with real ones, natural and man-made. Each simulated time series corresponds to one risk scenario, and the analysis of multiple time series allows for the probabilistic assessment of losses and for the recognition of more or less probable risk paths, including extremes or low-probability–high-consequences chains of events. We find that extreme events can be captured by adding more knowledge on potential interaction processes using in a brick-by-brick approach. We introduce the concept of risk migration matrix to evaluate how multi-risk participates to the emergence of extremes, and we show that risk migration (i.e., clustering of losses) and risk amplification (i.e., loss amplification at higher losses) are the two main causes for their occurrence.     

Protracted,coeval crust and mantle melting during Variscan late-orogenic evolution: U–Pb dating in the eastern French Massif Central

International Journal of Earth Sciences - The late stages of the Variscan orogeny are characterized by middle to lower crustal melting and intrusion of voluminous granitoids throughout the belt,...     

UPb emplacement and 40Ar/39Ar cooling ages of the eastern Mont-Louis granite massif (Eastern Pyrenees,France)

Zircon UPb dating by SIMS of the Mont-Louis granite yields an age of 305±5 Ma, intrepreted to reflect the igneous emplacement age of the massif. It is in agreement with the Hercynian syntectonic character of Pyrenees granite. ⁴⁰Ar/³⁹Ar on hornblende, biotite and K-feldspar permit, to estimate the massif cooling. A rapid temperature decrease (≈30 °C/Ma) is revealed from Westphalian to Late Stephanian, coeval with the emplacement of a laccolithe in the upper crust. Then, the cooling rate decreases to ≈1 °C/Ma. This would be consistent with a long time residence for the pluton from the Late Palaeozoic to the Early Cainozoic at 6–8 km depth. To cite this article: O. Maurel et al., C. R. Geoscience 336 (2004).     

Influence of rainfall spatial variability on flood prediction

This paper deals with the sensitivity of distributed hydrological models to different patterns that account for the spatial distribution of rainfall: spatially averaged rainfall or rainfall field. The rainfall data come from a dense network of recording rain gauges that cover approximately 2000 km² around Mexico City. The reference rain sample accounts for the 50 most significant events, whose mean duration is about 10 h and maximal point depth 170 mm. Three models were tested using different runoff production models: storm-runoff coefficient, complete or partial interception. These models were then applied to four fictitious homogeneous basins, whose sizes range from 20 to 1500 km². For each test, the sensitivity of the model is expressed as the relative differences between the empirical distribution of the peak flows (and runoff volumes), calculated according to the two patterns of rainfall input: uniform or non-uniform. Differences in flows range from 10 to 80%, depending on the type of runoff production model used, the size of the basin and the return period of the event. The differences are generally moderate for extreme events. In the local context, this means that uniform design rainfall combining point rainfall distribution and the probabilistic concept of the areal reduction factor could be sufficient to estimate major flood probability. Differences are more significant for more frequent events. This can generate problems in calibrating the hydrological model when spatial rainfall localization is not taken into account: a bias in the estimation of parameters makes their physical interpretation difficult and leads to overestimation of extreme flows.