Morphed block-crack preferential sedimentation in a reservoir bed: a smart design and evolution in nature

Abstract

A pedological study of the reservoir bed of Al-Khoud Dam, Oman, revealed an unusual sedimentation pattern which evolved into an intricate composition of silt blocks surrounded by vertical cracks and horizontal layers filled with a “proppant” sand. The discovered soil morphology reflects the complex topology of water motion (infiltration–seepage–evaporation) through the sand-filled cracks/layers and blocks during both the rare flood events and ensuing periods of ponding, and the long, intervening dry periods. These naturally formed soils demonstrate an ability to preserve a large quantity of water inside the silty blocks at depths of 0.5 to 1.5 m, despite the high temperature and dryness of the topsoil. The hydrological optimality and “smartness” of these soils is attributed to the unique block-crack system. Natural, lush vegetation was found in adjacent zones of the reservoir bed, and acted as a footprint of the shallow “fractured perched aquifer”. Planted “ivy” (Convolvulaceae) in the vertical face of one pedon showed intensive growth without irrigation. Soil moisture content data confirmed the hydrological immobility of water in the blocks if not depleted by transpiration. The novel phenomena reported unveil the possible alteration of soil heterogeneity for optimization of the soil–water system in arid zone soils.

Editor D. Koutsoyiannis; Associate editor F.F. Hattermann

Citation Al-Ismaily, S.S., Al-Maktoumi, A.K., Kacimov, A.R., Al-Saqri, S.M., Al-Busaidi, H.A., and Al-Haddabi, M.H., 2013. Morphed block-crack preferential sedimentation in a reservoir bed: a smart design and evolution in nature. Hydrological Sciences Journal, 58 (8), 1779–1788.     

The influence of Caragana korshinskii shrub on soil and hydrological properties in a revegetation-stabilized desert ecosystem

Abstract

Shrub-induced spatial heterogeneity of soil and hydrological properties are common in arid and semi-arid ecosystems. To examine the influence of shrubs on spatial patterns of soil physical, chemical and hydrological properties, the typical sand-fixation species, Caragana korshinskii, was studied in the Shapotou area of the Tengger Desert, China. Miniature cylinder infiltrometers were used to quantify the spatial variations of infiltration rate in the soils, and were installed at 20-cm intervals around the shrubs. Meanwhile, soil samples were collected at 0–5 cm depth every 10 cm to analyse their physical and chemical properties and soil moisture content. The results indicate that the various measured parameters showed a gradational change from sub-canopy to open space. The establishment of shrubs formed obvious “fertile islands” where more soil nutrients collected. The total nitrogen (TN), soil organic matter (SOM), electrical conductivity (EC) and surface soil moisture content decreased gradually from around shrub stems to the interspace. The sand content around shrub stems was significantly higher (p < 0.05), and decreased gradually from the centre towards the outside microsites. The silt and clay contents showed opposite variability characteristics. The variation of soil bulk density was less within 140 cm distance from the stem, and no abrupt change was found at the shrub’s drip line. No significant tendency was found for the soil pH values. The steady infiltration rates declined with increasing stem distance and then tended to be stable, and no abrupt change occurred at the position of the overhead canopy margin. The increase of infiltration rate was rapid nearer to the stem; the variability trend can be fitted by a log-log (power function) model. This study indicated the gradational change in soil and hydrological properties, which was not consistent with the binary division of shrubs into “canopy” and “interspace” zones.

Editor Z.W. Kundzewicz     

Liquefaction of silts: the Adapazari criteria

Ground failures in the form of liquefaction, loss of bearing capacity and soil softening have been observed during the 1999 Marmara (Turkey) earthquake. Research to understand the failure phenomena has been carried out since the earthquakes. This paper attempts to provide explanations to the liquefaction failure of silts in seismic conditions. Findings from a large amount of data collected in the city of Adapazari on the physical and mechanical properties of soils is presented. A geomorphological map of the city has shown that there are surprising horizontal and vertical variations of the facies due to the activity of rivers in the past. Cases of liquefaction appear to have concentrated in former backswamp areas where silts and sandy silts were deposited by crevasse splays. Properties of the soils in zones of liquefaction and non-liquefaction have been determined down to a reasonable depth by measuring the average size, clay content and liquidity index as well as cone penetration resistances with porewater pressures to discover that there is significant discrepancy among those profiles susceptible to liquefaction and non-liquefying deposits. A set of “Adapazari Criteria” is proposed which is intended to improve over the “Chinese Criteria” and is simple enough to be universally applicable. This classification is similar to the existing criteria but emphasizes more on the clay content in addition to measuring the liquid limit and the liquidity index as well as the average size.     

MESURE DE LA PERMÉABILITÉ DES SOLS EN PLACE: UN NOUVEL APPAREIL POUR LA MÉTHODE MUNTZ UNE EXTENSION DE LA MÉTHODE PORCHET AUX SOLS HÉTÉROGÈNES

ABSTRACT

The authors make a general survey of the methods of measuring soil permeability in the field, looking at them from the point of view of the hydrologist.

In the first part they describe new equipment built by the ORSTOM Centre in Lome, Togo, to apply the Muntz method of vertical infiltration under constant head, equipment using a double cylinder and semi-automatic regulation of the discharge.

In the second part, the authors make a mathematical analysis of the Porchet method of horizontal infiltration in a borehole full of water: they show that results may be interpreted even with heterogeneous (multi-layered) soils and that the calculations are simplified using a computer programme.

In conclusion they deal with some results of measurements obtained by the two methods, and examine their accuracy and their correlation with the physical and hydrodynamics characteristics of the soils.     

ORGANISATION DE L'E. G. I. G. Avril 1957

Abstract

The “non-Darcian” behaviour of water in soils is discussed considering the coupling that can arise between the different transport phenomena. In general, it is shown that there is a reason to expect that a difference in osmotic pressure will arise when solutions flow through porous media, and that due to that difference the resulting flow be related in a non-linear fashion to the applied pressure difference. Several numerical examples are presented illustrating this behaviour.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

Some effects of the air-water interface on gravity waves

Abstract

New rates of decay are presented for temporally-attenuated gravity waves in deep water, allowance being made for the energy dissipated in the Stokes interfacial boundary layer in the air. This decay-rate, involving air drag, may then be used to deduce a new “free-surface” boundary condition for the problem of the mass transport velocity due to progressive waves; for shallow-water waves, two specific velocity profiles are calculated, and indicate large differences in comparison with the corresponding profiles of Longuet-Higgins (1953) for a vacuum-water interface.     

Interpretation of the behaviour of agrowell systems in Sri Lanka using radial flow models

Abstract

This paper describes a study of the use of large diameter “agrowells” used for irrigation in Sri Lanka. Detailed field information for a single location is presented and a numerical radial flow model was used to interpret the field response. Initially the model was used to represent a pumping test in which the pumping and subsequent recovery were monitored. It was then used to investigate the operation of the well during a dry season. Satisfactory agreement could only be obtained when account was taken of the withdrawal of water from below the water table by trees. From the simulation it was possible to show that skilful farming practice could be used in steadily exploiting the resources for irrigation over the whole of the dry season.     

A combined methodology for estimating the potential natural aquifer recharge in an arid environment

ABSTRACT

An innovative methodology that combines an indirect physiography-based method for determining the runoff coefficient at a sub-basin scale and a water balance model applied on a daily time scale was developed to calculate the natural groundwater recharge in three watersheds within the Oum Zessar arid area, Tunisia. The effective infiltration was calculated as part of the water surplus by considering the average available water content (AWC) of soil and an average runoff coefficient for each sub-basin. The model indicates that the sub-basins covered mainly by the “artificial” soils of tabias and jessour, characterized by average AWC values greater than 150 mm, did not contribute to natural groundwater recharge over the 10-year period (2003–2012) considered. The estimated volume for the Triassic aquifer amounted to about 4.5 hm³ year^?1, which is consistent with previous studies. For the Jurassic and Cretaceous aquifers, the estimated volumes amounted to about 200 dm³ year^?1.     

Critical analysis of thermal inertia approaches for surface soil water content retrieval

Abstract

The “thermal inertia” method to retrieve surface soil water content maps on bare or sparsely-vegetated soils is analysed. The study area is a small experimental watershed, where optical and thermal images (in day and night time) and in situ data were simultaneously acquired. The sensitivity of thermal inertia to the phase difference between incoming radiation and soil temperature is demonstrated. Thus, to obtain an accurate value of the phase difference, the temporal distance between thermographs using a three-temperature approach is evaluated. We highlight when a cosine correction of the temperature needs to be applied, depending on whether the thermal inertia formulation includes two generic acquisition times, or not. Finally, the deviation in soil water content retrieval is quantifies for given values of each parameter by performing a sensitivity analysis on the basic parameters of the thermal inertia method that are usually affected by calibration errors.

Citation Maltese, A., Bates, P.D., Capodici, F., Cannarozzo, M., Ciraolo, G., and La Loggia, G., 2013. Critical analysis of thermal inertia approaches for surface soil water content retrieval. Hydrological Sciences Journal, 58 (5), 1144–1161.

Editor D. Koutsoyiannis; Associate editor D. Hughes     

Modelling water infiltration in cracked paddy field soil

Periodic paddy field flooding is a major source of groundwater recharge. Many paddy fields thus are used as groundwater recharge ponds after harvesting the first crop of the summer. Following rice harvesting, paddy field surfaces may crack into fissures as a result of drainage and exposure to sunlight. Field observation indicates that applying precipitation to the paddy field can increase the rate of infiltration. To quantitatively evaluate the amount of infiltration in a cracked paddy field, this study sets up a simple soil crack model to simulate the field infiltration process. A three‐dimensional groundwater model FEMWATER is adopted to simulate water movement in the paddy field subjected to various crack conditions. Using the field and laboratory data of irrigation water requirements, soil physical properties, hydraulic conductivities and soil profiles obtained from Ten‐Chung, FEMWATER simulates the water movement in the dry cracked paddy. Simulation results show that if the cracks develop extensively and penetrate the ploughed soil, the infiltration rate may increase significantly. The infiltration fluxes of crack with depths of 80, 60 and 27·5 cm are 18·77, 14·50 and 8·06 times higher than that of 20 cm, respectively. The simulation results of cracks with 80 cm depth correlated closely with field observations. The results of the study elucidate the processes of unsaturated water movement in a dry cracked paddy field. Copyright © 2004 John Wiley & Sons, Ltd.     

Real-Time Geoelectric Monitoring of Seepage into Sand and Clay Layer

The migration of groundwater in rock and soil can appear as abnormalities in geoelectric fields. It is therefore important to study the characteristics of seepage in porous media by measuring the geoelectric field signatures. In this study, a physical model with layers of sand and clay was constructed and an electrical resistivity meter was used to examine the changes in the geoelectric field parameters during the infiltration process. The results show that the infiltration could be detected based on the geoelectric signatures including temporal changes through the spontaneous potential, excitation currents, and apparent resistivity. Specifically, the spontaneous potential was reduced by 100 to 200 mV when the water reaches an electrode. During the second water injection in the experiment, the measured spontaneous potential of all the electrodes recovered to the previous extreme values that range from −200 to −550 mV, thus indicating a “memory” effect. With stepwise changes in the excitation current, it was possible to determine the seepage velocity in sand and clay layer. The apparent resistivity is reduced to less than 400 Ωm when the infiltration reaches the electrodes. These results indicate the potential for real-time monitoring of water flow.     

Analysing tension infiltrometer data from sloped surface using two‐dimensional approximation

Tension infiltrometers allow water to infiltrate into the soils at various specified pressure heads. The resulting infiltration rates can then be analysed for soil hydraulic properties by either analytical or inverse numerical methods. Tension infiltrometers however are primarily designed to be deployed on horizontal land surfaces, and their applications have been studied widely using an inverse numerical tool HYDRUS‐2D. However, natural landscapes are often nonhorizontal, and infiltration through tension infiltrometers on sloped surfaces is no longer an axisymmetrical two‐dimensional (2D) process but a fully three‐dimensional (3D) one. In addition, minimal research has examined the effect of simplifying the 3D problem to a 2D one on the hydraulic conductivity estimated using tension infiltrometer data from different land slopes of various soil types. Therefore, in this study, tension infiltrometer data on different slopes have been obtained from a catchment located at National University of Singapore. In addition, tension infiltrometer data of six soil types on different slopes and with different initial water content were simulated using HYDRUS‐3D. Combining field measurements, forward and inverse modelling, the influence of applying a 2D approximation on hydraulic property estimations using tension infiltrometer data was examined. The results show that the estimations for soils with high infiltration rates are more sensitive to application of the 2D approximation. The maximum allowable slopes for employing 2D approximation on clayey and sandy soils are 25° and less than 3°, respectively. Furthermore, the maximum allowable slope decreases with decreasing initial water content. Copyright © 2012 John Wiley & Sons, Ltd.     

REPORT ON THE ACTIVITIES OF THE INTERNATIONAL COMMISSION ON SNOW AND ICE FOR THE PERIOD 1971–1975

ABSTRACT

An infiltration model based on the Green-Ampt assumptions is developed for a class of non-uniform soils in which saturated hydraulic conductivity decreases as an exponential function of depth, and the storage-suction factor C = ΔθΔΨ is a constant. An analysis of measured hydraulic conductivity and porosity data suggests that the model may be useful in some soils where the changes of hydraulic characteristics with depth are of this form. A method for using the model with time variable rainfall rates is given. The model is applied to simulate the experimental results of Childs &; Bybordi (1969) who measured infiltration into layered sand profiles. Agreement was satisfactory even for this case.     

Quantifying potential recharge in mantled sinkholes using ERT

Potential recharge through thick soils in mantled sinkholes was quantified using differential electrical resistivity tomography (ERT). Conversion of time series two-dimensional (2D) ERT profiles into 2D volumetric water content profiles using a numerically optimized form of Archie's law allowed us to monitor temporal changes in water content in soil profiles up to 9 m in depth. Combining Penman-Monteith daily potential evapotranspiration (PET) and daily precipitation data with potential recharge calculations for three sinkhole transects indicates that potential recharge occurred only during brief intervals over the study period and ranged from 19% to 31% of cumulative precipitation. Spatial analysis of ERT-derived water content showed that infiltration occurred both on sinkhole flanks and in sinkhole bottoms. Results also demonstrate that mantled sinkholes can act as regions of both rapid and slow recharge. Rapid recharge is likely the result of flow through macropores (such as root casts and thin gravel layers), while slow recharge is the result of unsaturated flow through fine-grained sediments. In addition to developing a new method for quantifying potential recharge at the field scale in unsaturated conditions, we show that mantled sinkholes are an important component of storage in a karst system.     

Cyclonic-anticyclonic asymmetry and a new soliton concept for rossby vortices in the laboratory,oceans and the atmospheres of giant planets

Abstract

It is demonstrated in laboratory experiments with rotating shallow water that large scale Rossby vortices, greater than the Rossby-Obukhov radius in size, have dispersive and non-linear properties that are fundamentally different for the two possible polarities. We call this “cyclonic-anticyclonic asymmetry”. This asymmetry manifests itself in the following way: first, anticylones, unlike cyclones, do not undergo the dispersive spreading inherent in a linear wave packet. and therefore, having a considerably longer natural lifetime, are obvious candidates for Rossby solitons; second, dipolar vortices are, because of the comparatively rapid decay of a cyclone, transformed into anticyclonic solitons; third, anticyclones are much more readily generated by zonal flows of the type existing in planetary atmospheres. The evident dominance of anticyclones amongst the long-lived vortices in the atmospheres of giant planets strongly suggests that the cyclonic-anticyclonic symmetry plays a decisive role in the atmospheric cyclogenesis of large planets.

According to our concept, the Rossby soliton is a “real” vortex; unlike a wave, it retains some fluid particles within it throughout its lifetime. Two similar solitons can merge by mutual collisions. This picture of a “vortical” soliton differs in an essential way from the earlier idea due to Maxworthy and Redekopp (1976) of purely “wave-like” Rossby solitons that can freely pass through one another.

Laboratory experiments were performed by us to simulate the new Rossby solition, with special reference to naturally-occurring vortices of the same general type as Jupiter's great red spot. The experimental data presented contradict the idea of “pure wave solitons” but confirm our concept of “vortical solutions”.     

The effects of coal gangue and fly ash on the hydraulic properties and water content distribution in reconstructed soil profiles of coal‐mined land with a high groundwater table

Soil water systems have been severely degraded in coal‐mined and subsiding land, where a shallow groundwater table is also present. The present paper discussed the effects of fly ash (FA) and coal gangue (CG) as filling materials on the hydraulic properties and water content distribution in a profile for the purpose of rehabilitating subsided lands. The saturated water content, water characteristic curve, and water diffusivity of local soil, FA, CG, and a mixture of FA and CG (“mixed filling”) were characterized. A column experiment was conducted to investigate the changes of water content in profiles reconstructed from the combination of soil and filling materials, including soil only, FA, CG, and a mixture of FA and CG, which were used to fill the lower part of the reconstructed profile. The mixture of FA and CG was found to possess similar hydraulic properties to those of the soil, particularly high water‐holding capacity and permeability. Moreover, the volumetric water contents in the whole profile containing the mixture of FA and CG were consistent with those of the profile reconstructed with soil only. As a result, it is recommended to adopt the mixture of FA and CG for reconstructing the lower profile of the land to alleviate or rehabilitate subsided land in coal mines.     

A mathematical model of polythermal glaciers and ice sheets

Abstract

A mathematical model of the flow and temperature distribution of polythermal glaciers or ice sheets is deduced. Cold ice is treated as a non-linear viscous heat conducting fluid, while temperate ice is regarded as a binary mixture of ice and water. The simplest mixture concept with two balance laws of mass but only one balance law of momentum and energy is proposed. The field equations for the ice and water content and the boundary conditions which must hold at the free surface, at the ice-water interface, at the cold-temperature transition surface and at the rock-bed are deduced. In particular it is shown that an earlier formulation of polythermal ice due to Fowler and Larson (1978) is inconsistent. No boundary value problems are solved as the emphasis is on the physical motivation and justification of the principles.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.