Showing papers by "Philippe Davy published in 2005"

Discharge, discharge variability, and the bedrock channel profile

Abstract: Long-term bedrock incision is driven by daily discharge events of variable magnitude and frequency, with ineffective events below an incision threshold. We explore theoretically how this short-term stochastic behavior controls long-term steady state incision rates and bedrock channel profiles, combining a realistic frequency-magnitude distribution of discharge with a deterministic, detachment-limited incision model in which incision rate is a power function of basal shear stress above a critical shear stress. Our model predicts a power law relationship between steady state slope and drainage area consistent with observations. The exponent of this power law is independent of discharge mean and variability, while the amplitude factor, which controls mountain belt relief, is a power law function of mean runoff (with an exponent of -0.5) and a complex function of runoff variability. In accordance with evidence that incision occurs between 6 and 20% of time in rapidly incising rivers (>1 mm/yr) our model predicts that channel steepness is virtually insensitive to runoff variability. Runoff variability can only decrease channel steepness for very slow incision rates and/or weak lithologies. The relationship between channel steepness and incision rate is always a power law whose exponent depends on the channel cross-sectional geometry and runoff variability. This contradicts models neglecting discharge stochasticity in which the steepness-incision scaling is set by the incision law exponent. Our results suggest that changes in climate variability cannot explain an increase in bedrock incision rates during the Late Cenozoic within the context of a detachment limited model.

How did the Messinian Salinity Crisis end

Abstract: The cause of the desiccation of the Mediterranean Sea during the Messinian Salinity Crisis has been widely debated, but its re-flooding remains poorly investigated. Interpretations generally involve tectonic collapse of the Strait of Gibraltar or global sea-level rise, or even a combination of both. The dramatic sea-level fall in the Mediterranean has induced deep fluvial incision all around the desiccated basin. We investigate erosion dynamics related to this base level drop by using the numerical simulator EROS. We show that intense regressive erosion develops inevitably in the Gibraltar area eventually inducing the piracy of the Atlantic waters by an eastward-flowing stream and the subsequent re-flooding of the Mediterranean.

How to model shallow water‐table depth variations: the case of the Kervidy‐Naizin catchment, France

Abstract: The aim of this work is threefold: (1) to identify the main characteristics of water-table variations from observations in the Kervidy-Naizin catchment, a small catchment located in western France; (2) to confront these characteristics with the assumptions of the Topmodel concepts; and (3) to analyse how relaxation of the assumptions could improve the simulation of distributed water-table depth. A network of piezometers was installed in the Kervidy-Naizin catchment and the water-table depth was recorded every 15 min in each piezometer from 1997 to 2000. From these observations, the Kervidy-Naizin groundwater appears to be characteristic of shallow groundwaters of catchments underlain by crystalline bedrock, in view of the strong relation between water distribution and topography in the bottom land of the hillslopes. However, from midslope to summit, the water table can attain a depth of many metres, it does not parallel the topographic surface and it remains very responsive to rainfall. In particular, hydraulic gradients vary with time and are not equivalent to the soil surface slope. These characteristics call into question some assumptions that are used to model shallow lateral subsurface flow in saturated conditions. We investigate the performance of three models (Topmodel, a kinematic model and a diffusive model) in simulating the hourly distributed water-table depths along one of the hillslope transects, as well as the hourly stream discharge. For each model, two sets of parameters are identified following a Monte Carlo procedure applied to a simulation period of 2649 h. The performance of each model with each of the two parameter sets is evaluated over a test period of 2158 h. All three models, and hence their underlying assumptions, appear to reproduce adequately the stream discharge variations and water-table depths in bottom lands at the foot of the hillslope. To simulate the groundwater depth distribution over the whole hillslope, the steady-state assumption (Topmodel) is quite constraining and leads to unacceptable water-table depths in midslope and summit areas. Once this assumption is relaxed (kinematic model), the water-table simulation is improved. A subsequent relaxation of the hydraulic gradient (diffusive model) further improves water-table simulations in the summit area, while still yielding realistic water-table depths in the bottom land.

Brittle-ductile coupling: Role of ductile viscosity on brittle fracturing

Abstract: Localized or distributed deformations in continental lithosphere are supposed to be triggered by rheological contrasts, and particularly by brittle-ductile coupling. A plane-strain 2D finite-element model is used to investigate the mechanical role of a ductile layer in defining the transition from localized to distributed fracturing in a brittle layer. The coupling is performed through the shortening of a Von Mises elasto-visco-plastic layer rimed by two ductile layers. By increasing the viscosity of the ductile layers by only one order of magnitude, the fracturing mode in the brittle layer evolves from localized (few faults) to distributed (numerous faults), defining a viscosity-dependent fracturing mode. This brittle-ductile coupling can be explained by the viscous resistance of the ductile layer to fault motion, which limits the maximum displacement rate along any fault connected to the ductile interface. An increase of the viscosity will thus make necessary new faults nucleation to accommodate the boundary shortening rate.