Assessing cave internal aerology in understanding carbon dioxide (CO<Subscript>2</Subscript>) dynamics: implications on calcite mass variation on the wall of Lascaux Cave (France)

Carbon dioxide gas is a key component in dissolution and precipitation of carbonates in karst and cave systems. Therefore, characterizing the internal aerology of a cave is essential to obtain the spatiotemporal distribution of temperature and CO₂ level. In this research, Lascaux Cave (France), an important adorned cavity, was studied. First, the spatiotemporal distribution of CO₂ and temperatures were examined using continuous monitoring at a per minute basis. High-resolution spatial measurements (14 PCO₂ locations and 27 locations for temperature) were carried out for a year in the epikarst and the cave (February 2015 to February 2016). The spatiotemporal analysis presents that air and rock temperatures vary for less than a degree Celsius (12.4–12.9 °C). These are controlled by the conduction of the external thermal waves through the overlying calcarenite massif. As a consequence, two seasonal internal aerologic regimes were identified: stratification and convection. These regimes govern the spatiotemporal distribution of the CO₂ levels (1.1–3.7%), showing that this parameter is a good natural marker of the internal air movements. Second, a method was proposed to estimate the calcite mass potentially affected by condensation water (dissolution process) and exfiltration water (precipitation process). This method, based on numerical simulations, relies on CO₂ and air and rock temperature spatiotemporal distributions in the cave. Third, the method was applied using the case of the left wall of the Hall of the Bulls (one of the most adorned part of the cave). Results showed that the calcite mass, possibly dissolved, varies from 0.0002 to 0.006 g when the mass potentially precipitated is higher (from 0.013 to 0.067 g) depending on the aerologic conditions. This method allows determining which alteration process (e.g., precipitation or dissolution) could eventually lead to the largest variation of calcite on the wall. The results can serve as useful data to the cave experts of the French Ministry of Culture and Communication in terms of Lascaux Cave management policies.     

An Experimental High-Resolution Forecast System During the Vancouver 2010 Winter Olympic and Paralympic Games

Environment Canada ran an experimental numerical weather prediction (NWP) system during the Vancouver 2010 Winter Olympic and Paralympic Games, consisting of nested high-resolution (down to 1-km horizontal grid-spacing) configurations of the GEM–LAM model, with improved geophysical fields, cloud microphysics and radiative transfer schemes, and several new diagnostic products such as density of falling snow, visibility, and peak wind gust strength. The performance of this experimental NWP system has been evaluated in these winter conditions over complex terrain using the enhanced mesoscale observing network in place during the Olympics. As compared to the forecasts from the operational regional 15-km GEM model, objective verification generally indicated significant added value of the higher-resolution models for near-surface meteorological variables (wind speed, air temperature, and dewpoint temperature) with the 1-km model providing the best forecast accuracy. Appreciable errors were noted in all models for the forecasts of wind direction and humidity near the surface. Subjective assessment of several cases also indicated that the experimental Olympic system was skillful at forecasting meteorological phenomena at high-resolution, both spatially and temporally, and provided enhanced guidance to the Olympic forecasters in terms of better timing of precipitation phase change, squall line passage, wind flow channeling, and visibility reduction due to fog and snow.     

Corrective pattern-matching simulation with controlled local-mean histogram

This paper presents a new stochastic simulation method that improves on our unilateral patchwork simulation method (Faucher et al., Stoch Environ Res Risk Assess, 27:253–273, 2012) by eliminating anisotropy biases. As in our unilateral method, images are built by assembling patterns together while controlling the local-mean histogram. The patterns, which are square image-pieces, are picked in a reference image. The reference image is used as a data bank holding the statistical informations about the random field to simulate. In contrast with the unilateral method, the path followed by our new simulation is random and guided by local errors. The new method, called corrective pattern-matching simulation, proceeds iteratively by making local corrections to the simulated image. For several types of images, it is shown that our simulations respect conditioning data and reproduce faithfully the reference image visual appearance. It is shown that the control of the local-mean histogram allows to control one-point statistics and multi-point statistics.     

FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

An investigation of the SO2 content of the venusian mesosphere using SPICAV-UV in nadir mode

Using the SPICAV-UV spectrometer aboard Venus Express in nadir mode, we were able to derive spectral radiance factors in the middle atmosphere of Venus in the 170-320 nm range at a spectral resolution of R ? 200 during 2006 and 2007 in the northern hemisphere. By comparison with a radiative transfer model of the upper atmosphere of Venus, we could derive column abundance above the visible cloud top for SO₂ using its spectral absorption bands near 280 and 220 nm. SO₂ column densities show large temporal and spatial variations on a horizontal scale of a few hundred kilometers. Typical SO₂ column densities at low latitudes (up to 50°N) were found between 5 and 50 μm-atm, whereas in the northern polar region SO₂ content was usually below 5 μm-atm. The observed latitudinal variations follow closely the cloud top altitude derived by SPICAV-IR and are thought to be of dynamical origin. Also, a sudden increase of SO₂ column density in the whole northern hemisphere has been observed in early 2007, possibly related to a convective episode advecting some deep SO₂ into the upper atmosphere.     

Contamination of port zone sediments by metals from Large Marine Ecosystems of Brazil

Sediment contamination by metals poses risks to coastal ecosystems and is considered to be problematic to dredging operations. In Brazil, there are differences in sedimentology along the Large Marine Ecosystems in relation to the metal distributions. We aimed to assess the extent of Al, Fe, Hg, Cd, Cr, Cu, Ni, Pb and Zn contamination in sediments from port zones in northeast (Mucuripe and Pecém) and southeast (Santos) Brazil through geochemical analyses and sediment quality ratings. The metal concentrations found in these port zones were higher than those observed in the continental shelf or the background values in both regions. In the northeast, metals were associated with carbonate, while in Santos, they were associated with mud. Geochemical analyses showed enrichments in Hg, Cd, Cu, Ni and Zn, and a simple application of international sediment quality guidelines failed to predict their impacts, whereas the use of site-specific values that were derived by geochemical and ecotoxicological approaches seemed to be more appropriate in the management of the dredged sediments.     

The Influence of Topology on Hydraulic Conductivity in a Sand-and-Gravel Aquifer

A field experiment consisting of geophysical logging and tracer testing was conducted in a single well that penetrated a sand-and-gravel aquifer at the U.S. Geological Survey Toxic Substances Hydrology research site on Cape Cod, Massachusetts. Geophysical logs and flowmeter/pumping measurements were obtained to estimate vertical profiles of porosity ϕ, hydraulic conductivity K, temperature, and bulk electrical conductivity under background, freshwater conditions. Saline-tracer fluid was then injected into the well for 2 h and its radial migration into the surrounding deposits was monitored by recording an electromagnetic-induction log every 10 min. The field data are analyzed and interpreted primarily through the use of Archie's (1942) law to investigate the role of topological factors such as pore geometry and connectivity, and grain size and packing configuration in regulating fluid flow through these coarse-grained materials. The logs reveal no significant correlation between K and ϕ, and imply that groundwater models that link these two properties may not be useful at this site. Rather, it is the distribution and connectivity of the fluid phase as defined by formation factor F, cementation index m, and tortuosity α that primarily control the hydraulic conductivity. Results show that F correlates well with K, thereby indicating that induction logs provide qualitative information on the distribution of hydraulic conductivity. A comparison of α, which incorporates porosity data, with K produces only a slightly better correlation and further emphasizes the weak influence of the bulk value of ϕ on K.     

Simulation of the subsurface mobility of carbon nanoparticles at the field scale

The production and use of nanomaterials will inevitably lead to their disposal in the natural environment. To assess the risk that these materials pose to human and ecosystem health an understanding of their mobility and ultimate fate is essential. To date, however, relatively little research has been conducted on the fate of nanoparticles in subsurface systems. In this study the subsurface mobility of two carbon nanoparticles: nano-fullerenes (nC₆₀) and multi-walled carbon nanotubes (MWCNTs) is assessed. A two-dimensional finite element model was used to simulate the movement of these nanoparticles under a range of hydrologic and geological conditions, including a heterogeneous permeability field. The numerical model is based on colloid filtration theory (CFT) with a maximum retention capacity term. For the conditions evaluated the carbon nanotubes are much more mobile than nC₆₀ due to the smaller collector efficiency associated with carbon nanotubes. However, the mobility of nC₆₀ increased significantly when a maximum retention capacity term was included in the model. Model results also demonstrate that, for the systems examined, nanoparticles were predicted to be less mobile in heterogeneous systems compared to the homogeneous systems with the same average hydraulic properties.     

Morphology,structure and kinematics of a rainfall controlled slow‐moving Andean landslide,Peru

The large slow‐moving landslide of Maca is located in the upper Colca valley (southern Peru), a region characterized by a well pronounced rainy period, and intense and recurrent sustained seismicity. The landslide, developed in deep lacustrine deposits, has recently accelerated, threatening the Maca village. This work aims at understanding the rupture mechanism and the causes of the recent landslide reactivation/acceleration. We present a multidisciplinary characterization of the Maca landslide that includes: (i) geological and morphological mapping in the field; (ii) remote sensing analysis using an historical aerial photograph of 1955 and the Pléiades satellite images (2013); (iii) global positioning system (GPS) including time‐series of surveys over 13 years, and continuous measurements over 14 months; (iv) a geophysical campaign with deep electrical resistivity tomography profiles acquired across the landslide mass. Our study shows that this 60 Mm³ landslide, which can be classified as a clay/silt compound landslide, moved by 15 m between 2001 and 2014 with a large inter‐annual velocity variation (up to a factor of 500) depending on the rainfall intensity. We suggest that these dramatic changes in velocity are the result of the combination of a threshold mechanism and the short intense rainy season in Peru. This study reveals three main driving factors acting at different timescales: (i) over several decades, the river course has significantly changed, causing the Maca landslide reactivation in the 1980s due to the erosion of its toe; (ii) at the year scale, a minimum amount of rainfall is required to trigger the motion and this amount controls the landslide velocity; (iii) transient changes in slide velocity may occur anytime due to earthquakes. This study particularly highlights the non‐linear behaviour of the motion with rainfall. Copyright © 2016 John Wiley & Sons, Ltd.     

The Ereendavaa Range (north-eastern Mongolia): an additional argument for Mesozoic extension throughout eastern Asia

Mesozoic rift basins locally bounding metamorphic core complexes have been recognized in Transbaikalia and northern China. Numerous basement outcrops located between these two regions, in eastern Mongolia, are considered as pre-Palaeozoic in age. One of these, the Ereendavaa Range, appears as a gneissic core marked by amphibolite-facies metamorphic conditions. The range is overlain to the NW by the unmetamorphosed Mesozoic Onon Basin. Below the basin, the upper part of the range consists of a gently NW-dipping shear zone associated with top-to-the-NW motion. The structural pattern is consistent with syn-extensional exhumation of the range. Preliminary geochronological data indicate that the shear zone is late Jurassic to early Cretaceous in age, coeval with the Onon Basin. These new data from eastern Mongolia constitute a link between Transbaikalia and northern China, indicating that NW–SE extensional Mesozoic tectonics occurred throughout the entire region.