Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Modelling of post-irradiation events in polymer gel dosimeters

The nuclear magnetic resonance (NMR) spin-spin relaxation time (T-2) is related to the radiation-dependent concentration of polymer formed in polymer gel dosimeters manufactured from monomers in an aqueous gelatin matrix. Changes in T-2 with time post-irradiation have been reported in the literature but their nature is not fully understood. We investigated those changes with time after irradiation using FT-Raman spectroscopy and the precise determination of T-2 at high magnetic field in a polymer gel dosimeter, A model of fast exchange of magnetization taking into account ongoing gelation and strengthening of the gelatin matrix as well as the polymerization of the monomers with time is presented. Published data on the changes of T-2 in gelatin gels as a function of post-manufacture time are used and fitted closely by the model presented. The same set of parameters characterizing the variations of T-2 in gelatin gels and the increasing concentration of polymer determined from Fr-Raman spectroscopy are used successfully in the modelling of irradiated polymer gel dosimeters. Minimal variations in T-2 in an irradiated PAG dosimeter are observed after 13 h.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

BeachWin: modelling groundwater effects on swash sediment transport and beach profile changes

Field and laboratory observations have shown that a relatively low beach groundwater table enhances beach accretion. These observations have led to the beach dewatering technique (artificially lowering the beach water table) for combating beach erosion. Here we present a process-based numerical model that simulates the interacting wave motion on the beach. coastal groundwater flow, swash sediment transport and beach profile changes. Results of model simulations demonstrate that the model replicates accretionary effects of a low beach water table on beach profile changes and has the potential to become a tool for assessing the effectiveness of beach dewatering systems. (C) 2002 Elsevier Science Ltd. All rights reserved.

Diel vertical migration: Modelling light-mediated mechanisms

Light is generally regarded as the most likely cue used by zooplankton to regulate their vertical movements through the water column. However, the way in which light is used by zooplankton as a cue is not well understood. In this paper we present a mathematical model of diel vertical migration which produces vertical distributions of zooplankton that vary in space and time. The model is used to predict the patterns of vertical distribution which result when animals are assumed to adopt one of three commonly proposed mechanisms for vertical swimming. First, we assume zooplankton tend to swim towards a preferred intensity of light. We then assume zooplankton swim in response to either the rate of change in light intensity or the relative rate of change in light intensity. The model predicts that for all three mechanisms movement is fastest at sunset and sunrise and populations are primarily influenced by eddy diffusion at night in the absence of a light stimulus. Daytime patterns of vertical distribution differ between the three mechanisms and the reasons for the predicted differences are discussed. Swimming responses to properties of the light field are shown to be adequate for describing diel vertical migration where animals congregate in near surface waters during the evening and reside at deeper depths during the day. However, the model is unable to explain how some populations halt their ascent before reaching surface waters or how populations re-congregate in surface waters a few hours before sunrise, a phenomenon which is sometimes observed in the held. The model results indicate that other exogenous or endogenous factors besides light may play important roles in regulating vertical movement.

A formal representation of assumptions in process modelling

In this work, we present a systematic approach to the representation of modelling assumptions. Modelling assumptions form the fundamental basis for the mathematical description of a process system. These assumptions can be translated into either additional mathematical relationships or constraints between model variables, equations, balance volumes or parameters. In order to analyse the effect of modelling assumptions in a formal, rigorous way, a syntax of modelling assumptions has been defined. The smallest indivisible syntactical element, the so called assumption atom has been identified as a triplet. With this syntax a modelling assumption can be described as an elementary assumption, i.e. an assumption consisting of only an assumption atom or a composite assumption consisting of a conjunction of elementary assumptions. The above syntax of modelling assumptions enables us to represent modelling assumptions as transformations acting on the set of model equations. The notion of syntactical correctness and semantical consistency of sets of modelling assumptions is defined and necessary conditions for checking them are given. These transformations can be used in several ways and their implications can be analysed by formal methods. The modelling assumptions define model hierarchies. That is, a series of model families each belonging to a particular equivalence class. These model equivalence classes can be related to primal assumptions regarding the definition of mass, energy and momentum balance volumes and to secondary and tiertinary assumptions regarding the presence or absence and the form of mechanisms within the system. Within equivalence classes, there are many model members, these being related to algebraic model transformations for the particular model. We show how these model hierarchies are driven by the underlying assumption structure and indicate some implications on system dynamics and complexity issues. (C) 2001 Elsevier Science Ltd. All rights reserved.

Stay-green: A consequence of the balance between supply and demand for nitrogen during grain filling?

Retention of green leaf area in grain sorghum under post-anthesis drought, known as stay-green, is associated with greater biomass production, lodging resistance and yield. The stay-green phenomenon can be examined at a cell, leaf, or whole plant level. At a cell level, the retention of chloroplast proteins such as LHCP2, OEC33 and Rubisco until late in senescence has been reported in sorghum containing the KS19 source of stay-green, indicating that photosynthesis may be maintained for longer during senescence in these genotypes. At a leaf level, longevity of photosynthetic apparatus is intimately related to nitrogen (N) status. At a whole plant level, stay-green can be viewed as a consequence of the balance between N demand by the grain and N supply during grain filling. To examine some of these concepts, nine hybrids varying in the B35 and KS19 sources of stay-green were grown under a postanthesis water deficit. Genotypic variation in delayed onset and reduced rate of leaf senescence were explained by differences in specific leaf nitrogen (SLN) and N uptake during grain filling. Matching N supply from age-related senescence and N uptake during grain tilling with grain N demand found that the shortfall in N supply for grain filling was greater in the senescent than stay-green hybrids, resulting in more accelerated leaf senescence in the former. We hypothesise that increased N uptake by stay-green hybrids is a result of greater biomass accumulation during grain filling in response to increased sink demand (higher grain numbers) which, in turn, is the result of increased radiation use efficiency and transpiration efficiency due to higher SLN. Delayed leaf senescence resulting from higher SLN should, in turn, allow snore carbon and nitrogen to be allocated to the roots of stay-green hybrids during grain filling, thereby maintaining a greater capacity to extract N from the soil compared with senescent hybrids.

The modelling of froth zone recovery in batch and continuously operated laboratory flotation cells

This paper proposes an integrated methodology for modelling froth zone performance in batch and continuously operated laboratory flotation cells. The methodology is based on a semi-empirical approach which relates the overall flotation rate constant to the froth depth (FD) in the flotation cell; from this relationship, a froth zone recovery (R,) can be extracted. Froth zone recovery, in turn, may be related to the froth retention time (FRT), defined as the ratio of froth volume to the volumetric flow rate of concentrate from the cell. An expansion of this relationship to account for particles recovered both by true flotation and entrainment provides a simple model that may be used to predict the froth performance in continuous tests from the results of laboratory batch experiments. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Combined hydraulic and biological modelling and full-scale validation of SBR process

The biological reactions during the settling and decant periods of Sequencing Batch Reactors (SBRs) are generally ignored as they are not easily measured or described by modelling approaches. However, important processes are taking place, and in particular when the influent is fed into the bottom of the reactor at the same time (one of the main features of the UniFed process), the inclusion of these stages is crucial for accurate process predictions. Due to the vertical stratification of both liquid and solid components, a one-dimensional hydraulic model is combined with a modified ASM2d biological model to allow the prediction of settling velocity, sludge concentration, soluble components and biological processes during the non-mixed periods of the SBR. The model is calibrated on a full-scale UniFed SBR system with tracer breakthrough tests, depth profiles of particulate and soluble compounds and measurements of the key components during the mixed aerobic period. This model is then validated against results from an independent experimental period with considerably different operating parameters. In both cases, the model is able to accurately predict the stratification and most of the biological reactions occurring in the sludge blanket and the supernatant during the non-mixed periods. Together with a correct description of the mixed aerobic period, a good prediction of the overall SBR performance can be achieved.

Shear stress and sediment transport calculations for swash zone modelling

It is shown that the observed difference in sediment transporting efficiency by the swash uprush, compared with the downrush, could be mainly due to greater bed shear stress for a given velocity in the more abruptly accelerated uprush. The bed shear stress generated by an arbitrary free stream velocity time series is modelled in terms of usual wave boundary layer models plus a phase lead (phi(tau) of the bed shear stress compared with the free stream velocity at the peak frequency. With this approach, the total transport amounts in uprush and downrush can be modelled satisfactorily with the same sediment transport formula, without the need for different uprush and downrush coefficients. While the adaptation of sediment transport formulae from steady flow can thus lead to the right total amounts of sediment moved by this method, the timing of the instantaneous sediment transport rates are probably not accurately modelled due to the highly unsteady nature of the swash and the presence of pre-suspended sediment in the uprush. Nevertheless, the proposed method is a useful intermediate step before we have a complete understanding of sediment transport under very rapid accelerations and of the relative contribution of pre-suspended sediment to the onshore sediment transport in swash zones. (C) 2002 Published by Elsevier Science B.V.

Review and future directions in the modelling and control of continuous drum granulation

Many granulation plants operate well below design capacity, suffering from high recycle rates and even periodic instabilities. This behaviour cannot be fully predicted using the present models. The main objective of the paper is to provide an overview of the current status of model development for granulation processes and suggest future directions for research and development. The end-use of the models is focused on the optimal design and control of granulation plants using the improved predictions of process dynamics. The development of novel models involving mechanistically based structural switching methods is proposed in the paper. A number of guidelines are proposed for the selection of control relevant model structures. (C) 2002 Published by Elsevier Science B.V.

Numerical modelling of contaminant transport in coastal aquifers

This paper employs a two-dimensional variable density flow and transport model to investigate the transport of a dense contaminant plume in an unconfined coastal aquifer. Experimental results are also presented to show the contaminant plume in a freshwater-seawater flow system. Both the numerical and experimental results suggest that the neglect of the seawater interface does not noticeably affect the horizontal migration rate of the plume before it reaches the interface. However, the contaminant will travel further seaward and part of the solute mass will exit under the sea if the higher seawater density is not included. If the seawater density is included, the contaminant will travel upwards towards the beach along the freshwater-saltwater interface as shown experimentally. Neglect of seawater density, therefore, will result in an underestimate of solute mass rate exiting around the coastline. (C) 2002 IMACS. Published by Elsevier Science B.V. All rights reserved.

Using crop simulation to generate genotype by environment interaction effects for sorghum in water-limited environments

Multi-environment trials (METs) used to evaluate breeding lines vary in the number of years that they sample. We used a cropping systems model to simulate the target population of environments (TPE) for 6 locations over 108 years for 54 'near-isolines' of sorghum in north-eastern Australia. For a single reference genotype, each of 547 trials was clustered into 1 of 3 'drought environment types' (DETs) based on a seasonal water stress index. Within sequential METs of 2 years duration, the frequencies of these drought patterns often differed substantially from those derived for the entire TPE. This was reflected in variation in the mean yield of the reference genotype. For the TPE and for 2-year METs, restricted maximum likelihood methods were used to estimate components of genotypic and genotype by environment variance. These also varied substantially, although not in direct correlation with frequency of occurrence of different DETs over a 2-year period. Combined analysis over different numbers of seasons demonstrated the expected improvement in the correlation between MET estimates of genotype performance and the overall genotype averages as the number of seasons in the MET was increased.

Predicting plant leaf area production: shoot assimilate accumulation and partitioning, and leaf area ratio, are stable for a wide range of sorghum population densities

Predicting plant leaf area production is required for modelling carbon balance and tiller dynamics in plant canopies. Plant leaf area production can be studied using a framework based on radiation intercepted, radiation use efficiency (RUE) and leaf area ratio (LAR) (ratio of leaf area to net above-ground biomass). The objective of this study was to test this framework for predicting leaf area production of sorghum during vegetative development by examining the stability of the contributing components over a large range of plant density. Four densities, varying from 2 to 16 plants m(-2), were implemented in a field experiment. Plants were either allowed to tiller or were maintained as uniculm by systematic tiller removal. In all cases, intercepted radiation was recorded daily and leaf area and shoot dry matter partitioning were quantified weekly at individual culm level. Up to anthesis, a unique relationship applied between fraction of intercepted radiation and leaf area index, and between shoot dry weight accumulation and amount of intercepted radiation, regardless of plant density. Partitioning of shoot assimilate between leaf, stem and head was also common across treatments up to anthesis, at both plant and culm levels. The relationship with thermal time (TT) from emergence of specific leaf area (SLA) and LAR of tillering plants did not change with plant density. In contrast, SLA of uniculm plants was appreciably lower under low-density conditions at any given TT from emergence. This was interpreted as a consequence of assimilate surplus arising from the inability of the plant to compensate by increasing the leaf area a culm could produce. It is argued that the stability of the extinction coefficient, RUE and plant LAR of tillering plants observed in these conditions provides a reliable way to predict leaf area production regardless of plant density. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.