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.     

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.     

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.     

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.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

A 2D finite volume model for bebris flow and its application to events occurred in the Eastern Pyrenees

FLATModel is a 2D finite volume code that contains several original approaches to improve debris-flow simulation. Firstly, FLATModel incorporates a ＂stop-and-go＂ technique in each cell to allow continuous collapses and remobilizations of the debris-flow mass. Secondly, flow velocity and consequently yield stress is directly associated with the type of rheology to improve boundary accuracy. Thirdly, a simple approach for entrainment is also included in the model to analyse the effect of basal erosion of debris flows. FLATMODEL was tested at several events that occurred in the Eastern Pyrenees and simulation results indicated that the model can represent rather well the different characteristics observed in the field.     

甘青地区新石器-青铜时代居民颅面部形态变异及其在人群形成上的意义

通过对考古遗址出土古代居民遗骸的研究来探寻人群形成、适应、迁徙、融合历史是生物人类学研究的重要内容。甘青地区地处我国西北,位于新疆与中原之间,是人群迁徙与交流的关键地区。以往研究发现该地区新石器-青铜时代居民具有的一些独特的头面部特征,如具有偏狭的中颅型、高颅结合狭颅、中等偏狭的面宽、高而狭的面型、中等的面部扁平度、中眶型、狭鼻型和正颌型。一些学者就来自周边地区人群迁徙或基因交流对甘青地区古代人群形成的影响做过探讨。由于迄今开展的相关研究及掌握的数据资料有限,学术界对甘青地区古代人群内部的变异情况、时代变化及其与周边地区人群的差异程度还缺乏深入的了解,对新石器时代以来甘青地区人群的形成过程以及与东、西方人群的关系还不是很清楚。本文通过对243例甘青地区新石器-青铜时代人群颅骨测量数据的统计和分析,结合相关考古学和古环境学资料,对该地区古代居民颅面部形态特征的总体表现特点、时代变化、内部变异及与周边地区人群的关系进行了探讨,发现：1）从整体颅面部特征上看,甘青地区新石器-青铜时代人群介于典型西方人群与中原人群之间,但与中原人群更为接近;2）甘青地区新石器-青铜时代人群颅骨测量特征呈现一定的时代变化。从新石器时代到青铜时代,颅型趋向圆颅和阔颅化、面型趋狭、鼻型趋狭、眶型趋高;齿槽部的突颌水平变小;3）甘青地区古代人群颅面部测量特征的内部变异程度介于新疆与中原人群之间。从新石器时代到青铜时代,甘青地区人群颅面部测量特征的内部变异呈减小趋势;4）从新石器时代到青铜时代,甘青地区古代人群部分颅面部测量特征呈现出与西方人群接近的趋势。基于本文的分析,认为气候环境、生业方式以及人群之间的基因交流共同造成了甘青地区古代人群颅面部形态特征的变化,甘青地区古代人群形成以受到中原人群影响为主。

     

Runoff quality prediction from small urban catchments using SWMM

The RUNOFF block of EPA's storm water management model (SWMM) was used to simulate the quantity and quality of urban storm water runoff from four relatively small sites (i.e. 5·97–23·56 ha) in South Florida, each with a specific predominant land use (i.e. low density residential, high density residential, highway and commercial). The objectives of the study were to test the applicability of this model in small subtropical urban catchments and provide modellers with a way to select appropriate input parameters to be used in planning studies. A total of 58 storm events, measured by the US Geological Survey (USGS), provided hyetographs, hydrographs and pollutant loadings for biological oxygen demand (BOD₅), total suspended solids (TSS), total Kjeldahl nitrogen (TKN) and lead (Pb), and were used for calibration of the model. Several other catchment characteristics, also measured or estimated by USGS, were used in model input preparation. Application of the model was done using the Green–Ampt equation for infiltration loss computation, a pollutant accumulation equation using a power build-up equation dependent on the number of dry days, and a power wash-off equation dependent on the predicted runoff rate. Calibrated quantity input parameters are presented and compared with suggested values in the literature. The impervious depression storage was generally found to be the most sensitive calibration parameter, followed by the Manning's roughness coefficients of conduit and overland flow, the Green–Ampt infiltration parameters and, finally, the pervious depression storage. Calibrated quality input parameters are presented in the form of regression equations, as a function of rainfall depth and the number of antecedent dry days. A total of 16 independent rainfall events were used for verification of the model, which showed a good comparison with observed data for both hydrographs and pollutant loadings. Average model predictions for the four constituent concentrations from the verification runs also showed good agreement with NURP published values in Florida and US sites. © 1998 John Wiley & Sons, Ltd.