Simulations of water movement and solute transport through different soil texture configurations under negative‐pressure irrigation

This study examined the effects of different soil texture configurations on water movement and solute transport to provide a reliable scientific basis for the application of negative‐pressure irrigation (NPI) technology. HYDRUS‐2D was used to analyse water movement and solute transport under NPI. The main results are as follows: (a) HYDRUS‐2D can be used to simulate water movement and solute transport under NPI, as there was good agreement between the simulated and measured values for water contents, NaCl concentrations, cumulative water infiltration, and wetting distances in the horizontal and vertical directions; the Nash–Sutcliffe efficiency coefficients were in the range of 0.94–0.97. (b) Layered soils have obvious effects on water movement under NPI. With the emitter position in the loam layer, when a coarse texture of loamy sand was present below the loam layer (namely, L‐LS), irrigation water accumulated in the topsoil, and this led to an increase in evaporation compared with the homogeneous loam profile. However, fine texture silty loam or silty clay loam layers beneath the loam layer (namely, L‐SiL or L‐SiCL, respectively) was more conducive to water infiltration into the lower layer, and this increased the amount of water infiltration and simultaneously reduced the surface evaporation effectively. (c) Layered soils have obvious effects on solute transport under NPI, and salt accumulation will readily occur in the clay‐rich soil layer at the interface. The maximum soil salt accumulation of L‐LS occurred above the soil interface between the two soil layers with a value of 21.80 g/kg; however, for L‐SiCL and L‐SiL, the maximum salt accumulation occurred below the soil interface between the two soil layers, with values of 23.80 g/kg and 20.08 g/kg, respectively. (d) Interlayered soils showed remarkable changes in the water infiltration characteristics and salt‐leaching intensities under NPI, and the properties for the soil profile with a silty loam interlayer were better than those for the soil profile with a silty clay loam interlayer. The soil profile with a loamy sand interlayer had the lowest amount of water infiltration, which resulted in reductions of the salt‐leaching intensities. Thus, NPI is clearly not suitable for loamy sand soil. Overall, the results demonstrated that soil texture configurations affected water movement and solute transport under NPI. Therefore, careful consideration should be given to the use of NPI to achieve target soil water and solution conditions and reduce water loss.     

Investigation of liquefaction and pore water pressure development in layered sands

This paper presents the results of shaking table model tests which were carried out to investigate the pore water pressure generation and related liquefaction mechanism in layered sand deposits. The experiments were performed on uniform sand columns, silt interlayered sand columns and two layered sand columns deposited at various relative densities and subjected to different input excitations. During the experiments excess pore water pressures were measured by pore pressure transducers installed at three different depths and, surface settlements and thickness of water film developed under less permeable inclusions were measured by a digital camera. The experimental results are discussed and compared to demonstrate the effects of relative density, input acceleration and presence of a silt seam on the generation of excess pore water pressure in sand deposits subjected to dynamic loading. The results showed that the presence of a less permeable silt interlayer within the sand deposit and existence of a loose sand layer underlying dense sand deposits can have significant effect on the pore water pressure generation mechanism.     

Field Test of a Hybrid Finite‐Difference and Analytic Element Regional Model

Regional finite‐difference models often have cell sizes that are too large to sufficiently model well‐stream interactions. Here, a steady‐state hybrid model is applied whereby the upper layer or layers of a coarse MODFLOW model are replaced by the analytic element model GFLOW, which represents surface waters and wells as line and point sinks. The two models are coupled by transferring cell‐by‐cell leakage obtained from the original MODFLOW model to the bottom of the GFLOW model. A real‐world test of the hybrid model approach is applied on a subdomain of an existing model of the Lake Michigan Basin. The original (coarse) MODFLOW model consists of six layers, the top four of which are aggregated into GFLOW as a single layer, while the bottom two layers remain part of MODFLOW in the hybrid model. The hybrid model and a refined “benchmark” MODFLOW model simulate similar baseflows. The hybrid and benchmark models also simulate similar baseflow reductions due to nearby pumping when the well is located within the layers represented by GFLOW. However, the benchmark model requires refinement of the model grid in the local area of interest, while the hybrid approach uses a gridless top layer and is thus unaffected by grid discretization errors. The hybrid approach is well suited to facilitate cost‐effective retrofitting of existing coarse grid MODFLOW models commonly used for regional studies because it leverages the strengths of both finite‐difference and analytic element methods for predictions in mildly heterogeneous systems that can be simulated with steady‐state conditions.     

A study on water film in saturated sand

Water film can serve as a sliding surface and cause landslides on gentle slopes. The development of "water film" in saturated sand is analyzed numerically and theoretically based on a quasi-three-phase model. It is shown that stable water films initiate and grow if the choking state (where the fluid velocity decreases to near zero) remains steady in a liquefied sand column. Discontinuity can occur in pore water velocity, grain velocity and pore pressure after the initiation of a water film. However, the discontinuity and water film can disappear once the choking state is changed. The key to the formation of water film is the choking in the sand column caused by eroded fine grains.     

煤层上覆地层含水不均匀性电法探测的可能性

煤田上覆岩层中含水直接对煤层的开采构成威胁，本文先从理论上阐述了电阻率勘探可以探测煤层上覆地层含水不均匀性的可能性，然后用简单的断层型薄层含水模型的CSAMT正演模拟结果展示出：地下900m深处，尽管只有20m厚的的含水岩层，但在电阻率剖面上该薄层的反映具有非常明显的特征，从理论上和数值模拟实例说明，一定深度内一定厚度的薄层含水性的电阻率法探测是可能的。     

循环荷载下液化对土层水平往返变形的影响

采用多工况振动台实验研究液化对土层水平往返变形的影响.以干砂实验为参照,分析孔压增长与土层加速度和土层往返变形之间的关系.结果表明:液化将引起土表加速度显著降低,减小惯性力传递,但同时会引起土层往返剪应变明显增大.对往返变形而言,液化土层往返剪应变就可达到1%～5%的大变形状态,且液化土层往返剪应变沿深度呈下大上小分布.土层中孔压比0.4～0.8是往返变形出现放大的敏感段,在孔压比0.8左右而不是在1.0达到最大.作为其结果,土层液化将对刚性上部结构振动起减震作用,但同时增大的往返剪应变也易导致基础和地下结构破坏,特别是对液化层与下部非液化层交界处的构件更敏感.     

Removal of Perchloroethylene from a Layered Soil System by Steam Flushing

Steam flushing experiments were concluded in a two-dimensional chamber containing two layers of F75 silica sand separated by a layer of finer F10 silica sand. Perchloroelhylene (PCE), which had spilled into the chamber under water-saturated conditions, formed a pool on the F110 sand layer. Steam was injected above the F110 sand layer. Temperatures, moisture content, and PCE concentrations in the chamber were monitored. Samples, taken from the various locations in the sand chamber, indicated that complete removal of PCE from the steam zone was achieved, with an 84% overall recovery. Some downward displacement of PCE-contaminated water through the F110 sand layer was observed and a small amount of gravity override occurred. Channeling of steam was minimal. The experiment indicates that steam flushing may be used successfully for removal of PCE from relatively homogeneous soils. Issues of gravity override and downward mobilization of contaminants must be considered in applying steam flushing at the field scale.     

Physical modeling of lateral spreading induced by inclined sandy foundation in the state of zero effective stress

The state of zero effective stress is a situation at which the effective stress of saturated sand decreases to zero in the process of liquefaction. In the state of zero effective stress, sand particles suspend in water and the foundation is vulnerable to much large lateral deformation. The state of zero effective stress can be achieved through dynamic loading tests, but the obtained state is difficult to sustain a steady situation. To simulate the suspended sand in water under fully liquefied condition, plastic sand, characterized by small specific gravity, is used instead of quartz sand to build an inclined foundation. Salt water with similar density is used to pass in slowly near bottom of the foundation. As observed in tests, the plastic sand is able to suspend in sodium chloride solution (salt water) of a specific density and thus this model can be used to simulate the lateral spreading of a foundation under zero effective stress state. Lateral deformation occurs within a certain depth beneath the ground and the magnitude increases from the bottom up, showing nonlinear behaviors. This paper presents a physical modeling approach for achieving the state of zero effective stress under static laboratory condition.     

基于UPFs的珊瑚岛礁场地地震反应分析方法研究

合理且高效地模拟珊瑚砂非线性动力特性、远场无限地基辐射阻尼以及海域岛礁动水压力的影响是进行珊瑚岛礁抗震安全分析的关键技术问题.以通用有限元软件ANSYS为研究工具,基于UPFs二次开发灵活性的特点,建立了适用于珊瑚砂地基条件下的岛礁场地地震反应分析时域计算模型.通过创建一种新的珊瑚砂等价线性单元描述岛礁场地的非线性动力...     

Hydrological cycle and water balance estimates for the megadune–lake region of the Badain Jaran Desert,China

The hydrology and water balance of megadunes and lakes have been investigated in the Badain Jaran Desert of China. Field observations and analyses of sand layer water content, field capacity, secondary salt content, and grain size reveal 3 types of important natural phenomenon: (a) vegetation bands on the leeward slope of the megadunes reflect the hydrological regime within the sandy vadose zone; (b) seepage, wet sand deposits, and secondary salt deposits indicate the pattern of water movement within the sandy vadose zone; (c) zones of groundwater seeps and descending springs around the lakes reflect the influence of the local topography on the hydrological regime of the megadunes. The seepage exposed on the sloping surface of the megadunes and gravity water contained within the sand layer confirm the occurrence of preferential flow within the vadose zone of the megadunes. Alternating layers of coarse and fine sand create the conditions for the formation of preferential flows. The preferential flows promote movement of water within the sand layer water that leads to deep penetration of water within the megadunes and ultimately to the recharging of groundwater and lake water. Our results indicate that a positive water balance promotes recharge of the megadunes, which depends on the high permeability of the megadune material, the shallow depth of the surface sand layer affected by evaporation, the occurrence of rainfall events exceeding 15 mm, and the sparse vegetation cover. Water balance estimates indicate that the annual water storage of the megadunes is about 7.5 mm, accounting for only 8% of annual precipitation; however, the shallow groundwater per unit area under the megadunes receives only 3.6% of annual precipitation, but it is still able to maintain a dynamic balance of the lake water. From a water budget perspective, the annual water storage in the megadunes is sufficient to serve as a recharge source for lake water, thereby enabling the long‐term persistence of the lakes. Overall, our findings demonstrate that precipitation is a significant component of the hydrological cycle in arid deserts.     

MODEL STUDIES ON SOME ASPECTS OF RESISTIVITY AND MEMBRANE POLARIZATION BEHAVIOUR OVER A LAYERED EARTH*

Electrolytic model tank experiments to study resistivity and time domain induced polarization (IP) response over layered earth models were initiated primarily to facilitate the understanding of field results. Alternate layers of clay and sand (or clay-coated sand) with, in some cases, a surficial layer of water were assembled in the tank and resistivity and IP measurements made for a range of electrode spacings using the Wenner configuration. Graphite and silver-silver chloride electrodes were used as current and potential electrodes respectively. Clay-coated (3% by weight) sand was found to generate stronger polarization than either clay or sand alone. Apparent chargeability m_a was observed to be positive for a nonpolarizable surface layer. For a polarizable surface layer, the sign of IP was controlled by the polarizability, the thickness of the second layer, and the spacing of the electrode spreads. The apparent chargeability m_a can theoretically change sign from positive to negative and vice versa with a gradual increase in electrode spacing, and such negative IP effects were obtained in a few observations. A simultaneous decrease in IP and an increase in resistivity, which is a qualitative diagnostic feature for the occurrence of clean freshwater sand aquifers, could also be generated in the model tank experiment. Combined resistivity and IP soundings were carried out near Fredericton Junction and Tracy, New Brunswick, Canada. Field curves are presented along with the model curves for qualitative comparison and understanding of IP behaviour over a layered earth. Twenty-five out of twenty-seven soundings show only positive apparent chargeabilities, whereas two show chargeability sign changes (positive/negative/positive). The model study gives reason to believe that surface soils and Quaternary gravel boulder deposits near Fredericton Junction are relatively non-polarizable. As an auxiliary experiment, sand and clay were taken in different proportions by weight and mixed thoroughly with water in a cement mixer. The mixtures were then compressed with a suitable die and plunger under 3.6 Pa pressure to prepare cylindrical samples of height 18 cm and diameter 15.5 cm. IP measurements were done on the flat faces using the Wenner configuration with a= 2 cm. Chargeability was found to be negative for 100 and 90% clay mixtures. It reached a positive maximum for an 80% clay-20% sand mixture and then decreased gradually with increasing sand and decreasing clay content.     

Centrifuge modeling of the seismic responses of sand deposits with an intra-silt layer

Three dynamic centrifuge model tests were conducted at an acceleration of 80g to simulate the seismic responses of level sand deposits: an intra-silt layer was embedded in two of these sand deposits at different depths. The effects of a low-permeability intra-silt layer on the build-up and dissipation of excess pore-water pressure, surface settlement, and the related liquefaction mechanism were investigated. An intra-silt layer modifies the seismic response of the sand deposit, reduces the extent of liquefaction, and thus decreases surface settlement. The depth of the intra-silt layer is one of the factors influencing the seismic responses of the sand deposits. The magnitude of the surface settlement is proportional to the degree of liquefaction in the sand deposit. The high positive hydraulic gradients appearing in both the intra-silt layer and in the sand deposit lying on the intra-silt layer can break a thinner or weaker top layer and result in sand boiling. Our visual animation of the ratio of the excess pore-water pressure and the lateral displacement revealed that the liquefaction front travels upward during shaking and the solidification front travels upward after shaking.     

Study of the Santa Catarina aquifer system (Mexico Basin) using magnetotelluric soundings

A tensor magnetotelluric test survey was carried out in the region of Santa Catarina, located in the Chalco sub-basin of the Mexico Basin. The objective was to define the stratification at depth with an emphasis on the geometry of the main aquifer of that region which is partially known from DC resistivity soundings and drilling. High-quality magnetotelluric soundings could be recorded in the immediate vicinity of large urban zones because the sub-surface is very conductive. Interpretation shows that the solid bedrock is located at a depth of at least 800 m to the south and 1300 m to the north; it could, however, be much deeper. Using complementary DC resistivity sounding and well-logging data, three main layers have been defined overlying the bedrock. These layers are, from surface to bottom, an unsaturated zone of sand, volcanic ash and clay about 10 m thick, followed by a very conductive (1.5 ohm·m) 200 m thick layer of sand and ash with intercalated clay, saturated with highly mineralized water, and finally a zone with resistivity increasing gradually to 60 ohm·m. The investigated deep aquifer constitutes most of this third layer. It consists of a sequence of sand, gravel, pyroclastites and mainly fractured basalts. MT resistivity soundings and magnetic transfer functions also indicate that a shallow resistive structure is dipping, from the northwest, into the lacustrine deposits of the basin. This geologic feature is likely to be highly permeable fractured basaltic flows, which provide a channel by which water contaminated by the Santa Catarina landfill may leak into the basin.     

Inelastic response analysis of reinforced concrete structures using modified force analogy method

A modified force analogy method (MFAM) is developed to simulate the nonlinear inelastic response of reinforced concrete (RC) structures. Beam–column elements with three different plastic mechanisms are utilized to simulate inelastic response caused by moment and shear force. A multi‐linear hysteretic model is implemented to simulate the nonlinear inelastic response of RC member. The P‐Δ effect of the structure is also addressed in MFAM. Static and dynamic inelastic response of structure, damage condition and failure type for structural element, structural limit state and collapse time can also be simulated using MFAM. Compared with the general algorithm, the MFAM provides less computational time especially in the case of large structural system. It is also easier to be written as computer program. Three test data groups, which include cyclic loading test data of a non‐ductile RC bridge column, a two‐storey RC frame, and dynamic collapse test data of a non‐ductile RC portal frame, are selected to confirm the effectiveness of applying MFAM to simulate the inelastic behaviour of structures. Copyright © 2007 John Wiley & Sons, Ltd.     

Transboundary flow modeling: the Zohor depression of Austria and the Slovak Republic

Estimation of available ground water is a basic aspect of ground water management. Mathematical modeling is one of the methods that can be effectively used to obtain such estimates. A numerical model was used to calculate available ground water in the Zohor depression-an aquifer transcending national boundaries between the Slovak Republic and Austria. The aquifer, formed by Quaternary sediments overlying Neogene sequences, is composed of various clays interbedded with layers of sand, gravel, sandstones, and conglomerates. The AQUA computer model package was used to simulate flow in the aquifer. For model compilation, the following data were used: (1) effective precipitation; (2) surface water levels in surface water gauging profiles; and (3) withdrawal amounts. Hydraulic parameters of the aquifer were estimated based on information from 86 wells located in the area. The model was verified using data on ground water levels from a monitoring network. The simulation of the aquifer system permitted the estimation of the available ground water in the study area, showing that an additional 587 L/s can be abstracted. Ground water inflows to the Morava River, which flows through the region, range from 745 to 3100 L/s.     

Centrifuge modelling of reverse fault–foundation interaction

The propagation of reverse faults through soil to the ground surface has been observed to cause damage to surface infrastructure. However, the interaction between a fault propagating through a sand layer and a shallow foundation can be beneficial for heavily loaded foundations by causing deviation of the fault away from the foundation. This was studied in a series of centrifuge model tests in which reverse faults of dip angle 60° (at bedrock level) were initiated through a sand layer, close to shallow foundations. The tests revealed subtle interaction between the fault and the shallow foundation so that the foundation and soil response depend on the foundation loading, position, breadth and flexibility. Heavily loaded rigid foundations appeared best able to deviate fault rupture away from the foundation but this deviation could be associated with significant foundation rotations. However, a lightly loaded foundation was unable to deviate a reverse fault and the fault emerged beneath the foundation. This led to gapping beneath the foundation as well as significant rotations and may cause severe structural distress. As well as providing insight into the mechanisms of behaviour, the data from the tests is used to validate finite element analyses in a separate article.     

饱和弱化与液化土层的水平极限抗力

通过给饱和砂土层施加反压,模拟地震荷载作用下具有残余孔压的饱和弱化、液化土层。选择粉质细砂与细砂,进行了18组水平荷载作用下桩与饱和弱化、液化土层相互作用的模型试验,研究了饱和弱化、液化土层水平极限抗力随土层残余孔压增加的变化规律。结果表明,随土层中残余孔压增加,水平极限抗力逐渐降低,土层液化后的水平极限抗力大约降低80%~90%。通过定义饱和弱化、液化土层的强度,定量分析了饱和弱化、液化砂土的强度参数与水平极限抗力之间关系。又通过引入土层的残余孔压比折减系数,建立了确定饱和弱化、液化土层等效强度的关系式,进而提出了一种按等效强度确定饱和弱化、液化土层水平极限抗力的方法。     

Mathematical modelling of sea water intrusion,Nauru Island

A two-dimensional finite element model for density dependent groundwater flow was calibrated to simulate sea water intrusion in Nauru Island in the Central Pacific Ocean. Nauru Island occupies an area of 22 km² and supports a population of 8500. The island has been mined for its phosphate deposits and current reserves indicate that the mine has about eight years life remaining. The water supply of the island is about one third dependent on imported water which is also used as ballast on the phosphate ships. Imported water will not be available in the future, and a hydrogeological investigation shows that the island is underlain by a fresh water layer, less than 5 m thick. The freshwater layer overlies a thick transition zone of brackish water which in turn overlies sea water. Simulation of several management options shows that it is possible to substitute current importation of fresh water by careful extraction from the groundwater resources of the island.     

Analysis of saturated dense sand-structure system and comparison with results from shaking table test

An effective stress method is presented for the analysis of liquefaction of ground including soil-structure interaction, based on an explicit-implicit finite element method. A simple constitutive model is developed to be incorporated in the effective stress method. The constitutive model consists of the Ramberg-Osgood model extended to two-dimensional problems and a new dilatancy model. The effectiveness of the constitutive model is examined with results of a simple shear test. Besides, the effective stress method is verified by comparing its numerical results with results of a shaking table test. It is found that the present method can simulate well the response of a saturated dense sand-structure system. The difference of the response computed by the effective stress method and the total stress method is discussed. It is found that the total stress method can simulate the response of the saturated sand within an accumulating excess pore water pressure of less than 70 per cent of the initial overburden stress.     

Formulation and application of the unsaturated/saturated catchment models SUSCAT and WEC‐C

This paper describes the formulation and application of a coupled unsaturated/saturated model framework developed to investigate the impact of mining on catchment water yield and groundwater dynamics. The model conceptualization was implemented in both a finite‐element (SUSCAT) and finite‐difference (WEC‐C) solution scheme and found to give similar results. The model framework simulates a coupled surface‐water and groundwater system in which a physically based solution scheme was used to simulate one‐dimensional movement through the unsaturated zone, and a distributed model was used to simulate two‐dimensional saturated groundwater flow. Each soil column comprises a series of layers, each layer being connected to adjacent cells. Subsurface lateral flow is considered when any cell within a layer develops a saturated thickness. Simulation results presented are based on a catchment in the Darling Range, Western Australia that was progressively mined and subsequently rehabilitated. The results predicted the groundwater system beneath the mine areas to have a peak rise owing to mining of between 2 and 4 m. Six years after mining, and following vegetation rehabilitation, the groundwater rise had reduced to 1 m above simulated unmined levels. The corresponding streamflow increase as a result of mining was estimated to peak at 21 mm/year and declined to 7·4 mm/year eight years after revegetation of the mined areas. The simulated groundwater response and streamflow results derived from both models were found to be consistent with observed data. Copyright © 2002 John Wiley & Sons, Ltd.