Discrete modelling of rock-ageing in rockfill dams

The aim of this paper is to obtain, from a numerical study using the discrete element method, quantitative information on the behaviour of a rockfill dam, particularly the effect of rock ageing on the movements of the dam. A numerical model is defined at three different scales: the block scale, the elementary representative volume scale and the dam scale. The model uses the discrete element code PCF^2D, considering breakable clusters of 2D balls. Three rock-ageing laws, taking into account the time-evolving resistance for blocks, are proposed. The different parameters are determined from experimental data of laboratory tests performed on rock blocks. The relevance of the model seems acceptable, considering that the displacements obtained in the simulations are in the same order of magnitude as the displacements measured on the actual dam.     

Accurate and High-Precision Measurement of Lithium Isotopes in Two Reference Materials by MC-ICP-MS

We report here a newly developed method for measurement of Li isotopes by use of multi-collector ICP-MS (Neptune) allowing rapid and high precision determination of Li isotope ratios at low levels of lithium (15–20 ng). The lithium reference sample solution IRMM-016 was analysed over a period of ten months with an external reproducibility of 0.24% (2s, n = 52). Chemical separation of Li from matrix was performed on the seawater sample IRMM BCR-403, for which a mean δ⁷Li value of + 31.0 ± 0.1 % (2s/√n, n = 31) was obtained. This mean value is in good agreement with those previously published for other seawater samples. BCR-403 seawater being readily available, we propose that this seawater sample be used as a reference sample for Li isotope measurements.     

A statistical�Cdynamical scheme for reconstructing ocean forcing in the Atlantic. Part II: methodology, validation and application to high-resolution ocean models

A novel statistical?Cdynamical scheme has been developed to reconstruct the sea surface atmospheric variables necessary to force an ocean model. Multiple linear regressions are first built over a so-called learning period and over the entire Atlantic basin from the observed relationship between the surface wind conditions, or predictands, and the anomalous large scale atmospheric circulations, or predictors. The latter are estimated in the extratropics by 500?hPa geopotential height weather regimes and in the tropics by low-level wind classes. The transfer function further combined to an analog step is then used to reconstruct all the surface variables fields over 1958?C2002. We show that the proposed hybrid scheme is very skillful in reproducing the mean state, the seasonal cycle and the temporal evolution of all the surface ocean variables at interannual timescale. Deficiencies are found in the level of variance especially in the tropics. It is underestimated for 2-m temperature and humidity as well as for surface radiative fluxes in the interannual frequency band while it is slightly overestimated at higher frequency. Decomposition in empirical orthogonal function (EOF) shows that the spatial and temporal coherence of the forcing fields is however very well captured by the reconstruction method. For dynamical downscaling purposes, reconstructed fields are then interpolated and used to carry out a high-resolution oceanic simulation using the NATL4 (1/4°) model integrated over 1979?C2001. This simulation is compared to a reference experiment where the original observed forcing fields are prescribed instead. Mean states between the two experiments are virtually undistinguishable both in terms of surface fluxes and ocean dynamics estimated by the barotropic and the meridional overturning streamfunctions. The 3-dimensional variance of the simulated ocean is well preserved at interannual timescale both for temperature and salinity except in the tropics where it is underestimated. The main modes of interannual variability assessed through EOF are correctly reproduced for sea surface temperature, barotropic streamfunction and mixed layer depth both in terms of spatial structure and temporal evolution. Collectively, our results provide evidence that the statistical?Cdynamical scheme presented in this two-part study is an efficient and promising tool to infer oceanic changes (in particular those related to the wind-driven circulation) due to modifications in the large-scale atmospheric circulation. As a prerequisite, we have here validated the method for present-day climate; we encourage its use for climate change studies with some adaptations though.     

The Minorca Basin: a buffer zone between the Valencia and Liguro‐Provençal Basins (NW Mediterranean Sea)

Detailed analyses of seismic profiles and boreholes in the Valencia Basin (VB) reveal a differentiated basin, the Minorca Basin (MB), lying between the old Mesozoic Valencia Basin sensu strico (VBss) and the young Oligocene Liguro‐Provençal Basin (LPB). The relationship between these basins is shown through the correlation of four Miocene‐to‐present‐day megasequences. The Central and North Balearic Fracture Zones (CFZ and NBFZ) that border the MB represent two morphological and geodynamical thresholds that created an accommodation in steps between the three domains. Little to no horizontal Neogene movements are found for the Ibiza and Majorca Islands. In contrast, the counterclockwise movement of the Corso‐Sardinian blocks induced a counterclockwise movement of the Minorca block towards the SE along the CFZ and NBFZ, during the exhumation of lower continental crust in the LPB. This new understanding implies pure Neogene vertical subsidence in the VBss and places the AlKaPeCa northeastward of the present‐day Alboran Area.     

Exploration of tektite formation processes through water and metal content measurements

Abstract– Impact events are a significant surface‐modifying process on solar system objects lacking frequent resurfacing by atmospheric or igneous processes. To better explore the effects of impacts on surface materials, we measured the water and trace element compositions of tektites from Da Lat, Vietnam. Fourier transform infrared spectroscopy was used for water measurements and laser ablation inductively coupled plasma mass spectrometry was used for trace element measurements. Consistent with previous investigations of tektites, we found that the samples are depleted in volatile metals (e.g., Zn, Pb) as well as in water compared with the average continental crust, although water contents are still extremely high for rocks melted at atmospheric pressure. While Zn and Pb concentrations are linearly correlated with each other, there is no correlation between H₂O and Zn or Pb contents. Combined with water contents of other tektites in the Australasian strewn field, our results demonstrate that the source impact occurred at a wet site in or off the coast of the Indochina peninsula.     

Evaluation of aDcp processing options for secondary flow identification at river junctions

Secondary circulation in river confluences results in a spatial and temporal variation of fluid motion and a relatively high level of morphodynamic change. Acoustic Doppler current profiler (aDcp) vessel-mounted flow measurements are now commonly used to quantify such circulation in shallow water fluvial environments. It is well established that such quantification using vessel-mounted aDcps requires repeated survey of the same cross-section. However, less attention has been given to how to process these data. Most aDcp data processing techniques make the assumption of homogeneity between the measured radial components of velocity. As acoustic beams diverge with distance from the aDcp probe, the volume of the flow that must be assumed to be homogeneous between the beams increases. In the presence of secondary circulation cells, and where there are strong rates of shear in the flow, the homogeneity assumption may not apply, especially deeper in the water column and close to the bed. To reduce dependence on this assumption, we apply a newly-established method to aDcp data obtained for two medium-sized (~60–80 m wide) gravel-bed river confluences and compare the results with those from more conventional data processing approaches. The comparison confirms that in the presence of strong shear our method produces different results to more conventional approaches. In the absence of a third set of fully independent data, we cannot demonstrate conclusively which method is best, but our method involves less averaging and so in the presence of strong shear is likely to be more reliable. We conclude that it is wise to apply both our method and more conventional methods to identify where data analysis might be impacted upon by strong shear and where inferences of secondary circulation may need to be made more cautiously. © 2019 John Wiley & Sons, Ltd.     

The application of an innovative inverse model for understanding and predicting landslide movements (Salazie cirque landslides,Reunion Island)

The prediction of landslide movement acceleration is a complex problem, among others identified for deep-seated landslides, and represents a crucial step for risk assessment. Within the scope of this problem, the objective of this paper is to explore a modelling method that enables the study of landslide function and facilitates displacement predictions based on a limited data set. An inverse modelling approach is proposed for predicting the temporal evolution of landslide movement based on rainfall and displacement velocities. Initially, the hydrogeology of the studied landslides was conceptualised based on correlative analyses. Subsequently, we applied an inverse model with a Gaussian-exponential transfer function to reproduce the displacements. This method was tested on the Grand Ilet (GI) and Mare-à-Poule-d’Eau (HB) landslides on Reunion Island in the Indian Ocean. We show that the behaviour of landslides can be modelled by inverse models with a bimodal transfer function using a Gaussian-exponential impulse response. The cumulative displacements over 7 years of modelling (2 years of calibration period for GI, and 4 years for HB) were reproduced with an RMSE above 0.9. The characteristics of the bimodal transfer function are directly related to the hydrogeological functioning demonstrated by the correlative analyses: the rapid reaction of a landslide can be associated with the effect of a preferential flow path on groundwater level variations. Thus, this study shows that the inverse model using a Gaussian-exponential transfer function is a powerful tool for predicting deep-seated landslide movements and for studying how they function. Beyond modelling displacements, our approach effectively demonstrates its ability to contribute relevant data for conceptualising the sliding mechanisms and hydrogeology of landslides.     

Voronoi grids conforming to 3D structural features

Flow simulation in a reservoir can be highly impacted by upscaling errors. These errors can be reduced by using simulation grids with cells as homogeneous as possible, hence conformable to horizons and faults. In this paper, the coordinates of 3D Voronoi seeds are optimized so that Voronoi cell facets honor the structural features. These features are modeled by piecewise linear complex (PLC). The optimization minimizes a function made of two parts: (1) a barycentric function, which ensures that the cells will be of good quality by maximizing their compactness; and (2) a conformity function, which allows to minimize the volume of cells that is isolated from the Voronoi seed w.r.t., a structural feature. To determine the isolated volume, a local approximation of the structural feature inside the Voronoi cells is used to cut the cells. It improves the algorithm efficiency and robustness compared to an exact cutting procedure. This method, used jointly with an adaptive gradient solver to minimize the function, allows dealing with complex 3D geological cases. It always produces a Voronoi simulation grid with the desired number of cells.     

A combined modelling and geochemical study of the fate of terrigenous inputs from mixed natural and mining sources in a coral reef lagoon (New Caledonia)

Open-cast mining for Ni, Cr and Co was conducted in the south-west part of New Caledonia during the 20th century. Abandoned mining and prospecting sites were severely affected by erosion, resulting in an increase in the load of terrigenous particles transported to the coral reef lagoon. This article assesses the impact of a typical small catchment area (La Coulée River, 85 km2 watershed) on two bays (Boulari and Sainte Marie) located near Noumea, New Caledonia's main city. This multi-disciplinary study combines geochemical, sedimentological, and hydrodynamic approaches. Ni and Cr concentrations contained in the geochemical matrix phase of the pelitic fraction were determined. The study of the geochemical signatures together with sedimentological data and 3D numerical simulations of the transport of non-settling particles throughout the lagoon demonstrated that terrigenous inputs from the Coulée River were mainly transported and deposited along the shoreline, reaching areas as distant as Sainte Marie Bay. Although quantitatively low (about 3% of the pelite mass of the bay sediments), the terrigenous inputs in Sainte Marie Bay originating from the Coulée River could be traced. The metal content in suspended matter was over 7000 mg kg(-1) for Ni and 4200 mg kg(-1) for Cr.