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.     

A simulation tool for high data-rate acoustic communication in ashallow-water, time-varying channel

The paper discusses the development of a simulation tool to model high data-rate acoustic communication in shallow water. The simulation tool is able to generate synthetic time series of signals received at a transducer array after transmission across a shallow-water communication channel. The simulation tool is suitable for testing advanced signal processing techniques for message recovery. A channel model has been developed based on the physical aspects of the acoustic channel. Special emphasis has been given to fluctuations of the signal transmission caused by time-varying multipath effects. At shorter ranges, the temporal variations are dominated by acoustic scattering from the moving sea surface. Therefore, the channel model produces a coherence function which may be interpreted as a time-varying reflection coefficient for the surface scattered acoustical path. A static, range-independent ray model identifies the significant multipaths, and the surface path is modulated with the time-varying reflection coefficient. The advantages and limitations of the channel model are discussed and assumptions necessary to overcome the limitations are emphasised. Based on the assumptions, an algorithm has been developed and implemented to model how a binary message will be modulated when transmitted by a transducer, is distorted in the channel and finally is received by a transducer array     

OSCR wave measurements-some preliminary results

OSCR is an HF radar system that has been developed for high spatial resolution coastal surface current measurement. This paper describes preliminary results that demonstrate that wave measurement can be successfully obtained from suitably processed OSCR data. Comparisons with data from a WAVEC directional buoy are presented and show encouraging agreement. Some of the limitations to the measurement process are discussed and indicate a maximum range of about 20 km. Surface current variability on short time scales presents the most serious obstacle to wave measurement. This appears to be more of a problem when the mean currents are large, in that in these circumstances the data fail initial quality control criteria. However, in lower mean currents, the effect is often still present and leads to errors in long wave measurement     

Application of an improved self-starter to geoacoustic inversion

The self-starter is improved using the operator of the split-step Pade solution. In addition to providing greater stability and being applicable closer to the source, the improved self-starter is an efficient forward model for geoacoustic inversion. It is necessary to solve only O(10) tridiagonal systems of equations to obtain the acoustic field on a vertical array located O(10) wavelengths from a source. This experimental configuration is effective for geoacoustic inverse problems involving unknown parameters deep in the ocean bottom. For problems involving depth-dependent acoustic parameters, the improved self-starter can be used to solve nonlinear inverse problems involving O(10) unknown sediment parameters in less than a minute on the current generation of workstations     

A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. 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.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

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.     

Rapid vegetation change during the early Holocene in the Faroe Islands detected in terrestrial and aquatic ecosystems

High‐resolution pollen, plant macrofossil and sedimentary analyses from early Holocene lacustrine sediments on the Faroe Islands have detected a significant vegetation perturbation suggesting a rapid change in climate between ca. 10 380 cal. yr BP and the Saksunarvatn ash (10 240±60 cal. yr BP). This episode may be synchronous with the decline in δ¹⁸O values in the Greenland ice‐cores. It also correlates with a short, cold event detected in marine cores from the North Atlantic that has been ascribed to a weakening of thermohaline circulation associated with the sudden drainage of Lake Agassiz into the northwest Atlantic, or, alternatively, a period with distinctly decreased solar forcing. The vegetation sequence begins at ca. 10 500 cal. yr BP with a succession from tundra to shrub‐tundra and increasing lake productivity. Rapid population increases of aquatic plants suggest high summer temperatures between 10 450 and 10 380 cal. yr BP. High pollen percentages, concentrations and influx of Betula, Juniperus and Salix together with macrofossil leaves indicate shrub growth around the site during the initial phases of vegetation colonisation. Unstable conditions followed ca. 10 380 cal. yr BP that changed both the upland vegetation and the aquatic plant communities. A decrease in percentage values of shrub pollen is recorded, with replacement of both aquatics and herbaceous plants by pioneer plant communities. An increase in total pollen accumulation rates not seen in the concentration data suggests increased sediment delivery. The catchment changes are consistent with less seasonal, moister conditions. Subsequent climatic amelioration reinitiated a warmth‐driven succession and catchment stabilisation, but retained high precipitation levels influencing the composition of the post‐event communities. Copyright © 2003 John Wiley & Sons, Ltd.