Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

An agent-based approach to modelling driver route choice behaviour under the influence of real-time information

This paper presents an agent-based approach to modelling individual driver behaviour under the influence of real-time traffic information. The driver behaviour models developed in this study are based on a behavioural survey of drivers which was conducted on a congested commuting corridor in Brisbane, Australia. Commuters' responses to travel information were analysed and a number of discrete choice models were developed to determine the factors influencing drivers' behaviour and their propensity to change route and adjust travel patterns. Based on the results obtained from the behavioural survey, the agent behaviour parameters which define driver characteristics, knowledge and preferences were identified and their values determined. A case study implementing a simple agent-based route choice decision model within a microscopic traffic simulation tool is also presented. Driver-vehicle units (DVUs) were modelled as autonomous software components that can each be assigned a set of goals to achieve and a database of knowledge comprising certain beliefs, intentions and preferences concerning the driving task. Each DVU provided route choice decision-making capabilities, based on perception of its environment, that were similar to the described intentions of the driver it represented. The case study clearly demonstrated the feasibility of the approach and the potential to develop more complex driver behavioural dynamics based on the belief-desire-intention agent architecture. (C) 2002 Elsevier Science Ltd. All rights reserved.

Prediction of carbon monoxide in fires by conditional moment closure

Carbon monoxide is the chief killer in fires. Dangerous levels of CO can occur when reacting combustion gases are quenched by heat transfer, or by mixing of the fire plume in a cooled under- or overventilated upper layer. In this paper, carbon monoxide predictions for enclosure fires are modeled by the conditional moment closure (CMC) method and are compared with laboratory data. The modeled fire situation is a buoyant, turbulent, diffusion flame burning under a hood. The fire plume entrains fresh air, and the postflame gases are cooled considerably under the hood by conduction and radiation, emulating conditions which occur in enclosure fires and lead to the freezing of CO burnout. Predictions of CO in the cooled layer are presented in the context of a complete computational fluid dynamics solution of velocity, temperature, and major species concentrations. A range of underhood equivalence ratios, from rich to lean, are investigated. The CMC method predicts CO in very good agreement with data. In particular, CMC is able to correctly predict CO concentrations in lean cooled gases, showing its capability in conditions where reaction rates change considerably.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Simultaneous X-ray imaging of plant root growth and water uptake in thin-slab systems

Direct and simultaneous observation of root growth and plant water uptake is difficult because soils are opaque. X-ray imaging techniques such as projection radiography or Computer Tomography (CT) offer a partial alternative to such limitations. Nevertheless, there is a trade-off between resolution, large field-of-view and 3-dimensionality: With the current state of the technology, it is possible to have any two. In this study, we used X-ray transmission through thin-slab systems to monitor transient saturation fields that develop around roots as plants grow. Although restricted to 2-dimensions, this approach offers a large field-of-view together with high spatial and dynamic resolutions. To illustrate the potential of this technology, we grew peas in 1 cm thick containers filled with soil and imaged them at regular intervals. The dynamics of both the root growth and the water content field that developed around the roots could be conveniently monitored. Compared to other techniques such as X-ray CT, our system is relatively inexpensive and easy to implement. It can potentially be applied to study many agronomic problems, such as issues related to the impact of soil constraints (physical, chemical or biological) on root development.

Measuring the decoherence rate in a semiconductor charge qubit

We describe a method by which the decoherence time of a solid-state qubit may be measured. The qubit is coded in the orbital degree of freedom of a single electron bound to a pair of donor impurities in a semiconductor host. The qubit is manipulated by adiabatically varying an external electric field. We show that by measuring the total probability of a successful qubit rotation as a function of the control field parameters, the decoherence rate may be determined. We estimate various system parameters, including the decoherence rates due to electromagnetic fluctuations and acoustic phonons. We find that, for reasonable physical parameters, the experiment is possible with existing technology. In particular, the use of adiabatic control fields implies that the experiment can be performed with control electronics with a time resolution of tens of nanoseconds.

Boys performing English

This article explores the gender dynamics of boys' responses to one Particular aspect of English teaching. oral performance work. It focuses on the possibility that the requirement to perform publicly in dramatic and other oral tasks may be an important factor in the rejection of English by many boys, and may contribute to boys' relatively poor achievement in English. The article provides a study of boys' engagement in English oral activities in two classrooms, and identifies a number of factors influencing boys' English learning. In particular, it shows that there is no simple relation between the performance requirements of English learning activities and boys' disengagement with English.

Controlling chaos in a Lorenz-like system using feedback

We demonstrate that the dynamics of an autonomous chaotic laser can be controlled to a periodic or steady state under self-synchronization. In general, past the chaos threshold the dependence of the laser output on feedback applied to the pump is submerged in the Lorenz-like chaotic pulsation. However there exist specific feedback delays that stabilize the chaos to periodic behavior or even steady state. The range of control depends critically on the feedback delay time and amplitude. Our experimental results are compared with the complex Lorenz equations which show good agreement.

The effect of nutrient enrichment on growth, photosynthesis and hydraulic conductance of dwarf mangroves in Panama

1. Dwarf stands of the mangrove Rhizophora mangle L. are extensive in the Caribbean. We fertilized dwarf trees in Almirante Bay, Bocas del Toro Province, north-eastern Panama with nitrogen (N) and phosphorus (P) to determine (1) if growth limitations are due to nutrient deficiency; and (2) what morphological and/or physiological factors underlie nutrient limitations to growth. 2. Shoot growth was 10-fold when fertilized with P and twofold with N fertilization, indicating that stunted growth of these mangroves is partially due to nutrient deficiency. 3. Growth enhancements caused by N or P enrichment could not be attributed to increases in photosynthesis on a leaf area basis, although photosynthetic nutrient-use efficiency was improved. The most dramatic effect was on stem hydraulic conductance, which was increased sixfold by P and 2.5-fold with N enrichment. Fertilization with P enhanced leaf and stem P concentrations and reduced C : N ratio, but did not alter leaf damage by herbivores. 4. Our findings indicate that addition of N and P significantly alter tree growth and internal nutrient dynamics of mangroves at Bocas del Toro, but also that the magnitude, pattern and mechanisms of change will be differentially affected by each nutrient.

Use of bulk density profiles from X-radiography to examine structural crust models

Exponential and sigmoidal functions have been suggested to describe the bulk density profiles of crusts. The present work aims to evaluate these conceptual models using high resolution X-radiography. Repacked seedbeds from two soil materials, air-dried or prewetted by capillary rise, were subjected to simulated rain, which resulted in three types of structural crusts, namely, slaking, infilling, and coalescing. Bulk density distributions with depth were generated using high-resolution (70 mum), calibrated X-ray images of slices from the resin-impregnated crusted seedbeds. The bulk density decreased progressively with depth, which supports the suggestion that a crust should be considered as a nonuniform layer. For the slaking and the coalescing crusts, the exponential function underestimated the strong change in bulk density across the morphologically defined transition between the crust and the underlying material; the sigmoidal function provided a better description. Neither of these crust models effectively described the shape of the bulk density profiles through the whole seedbed. Below the infilling and slaking crusts, bulk density increased linearly with depth as a result of slumping. In the coalescing crusted seedbed, the whole seedbed uniformly collapsed and most of the bulk density change within the crust could be ascribed to slumping (0.33 g cm(-3)) rather than to crusting (0.12 g cm(-3)). Finally, (i) X-radiography appears as a unique tool to generate high resolution bulk density profiles and (ii) in structural crusts, bulk density profiles could be modeled using the existing exponential and sigmoidal crusting models, provided a slumping model would be coupled.

Comparing predictions of extinction risk using models and subjective judgement

Models of population dynamics are commonly used to predict risks in ecology, particularly risks of population decline. There is often considerable uncertainty associated with these predictions. However, alternatives to predictions based on population models have not been assessed. We used simulation models of hypothetical species to generate the kinds of data that might typically be available to ecologists and then invited other researchers to predict risks of population declines using these data. The accuracy of the predictions was assessed by comparison with the forecasts of the original model. The researchers used either population models or subjective judgement to make their predictions. Predictions made using models were only slightly more accurate than subjective judgements of risk. However, predictions using models tended to be unbiased, while subjective judgements were biased towards over-estimation. Psychology literature suggests that the bias of subjective judgements is likely to vary somewhat unpredictably among people, depending on their stake in the outcome. This will make subjective predictions more uncertain and less transparent than those based on models. (C) 2004 Elsevier SAS. All rights reserved.

Variable effects of larval size on post-metamorphic performance in the field

Larval quality may be capable of explaining much of the variation in the recruitment and subsequent population dynamics of benthic marine invertebrates. Whilst the effects of larval nutritional condition on adult performance have received the most attention, recent work has shown that larval size may also be an important and ubiquitous source of variation in larval quality. We examined the effects of variation in larval size on the post-metamorphic survival and growth of Watersipora subtorquata in 2 very different habitats - experimental substrata and pier pilings. We found strong effects of larval size on colony performance, although these varied among experiments. For colonies on experimental substrata, larval size positively affected adult survival and, initially, growth. However, after 3 wk in the field, there was no relationship between larval size and colony size, possibly because colonies were completely surrounded by newly settled organisms. Larval size also positively affected post-metamorphic growth of colonies on pier pilings, but, surprisingly, colonies that came from larger larvae had lower survival than colonies from smaller larvae. Overall, variation in larval size will strongly affect the recruitment and subsequent performance of adults in this species, although this may vary among different habitats. This study highlights the importance of examining the effects of larval quality on adult performance in as realistic conditions as possible, because of the strong interaction between larval size effects and the environment.

Quantum phase-space analysis of the pendular cavity

We perform a quantum-mechanical analysis of the pendular cavity, using the positive-P representation, showing that the quantum state of the moving mirror, a macroscopic object, has noticeable effects on the dynamics. This system has previously been proposed as a candidate for the quantum-limited measurement of small displacements of the mirror due to radiation pressure, for the production of states with entanglement between the mirror and the field, and even for superposition states of the mirror. However, when we treat the oscillating mirror quantum mechanically, we find that it always oscillates, has no stationary steady state, and exhibits uncertainties in position and momentum which are typically larger than the mean values. This means that previous linearized fluctuation analyses which have been used to predict these highly quantum states are of limited use. We find that the achievable accuracy in measurement is fat, worse than the standard quantum limit due to thermal noise, which, for typical experimental parameters, is overwhelming even at 2 mK