Catherine Mériaux

DYKE PROPAGATION WITH DISTRIBUTED DAMAGE OF THE HOST ROCK

Observations of off-plane inelastic deformation around dykes motivate consideration of models of fluid-driven crack propagation in a solid which can undergo material degradation, or damage. The application to dyke propagation of a recently proposed damage rheology [Lyakhovsky et al., J. Geophys. Res. 102 (1997) 27635‐27649] based on thermodynamical principles and experimental measurements is discussed. The rate of accumulation of damage in this rheology is the product of a material-dependent parameter cd and the square of the strain. For geological values, a dimensionless parameter cd=DP characterizing the ratio of a damage timescale to a flow timescale is very small, where is the magmatic viscosity and 1P the driving pressure. As a result, significant rates of damage are confined to a small region near the dyke tip, where the strain is large. Consideration of possible singularities in near-tip solutions, shows that the rate of propagation is governed by the viscous fluid mechanics. To a good approximation, the rate has a value equal to that given by the zero-stress-intensity solutions of previous models based on linear elastic fracture mechanics. Predictions from the damage rheology both of a narrow damage zone and of the rate of propagation are in good agreement with observations.

Two dimensional fall of granular columns controlled by slow horizontal withdrawal of a retaining wall

This paper describes a series of experiments designed to investigate the fall of granular columns in a quasi-static regime. Columns made of alternatively green and red sand layers were initially laid out in a box and then released when a retaining wall was set in slow motion with constant speed. The dependence of the dynamics of the fall on the initial aspect ratio of the columns, the velocity of the wall, and the material properties was investigated within the quasi-static regime. A change in the behavior of the columns was identified to be a function of the aspect ratio (height/length) of the initial sand column. Columns of high aspect ratio first subsided before sliding along failure planes, while columns of small aspect ratio were only observed to slide along failure planes. The transition between these two characteristic falls occurred regardless of the material and the velocity of the wall in the context of the quasi-static regime. When the final height and length of the piles were analyzed, we foun...

A two-way interaction between the Hainan plume and the Manila subduction zone

The interaction between mantle plumes and subducting slabs is well accepted, but the influence of slabs on plumes has more often been portrayed than the reverse. Here we present three-dimensional upper mantle laboratory models in which a compositional plume rises underneath a subducting plate. Slab/plume buoyancy flux ratios ranged between 7 and 18. The models exhibit a two-way interaction. While the plume conduit increasingly tilts away from the trench as a result of slab rollback-induced toroidal mantle flow, the slab subduction rate decreases as a function of the amount of plume buoyancy opposing that of the slab, which gets subducted beneath the slab. We propose that our models apply to the Hainan/Manila system and explain the recently imaged tilt of the Hainan plume by the Manila slab-induced mantle return flow. The Hainan plume could lessen the Manila subduction rate from 8 Ma into the future.

Scaling the final deposits of dry cohesive granular columns after collapse and quasi-static fall

This paper reports on laboratory experiments that were designed to investigate the collapse and quasi-static fall of dry cohesive granular columns. These experiments were compared with similar experiments that were performed with non-cohesive dry sand columns. A powder of gypsum (calcium sulphate dihydrate) was used to represent cohesive granular material. In all the experiments, the cohesive granular columns fractured and flowed in coherent blocks but, while faults remained steep in the quasi-static fall experiments, they flattened in the collapse experiments as the initial aspect ratio of the columns increased. Dilation was seen in the quasi-static fall experiments, while some air entrapment within the columns occurred in the collapse experiments. The final deposits of the cohesive granular columns were found to satisfy power law relationships as a function of the initial aspect ratio of the columns. Two asymptotes were found for the lower and higher range of initial aspect ratios, which varied between ...

Effects of Recent Minimum Temperature and Water Deficit Increases on Pinus pinaster Radial Growth and Wood Density in Southern Portugal

Western Iberia has recently shown increasing frequency of drought conditions coupled with heatwave events, leading to exacerbated limiting climatic conditions for plant growth. It is not clear to what extent wood growth and density of agroforestry species have suffered from such changes or recent extreme climate events. To address this question, tree-ring width and density chronologies were built for a Pinus pinaster stand in southern Portugal and correlated with climate variables, including the minimum, mean and maximum temperatures and the number of cold days. Monthly and maximum daily precipitations were also analyzed as well as dry spells. The drought effect was assessed using the standardized precipitation-evapotranspiration (SPEI) multi-scalar drought index, between 1 to 24-months. The climate-growth/density relationships were evaluated for the period 1958-2011. We show that both wood radial growth and density highly benefit from the strong decay of cold days and the increase of minimum temperature. Yet the benefits are hindered by long-term water deficit, which results in different levels of impact on wood radial growth and density. Despite of the intensification of long-term water deficit, tree-ring width appears to benefit from the minimum temperature increase, whereas the effects of long-term droughts significantly prevail on tree-ring density. Our results further highlight the dependency of the species on deep water sources after the juvenile stage. The impact of climate changes on long-term droughts and their repercussion on the shallow groundwater table and P. pinaster’s vulnerability are also discussed. This work provides relevant information for forest management in the semi-arid area of the Alentejo region of Portugal. It should ease the elaboration of mitigation strategies to assure P. pinaster’s production capacity and quality in response to more arid conditions in the near future in the region.

The propagation of particulate gravity currents in a V-shaped triangular cross section channel: Lock-release experiments and shallow-water numerical simulations

In this study, we investigate the motion of particulate gravity currents in a horizontal V-shaped channel. The particulate currents consisted of particles whose size varied between 0 and 100 μm but whose mean size increased. Particles were poorly sorted as the variance of the grain size distributions varied between 50 and 200. While the phases of propagation of homogeneous currents in such a geometry have been studied in the literature, this study considers the effects of the grain size on the propagation. The distance of propagation and front velocity of full-depth high-Reynolds-number lock-release experiments and shallow-water equation simulations were analyzed as the mean grain size of the initial particle distributions, defined by mass, was increased from 19 to 58 μm. Similar to the homogeneous currents, three consecutive phases of the front velocity could be identified but their characteristics and extent depend on the particle size. The initial phase, in particular, depends on a dimensionless settling number β that is defined as the ratio of two characteristic time scales, the propagation time x0/U, where U is the scale for the front speed and x0 the lock length, and the settling time h0/vs, where vs is the scale for the settling velocity and h0 the initial height of the current. For dimensionless settling numbers less than 0.001, the initial phase is characterized by a constant velocity for over about 6-7 lock lengths that is alike the initial slumping phase of perfectly constant velocity of the homogeneous currents. For dimensionless settling numbers greater than 0.001 and less than 0.015, the initial phase is no longer characterized by a constant velocity but an almost constant velocity for over about a similar 6-7 lock lengths. For dimensionless settling numbers greater than 0.015, however, as such, this phase is no longer seen. This initial phase is followed by a continuous decrease of the front advance, which results from the sedimentation of the particles. Unlike the homogeneous currents, this phase is a non-self-similar propagation. This phase is ended by a viscosity-dominated phase appearing to vary as ∼t1/7. The good agreement between the front advance of the experiments and shallow-water equation simulations demonstrates that the mean size by mass is a fairly good proxy of poorly sorted particles.

The thermal signature of subducted lithospheric slabs at the core–mantle boundary

Abstract We study the thermal structure around a cold deformable lithospheric slab as it sinks to the core–mantle boundary and migrates along it. We present analytical results for the steady thermal structure established by a steady but spatially varying motion. The analysis gives a time-like criterion for the thermal signature of a cold slab to persist by the time that the slab moves along the core–mantle boundary. The model is used to assess the feasibility of a purely thermal origin for some of the observed seismic reflectors near the core–mantle boundary. Calculations of the time-like criterion show that the dynamical conditions in our model, namely the velocity and the thickness of the descending slab, are hard to reconcile with observations of subduction and seismic features. Seismic reflections and refractions from anomalously fast regions above the core–mantle boundary could be explained as thermal slabs if the thickness of the slab at subduction was larger than 200 km or somewhat less if the slab did not split at the core–mantle boundary. A simple thermal model also predicts from mineral physics a certain correlation between S- and P-wave velocity anomalies, which is not observed. However, a purely thermal origin cannot be ruled out if the slab is buckling. This process could be in agreement with the observations: the amplitude of the seismic anomalies, the vertical extent of high-gradient zones and the P versus S comparisons. Chemical heterogeneities and phase transformations remain alternative or complementary explanations.