Selection of sand models and identification of parameters using an enhanced genetic algorithm

Numerous constitutive models of granular soils have been developed during the last few decades. As a consequence, how to select an appropriate model with the necessary features based on conventional tests and with an easy way of identifying parameters for geotechnical applications has become a major issue. This paper aims to discuss the selection of sand models and parameters identification by using genetic algorithm. A real‐coded genetic algorithm is enhanced for the optimization with high efficiency. Models with gradually varying features (elastic‐perfectly plastic modelling, nonlinear stress–strain hardening, critical state concept and two‐surface concept) are selected from numerous sand models as examples for optimization. Conventional triaxial tests on Hostun sand are selected as the objectives in the optimization. Four key points are then discussed in turn: (i) which features are necessary to be accounted for in constitutive modelling of sand; (ii) which type of tests (drained and/or undrained) should be selected for an optimal identification of parameters; (iii) what is the minimum number of tests that should be selected for parameter identification; and (iv) what is the suitable and least strain level of objective tests to obtain reliable and reasonable parameters. Finally, a useful guide, based on all comparisons, is provided at the end of the discussion. Copyright © 2015 John Wiley & Sons, Ltd.     

TWO REFERENCE MATERIALS, TRACHYTES MDO-G AND ISH-G, FOR ARGON DATING (K-AR AND 40AR/39AR) OF PLEISTOCENE AND HOLOCENE ROCKS

With the recent extension of K-Ar dating methods within the Upper Pleistocene time, the use of new standards is recommended to allow calibration and inter-laboratory comparisons. For this purpose, two pure microlitic groundmass reference materials were prepared from trachy-basaltic lavas: Trachyte (MDO-G) from Mont Dore, Puy de Sancy, Massif Central is around 250 000 years old; the second Trachyte (ISH-G) from the Isle of Ischia, Gulf of Naples (Italy), is historical: Arso flow erupted on Christmas 1301 - January 1302 AD.     

On the age calibration of the geomagnetic polarity timescale

Knowledge of the age of undated events is not null if a time-order relationship can be found among these events. The knowledge of such a time-ordered sequence can be formalized by using non-informative (uniform) prior probability densities for the ages of undated events and Bayes' theorem to introduce the time-order relationship condition. We show that the conditional probability densities of the ages of events of unknown age are given by various forms of Euler's beta distribution. These distributions yield an estimate of the probability for an undated event to occur in a given age interval.
  We use this method to propose appropriate probabilistic representations of our actual knowledge of the dating of the magnetic polarity reversals during the Cenozoic. These representations take into account the uncertainties arising from irregularities in accretion process and from the quality of a few calibration points. Both types of uncertainties generate large ambiguities in the age of magnetic reversals, which should be taken into consideration when the geomagnetic polarity timescale is used for dating purposes. We propose to use the entropy function to quantify these ambiguities.     

Detection of turbulent coherent motions in a forest canopy part I: Wavelet analysis

Turbulence measurements performed at high frequencies yield data revealing intermittent and multi-scale processes. Analysing time series of turbulent variables thus requires extensive numerical treatment capable, for instance, of performing pattern recognition. This is particularly important in the case of the atmospheric surface layer and specifically in the vicinity of plant canopies, where largescale coherent motions play a major role in the dynamics of turbulent transport processes. In this paper, we examine the ability of the recently developedwavelet transform to extract information on turbulence structure from time series of wind velocities and scalars. It is introduced as a local transform performing a time-frequency representation of a given signal by a specific wavelet function; unlike the Fourier transform, it is well adapted to studying non-stationary signals. After the principles and the most relevant mathematical properties of wavelet functions and transform are given, we present various applications of relevance for our purpose: determination of time-scales, data reconstruction and filtering, and jump detection. Several wavelet functions are inter-compared, using simple artificially generated data presenting large-scale features similar to those observed over plant canopies. Their respective behaviour in the time-frequency domain leads us to assign a specific range of applications for each.     

Perspectives on present-day sea level change: a tribute to Christian le Provost

In this paper, we first discuss the controversial result of the work by Cabanes et al. (Science 294:840–842, 2001), who suggested that the rate of past century sea level rise may have been overestimated, considering the limited and heterogeneous location of historical tide gauges and the high regional variability of thermal expansion which was supposed to dominate the observed sea level. If correct, this conclusion would have solved the problem raised by the IPCC third assessment report [Church et al, Cambridge University Press, Cambridge, pp 881, 2001], namely, the factor two difference between the 20th century observed sea level rise and the computed climatic contributions. However, recent investigations based on new ocean temperature data sets indicate that thermal expansion only explains part (about 0.4 mm/year) of the 1.8 mm/year observed sea level rise of the past few decades. In fact, the Cabanes et al.’s conclusion was incorrect due to a contamination of abnormally high ocean temperature data in the Gulf Stream area that led to an overestimate of thermal expansion in this region. In this paper, we also estimate thermal expansion over the last decade (1993–2003), using a new ocean temperature and salinity database. We compare our result with three other estimates, two being based on global gridded data sets, and one based on an approach similar to that developed here. It is found that the mean rate of thermosteric sea level rise over the past decade is 1.5±0.3 mm/year, i.e. 50% of the observed 3 mm/year by satellite altimetry. For both time spans, past few decades and last decade, a contribution of 1.4 mm/year is not explained by thermal expansion, thus needs to be of water mass origin. Direct estimates of land ice melt for the recent years account for about 1 mm/year sea level rise. Thus, at least for the last decade, we have moved closer to explaining the observed rate of sea level rise than the IPCC third assessment report.     

Spectral and geological study of the sulfate-rich region of West Candor Chasma, Mars

Sulfates have been discovered by the OMEGA spectrometer in different locations of the planet Mars. They are strongly correlated to light toned layered deposits in the equatorial regions. West Candor Chasma is the canyon with the thickest stack of layers and one with the largest area covered by sulfates. A detailed study coupling mineralogy derived from OMEGA spectral data and geology derived from HRSC imager and other datasets leads to some straightforward issues. The monohydrated sulfate kieserite is found mainly over heavily eroded scarps of light toned material. It likely corresponds to a mineral present in the initial rock formed either during formation and diagenesis of sediments, or during hydrothermal alteration at depth, because it is typically found on outcrops that are eroded and steep. Polyhydrated sulfates, that match any Ca-, Na-, Fe-, or Mg-sulfates with more than one water molecule, are preferentially present on less eroded and darker outcrops than outcrops of kieserite. These variations can be the result of a diversity in the composition and/or of the rehydration of kieserite on surfaces with longer exposure. The latter possibility of rehydration in the current, or recent, atmosphere suggests the low surface temperatures preserve sulfates from desiccation, and, also can rehydrate part of them. Strong signatures of iron oxides are present on sulfate-rich scarps and at the base of layered deposits scarps. They are correlated with TES gray hematite signature and might correspond to iron oxides present in the rock as sand-size grains, or possibly larger concretions, that are eroded and transported down by gravity at the base of the scarp. Pyroxenes are present mainly on sand dunes in the low lying terrains. Pyroxene is strongly depleted or absent in the layered deposits. When mixed with kieserite, local observations favor a spatial mixing with dunes over layered deposits. Sulfates such as those detected in the studied area require the presence of liquid water to form by precipitation, either in an intermittent lacustrine environment or by hydrothermal fluid circulation. Both possibilities require the presence of sulfur-rich groundwater to explain fluid circulation. The elevation of the uppermost sulfate signatures suggests the presence of aquifers up to 2.5 km above datum, only 1 km below the plateau surface.     

Structure cristalline d'une ménéghinite naturelle pauvre en cuivre,Cu0,58Pb12,72(Sb7,04Bi0,24)S24Crystal structure of a Cu-poor natural meneghinite,Cu0.58Pb12.72(Sb7.04Bi0.24)S24

Cu-poor meneghinite from La Lauzière Massif (Savoy, France) has the composition (electron microprobe) (in wt%): Pb 59.50, Sb 20.33, Bi 1.19, Cu 0.87, Ag 0.05, Fe 0.03, S 17.62, Se 0.05, Total 99.64. Its crystal structure (X-ray on a single crystal) was solved with R1=0.0506, wR2=0.1026, with an orthorhombic symmetry, space group Pnma, and a=24.080(5) Å, b=4.1276(8) Å, c=11.369(2) Å, V=1130.0(4) Å³, Z=4. Relatively to the model of Euler and Hellner (1960), this structure shows a significantly lower site occupancy factor for the tetrahedral Cu site (0.146 against 0.25). Among the five other metallic sites, Bi appears in the one with predominant Sb. Developed structural formula: Cu_0.15Pb₂(Pb_0.53Sb_0.47)(Pb_0.46Sb_0.54)(Sb_0.75Pb_0.19Bi_0.06)S₆; the reduced one: Cu_0.58Pb_12.72(Sb_7.04Bi_0.24)S₂₄. The formation of such a Cu-poor variety seems to be related to specific paragenetic conditions (absence of coexisting galena), or to crystallochemical constraints (minor Bi). To cite this article: Y. Moëlo et al., C. R. Geoscience 334 (2002) 529–536.     

La prévision du climat : de l'échelle saisonnière à l'échelle décennale

The atmosphere and the ocean are subject to many dynamical instabilities, which limit the time during which their behaviour can be deterministically forecasted. At longer timescales, the atmosphere can be predicted at best using statistical methods, as a response to external forcing linked to sea- and land-surface anomalies. Climate being defined as the mean of atmospheric states, it appears that it can be predicted up to a few months in advance, which is the characteristic time of the so-called slow components of the climate system. Forecasting can sometimes be extended to longer time ranges, especially when the coupled ocean–atmosphere system exhibits internal variability modes, with characteristic times of a few years. Seasonal climate forecasting is most often based upon Monte-Carlo simulations, where the various realisations correspond to slightly different initial conditions. The present sate-of-the-art in Europe (ECMWF) and/or in the USA (IRI) allows to forecast such major phenomena, as El Niño, up to six months in advance. Finally, some parameters may exhibit predictability at still longer time-ranges (inter-annual to decadal), but only for certain regions. The example of electricity production is used to underline the potentially large economical benefit of seasonal climate forecasting. To cite this article: J.-C. André et al., C. R. Geoscience 334 (2002) 1115–1127.

Résumé

L'atmosphère et l'océan sont le siège d'instabilités dynamiques, qui limitent la durée pendant laquelle il est possible d'en prévoir l'évolution de façon déterministe. Au-delà, l'atmosphère n'est plus prévisible, au mieux, que de façon statistique, en fonction du forçage externe qu'exerce(nt) sur elle l'océan et/ou la surface des continents. Le climat (au sens d'une moyenne des états atmosphériques) se révèle ainsi prévisible jusqu'à des échéances temporelles de quelques mois, échelle de temps caractéristique des composantes dites « lentes » du système climatique. La prévision peut s'étendre à des échéances parfois plus longues, dans le cas où le système couplé océan–atmosphère posséderait des modes de variabilité temporelle de périodes caractéristiques de quelques années. La prévision climatique saisonnière est très souvent construite à partir de simulations de type Monte-Carlo, avec des ensembles de réalisations utilisant des conditions initiales légèrement différentes. Dans l'état actuel de ces prévisions, qu'elles soient réalisées en Europe (CEPMMT) ou aux États-Unis (IRI), il est possible de prévoir environ six mois à l'avance un certain nombre de phénomènes climatiques, en particulier ceux liés aux épisodes dits « El Niño », pour lesquels l'amplitude des variations est suffisamment importante. Il existe, par ailleurs, une prévisibilité à encore plus longue échéance (inter-annuelle à décennale), mais seulement pour certains paramètres et certaines régions. L'exemple de la production d'électricité montre l'importance économique potentielle très grande de la prévision climatique saisonnière. Pour citer cet article : J.-C. André et al., C. R. Geoscience 334 (2002) 1115–1127.