Evaluating soil water routing approaches in watershed-scale,ecohydrologic modelling

Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.     

Modelling the poroelastoplastic behaviour of soils subjected to internal erosion by suffusion

Granular soils subjected to seepage flow may suffer suffusion, ie, a selective internal erosion. Extending the classical approach of poromechanics, we deduce a new form of the Clausius-Duhem inequality accounting for dissipation due to suffusion, and we deduce restrictions on the constitutive laws of the soil. We suggest (a) a possible coupling between the seepage forces and the suffusion kinetics and (b) an extension of an existing elastoplastic model for the skeleton mechanical behaviour. Numerical integrations of the elastoplastic model are carried out under drained axisymmetric triaxial and oedometric conditions. As a result, we prove that the extended model is able to qualitatively reproduce the suffusion induced strains and the strength reduction experimentally observed. Predictions on the oedometric behaviour of suffusive soils are also provided.     

自然和自然贡献情景模型研究综述

本文回顾了自然和自然贡献情景模型发展的背景、历史和内容,概括总结了自然和自然贡献情景模型的发展进程以及联合国生物多样性与生态系统服务政府间科学—政策平台(IPBES)情景模型的概念框架,讨论了自然和自然贡献情景模型存在的问题和发展方向。为了在全球层面解决现有综合集成模型存在的问题,根据地球表层建模基本定理和生态环境曲面建模基本定理,提出了具有中国原创特点的自然与自然贡献情景模型概念框架。     

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.     

Reef exploration in the East Sengkang Basin,Sulawesi, Indonesia

Reconnaissance seismic shot in 1971/72 showed a number of well defined seismic anomalies within the East Sengkang Basin which were interpreted as buried reefs. Subsequent fieldwork revealed that Upper Miocene reefs outcropped along the southern margin of the basin. A drilling programme in 1975 and 1976 proved the presence of shallow, gas-bearing, Upper Miocene reefs in the northern part of the basin. Seismic acquisition and drilling during 1981 confirmed the economic significance of these discoveries, with four separate accumulations containing about 750 × 10⁹ cubic feet of dry gas in place at an average depth of 700 m. Kampung Baru is the largest field and contains over half the total, both reservoir quality and gas deliverability are excellent. Deposition in the East Sengkang Basin probably started during the Early Miocene. A sequence of Lower Miocene mudstones and limestones unconformably overlies acoustic basement which consists of Eocene volcanics. During the tectonically active Middle Miocene, deposition was interrupted by two periods of deformation and erosion. Carbonate deposition became established in the Late Miocene with widespread development of platform limestones throughout the East Sengkang Basin. Thick pinnacle reef complexes developed in the areas where reef growth could keep pace with the relative rise in sea level. Most reef growth ceased at the end of the Miocene and subsequent renewed clastic sedimentation covered the irregular limestone surface. Late Pliocene regression culminated in the Holocene with erosion. The Walanae fault zone, part of a major regional sinistral strike-slip system, separates the East and West Sengkang Basins. Both normal and reverse faulting are inferred from seismic data and post Late Pliocene reverse faulting is seen in outcrop.     

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.