A thermo‐mechanical model for variably saturated soils based on hypoplasticity

The paper presents a mechanical model for non‐isothermal behaviour of unsaturated soils. The model is based on an incrementally non‐linear hypoplastic model for saturated clays and can therefore tackle the non‐linear behaviour of overconsolidated soils. A hypoplastic model for non‐isothermal behaviour of saturated soils was developed and combined with the existing hypoplastic model for unsaturated soils based on the effective stress principle. Features of the soil behaviour that are included into the model, and those that are not, are clearly distinguished. The number of model parameters is kept to a minimum, and they all have a clear physical interpretation, to facilitate the model usefulness for practical applications. The step‐by‐step procedure used for the parameter calibration is described. The model is finally evaluated using a comprehensive set of experimental data for the thermo‐mechanical behaviour of an unsaturated compacted silt. Copyright © 2011 John Wiley & Sons, Ltd.     

Factors and processes of permanent gully evolution in a Mediterranean marly environment (Cape Bon,Tunisia)

Abstract

Gully erosion is considered to be one of the most important soil erosion processes in Mediterranean marly environments, but its actual contribution to total soil loss is still under discussion. The objectives of this paper are: (a) to acquire the distributed value of erosion rate in a permanent gully developed on a marly substratum in a Mediterranean environment; and (b) to quantify the key factors responsible for the spatial and temporal differences in erosion rates observed within the gully. A permanent gully located in Cap Bon (northeastern Tunisia) has been intensively and regularly monitored over a 7-year period with electronic survey equipment (total station) to give five field topographic surveys, as well as hydrological measurements at the gully outlet. The net soil loss for the 7-year period comprised a denudation of 51 m³ of sediment on the gully bank slopes, which corresponds to a mean soil loss of 61 m³ ha^?1 year^?1 or 6.1 mm year^?1. Denudation was observed on bed units with a slope gradient greater than 20%, while the remainder showed deposition. By confirming the factors involved in gully evolution, and by refining the statistical link between factors and erosion rates within the gully, the results provide important information to predict gully erosion rates in Mediterranean marly environments.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation El Khalili, A., Raclot, D., Habaieb, H., and Lamachère, J.M., 2013. Factors and processes of permanent gully evolution in a Mediterranean marly environment (Cape Bon, Tunisia). Hydrological Sciences Journal, 58 (7), 1519–1531.     

Evaluation of Drought Impacts on Vegetation Changes Trend Using Hydrological Re-Analysis Models and Remote Sensing

Water Resources - Monitoring and reviewing vegetation changes are one of the important strategies for multi-year planning in an area regarding the environmental problems of the last decade. The...     

Probabilistic analysis of extreme regional meteorological droughts by L-moments in a semi-arid environment

Regional extreme value analyses of drought characteristics provide information on probabilistic nature of drought occurrence, viewed as an essential tool in drought mitigation and planning. In this paper, L-moments are used to investigate the regional characteristics and probabilistic behavior of drought severity levels, represented by the Standardized Precipitation Index (SPI) annual minima (the minimum monthly SPI value). Rainfall data of 3, 6, 12, and 24 month time scales are investigated. A regional watershed in southwestern Iran is used as a case study area. The semi-arid nature of the study area requires appropriate selection of rainfall data. The boxplot approach is used to select those months with adequate data time series for the SPI analysis. Appropriateness of the suggested data time series is discussed in the context of the research by Wu et al. (2007). Based on the results, all of the suggested time scales are found appropriate for SPI investigations. For each time scale of interest regional homogeneity is evaluated and the best regional/sub-regional probability distribution function is selected. Regional quantiles are estimated for different time scales and their variability with respect to return period is discussed.     

The effect of the North Sea-Caspian pattern (NCP) on winter temperatures in Iran

Summary This study attempts to find possible linkages between the NCP index and the winter temperature variability over Iran. The investigation is based on statistical analysis of simple, partial and multiple correlations and also evaluation of composite maps of the extreme NCP index and maps of correlation between atmospheric variables and the temperature time series. Our results show that the NCP has a strong negative correlation with the winter temperature in Iran. Furthermore, combination of both the NCP and the AO (Arctic Oscillation) indices improve the correlations in all stations, implying both NCP and AO can be considered as major patterns for explaining the Iranian winter temperature variability. The results show that the positive NCP is associated with enhanced precipitation and cloudy conditions, consequently causing below normal temperature over Iran. The anomalies of OLR in this phase are also negative, implying a cloudy sky. For the negative NCP phase these results are completely reversed. The correlation maps indicate that the NCP is negatively/positively correlated with winter Outgoing Long-wave Radiation/precipitation over Iran. The results also show that the SLP and GPH patterns are quite different for the positive and negative NCP phases over Iran. During the negative NCP a small cyclone is formed over the Arabian Sea causing a strong easterly towards Iran. During the positive NCP this cyclone is removed. Our results show that for the positive NCP years an upper-level trough is formed over northern Iran and the eastern Mediterranean. For the negative NCP years this trough becomes weak and is located over central European regions. This trough is closely linked with the winter temperature over Iran. This is expressed by a high correlation between 500-hPa geopotential height at this region and Iranian winter temperature. Authors’ addresses: A. R. Ghasemi, Climate Research Center, Water Engineering Department, Agricultural Faculty, Shiraz University, Shiraz, Iran; D. Khalili, Water Engineering Department, Agricultural Faculty, Shiraz University, Shiraz, Iran.     

Magnetized Kelvin-Helmholtz instability in the presence of a radiation field

The purpose of this study is to analyze the dynamical role of a radiation field on the growth rate of the unstable Kelvin-Helmholtz (KH) perturbations. As a first step toward this purpose, the analyze is done in a general way, irrespective of applying the model to a specific astronomical system. The transition zone between the two layers of the fluid is ignored. Then, we perform a linear analysis and by imposing suitable boundary conditions and considering a radiation field, we obtain appropriate dispersion relation. Unstable modes are studied by solving the dispersion equation numerically, and then growth rates of them are obtained. By analyzing our dispersion relation, we show that for a wide range of the input parameters, the radiation field has a destabilizing effect on KH instability. In eruptions of the galaxies or supermassive stars, the radiation field is dynamically important and because of the enhanced KH growth rates in the presence of the radiation; these eruptions can inject more momentum and energy into their environment and excite more turbulent motions.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

An axi-symmetric infinite element for transient radial flow problems

The formulation of an axi-symmetric infinite element for transient analysis of flow problems in unbounded domain is presented. The theoretical basis as well as the implementation of the element is discussed, and the element decay function is derived using the analytical solution of a one-dimensional axially symmetric configuration. The form of decay within the element is described as a function of both time and space, and thus the hydraulic head distribution in the far field is simulated rigorously. The accuracy and the efficiency of the proposed element are demonstrated through several numerical examples in infinite media. In general, it is shown that using the present infinite element transient flow problems in unbounded domains can be simulated effectively. Copyright © 1999 John Wiley & Sons, Ltd.