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.     

Physical model experiments for sediment supply to the old Rhine through induced bank erosion

A mobile-bed, undistorted physical model （1：40） has been used to investigate different sediment supply strategies to the Old Rhine through bank protection removal and modifications of groyne dimensions and configuration, which cause bank erosion. This trained channel was previously the main bed of the upper Rhine downstream of Basel （Switzerland）, but it has an artificially low flow regime since the construction of the Grand Canal d＇Alsace, a navigation canal, and a flow control dam at Kembs （France）. Training works and subsequent channel incision have also greatly reduced sediment transport rates and created a heavily armoured bed. The modelled pilot site has a groyne field on the left bank. Results show that the currently existing groynes at the site are not effective in creating high bank-side velocities conducive to bank erosion, even for flow rates significantly higher than the mean annual flow rate. The river bank has also proved to be more resistant than previously thought, allowing long stretches of bank protection to be safely removed. The physical model testing process has produced a new configuration for the groyne field, where two higher, larger island groynes are placed further apart than the three existing attached groynes. This innovative approach has proved effective, causing bank erosion for flow rates below the mean annual flow rate, with consistent erosion being observed. It has also been found that such a configuration does not pose a hazard for the Grand Canal d＇Alsace, which is situated next to the Old Rhine, through excessive bank erosion during high flow events.     

Evaluation of equivalence between different methods for enumeration of fecal indicator bacteria before and after adoption of the new Bathing Water Directive and risk assessment of pollution

The quality of bathing water is of considerable public importance due to the possibility of fecal contamination. In 2009, Croatia implemented the new European Bathing Water Directive (BWD, 2006/7/EC) establishing stricter microbiological standards for new parameters with new reference methods. This study aims to evaluate the equivalence of different methods according to the old and revised BWD and to provide the possibility of data comparison. Furthermore, the directive requires the establishment of the bathing water profile (BWP) for pollution risk assessment. The estimation of consistency of pollution risk assessment with obtained microbiological results was also performed.     

Finding Optimal Travel Routes with Uncertain Cost Data

Geospatial data analysis techniques are widely used to find optimal routes from specified starting points to specified destinations. Optimality is defined in terms of minimizing some impedance value over the length of the route – the value to be minimized might be distance, travel time, financial cost, or any other metric. Conventional analysis procedures assume that impedance values of all possible travel routes are known a priori, and when this assumption holds, efficient solution strategies exist that allow truly optimal solutions to be found for even very large problems. When impedance values are not known with certainty a priori, exact solution strategies do not exist and heuristics must be employed. This study evaluated how the quality of the solutions generated by one such heuristic were impacted by the nature of the uncertainty in the cost database, the nature of the costs themselves, and the parameters used in the heuristic algorithm. It was found that all of these factors influenced the qualities of the solutions produced by the heuristic, but encouragingly, an easily controlled parameter of the heuristic algorithm itself played the most important role in controlling solution quality.     

Structure and nuclear density distribution in the cheralite—CaTh(PO4)2: studies of its behaviour under high pressure (36?GPa)

The crystal structure of the cheralite—CaTh(PO₄)₂—has been revisited by neutron diffraction and its behaviour under high pressure investigated by X-ray diffraction up to 36?GPa. The neutron diffraction data at ambient pressure gave a more accurate determination of the Ca/Th cation position than previous XRD data, taking advantage that the neutron scattering lengths of calcium and thorium are of same order of magnitude. The nuclear density distribution was also determined using the maximum entropy method (MEM) confirming that the two cations are not located at the same position in the unit cell but are slightly displaced from one another along a specific direction in order to minimize the electrostatic repulsion with the surrounding phosphorus atoms. At high pressure, the compound did not show any phase transition or amorphization. From the evolution of the unit-cell volume as a function of the pressure, the zero-pressure bulk modulus B₀ and its pressure derivative B₀ ^′ have been determined by fitting the experimental compressibility curve to the Birch–Murnaghan equation of state. The results are B₀?=?140(2) GPa and B ₀ ^′ ?=?4.4(4) GPa.     

Boundary integral formulation and two-dimensional fundamental solutions for dynamic behavior analysis of unsaturated soils

In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–p_w–p_a theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case.     

StalAge – An algorithm designed for construction of speleothem age models

Here we present a new algorithm (StalAge), which is designed to construct speleothem age models. The algorithm uses U-series ages and their corresponding age uncertainty for modelling and also includes stratigraphic information in order to further constrain and improve the age model. StalAge is applicable to problematic datasets that include outliers, age inversions, hiatuses and large changes in growth rate. Manual selection of potentially inaccurate ages prior to application is not required. StalAge can be applied by the general, non-expert user and has no adjustable free parameters. This offers the highest degree of reproducibility and comparability of speleothem records from different studies. StalAge consists of three major steps. Firstly, major outliers are identified. Secondly, age data are screened for minor outliers and age inversions, and the uncertainty of potential outliers is increased using an iterative procedure. Finally, the age model and corresponding 95%-confidence limits are calculated by a Monte-Carlo simulation fitting ensembles of straight lines to sub-sets of the age data.We apply StalAge to a synthetic stalagmite ’sample’ including several problematic features in order to test its performance and robustness. The true age is mostly within the 95%-confidence age limits of StalAge showing that the calculated age models are accurate even for very difficult samples. We also apply StalAge to three published speleothem datasets. One of those is annually laminated, and the lamina counting chronology agrees with the age model calculated by StalAge. For the other two speleothems the resulting age models are similar to the published age models, which are both based on smoothing splines. Calculated uncertainties are in the range of those calculated by combined application of Bayesian chronological ordering and a spline, showing that StalAge is efficient in using stratigraphic information in order to reduce age model uncertainty.The algorithm is written in the open source statistical software R and available from the authors or as an electronic supplement of this paper.     

Timing of advance and basal condition of the Laurentide Ice Sheet during the last glacial maximum in the Richardson Mountains,NWT

This study presents new ages for the northwest section of the Laurentide Ice Sheet (LIS) glacial chronology from material recovered from two retrogressive thaw slumps exposed in the Richardson Mountains, Northwest Territories, Canada. One study site, located at the maximum glacial limit of the LIS in the Richardson Mountains, had calcite concretions recovered from aufeis buried by glacial till that were dated by U/Th disequilibrium to 18,500 cal yr BP. The second site, located on the Peel Plateau to the east yielded a fossil horse (Equus) mandible that was radiocarbon dated to ca. 19,700 cal yr BP. These ages indicate that the Peel Plateau on the eastern flanks of the Richardson Mountains was glaciated only after 18,500 cal yr BP, which is later than previous models for the global last glacial maximum (LGM). As the LIS retreated the Peel Plateau around 15,000 cal yr BP, following the age of the Tutsieta phase, we conclude that the presence of the northwestern margin of the LIS at its maximum limit was a very short event in the western Canadian Arctic.