Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Dust supply from river floodplains: the case of the lower Huang He (Yellow River) recorded in a loess–palaeosol sequence from the Mangshan Plateau

The Mangshan Plateau is located on the south bank of the Huang He (Yellow River) just west of the city of Zhengzhou, well outside the Loess Plateau in central China. Mixing models of the grain‐size data indicate that the loess deposits are mixtures of three loess components. Comparison of the mixing model with existing models established for a series of loess–palaeosol sequences from the Loess Plateau indicates that the Mangshan loess has been supplied from a proximal dust source, the Huang He floodplain, during major dust outbreaks. The high accumulation rates, the composition of the loess components, and especially the high proportions of a sandy loess component support this. Owing to the exceptionally high accumulation rates, the Mangshan grain size, magnetic susceptibility and carbonate records provide a high‐resolution archive of environmental and climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

Identification of physical processes inherent in artificial neural network rainfall runoff models

The emergence of artificial neural network (ANN) technology has provided many promising results in the field of hydrology and water resources simulation. However, one of the major criticisms of ANN hydrologic models is that they do not consider/explain the underlying physical processes in a watershed, resulting in them being labelled as black‐box models. This paper discusses a research study conducted in order to examine whether or not the physical processes in a watershed are inherent in a trained ANN rainfall‐runoff model. The investigation is based on analysing definite statistical measures of strength of relationship between the disintegrated hidden neuron responses of an ANN model and its input variables, as well as various deterministic components of a conceptual rainfall‐runoff model. The approach is illustrated by presenting a case study for the Kentucky River watershed. The results suggest that the distributed structure of the ANN is able to capture certain physical behaviour of the rainfall‐runoff process. The results demonstrate that the hidden neurons in the ANN rainfall‐runoff model approximate various components of the hydrologic system, such as infiltration, base flow, and delayed and quick surface flow, etc., and represent the rising limb and different portions of the falling limb of a flow hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

试论城市地貌环境与城市道路系统的关系*——以四川省几个城市为例

本文分析了城市地貌环境对城市道路系统的影响,认为地貌过程不仅影响了城市道路的建设和管理,而且是城市道路系统布局和发展的重要控制因素。以四川省几个城市为例,分析了山地城市、丘陵城市和平原城市道路系统的不同特点。     

Design of stable drainage networks in upland forestry plantations

Erosion of man-made, forestry drainage channels occurs when the plough cuts through the surface peat layer into the underlying erodible mineral soil. A procedure is developed, based on hydraulic considerations, which will allow the drainage engineer to design stable drainage networks in upland forestry plantations. An example design chart is given for an erodible sandy loam type soil.