The value of a physically based model versus an empirical approach in the prediction of ephemeral gully erosion for loess-derived soils

A data set on soil losses and controlling factors for 58 ephemeral gullies has been collected in the Belgian loess belt from March 1997 to March 1999. Of the observed ephemeral gullies, 32 developed at the end of winter or in early spring (winter gullies) and 26 ephemeral gullies developed during summer (summer gullies). The assessed data have been used to test the physically based Ephemeral Gully Erosion Model (EGEM) and to compare its performance with the value of simple topographical and morphological indices in the prediction of ephemeral gully erosion.Analysis shows that EGEM is not capable of predicting ephemeral gully cross-sections well. Although conditions for input parameter assessment were ideal, some parameters such as channel erodibility, critical flow shear stress and local rainfall depth showed great uncertainty. Rather than revealing EGEM's inability of predicting ephemeral gully erosion, this analysis stresses the problematic nature of physically based models, since they often require input parameters that are not available or can hardly be obtained.With respect to the value of simple topographical and morphological indices in predicting ephemeral gully erosion, this study shows that for winter gullies and summer gullies, respectively, over 80% and about 75% of the variation in ephemeral gully volume can be explained when ephemeral gully length is known. Moreover, when previously collected data for ephemeral gullies in two Mediterranean study areas and the data for summer gullies formed in the Belgian loess belt are pooled, it appears that one single length (L)–volume (V) relation exists (V=0.048 L^1.29; R²=0.91). These findings imply that predicting ephemeral gully length is a valuable alternative for the prediction of ephemeral gully volume. A simple procedure to predict ephemeral gully length based on topographical thresholds is presented here. Secondly, the empirical length–volume relation can also be used to convert ephemeral gully length data extracted from aerial photos into ephemeral gully volumes.     

A note on magnetic modelling with remanence

The general problem of magnetic modelling involves accounting for the effect of both remanent magnetization and the application of an external magnetic field. However, as far as the disturbing field of a magnetic body in a magnetic environment is concerned, there is an equivalence between the effects of these two causations that allows the remanence to be represented in terms of an equivalent primary magnetic H field. Moreover, due to the linearity of the magnetic field in terms of its causations, the general modelling problem involving an applied magnetic field in the presence of remanence can be simply and more efficiently dealt with in terms of an equivalent primary field acting in the absence of any remanent magnetization.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

Distorted Froude‐scaled flume analysis of large woody debris

This paper presents the results of a movable‐boundary, distorted, Froude‐scaled hydraulic model based on Abiaca Creek, a sand‐bedded channel in northern Mississippi. The model was used to examine the geomorphic and hydraulic impact of simplified large woody debris (LWD) elements. The theory of physical scale models is discussed and the method used to construct the LWD test channel is developed. The channel model had bed and banks moulded from 0·8 mm sand, and flow conditions were just below the threshold of motion so that any sediment transport and channel adjustment were the result of the debris element. Dimensions and positions of LWD elements were determined using a debris jam classification model. Elements were attached to a dynamometer to measure element drag forces, and channel adjustment was determined through detailed topographic surveys. The fluid drag force on the elements decreased asymptotically over time as the channel boundary eroded around the elements due to locally increased boundary shear stress. Total time for geomorphic adjustment computed for the prototype channel at the Q₂ discharge (discharge occurring once every two years on average) was as short as 45 hours. The size, depth and position of scour holes, bank erosion and bars created by flow acceleration past the elements were found to be related to element length and position within the channel cross‐section. Morphologies created by each debris element in the model channel were comparable with similar jams observed in the prototype channel. Published in 2001 John Wiley & Sons, Ltd.     

Factors affecting sediment trapping in vegetated filter strips: simulation study using VFSMOD

Soil and water conservation practices have been promoted for a long time, in order to sustain agricultural activities and prevent environmental pollution. Vegetated filter strips (VFS) have been used to reduce sediment pollution into water bodies at or near the pollutant source. However, factors effecting VFS performance under natural conditions have not been well understood owing to the physical, time and financial limitations of field experiments. The use of well‐validated simulation models to understand the performance of VFS and factors affecting sediment deposition is highly justified. The objective of this research is to investigate sediment trapping in VFS and to study various factors affecting VFS performance using the simulation model VFSMOD, which was developed by researchers at University of North Carolina. Recently, VFSMOD has been validated successfully by using 21 filters with varying length, slope and vegetated cover. A wide range of five parameters was selected for the simulations, namely filter length, filter slope, manning roughness coefficient, soil type and characteristics of incoming sediment from adjacent fields. Computer simulations revealed that the length of filter is the most significant factor affecting sediment trapping in VFS. The relative increase in trapping efficiencies was not linearly related to an increase in filter length. Inflow sediment class also has a major influence on sediment trapping in VFS. The trapping efficiency of clay sediments in a 15 m length VFS was 47% compared with 92% for silt from incoming sediment. Manning roughness coefficient had a moderate effect on sediment trapping and was more significant in short filters. Land slope and soil type of VFS had a minor influence on the performance of VFS. Copyright © 2001 John Wiley & Sons, Ltd.     

A New Methodology for Singularity Modelling in Flow Simulations in Reservoir Engineering

In this paper, numerical methods for flow simulation in the well vicinity are discussed. Flux truncation errors are analyzed on gridblocks near the well. Due to the singularity of the well, standard numerical schemes are not efficient for near well flow simulations. To improve simulation accuracy, a new methodology, using a change of coordinates to make near well pressure linear, is presented for the singularity modelling. Based on this new approach, two numerical schemes, a two-point flux approximation scheme and a multi-point one, are proposed. These schemes are accurate for the near well modelling and are suitable for any kinds of gridblocks.     

复杂地震系统的数字模拟

借用细胞自动机的计算处理方法，对地震孕育过程中整个系统能量输入，积累，释放，传播等进行了数字模拟，并研究了复杂地震系统的整体特征，结果表明，整体调整系统，地震时能量的损耗系数，外部构造应力场的驱动力大小会影响到地震系统中能量积累和释放的状况，影响到Ｇ－Ｒ关系中ｂ值的大小。     

Processing of oxidised nitrogen compounds by passage through winter-time orographic cloud

Four case studies are described, from a three-site field experiment in October/November 1991 using the Great Dun Fell flow-through reactor hill cap cloud in rural Northern England. Measurements of total odd-nitrogen nitrogen oxides (NO _y ) made on either side of the hill, before and after the air flowed through the cloud, showed that 10 to 50% of the NO _y , called NO _z , was neither NO nor NO₂. This NO _z failed to exhibit a diurnal variation and was often higher after passage through cloud than before. No evidence of conversion of NO _z to NO₃ ^- in cloud was found. A simple box model of gas-phase chemistry in air before it reached the cloud, including scavenging of NO₃ and N₂O₅ by aerosol of surface area proportional to the NO₂ mixing ratio, shows that NO₃ and N₂O₅ may build up in the boundary layer by night only if stable stratification insulates the air from emissions of NO. This may explain the lack of evidence for N₂O₅ forming NO₃ ^- in cloud under well-mixed conditions in 1991, in contrast with observations under stably stratified conditions during previous experiments when evidence of N₂O₅ was found. Inside the cloud, some variations in the calculated total atmospheric loading of HNO₂ and the cloud liquid water content were related to each other. Also, indications of conversion of NO _x to NO _z were found. To explain these observations, scavenging of NO _x and HNO₂ by cloud droplets and/or aqueous-phase oxidation of NO₂ ^- by nitrate radicals are considered. When cloud acidity was being produced by aqueous-phase oxidation of NO _x or SO₂, NO₃ ^- which had entered the cloud as aerosol particles was liberated as HNO₃ vapour. When no aqueous-phase production of acidity was occurring, the reverse, conversion of scavenged HNO₃ to particulate NO₃ ^-, was observed.     

Modelling of ground motion in the vicinity of massive structures

A two-dimensional elastic Chebyshev spectral element method (SPEM) is used to model the seismic wavefield within a massive structure and in its vicinity. We consider 2-D models where a linear elastic structure, with quadrangular cross-section, resting on an elastic homogeneous half-space, is impinged upon by the waves generated by a surface impulse at some distance. The scattering of Rayleigh waves and the response of the structure are extensively analysed in a parametric way, varying size, mechanical parameters and shape of the load. Some of the models considered are representative of embankments and earth dams. The simulation shows that some models resonate, storing part of the incoming energy. With realistic parameters, the lowest resonance frequency is due to pure shear deformation and is controlled by the shear velocity and height of the load. Flexural modes are excited only at higher frequencies. The acceleration at the top of the structure may be five/seven times higher than at the base, depending on the mass of the structure. The gradual release of trapped energy produces a ground roll lasting several seconds after the wave front has passed. The ground-roll amplitude depends on the sturcture's mass and can be as large as 30% of the peak acceleration. Outside resonance conditions, the ground motion is almost unaffected by the presence of the artefact; the horizontal motion on top of it is nearly twice the motion at ground level. Similar results should be expected when the incident field is an upcoming shear wave. A qualitative discussion shows that the presence of anelastic attenuation in the embankment does not significantly alter the preceding conclusions, unless it is of very low values (e.g. Q < 15).The modelling results that we discuss indicate that the soil-structure interaction may substantially alter the ‘free-field’ ground motion. From a practical point of view, the main conclusions are: (1) careful analysis is necessary when interpreting seismic records collected in the vicinity of large artefacts; (2) seismic hazard at a site may depend on the presence of man-made structures such as embankments, dams, tall and massive buildings.