The stratospheric version of LMDz: dynamical climatologies, arctic oscillation, and impact on the surface climate

A climatology of the stratosphere is determined from a 20-year integration with the stratospheric version of the Atmospheric General Circulation Model LMDz. The model has an upper boundary at near 65 km, uses a Doppler spread non-orographic gravity waves drag parameterization and a subgrid-scale orography parameterization. It also has a Rayleigh damping layer for resolved waves only (not the zonal mean flow) over the top 5 km. This paper describes the basic features of the model and some aspects of its radiative-dynamical climatology. Standard first order diagnostics are presented but some emphasis is given to the model’s ability to reproduce the low frequency variability of the stratosphere in the winter northern hemisphere. In this model, the stratospheric variability is dominated at each altitudes by patterns which have some similarities with the arctic oscillation (AO). For those patterns, the signal sometimes descends from the stratosphere to the troposphere. In an experiment where the parameterized orographic gravity waves that reach the stratosphere are exaggerated, the model stratosphere in the NH presents much less variability. Although the stratospheric variability is still dominated by patterns that resemble to the AO, the downward influence of the stratosphere along these patterns is near entirely lost. In the same time, the persistence of the surface AO decreases, which is consistent with the picture that this persistence is linked to the descent of the AO signal from the stratosphere to the troposphere. A comparison between the stratospheric version of the model, and its routinely used tropospheric version is also done. It shows that the introduction of the stratosphere in a model that already has a realistic AO persistence can lead to overestimate the actual influence of the stratospheric dynamics onto the surface AO. Although this result is certainly model dependent, it suggests that the introduction of the stratosphere in a GCM also call for a new adjustment of the model parameters that affect the tropospheric variability.     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Crack nucleation in saturated porous media

Crack nucleation has been the subject of important contributions in the last two last decades. However, it seems that few attention has been granted to the case of saturated porous media. This is the question addressed in the present paper which is devoted to nucleation in traction mode. From a physical point of view, nucleation is a sudden phenomenon, so that the material response is both adiabatic and undrained. In the spirit of the variational approach, the nucleated crack is viewed as the final state of a region of space in which the material undergoes a full damage process. In traction mode, the opening of a saturated crack in undrained condition induces a drop of fluid pressure. In case of low fluid compressibility, the presence of the fluid delays the brittle failure usually associated with nucleation, as long as the fluid pressure remains above the saturation vapor pressure. Nucleation is therefore possible only if a partial vaporization of the fluid takes place.     

Magnitude of latent heat in thermally loaded clays

Temperature changes are known to induce specific couplings in clay, in particular, an anomalously high thermal pressurization in undrained conditions or a thermal compaction in drained conditions, both of which are potential threats for the mechanical stability and sealing capacity of the geomaterials. Thermodynamical analysis of those peculiar thermomechanical couplings points to a potentially important latent energy, which in turn could limit the temperature change upon heating or cooling. The direct measurement of latent energy developed during a laboratory geomechanical test is challenging. Instead, proper identification of thermal hardening in conventional experiments with temperature changes provides an alternative route to estimate latent energy. In this work, existing laboratory thermomechanical tests of clays are analyzed with a rigorous thermodynamic framework to quantify the magnitude of latent energy in thermomechanically loaded clays. A thermodynamically consistent constitutive model for fully saturated clays that combines two key features, (a) the temperature dependence of the blocked energy and (b) the framework of bounding plasticity, is proposed. The performance of the model is validated by reproducing results obtained in laboratory tests for Boom and Opalinus clays. The thermomechanical loads considered to validate the model performance were then used to estimate the percentage of work that remains latent in the clayey material during plastic yielding. We find that the magnitude of latent energy is quite significant, typically a few tens of percent of the total dissipated energy, and increases significantly with temperature. Accordingly, it is expected to play an important role in the thermomechanical response of clays.     

Provenance discrimination of fine sediments by mid-infrared spectroscopy: Calibration and application to fluvial palaeo-environmental reconstruction

Determining sediment provenance allows a better understanding of fluvial palaeo-dynamics, and identifying involved watersheds, at broad spatio-temporal scales. Conventional approaches for source identification are usually based on the physical, mineralogical, geochemical, magnetic or isotopic properties of sediments. Rapid, non-destructive and, in well-established contexts, highly accurate, mid-infrared spectroscopy is an alternative method for investigating sediment sources. The present research objectives are: (i) to use the mid-infrared spectroscopy method to discriminate the provenance of fine sediments, by applying discriminant analysis on a large set of reference samples from three different watersheds in the Upper Rhine area (associated with the Rhine, Ill and Vosges tributaries); (ii) to clarify whether the provenance spectra signatures are influenced by riverine depositional contexts (bars versus banks) and, to some extent, by grain size and/or high organic matter content; and (iii) to apply the mid-infrared spectroscopy – discriminant analysis method to a study of fluvial palaeo-dynamics and determine the provenance of palaeo-channel infillings. Three main sedimentary sources, divided into eight sub-categories, have been characterized by 196 modern reference samples from 78 collecting sites. Discriminant analysis displayed a strong separating power by classifying correctly the origin of samples without any inter-group overlap, independently from the geomorphological context (bar or bank) and associated slight changes in organic matter contents or grain size. Mid-infrared spectroscopy – discriminant analysis investigations of the palaeo-channel infill, complemented by radiocarbon dates and mineralogical data, allowed reconstructing general trends for the local morpho-sedimentary dynamics over the last ca 12 millennia.     

The junction of Fensch and Moselle rivers,France; mineralogy and composition of river materials

The Moselle river flows in the north-east of France, from Vosges Mountains to neighboring countries Luxembourg and Germany. One of its tributaries, the Fensch river, drains a highly industrialized watershed, strongly impacted by mining, smelting and surfacing activities. The objective of this work, part of a general research program on Moselle watershed (Zone Atelier Moselle) was to assess the impact of the polluted Fensch river on the global quality of the Moselle river. For that purpose, water, sediments and suspended particulate matter were sampled in both rivers, upstream and downstream the junction. Four main sampling campaigns were carried out, in winter during a flood event and in spring at low water level. On a first step, mineralogical analyses (XRD and FTIR) and chemical analyses (ICP-MS, ICP-AES), were performed on sediments, suspended particulate matters and filtered waters. Major and trace elements concentrations were obtained on two different granulometric fractions (0–2 mm and 0–50 μm) revealing the enrichment in heavy metals of fine particles. From one collecting campaign to another, seasonal variations could be evidenced on suspended matter composition even though major minerals (calcite, quartz and kaolinite) were always present. Furthermore, spatial variations were evidenced for Fensch and Moselle downstream sediments. Thus, very fine-grained sediments, poorly crystallized, displaying at the same time higher metal concentrations and higher organic matter content than in Fensch river material, were collected downstream, in a low hydrodynamic conditions zone, assumed as a preferential sedimentary zone or even as a placer. Strong correlations could be revealed between iron content and contaminant concentrations, confirming the origin of polluted material.     

Garnet reequilibration and growth in the eclogite facies and geodynamical evolution near peak metamorphic conditions

Caledonian eclogite-facies metamorphism partially reworking Grenvillian granulite-facies anorthosite allows us to study the processes of garnet reequilibration at high pressure and to reconstruct the evolution of the unit near metamorphic peak conditions. Our results indicate that eclogite-facies metamorphism happened in two successive phases: first, inherited granulitic garnet was fractured and reequilibrated from their boundaries (crystal or fracture rims); then eclogite-facies minerals were crystallised in the fractures as overgrowths on inherited garnets. The reequilibration of inherited garnets is achieved through Fe²⁺Mg⁻¹ exchange, whereas eclogite-facies garnets crystallised during the subsequent phase are notably richer in Ca than un- and re-equilibrated granulitic garnet. Pseudosection construction shows that this lack in Ca reequilibration cannot be related to variations in thermodynamic conditions (a _H2O, reacting system composition) between the two phases. From the compilation of the available data, the reequilibration of granulitic garnet seems to be controlled by the inefficient intra- and inter-granular transport properties of Ca compared to Fe²⁺ and Mg. While these kinetic factors confine garnet reequilibration to Fe²⁺Mg⁻¹ exchange, the extent of reequilibration along this exchange vector is controlled by partitioning with adjacent omphacite. On the contrary to the diffusional reequilibration of granulitic garnet that lasted for several My according to our modelling of the diffusional relaxation, the strong compositional gradients between eclogite-facies and reequilibrated garnets, which are almost unaffected by diffusional reequilibration, provide evidence that rapid exhumation followed the crystallisation of eclogite-facies minerals. We propose that the movement reversal itself, from burial to exhumation, and associated deformation and fluid flow, triggered this crystallisation event. The resulting evolution near metamorphic peak conditions is therefore strongly asymmetrical: on the one hand, the prograde diffusional relaxation profiles indicate slow movement during the last stages of burial, whereas the unaffected retrograde overgrowth indicates fast exhumation rates.     

Grain-size-sensitive flow and shear-stress enhancement at the brittle–ductile transition of the continental crust

Localized shear zones along low-angle normal faults have been identified in regions of extension at the brittle-ductile transition of the continental crust. The possibility of the strain localizing at a depth of 10 km is interpreted here as a consequence of an increase in the equivalent shear stress applied to the flow of the lower crust. This enhancement of the flow stress is seen as a prerequisite for the triggering of brittle deformation mechanisms leading to strain localization. The lower crust rheology used to examine this stress increase is strain-rate, temperature and grain-size dependent, due to the coupling of dislocation and diffusion creep. The model structure proposed consists of a top layer, the upper crust, gliding rigidly above a bottom layer, the lower crust, which deforms in simple shear. During a short time interval (1400 years), the equivalent shear stress is found to increase by a factor of up to 3 (67 MPa for anorthite and 17 MPa for quartz). For anorthite, this stress could explain the activation of a Mohr-Coulomb failure with a friction coefficient of 0.2, which is reasonable at the depth of 10 km. Dislocation creep is activated during a rapid change in the prescribed velocity, whereas diffusion creep dominates if the velocity is held constant, highlighting the importance of grain-size sensitivity for lower crustal rheology.     

Construction and analysis of higher order finite difference schemes for the 1D wave equation

In this article, we propose to: 1. Establish most of the properties conjectured in [2] about the higher order finite difference approximation of the 1D Laplace operator. 2. Generalize to any order the fourth-order accurate scheme in space and time of Shubin and Bell [20] and Cohen [6]. For this new family of 2m–2m schemes, we establish, via elementary mathematics, various stability and dispersion results that are helpful to compare these schemes to the 2–2m schemes of Anné et al. [2]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Influence of uncertain boundary conditions and model structure on flood inundation predictions

In this study, the GLUE methodology is applied to establish the sensitivity of flood inundation predictions to uncertainty of the upstream boundary condition and bridges within the modelled region. An understanding of such uncertainties is essential to improve flood forecasting and floodplain mapping. The model has been evaluated on a large data set. This paper shows uncertainty of the upstream boundary can have significant impact on the model results, exceeding the importance of model parameter uncertainty in some areas. However, this depends on the hydraulic conditions in the reach e.g. internal boundary conditions and, for example, the amount of backwater within the modelled region. The type of bridge implementation can have local effects, which is strongly influenced by the bridge geometry (in this case the area of the culvert). However, the type of bridge will not merely influence the model performance within the region of the structure, but also other evaluation criteria such as the travel time. This also highlights the difficulties in establishing which parameters have to be more closely examined in order to achieve better fits. In this study no parameter set or model implementation that fulfils all evaluation criteria could be established. We propose four different approaches to this problem: closer investigation of anomalies; introduction of local parameters; increasing the size of acceptable error bounds; and resorting to local model evaluation. Moreover, we show that it can be advantageous to decouple the classification into behavioural and non-behavioural model data/parameter sets from the calculation of uncertainty bounds.