Use of remote sensing, GIS and groundwater monitoring to estimate artificial groundwater recharge in Riyadh, Saudi Arabia

This study aims to estimate artificial recharge of groundwater by using remote sensing technology, geographical information systems, and groundwater surveys. This study is part of the King Fahd project for rainfall and runoff water harvesting, within the premises of Alilb Dam in Diriyah to the west of Riyadh. Digital elevation models were obtained with the help of aerial photography from the year 2007. These models were used to delineate watershed. Average rainfall was calculated using isoheytal method, and the area of each of the storage lakes was measured using SPOT 5 satellite images from 2007. Fluctuations in groundwater levels, evaporation, and infiltration rates were used to determine the water balance for the purpose of estimating of artificial recharge. Artificial recharge rates were found to surpass natural recharge from rainfall. Recharge wells caused a reduction in the effect of evaporation on storage lakes and helped in supplying water to the groundwater reservoir. Moreover, 80% and 86% of the rainwater was found to be available for artificial recharge in Alilb at 2005 and 2007, respectively. The study recommends the establishment of strategic projects for water storage using artificial recharge wells, an increase in the number of monitoring wells around the dams, and the monitoring of hydrochemical changes in groundwater both before and after the artificial recharge. It also recommends the erection of a weather station in the northwest of Wadi Hanifa.     

G23FM: a tool for meshing complex geological media

We present G23FM, a mesh generation tool for discretizing two- and three-dimensional complex fractured geological media. G23FM includes different techniques to generate finite element grids that maintain the geometric integrity of input surfaces, and geologic data and produce optimal triangular/tetrahedral grids for flow and transport simulations. G23FM generates grid for two-dimensional cross-sections, represents faults and fractures, for three-dimensional fractured media, and has the capability of including finer grids. Different examples are presented to illustrate some of the main features of G23FM.     

A 2-year study of soil tillage and cattle manure application effects on soil fauna populations under Zea Mays cultivation, in western Iran (Sanandaj)

The absence of environmentally sensitive soil management systems can be considered as one of the major risks to sustainability of agricultural soils in Iran. Tillage is the most critical operation in soil management designed to achieve high crop yield, but it can adversely affect the soil fauna in several ways. In the present study, assessment of soil fauna was carried out in Western Iran in 2008 and 2009 in soil subjected to conventional (CT), minimum (MT) and no (NT)-tillage systems and amended with three levels of cattle manure (CM). Earthworm, mite, springtail and nematode populations were measured as indicators of macro, meso and micro fauna groups, respectively. Soil moisture and bulk density were also determined. Generally, low populations of soil fauna were observed consistent with expectations under similar conditions for this region. Earthworm populations were low and had a patchy distribution. Tillage and CM were found to have no effects on soil mites in both years. Soil springtails were reduced by soil tillage, indicating their sensivity to soil disturbance induced by tillage. In 2008, the nematode population was greater with application of 40 ton ha^?1 CM applications (113 N.100 g soil^?1). Soil tillage-induced disturbance reduced nematode population in 2009 (214 N.100 g soil^?1 at CT). Minimum seedbed preparation besides less soil disturbance makes MT a proper tillage system for Zea mays cultivation. Cattle manure application increased Z. mays’ biomass, but according to our results its annual application is not recommended. There were no changes in BD in both years. We conclude that in short-term studies, soil nematode populations are suitable biological indices (under similar soil and climatic conditions) for the ecological comparison of agricultural management systems in Iran.     

Experimental Study Optimizing Hole Cleaning-Cuttings Transport in Oil Drilling Engineering

With a clear understanding of the drilling fluid techniques and the cutting-taking mechanism, a new advanced model is set up for analyzing field data and quantitative forecast of cutting-taking mechanism. Therefore, a number of values affecting the drilling rate and the hole cleaning are studied over a wide range of parameters. Drilling data obtained under high-borehole-pressure conditions are analyzed to determine the causes of the reduction in rate of penetration (ROP) as the borehole pressure increases, which in some cases is caused by the buildup of rock debris under the bit. The theoretical achievement and testing conclusions can be very instructional for horizontal well drilling. Much higher annular velocities are required for effective hole cleaning in directional wells than in vertical wells. Highviscosity muds are observed to provide better transport than low-viscosity muds.     

Magmatic processes and mixing origin of andesite: Miocene Karamağara volcanics,Central Anatolia,Turkey

The Miocene Karamağara volcanics (KMV) crop out in the Saraykent region (Yozgat) of Central Anatolia. The KMV include four principal magmatic components based on their petrography and compositional features: basaltic andesites (KMB); enclaves (KME); andesites (KMA); and dacites (KMD). Rounded and ellipsoidal enclaves occur in the andesites, ranging in diameter from a few millimetres to ten centimetres. A non‐cognate origin for the enclaves is suggested due to their mineralogical dissimilarity to the enclosing andesites. The enclaves range in composition from basaltic andesite to andesite. Major and trace element data and primitive mantle‐normalized rare‐earth element (REE) patterns of the KMV exhibit the effects of fractional crystallization on the evolution of the KME which are the product of mantle‐derived magma. The KMA contain a wide variety of phenocrysts, including plagioclase, clinopyroxene, orthopyroxene, hornblende and opaque minerals. Comparison of textures indicates that many of the hornblende phenocrysts within the KMA were derived from basaltic andesites (KMB) and are not primary crystallization products of the KMA. Evidence of disequilibrium in the hybrid andesite includes the presence of reacted hornblendes, clinopyroxene mantled by orthopyroxene and vice versa, and sieve‐texture and inclusion zones within plagioclase. The KMV exhibit a complex history, including fractional crystallization, magma mixing and mingling processes between mantle and crust‐derived melts. Textural and geochemical characteristics of the enclaves and their hosts require that mantle‐derived basic magma intruded the deep continental crust followed by fractional crystallization and generation of silicic melts from the continental material. Hybridization between basic and silicic melts subsequently occurred in a shallow magma chamber. Modelling of major element geochemistry suggests that the hybrid andesite represents a 62:38 mix of dacite and basaltic andesite. The implication of this process is that calc‐alkaline intermediate volcanic rocks in the Saraykent region represent hybrids resulting from mixing between basic magma derived from the mantle and silicic magma derived from the continental crust. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures

Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm.     

Application of a revised hydrostratigraphical classification and nomenclature to the Mesozoic and Cenozoic succession of Saudi Arabia

At the dawn of the 21st century many parts of the world are suffering a shortage of water resources. Arid and semi-arid areas in particular are facing challenges and increasing pressure is being put on their groundwater management plans. Such a situation created an urgent need to put forward a design for classification and nomenclature that can help to differentiate highly productive zones from small locally exploited productive zones. A previously proposed hydrostratigraphical classification and nomenclature scheme is revised, modified, and applied to the Mesozoic-Cenozoic sedimentary succession of Saudi Arabia. The scheme utilised in the present paper is modified to include, in addition, productivity, water quality criteria, usability, and recharge as criteria reflecting the hydrogeological importance of the productive zones.The Triassic Sudair Mega-aquitard forms the lower hydrostratigraphic boundary for the Ad-Dahna'a Aquasystem. The Ad-Dahna'a Aquasystem comprises two aquagroups: the Riyadh and Rub'Al-Khali Aquagroups. Each aquagroup is in turn subdivided into superaquifers, aquifers and possibly subaquifers, separated by aquitards of different ranks. The physical and hydrogeological characteristics of each unit are also discussed.     

Oil spill sorption using carbonized pith bagasse: trial for practical application

In the present work, an attempt was made to provide an efficient, easily deployable method of cleaning up oil spills and recovering of the oil. Carbonized pith bagasse, a relatively abundant and inexpensive material is currently being investigated as an adsorbent to remove contaminants “oil” from water. Fibers extracted from bagasse and carbonized at 300 °C were found to have a high performance for sorption and recovery of light, heavy oils and even the viscous ones. The physical properties of pith bagasse were investigated using scanning electronic microscope to show the inner and the outer surface and the cross section area of the pith bagasse and thermo gravimetric analyzer to investigate the degradation profile of the pith bagasse. The carbonized pith bagasse was packed into a polypropylene bag and its sorption behavior was studied. A comparison was made between the prepared pad and the commercial sorbents show that the pad containing carbonized pith bagasse has higher sorption capacity in comparison to the commercial sorbents. The pad exhibited high oil retention ability and a high selectivity for the oils over the water. The pad showed a possibility of reuse for eight times. The sorption capacity of the pads containing carbonized pith bagasse was found to increase with increasing the time of sorption till it reaches the maximum value at the time of sorption equal to 60 min.     

Spatial variability of soil erodibility and its correlation with soil properties in semi-arid mountainous watershed,Saudi Arabia

Soil erodibility values are best estimated from long-term direct measurements on runoff-plots; however, in lack of field tests, these values can be estimated using relationships based on physico-chemical soil properties. The study objective was to assess the erodibility and its correlation with soil properties. The average erodibility value was estimated 0.043 t ha h ha^?1 MJ^?1 mm^?1. The areas with heavy textured soil and low organic matter content had the lowest values of erodibility. The erodibility decreases as the sand content increases, whereas silt showed a positive correlation. The erodibility factors and its relation to soil properties were evaluated using multiple regression analysis. Results revealed that sand and organic matter content of soil combinedly explained 78% of variation. Altitudinal increases also seem to affect the soil texture. This study has demonstrated that soil properties and erodibility values can be used as assistance for soil conservation practices and modelling of landscape processes.     

Influence of diaphragm wall installation on the numerical analysis of deep excavation

This paper presents a numerical analysis of the influence of initial stress state on the response of deep excavation supported by retaining wall. Indeed, the influence of diaphragm wall installation prior to excavation works may affect the soil response and lateral wall deflection induced by excavation process. The first part of this paper gives a short review of the numerical methods aimed to reproduce the retaining wall installation. Numerical analysis of a deep excavation in two‐dimensional and three‐dimensional conditions is then performed according to the methods previously presented. In three‐dimensional conditions, diaphragm wall installation is performed considering a sequence of panels, described by their number and length. Results of three‐dimensional calculations confirm that stress state is disturbed by wall installation, which has a sensitive effect on the ground response induced by soil excavation. It is also noted that these results are not easily reproduced in two‐dimensional conditions. Copyright © 2012 John Wiley & Sons, Ltd.