Soil cracking effects on hydrological and erosive processes: a study case in Mediterranean cultivated vertisols

Shrink–swell soils, such as those in a Mediterranean climate regime, can cause changes in terms of hydrological and erosive responses because of the changing soil water storage conditions. Only a limited number of long‐term studies have focused on the impacts on both hydrological and erosive responses and their interactions in an agricultural environment. In this context, this study aims to document the dynamics of cracks, runoff and soil erosion within a small Mediterranean cultivated catchment and to quantify the influence of crack processes on the water and sediment supplied to a reservoir located at the catchment outlet using water and sediment measurements at a cultivated field outlet as baseline. Detailed monitoring of the presence of topsoil cracks was conducted within the Kamech catchment (ORE OMERE, Tunisia), and runoff and suspended sediment loads were continuously measured over a long period of time (2005–2012) at the outlets of a field (1.32 ha) and a catchment (263 ha). Analysis of the data showed that topsoil cracks were open approximately half of the year and that the rainfall regime and water table level conditions locally control the seasonal cracking dynamics. Topsoil cracks appeared to seriously affect the generation of runoff and sediment concentrations and, consequently, sediment yields, with similar dynamics observed at the field and catchment outlets. A similar time lag in the seasonality between water and sediment delivery was observed at these two scales: although the runoff rates were globally low during the presence of topsoil cracks, most sediment transport occurred during this period associated with very high sediment concentrations. This study underlines the importance of a good prediction of runoff during the presence of cracks for reservoir siltation considerations. In this context, the prediction of cracking effects on runoff and soil erosion is a key factor for the development of effective soil and water management strategies and downstream reservoir preservation. Copyright © 2016 John Wiley & Sons, Ltd.     

A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

A robust debris-flow and GLOF risk management strategy for a data-scarce catchment in Santa Teresa,Peru

The town of Santa Teresa (Cusco Region, Peru) has been affected by several large debris-flow events in the recent past, which destroyed parts of the town and resulted in a resettlement of the municipality. Here, we present a risk analysis and a risk management strategy for debris-flows and glacier lake outbursts in the Sacsara catchment. Data scarcity and limited understanding of both physical and social processes impede a full quantitative risk assessment. Therefore, a bottom-up approach is chosen in order to establish an integrated risk management strategy that is robust against uncertainties in the risk analysis. With the Rapid Mass Movement Simulation (RAMMS) model, a reconstruction of a major event from 1998 in the Sacsara catchment is calculated, including a sensitivity analysis for various model parameters. Based on the simulation results, potential future debris-flows scenarios of different magnitudes, including outbursts of two glacier lakes, are modeled for assessing the hazard. For the local communities in the catchment, the hazard assessment is complemented by the analysis of high-resolution satellite imagery and fieldwork. Physical, social, economic, and institutional vulnerability are considered for the vulnerability assessment, and risk is eventually evaluated by crossing the local hazard maps with the vulnerability. Based on this risk analysis, a risk management strategy is developed, consisting of three complementing elements: (i) standardized risk sheets for the communities; (ii) activities with the local population and authorities to increase social and institutional preparedness; and (iii) a simple Early Warning System. By combining scientific, technical, and social aspects, this work is an example of a framework for an integrated risk management strategy in a data scarce, remote mountain catchment in a developing country.     

Nonlinear Mechanical Behavior Analysis of Flexible Lateritic Pavements of Senegal (West Africa) by FEM for M.-E. Pavement Design

Pavement design in Senegal is based on a linear elastic behavior of pavement materials and the hypothesis of a static loading. However, previous works on the mechanical behavior of road materials showed that this one is reversible after several loading cycles and depends on the applied stress. The described behavior is from then on, purely nonlinear. One of the objectives of this research is to determine the parameters of response of lateritic pavement materials submitted to road traffic by using FEM. Therefore, experiments were made on gravel lateritic soils from Dougar, Sébikotane, Mont-Rolland, Pâ Lo and Ngoundiane. The Young’s modulus of the materials was defined in unconfined compression test while repeated load triaxial test was performed to determine the resilient modulus of the gravels and the appropriate model (Uzan model). An implementation was realized with Cast3M©. The importance of the nonlinearity was revealed in a very clear way and was crucial in the construction of the calculation algorithm. The observations for certain conditions showed that the values of the critical responses are more important for the linear model than for the nonlinear model. However, this trend should be validated by further studies.     

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.     

The position of the lowest levels in the boundary layer of atmospheric circulation models

Abstract

The error associated with the position of the lowest wind level in atmospheric boundary‐layer modelling is studied in connection with vertical resolutions typical of parametrization schemes for atmospheric circulation models. The test case is the neutrally stratified steady‐state boundary layer. Finite‐difference and finite‐element schemes of two types are used: in one case the lowest wind level is centred with respect to the surface and the lowest internal level where the shear stress is calculated, in the other case the lowest wind level is set very close to the surface (5 m). It is found that schemes of the latter type underestimate the friction force at the lowest level and consequently overestimate the wind and the surface stress. This error is largest at low resolution since it is due to the uncentring of the lowest wind level with respect to the stress levels. The error in schemes of the former type is different, and is associated with the determination of the surface stress from a wind at a height that may exceed the extent of the surface layer. For all neutral cases, this error can be made small by adding a corrective term to the traditional logarithmic formulation. This paper shows that considerably more error is generated by uncentred differencing than by deepening the surface layer.     

Internal erosion of chemically reinforced granular materials: a mathematical modeling approach

Internal erosion (IE) affects the stability of natural and reinforced materials by causing instability within their granular structure. The dislodgement and transport of eroded particles affect both the particulate concentration of eroding fluid and the pore network of eroded material. In this study, we examined these modifications using a transport model with a finite element code. First, IE tests on chemically reinforced sand columns were performed to obtain information about eroded material loss of mass, particulate concentration of effluent, porosity and permeability modifications, and existing IE stages. Second, based on experimental results, a mathematical one‐dimensional model has been formulated to monitor the evolution and spatial distribution of erodible solids, fluidized particles, porosity, permeability, and seepage stresses. The model consists of a set of coupled nonlinear differential equations solved in sequence. It provides valuable information about the extent and the dynamics of structural changes, which can be used to estimate an IE time for the hydraulic work to reach failure. Copyright © 2011 John Wiley & Sons, Ltd.     

Distribution and abundance of ichthyoplankton in the Manicouagan River estuary,a tributary of the lower St. Lawrence estuary

The species composition and relative abundance of ichthyoplankton were investigated during summer 1986 at four stations along the salinity gradient in the Manicouagan River estuary, a tributary of the lower St. Lawrence estuary. Physical characteristics of water masses indicated the presence of a strong saline front (>10‰ per km) delineating the freshwater and marine section of the Manicouagan estuary. The estuary supports a depauperate ichthyoplankton community, including four species of pelagic fish eggs and eight species of fish larvae. Species richness increased with salinity. The ichthyoplankton fauna can be divided into two distinct groups: freshwater and marine. These two groups result initially from spawning preferences exhibited by the different species abundance of freshwater larvae was maximal at the head of the estuary and marine larvae were most abundant at the most saline station. The length frequency distribution suggests that marine larvae are not effectively retained within the estuary. The Manicouagan estuary cannot be considered as a major spawning site nor an important nursery zone for any fish found in this area.     

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.