A 17,900-year multi-proxy lacustrine record of Lago Puyehue (Chilean Lake District): introduction

This paper introduces the background and main results of a research project aimed at unravelling the paleolimnological and paleoclimatological history of Lago Puyehue (40° S, Lake District, Chile) since the Last Glacial Maximum (LGM), based on the study of several sediment cores from the lake and on extensive fieldwork in the lake catchment. The longest record was obtained in an 11-m-long piston core. An age-depth model was established by AMS ¹⁴C dating, ²¹⁰Pb and ²³⁷Cs measurements, identification of event-deposits, and varve-counting for the past 600 years. The core extends back to 17,915 cal. yr. BP, and the seismic data indicate that an open-lake sedimentary environment already existed several thousands of years before that. The core was submitted to a multi-proxy analysis, including sedimentology, mineralogy, grain-size, major geochemistry and organic geochemistry (C/N ratio, δ¹³C), loss-on-ignition, magnetic susceptibility, diatom analysis and palynology. Along-core variations in sediment composition reveal that the area of Lago Puyehue was characterized since the LGM by a series of rapid climate fluctuations superimposed on a long-term warming trend. Identified climate fluctuations confirm a.o. the existence of a Late-Glacial cold reversal predating the northern-hemisphere Younger Dryas cold period by 500–1,000 years, as well as the existence of an early southern-hemisphere Holocene climatic optimum. Varve-thickness analyses over the past 600 years reveal periodicities similar to those associated with the El Niño Southern Oscillation and the Pacific Decadal Oscillation, as well as intervals with increased precipitation, related to an intensification of the El Niño impact during the southern-hemisphere equivalent of the Little Ice Age.     

Mineral Paragenesis and Ore‐Forming Processes of the Dongping Gold Deposit,Hebei Province,China

The Dongping gold deposit is located near the center of the northern margin of the North China Craton. It is hosted in the Shuiquangou syenite and characterized by large amounts of tellurides. Numerous studies have addressed this deposit; the mineral paragenesis and ore‐forming processes, however, are still poorly studied. In this contribution, a new mineral paragenesis has been evaluated to further understand ore formation, including sulfides (pyrite, chalcopyrite, galena, sphalerite, molybdenite, and bornite), tellurides (altaite, calaverite, hessite, muthmannite, petzite, rucklidgeite, sylvanite, tellurobismuthite, tetradymite, and volynskite), and native elements (tellurium and gold). Molybdenite, muthmannite, rucklidgeite, and volynskite are reported for the first time in this deposit. We consider the Dongping gold deposit mainly formed in the Devonian, and the ore‐forming processes and the physicochemical conditions for ore formation can be reconstructed based on our newly identified ore paragenesis, that is, iron oxides → (CO₂ effervescence) → sulfides → (fTe₂/fS₂ ratio increase) → Pb‐Bi‐tellurides → (condensation of H₂Te vapor) → Au‐Ag‐tellurides → (mixing with oxidizing water) → carbonate and microporous gold → secondary minerals → secondary minerals. The logfO₂ values increase from the early to late stages, while the fH₂S and logfS₂ values increase initially and then decrease. CO₂ effervescence is the main mechanism of sulfides precipitation; this sulfidation and condensation of H₂Te vapor lead to deposition of tellurides. The development of microporous gold indicates that the deposit might experience overprint after mineralization. The Dongping gold deposit has a close genetic relationship with the Shuiquangou syenite, and tellurium likely originated from Shuiquangou alkaline magmatic degassing.     

Giant submarine collapse of a carbonate platform at the Turonian–Coniacian transition: The Ayabacas Formation, southern Peru

The Ayabacas Formation of southern Peru is an impressive unit formed by the giant submarine collapse of the mid‐Cretaceous carbonate platform of the western Peru back‐arc basin (WPBAB), near the Turonian–Coniacian transition (～90–89 Ma). It extends along the southwestern edge of the Cordillera Oriental and throughout the Altiplano and Cordillera Occidental over >80 000 km² in map view, and represents a volume of displaced sediments of >10 000 km³. The collapse occurred down the basin slope, i.e. toward the SW. Six zones are characterised on the basis of deformational facies, and a seventh corresponds to the northeastern ‘stable’ area (Zone 0). Zones 1–3 display increasing fragmentation from NE to SW, and are composed of limestone rafts and sheets embedded in a matrix of mainly red, partly calcareous and locally sandy, mudstones to siltstones. In contrast, in Zones 4 and 5 the unit consists only of displaced and stacked limestone masses forming a ‘sedimentary thrust and fold system’, with sizes increasing to the southwest. In Zone 6, the upper part of the limestone succession consists of rafts and sheets stacked over the regularly bedded lower part. The triggering of this extremely large mass wasting clearly ensued from slope creation, oversteepening and seismicity produced by extensional tectonic activity, as demonstrated by the observation of synsedimentary normal faults and related thickness variations. Other factors, such as pore pressure increases or lithification contrasts probably facilitated sliding. The key role of tectonics is strengthened by the specific relationships between the basin and collapse histories and two major fault systems that cross the study area. The Ayabacas collapse occurred at a turning point in the Central Andean evolution. Before the event, the back‐arc basin had been essentially marine and deepened to the west, with little volcanic activity taking place at the arc. After the event, the back‐arc was occupied by continental to near‐continental environments, and was bounded to the southwest by a massive volcanic arc shedding debris and tuffs into the basin.     

Upscaling of CO2 vertical migration through a periodic layered porous medium: The capillary-free and capillary-dominant cases

We present an upscaled model for the vertical migration of a CO₂ plume through a vertical column filled with a periodic layered porous medium. This model may describe the vertical migration of a CO₂ plume in a perfectly layered horizontal aquifer. Capillarity and buoyancy are taken into account and semi-explicit upscaled flux functions are proposed in the two following cases: (i) capillarity is the main driving force and (ii) buoyancy is the only driving force. In both cases, we show that the upscaled buoyant flux is a bell-shaped function of the saturation, as in the case of a homogeneous porous medium. In the capillary-dominant case, we show that the upscaled buoyant flux is the harmonic mean of the buoyant fluxes in each layer. The upscaled saturation is governed by the continuity of the capillary pressure at the interface between layers. In the capillary-free case, the upscaled buoyant flux and upscaled saturation are determined by the flux continuity condition at the interface. As the flux is not continuous over the entire range of saturation, the upscaled saturation is only defined where continuity is verified, i.e. in two saturation domains. As a consequence, the upscaled buoyant flux is described by a piecewise continuous function. Two analytical approximations of this flux are proposed and this capillary-free upscaled model is validated for two cases of heterogeneity. Upscaled and cell averaged saturations are in good agreement. Furthermore, the proposed analytical upscaled fluxes provide satisfactory approximations as long as the saturation set at the inlet of the column is in a range where analytical and numerical upscaled fluxes are close.     

A Bayesian methodology to update the probabilistic seismic hazard assessment

This paper presents a Bayesian methodology for updating the seismic hazard curves. The methodology is based on the comparison of predictive exceedance rates of a fixed acceleration level (given by the seismic hazard curves) and the observed exceedance rates in some selected sites. The application of the methodology needs, firstly, the definition of a prior probabilistic seismic hazard assessment based in a logic tree. Each main branch corresponds to a probabilistic model of calculus of seismic hazard. The method considers that, initially (or a priori), the weights of all branches of the logic tree are equivalent. Secondly, the method needs to compile the observations in the region. They are introduced in a database containing the recorded acceleration data (during the instrumental period). Nevertheless, the instrumental period in stable zones (as France) shows only very low acceleration levels recorded during a short observation period. Then, a method to enlarge the REX (number of observations) is presented taking into account the historical data and defining “synthetic” accelerations in the sites of observation. The synthetic REX allows to expand the period of observation and to increase the acceleration thresholds used in the Bayesian updating process. The application of the Bayesian approach leads to a new and more objective definition of the weights of each branch of the logic tree and, therefore, to new seismic hazard curves (mean and centiles). The Bayesian approach doesn’t change the probabilistic models (seismic hazard curves). It only modifies the weights of each branch of the logic tree.     

Comparing the performances of four stochastic optimisation methods using analytic objective functions, 1D elastic full‐waveform inversion,and residual static computation

We compare the performances of four stochastic optimisation methods using four analytic objective functions and two highly non‐linear geophysical optimisation problems: one‐dimensional elastic full‐waveform inversion and residual static computation. The four methods we consider, namely, adaptive simulated annealing, genetic algorithm, neighbourhood algorithm, and particle swarm optimisation, are frequently employed for solving geophysical inverse problems. Because geophysical optimisations typically involve many unknown model parameters, we are particularly interested in comparing the performances of these stochastic methods as the number of unknown parameters increases. The four analytic functions we choose simulate common types of objective functions encountered in solving geophysical optimisations: a convex function, two multi‐minima functions that differ in the distribution of minima, and a nearly flat function. Similar to the analytic tests, the two seismic optimisation problems we analyse are characterised by very different objective functions. The first problem is a one‐dimensional elastic full‐waveform inversion, which is strongly ill‐conditioned and exhibits a nearly flat objective function, with a valley of minima extended along the density direction. The second problem is the residual static computation, which is characterised by a multi‐minima objective function produced by the so‐called cycle‐skipping phenomenon. According to the tests on the analytic functions and on the seismic data, genetic algorithm generally displays the best scaling with the number of parameters. It encounters problems only in the case of irregular distribution of minima, that is, when the global minimum is at the border of the search space and a number of important local minima are distant from the global minimum. The adaptive simulated annealing method is often the best‐performing method for low‐dimensional model spaces, but its performance worsens as the number of unknowns increases. The particle swarm optimisation is effective in finding the global minimum in the case of low‐dimensional model spaces with few local minima or in the case of a narrow flat valley. Finally, the neighbourhood algorithm method is competitive with the other methods only for low‐dimensional model spaces; its performance sensibly worsens in the case of multi‐minima objective functions.     

Microbial deposits in the aftermath of the end‐Permian mass extinction: A diverging case from the Mineral Mountains (Utah,USA)

The Lower Triassic Mineral Mountains area (Utah, USA) preserves diversified Smithian and Spathian reefs and bioaccumulations that contain fenestral‐microbialites and various benthic and pelagic organisms. Ecological and environmental changes during the Early Triassic are commonly assumed to be associated with numerous perturbations (productivity changes, acidifica‐tion, redox changes, hypercapnia, eustatism and temperature changes) post‐dating the Permian–Triassic mass extinction. New data acquired in the Mineral Mountains sediments provide evidence to decipher the relationships between depositional environments and the growth and distribution of microbial structures. These data also help to understand better the controlling factors acting upon sedimentation and community turnovers through the Smithian–early Spathian. The studied section records a large‐scale depositional sequence during the Dienerian(?)–Spathian interval. During the transgression, depositional environments evolved from a coastal bay with continental deposits to intertidal fenestral–microbial limestones, shallow subtidal marine sponge–microbial reefs to deep subtidal mud‐dominated limestones. Storm‐induced deposits, microbialite–sponge reefs and shallow subtidal deposits indicate the regression. Three microbialite associations occur in ascending order: (i) a red beds microbialite association deposited in low‐energy hypersaline supratidal conditions where microbialites consist of microbial mats and poorly preserved microbially induced sedimentary structure; (ii) a Smithian microbialite association formed in moderate to high‐energy, tidal conditions where microbialites include stromatolites and associated carbonate grains (oncoids, ooids and peloids); and (iii) a Spathian microbialite association developed in low‐energy offshore conditions that is preserved as multiple decimetre thick isolated domes and coalescent domes. Data indicate that the morphologies of the three microbialite associations are controlled primarily by accommodation, hydrodynamics, bathymetry and grain supply. This study suggests that microbial constructions are controlled by changes between trapping and binding versus precipitation processes in variable hydrodynamic conditions. Due to the presence of numerous metazoans associated with microbialites throughout the Smithian increase in accommodation and Spathian decrease in accommodation, the commonly assumed anachronistic character of the Early Triassic microbialites and the traditional view of prolonged deleterious conditions during the Early Triassic time interval is questioned.     

Analysis of lateral rock spreading process initiation with a numerical modelling approach

The first erosional front of the south‐eastern Alps has been affected by several deep‐seated landslides. The movements affected a Jurassic limestone sequence underlain by an upper Trias unit (Keuper) mainly composed of marl and clays. Horst‐ and graben‐like structures, newly formed gravitational faults, outflow and other field evidence suggest that such movements are consistent with lateral rock spreading processes. An example of such processes is the ‘La Marbrière’ slope, where field investigations have revealed that three distinct zones corresponding to three evolutionary stages of movement exist. The question is: What are the structural and mechanical modifications undergone by the slope between those evolution stages? To tackle this question, numerical models based on a 2D cross section of the ‘La Marbrière’ slope have been performed. Results show that, to fit the structural evidence of the precollapse situation, elasto‐plastic behaviour of the Trias Keuper unit is needed. The transition between the precollapse and fresh post‐collapse situations requires a modification of the mechanical behaviour of the Trias Keuper unit to Burger type, resulting in the formation of antithetic and synthetic gravitational faults (producing a graben‐like structure).     

Seasonal variation of infiltration rates through pond bed in a managed aquifer recharge system in St-André, Belgium

In Belgium, IWVA uses managed aquifer recharge (MAR) to recharge the aquifer with treated wastewater generated from the communities to sustain the potable water supply on the Belgian coast. This MAR facility is faced with a challenge of reduced infiltration rates during the winter season when pond water temperatures near 4°C. This study involves the identification of the predominant factor influencing the rate of infiltration through the pond bed. Several factors, including pumping rates, natural recharge, tidal influences of the North Sea and pond-water temperature, were identified as potential causes for variation of the recharge rate. Correlation statistics and linear regression analysis were used to determine the sensitivity of the infiltration rate to the aforementioned factors. Two groundwater flow models were developed in visual MODFLOW to simulate the water movement under the pond bed and to obtain the differences in flux to track the effects of variation of hydraulic conductivity during the two seasons. A 32% reduction in vertical hydraulic gradient in the top portion of the aquifer was observed in winter, causing the recharge rates to fluctuate. Results showed that water temperature caused a 30% increase in hydraulic conductivity in summer as compared with winter and has the maximum impact on infiltration rate. Cyclic variations in water viscosity, occurring because of seasonal temperature changes, influence the saturated hydraulic conductivity of the pond bed. Results from the models confirm the impact on infiltration rate by temperature-influenced hydraulic conductivity.