TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel

Since the 1990s a large number of sinkholes have appeared in the Dead Sea (DS) coastal area. Sinkhole development was triggered by the lowering of the DS level. In the literature the relationship between the sinkholes and the DS level is explained by intrusion of relatively fresh water into the aquifer thereby dramatically accelerating the salt dissolution with creation of subsurface caverns, which in turn cause sinkholes. The main goal of our project was detection and localization of relatively fresh groundwater. During our study we used the transient electromagnetic method (TEM) to measure the electrical resistivity of the subsurface. As a test site we selected Nahal Hever South which is typical for the DS coast. Our results show that resistivity of the shallow subsurface reflects its vertical and lateral structure, e.g., its main hydrogeological elements explain the inter-relations between geology, hydrogeology, and sinkholes. The TEM method has allowed detailed differentiation of layers (clay, salt, etc.) in the subsurface based on their bulk resistivity. The 10 m-thick salt layer composed of idiomorphic crystals of halite deposited during the earlier Holocene period was extrapolated from borehole HS-2 through the study area. It was found that in Nahal Hever the typical value of the bulk resistivity of clay saturated with the DS brine varies between 0.2 and 0.3 Ωm, whereas saturated gravel and sandy sediments are characterized by resistivity between 0.4 and 0.6 Ωm. The high water salinity of the aquifer (enveloping the salt layer) expressed in terms of resistivity is also an important characterization of the sinkhole development mechanism. The electrical resistivity of the aquifer in the vicinity of the salt unit and its western border did not exceed 1 Ωm (in most cases aquifer resistivity was 0.2-0.6 Ωm) proving that, in accordance with existing criteria, the pores of the alluvial sediments are filled with highly mineralized DS brine. However, we suggest that the criterion of the aquifer resistivity responsible for the salt dissolution should be decreased from 1 Ωm to 0.6 Ωm corresponding to the chloride concentration of approximately 100 g/l (the chloride saturation condition reaches 224 g/l in the northern DS basin and 280 g/l in the southern one).Based on TEM results we can reliably conclude that in 2005, when most of sinkholes had appeared at the surface, salt was located within a very low resistivity environment inside sediments saturated with DS brine. Intrusion of relatively fresh groundwater into the aquifer through the 600 × 600 m²area affected by sinkholes has not been observed.     

Advanced Mathematical Models to Predict the Compaction Properties of Coarse-Grained Soils from Various Physical Properties

An essential task in the process of construction is the determination of compaction properties of soils. Many years of laboratory test experience strengthen our belief in the existence of predictive equations that govern the compaction characteristics of soils. An advanced mathematical model developed in this research in order to uncertain the governing equations. An advanced mathematical model developed in this research in order to uncertain the governing equations. Through a comparative study among a Multiple Linear Regression (MLR) model, an Artificial Neural Network (ANN) model, Extreme Learning Machine (ELM) and a Support Vector Machine (SVM) model, the best predicting model was determined. For this purpose, Six hundred and six (606) samples collected and split into a dataset used for training the models and another used for validation of the derived model. 8 neural networks with a varying number of hidden layers and a varying number of nodes in hidden layers were employed. In ELM 1 hidden layer with varying number of units were employed. It was found that the equations derived from the ELM models described the relationship with superiority over multiple regression, ANN and SVM models for Maximum Dry Density and MLR models described the relationship with superiority over ANN, ELM and SVM models for Optimum Moisture Content.     

Impact of the geomorphology and soil management on the development of waterlogging in closed drainage basins of Egypt and Saudi Arabia

This paper investigates the interplay of the catchment geomorphology, hydrology and soil properties on the development of waterlogging and land degradation within different dry land catchments in Egypt and Saudi Arabia. Multi-temporal remote sensing data of the Landsat Thematic Mapper and Enhanced Thematic Mapper were collected and processed to detect the land cover changes and development of cultivations within the two areas. The Shuttle Radar Topography Mission Digital Elevation Model were hydrologically processed to delineate the catchment morphometrical parameters and to examine the spatial distribution of cultivated fields and their relation to the extracted drainage networks. The fluvial channels of the Farafra Oasis have largely been obliterated by the prevailing aridity and often buried under aeolian deposits. The soil of these areas are mainly lithic with a high calcium carbonate content, thus limiting the downward percolation of excess irrigation water and therefore develop perched water table and seepage through the buried fluvial channels. On the other hand, the cultivations of Tabuk catchment (which has similar geomorphologic setting to the Farafra Oasis) have shown no signs for waterlogging. This situation is could be related to the different soil properties; the spatial distribution of cultivated areas and the adopted irrigation methods via pivots. The inactive alluvial channels, landforms and irrigation methods have to be considered when planning for a new cultivation in dry land catchments to better control waterlogging hazard. The ‘dry-drainage’ concept can be implemented as the drainage and seepage water, which can be conveyed into certain abandoned playas for evaporation.     

Interactions between the surface water and groundwater in the western shoreline of Lake Nasser, Upper Egypt

The present study indicates that the factors controlling the hydraulic relation between surface water and groundwater at the western lake shoreline change from one locality to another. This depends upon the lithological characteristics and the major structures. In the southern sectors, sedimentation at the bottom and sides of the lake prevents the water movement to the Nubian sandstone aquifer. The potentiometric map reveals that the water level altitudes range between 170 m in the vicinity of the lakeshore line and 110 m west of the lake. The groundwater flow lines show that the main recharge to the aquifer comes from the southwest direction, as well as from the lake inland to variable distances (about 30 Km). During the present study, Darcy’s law was applied to calculate the recharge from the western shoreline of Lake Nasser to the adjacent Nubian aquifer. The maximum value of seepage was at Garf Hussein (27.71?×?10⁶ m³/year), which may be related to high permeability and hydraulic gradient. Also, it may be related to the N–S strike faults that cut the area on both sides of the Lake, and the groundwater is expected to have free circulation through the faults of this trend. The minimum value was recorded in Adindan section (0.61?×?10⁶ m³/year). This may be related to the limited recharge from the lake to the aquifer, due to the sedimentation that dislocates this recharge.     

Soil and bedrock distribution estimated from gully form and frequency: A GIS-based decision-tree model for Lebanon

Torrential rainfall and relatively sparse vegetation in the Mediterranean region result in the development of gully systems and land degradation, notably on lands with specific types of soil and bedrock. This paper proposes a decision-tree model to predict the distribution of soil and bedrock susceptible to gully erosion (white Rendzinas and marly rocks) from the form and frequency of gullies. The study area is located in Lebanon and the model is linked to GIS. V-fold cross-validation of the pruned model indicates that gully features including cross-section size and shape, network frequency, types of meandering, and catchment area can explain 80% of variance in soil/rock properties. The overall accuracy of the soil/rock map was estimated to be ca. 87%. The proposed model is relatively simple, and may also be applied to other areas. It is particularly useful when information about soil and rock obtained from conventional field surveys is limited.     

Pressuremeter test parameters of a compacted illitic soil under thermal cycling

The incorporation of heat exchangers in geostructures changes the temperature of the adjacent soil, raising important issues concerning the effect of temperature variations on hydro-mechanical soil behaviour. The objective of this paper is to improve the understanding and quantification of the impact of temperature variation on the bearing capacity of thermo-active piles. Currently, the design of deep foundations is based on the results of in situ penetrometer or pressuremeter tests. However, there are no published data on the effect of temperature on in situ soil parameters, preventing the specific assessment of the behaviour of thermo-active piles. In this study, an experimental device is developed to perform mini-pressuremeter tests under controlled laboratory conditions. Mini-pressuremeter tests are performed on an illitic soil in a thermo-regulated metre-scale container subjected to temperatures from 1 to 40 °C. The results reveal a slight decrease in the pressuremeter modulus (E _p) and a significant decrease in the creep pressure (p _f) and limit pressure (p _l) with increasing temperature. The results also reveal the reversibility of this effect during a heating–cooling cycle throughout the investigated temperature range, whereas the effect of a cooling–heating cycle was only partially reversible. In the case of several thermal cycles, the effect of the first cycle on the soil parameters is decisive.     

Understanding groundwater chemistry in Mediterranean semi-arid system using multivariate statistics techniques and GIS methods: case of Manouba aquifer (Northeastern Tunisia)

The rapid urbanization and industrialization of the Manouba plain (Northeastern Tunisia), the extensive agricultural expansion and the succession of dry years during recent decades have exerted greatly load on the water needs and lead to groundwater quality degradation. The aim of this study is to evaluate the processes controlling the groundwater mineralization of the shallow aquifer for determining its suitability for drinking and agricultural purposes. For establishing that, we combine several geological, hydrological and hydrochemical data with geostatistical techniques. The samples were collected at 17 sites covering 230 km² of the study area and analyzed for major and trace components. The total dissolved solid (TDS) content ranges from 1372 to 3999 mg/l. The results of Piper diagram indicate that Na⁺/Cl^? and Ca²⁺?>?Na⁺/SO₄ ^2? were the main dominant water types localized in the sloping sides of the watershed and near the saline depression; the suitability for irrigation use was also evaluated. The high concentrations of nitrates and chlorides are indicators of anthropogenic pollution, like the agricultural over application of nitrogen fertilizers and the discharge of domestic and industrial wastewater. Saturation indexes calculated by using PHREEQC (USGS) program show that groundwaters are undersaturated with evaporitic minerals (halite, gypsum) and saturated with carbonates (calcite, aragonite). The use of principal component analysis and hierarchical cluster analysis has shown that two main factors accounting 67.13% of the information of variability within the dataset confirm the existence of dissolution of evaporitic minerals and the mechanisms of nitrate increasing the salinity of the Manouba groundwater.     

Parametric study using FEM for the stability of the RCC Tannur dam

Parametric study was carried out in this paper using the Finite Element Method (FEM) to study the effects of variations of the two parameters, i.e., variations in batter slope and foundation properties under static and dynamic loading conditions on the distribution of the peak stresses within the dam and foundation soils. The response spectra for the operating basis earthquake (OBE) and the maximum credible earthquake (MCE) were used in the dynamic analysis. It is shown that the distribution of peak stresses across base of dam and the extent of tensile zone on the foundation soils at the foundation level vary with the variations in the batter slope. Generally, as the batter slope increases the extent of the tensile zone at the foundation level decreases but not necessarily the maximum tensile stresses. However, an optimum batter slope could be achieved based on the results of this parametric study. It should be noted that the greatest tension is developed in the rock adjacent to the toe of the batter. For the foundation mechanical properties, it has been shown that the assessment of accurate soil mechanical properties has a great effect on the computed design stresses from finite element analysis. However, the greatest effect of the foundation properties was on the stresses at the base of dam while it has a negligible effect within the dam body.     

Dual-porosity modeling of groundwater recharge: testing a quick calibration using in situ moisture measurements,Areuse River Delta,Switzerland

A simple method for calibrating the dual-porosity MACRO model via in situ TDR measurements during a brief infiltration run (2.8 h) is proposed with the aim of estimating local groundwater recharge (GR). The recharge was modeled firstly by considering the entire 3 m of unsaturated soil, and secondly by considering only the topsoil to the zero-flux plane (0–0.70 m). The modeled recharge was compared against the GR obtained from field measurements. Measured GR was 313 mm during a 1-year period (15 October 1990–15 October 1991). The best simulation results were obtained when considering the entire unsaturated soil under equilibrium conditions excluding the macropore flow effect (330 mm), whereas under non-equilibrium conditions GR was overestimated (378 mm).Sensitivity analyses showed that the investigation of the topsoil is sufficient in estimating local GR in this case, since the water stored below this depth appears to be below the typical rooting depth of the vegetation and is not available for evapotranspiration. The modeled recharge under equilibrium conditions for the 0.7-m-topsoil layer was found to be 364 mm, which is in acceptable agreement with measurements.

---

Resumen Se propone un método sencillo para calibrar el modelo de doble porosidad "MACRO" mediante medidas in-situ obtenidas por TDR durante un breve ensayo de infiltración (2,8 horas), con el objetivo de estimar la recarga local al acuífero. Ésta ha sido modelada de dos formas: considerando los 3 m de suelo no saturado y empleando sólo desde la capa superior hasta el plano de flujo nulo (de 0 a 0,70 m). Se compara la recarga modelada con la recarga local medida en campo, la cual fue de 313 mm durante un ciclo anual (del 15 de octubre de 1990 al 15 de octubre de 1991). Las mejores simulaciones corresponden a la hipótesis de columna entera no saturada en condiciones de equilibrio, excluyendo el efecto de macroporos (valor de 330 mm), mientras que el resultado obtenido para condiciones de no equilibrio en la recarga local está sobreestimado (378 mm).Los análisis de sensibilidad muestran que la investigación del horizonte superior del suelo es suficiente para estimar la recarga local en este caso, ya que el agua almacenada por debajo de esta profundidad parece estar fuera del alcance típico de las raíces de la vegetación y no puede ser evapotranspirada. La recarga modelada en condiciones de equilibrio para la capa superior de 0,70 m de espesor es de 364 mm, valor aceptable respecto a las medidas.

---

Résumé Une méthode simple pour la calibration du modèle à double porosité MACRO par des mesures TDR in situ durant un bref essai d'infiltration (2.8 h) a été proposée pour l'estimation locale de la recharge de la nappe (RN). La RN a été d'abord simulée en tenant compte de toute la zone non saturée (3 m) et ensuite, en considérant uniquement la couverture du sol entre zéro et le plan du flux nul (0.70 m). La RN simulée a été comparée à la RN observée. La RN mesurée durant une année (15 octobre 1990–15 octobre 1991) était de 313 mm. Les meilleures simulations ont été obtenues en tenant compte de toute la zone non saturée sous les conditions d'équilibre excluant le flux préférentiel (330 mm). Sous les conditions de non équilibre, la RN a été surestimée (378 mm).Les analyses de sensitivité ont montré que l'investigation de la couverture du sol est suffisante pour l'estimation locale de la RN du fait que l'eau traversant le plan du flux nul se trouverait sous la zone des racines et échapperait à l'évapotranspiration. La RN simulée sur les 0.70 m du sol sous les conditions d'équilibre était de 364 mm, ce qui est comparable aux mesures.