Non-invasive imaging of roots with high resolution X-ray micro-tomography

ray micro-tomography is a well-established technique for non-invasive imaging and evaluation of heterogeneous materials. An inexpensive X-ray micro-tomography system has been designed and built for the specific purposes of examining root growth and root/soil interactions. The system uses a silver target X-ray source with a focal spot diameter of 80 mum, an X-ray image intensifier with a sampling aperture of about 100 mum, and a sample with a diameter of 25 mm. Pre-germinated wheat and rape seeds were grown for up to 8-10 days in plastic containers in a sandy loam soil sieved to < 250 μm, and imaged with the X-ray system at regular intervals. The quality of 3 D image obtained was good allowing the development and growth of both root axes and some first-order laterals to be observed. The satisfactory discrimination between soil and roots enabled measurements of root diameter (wheat values were 0.48-1.22 mm) in individual tomographic slices and, by tracking from slice to slice, root lengths were also measured. The measurements obtained were generally within 10% of those obtained from destructive samples measured manually and with a flat-bed scanner. Further developments of the system will allow more detailed examination of the root: soil interface.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Modelling impacts of pollution in river systems: a new dispersion model and a case study of mine discharges in the Abrud, Aries and Mures river system in Transylvania, Romania

A new model of dispersion has been developed to simulate the impact of pollutant discharges on river systems. The model accounts for the main dispersion processes operating in rivers as well as the dilution from incoming tributaries and first-order kinetic decay processes. The model is dynamic and simulates the hourly behaviour of river flow and pollutants along river systems. The model has been applied to the Aries and Mures River System in Romania and has been used to assess the impacts of potential dam releases from the Roia Montan Mine in Transylvania, Romania. The question of mine water release is investigated under a range of scenarios. The impacts on pollution levels downstream at key sites and at the border with Hungary are investigated.

Monte Carlo phase diagram for diblock copolymer melts

A partial phase diagram is constructed for diblock copolymer melts using lattice-based Monte Carlo simulations. This is done by locating the order-disorder transition (ODT) with the aid of a recently proposed order parameter and identifying the ordered phase over a wide range of copolymer compositions (0.2 <= f <= 0.8). Consistent with experiments, the disordered phase is found to exhibit direct first-order transitions to each of the ordered morphologies. This includes the spontaneous formation of a perforated-lamellar phase, which presumably forms in place of the gyroid morphology due to finite-size and/or nonequilibrium effects. Also included in our study is a detailed examination of disordered cylinder-forming (f=0.3) diblock copolymers, revealing a substantial degree of pretransitional chain stretching and short-range order that set in well before the ODT, as observed previously in analogous studies on lamellar-forming (f=0.5) molecules. (c) 2006 American Institute of Physics.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

LiDAR mapping of tidal marshes for ecogeomorphological modelling in the TIDE project

The European research project TIDE (Tidal Inlets Dynamics and Environment) is developing and validating coupled models describing the morphological, biological and ecological evolution of tidal environments. The interactions between the physical and biological processes occurring in these regions requires that the system be studied as a whole rather than as separate parts. Extensive use of remote sensing including LiDAR is being made to provide validation data for the modelling. This paper describes the different uses of LiDAR within the project and their relevance to the TIDE science objectives. LiDAR data have been acquired from three different environments, the Venice Lagoon in Italy, Morecambe Bay in England, and the Eden estuary in Scotland. LiDAR accuracy at each site has been evaluated using ground reference data acquired with differential GPS. A semi-automatic technique has been developed to extract tidal channel networks from LiDAR data either used alone or fused with aerial photography. While the resulting networks may require some correction, the procedure does allow network extraction over large areas using objective criteria and reduces fieldwork requirements. The networks extracted may subsequently be used in geomorphological analyses, for example to describe the drainage patterns induced by networks and to examine the rate of change of networks. Estimation of the heights of the low and sparse vegetation on marshes is being investigated by analysis of the statistical distribution of the measured LiDAR heights. Species having different mean heights may be separated using the first-order moments of the height distribution.

A pattern language for goals

Previous work has established the value of goal-oriented approaches to requirements engineering. Achieving clarity and agreement about stakeholders’ goals and assumptions is critical for building successful software systems and managing their subsequent evolution. In general, this decision-making process requires stakeholders to understand the implications of decisions outside the domains of their own expertise. Hence it is important to support goal negotiation and decision making with description languages that are both precise and expressive, yet easy to grasp. This paper presents work in progress to develop a pattern language for describing goal refinement graphs. The language has a simple graphical notation, which is supported by a prototype editor tool, and a symbolic notation based on modal logic.

The propagation mechanism of a vortex on the β plane

The propagation velocity and propagation mechanism for vortices on a β plane are determined for a reduced-gravity model by integrating the momentum equations over the β plane. Isolated vortices, vortices in a background current, and initial vortex propagation from rest are studied. The propagation mechanism for isolated anticyclones as well as cyclones, which has been lacking up to now, is presented. It is shown that, to first order, the vortex moves to generate a Coriolis force on the mass anomaly of the vortex to compensate for the force on the vortex due to the variation of the Coriolis parameter. Only the mass anomaly of the vortex is of importance, because the Coriolis force due to the motion of the bulk of the layer moving with the vortex is almost fully compensated by the Coriolis force on the motion of the exterior flow. Because the mass anomaly of a cyclone is negative the force and acceleration have opposite sign. The role of dipolar structures in steadily moving vortices is discussed, and it is shown that their overall structure is fixed by the steady westward motion of the mass anomaly. Furthermore, it is shown that reduced-gravity vortices are not advected with a background flow. The reason for this behavior is that the background flow changes the ambient vorticity gradient such that the vortex obtains an extra self-propagation term that exactly cancels the advection by the background flow. Last, it is shown that a vortex initially at rest will accelerate equatorward first, after which a westward motion is generated. This result is independent of the sign of the vortex.

Selection rules for rovibronic transitions in symmetric top molecules

A simple diagrammatic rule is presented for determining the rotational selection rules governing transitions between any pair of vibronic states in electric dipole spectra of symmetric top molecules. The rule is useful in cases where degenerate vibronic levels with first-order Coriolis splittings occur, because it gives immediately the selection rule for the (+l) and (-l) components in any degenerate state. The rule is also helpful in determining the symmetry species and the effective zeta constants in overtone and combination levels involving degenerate vibrations. Particular attention is devoted to the conventions concerning the signs of zeta constants.

Raman selection rules for vibration-rotation transitions in symmetric top molecules

Symmetry restrictions on Raman selection rules can be obtained, quite generally, by considering a Raman allowed transition as the result of two successive dipole allowed transitions, and imposing the usual symmetry restrictions on the dipole transitions. This leads to the same results as the more familiar polarizability theory, but the vibration-rotation selection rules are easier to obtain by this argument. The selection rules for symmetric top molecules involving the (+l) and (-l) components of a degenerate vibrational level with first-order Coriolis splitting are derived in this paper. It is shown that these selection rules depend on the order of the highest-fold symmetry axis Cn, being different for molecules with n=3, n=4, or n ≧ 5; moreover the selection rules are different again for molecules belonging to the point groups Dnd with n even, and Sm with 1/2m even, for which the highest-fold symmetry axes Cn and Sm are related by m=2n. Finally it is shown that an apparent anomaly between the observed Raman and infra-red vibration-rotation spectra of the allene molecule is resolved when the correct selection rules are used, and a value for the A rotational constant of allene is derived without making use of the zeta sum rule.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part II: Rockfall activity and talus deformation

During many lava dome-forming eruptions, persistent rockfalls and the concurrent development of a substantial talus apron around the foot of the dome are important aspects of the observed activity. An improved understanding of internal dome structure, including the shape and internal boundaries of the talus apron, is critical for determining when a lava dome is poised for a major collapse and how this collapse might ensue. We consider a period of lava dome growth at the Soufrière Hills Volcano, Montserrat, from August 2005 to May 2006, during which a 100 × 106 m3 lava dome developed that culminated in a major dome-collapse event on 20 May 2006. We use an axi-symmetrical Finite Element Method model to simulate the growth and evolution of the lava dome, including the development of the talus apron. We first test the generic behaviour of this continuum model, which has core lava and carapace/talus components. Our model describes the generation rate of talus, including its spatial and temporal variation, as well as its post-generation deformation, which is important for an improved understanding of the internal configuration and structure of the dome. We then use our model to simulate the 2005 to 2006 Soufrière Hills dome growth using measured dome volumes and extrusion rates to drive the model and generate the evolving configuration of the dome core and carapace/talus domains. The evolution of the model is compared with the observed rockfall seismicity using event counts and seismic energy parameters, which are used here as a measure of rockfall intensity and hence a first-order proxy for volumes. The range of model-derived volume increments of talus aggraded to the talus slope per recorded rockfall event, approximately 3 × 103–13 × 103 m3 per rockfall, is high with respect to estimates based on observed events. From this, it is inferred that some of the volumetric growth of the talus apron (perhaps up to 60–70%) might have occurred in the form of aseismic deformation of the talus, forced by an internal, laterally spreading core. Talus apron growth by this mechanism has not previously been identified, and this suggests that the core, hosting hot gas-rich lava, could have a greater lateral extent than previously considered.

Quartic anharmonic resonances in acetylene molecules

Formulas are derived for the quartic anharmonic resonance coefficients observed to be important between C–H stretching and the combination of one quantum of C≡C stretching and two quanta of H–C≡C bending in a number of acetylene molecules. Examples of this resonance are ν3 with ν2+ν4+ν5 in 12C2H2, ν1 with ν2+2ν5 in 13C2H2, and ν1 with ν2+2ν4 in monofluoroacetylene and monochloroacetylene. The coefficients characterizing the resonances in these examples, which we denote K3,245, K1,255, and K1,244, arise from cubic and quartic terms in the anharmonic force field, in the normal coordinate representation, through second order and first order perturbation treatments respectively, where the second order resonances are calculated by a Van Vleck resonance formalism. The experimentally determined values of these coefficients are compared with values calculated from model anharmonic force fields.

Arthropod community structure in pastures of an island archipelago (Azores): looking for local-regional species richness patterns at fine-scales

The arthropod species richness of pastures in three Azorean islands was used to examine the relationship between local and regional species richness over two years. Two groups of arthropods, spiders and sucking insects, representing two functionally different but common groups of pasture invertebrates were investigated. The local-regional species richness relationship was assessed over relatively fine scales: quadrats (= local scale) and within pastures (= regional scale). Mean plot species richness was used as a measure of local species richness (= alpha diversity) and regional species richness was estimated at the pasture level (= gamma diversity) with the 'first-order-Jackknife' estimator. Three related issues were addressed: (i) the role of estimated regional species richness and variables operating at the local scale (vegetation structure and diversity) in determining local species richness; (ii) quantification of the relative contributions of alpha and beta diversity to regional diversity using additive partitioning; and (iii) the occurrence of consistent patterns in different years by analysing independently between-year data. Species assemblages of spiders were saturated at the local scale (similar local species richness and increasing beta-diversity in richer regions) and were more dependent on vegetational structure than regional species richness. Sucking insect herbivores, by contrast, exhibited a linear relationship between local and regional species richness, consistent with the proportional sampling model. The patterns were consistent between years. These results imply that for spiders local processes are important, with assemblages in a particular patch being constrained by habitat structure. In contrast, for sucking insects, local processes may be insignificant in structuring communities.

Effect of breed, gender, housing system and dietary crude protein content on performance of finishing beef cattle fed maize-silage-based diets

Maize silage-based diets with three dietary crude protein (CP) supplements were offered to 96 finishing cattle of contrasting breed (Holstein Friesian (HF) v. Simmental x HF (SHF)) and gender (bull v. steer) housed in two types of feeding system (group fed v. individually fed). The three protein supplements differed either in CP or protein degradability (degradable (LUDP) v. rumen undegradable (HUDP)) and provided CP concentrations of 142 (Con), 175 (LUDP) and 179 (HUDP) g/kg dry matter (DM) respectively, with ratios of degradable to undegradable of 3.0, 1.4 and 0.9:1 for diets Con, LOP and HUDP respectively. DM intakes were marginally higher (P = 0. 102) for LOP when compared with Con and HOP Rates of daily live-weight gain (DLWG) were higher (P = 0.005) in LUDP and HOP when compared with Con. HF had higher DM intakes than SHF although this did not result in any improvement in HF DLWG. Bulls had significantly better DM intakes, DLWG and feed conversion efficiency than steers. Conformation scores were better in SHF than HF (P < 0.001) and fat scores lower in bulls than steers (p < 0.001). There was a number of first order interactions established between dietary treatment, breed, gender and housing system with respect to rates of gain and carcass fat scores.

Large supersaturated designs

A supersaturated design (SSD) is an experimental plan, useful for evaluating the main effects of m factors with n experimental units when m > n - 1, each factor has two levels and when the first-order effects of only a few factors are expected to have dominant effects on the response. Use of these plans can be extremely cost-effective when it is necessary to screen hundreds or thousands of factors with a limited amount of resources. In this article we describe how to use cyclic balanced incomplete block designs and regular graph designs to construct E (s(2)) optimal and near optimal SSDs when m is a multiple of n - 1. We also provide a table that can be used to construct these designs for screening thousands of factors. We also explain how to obtain SSDs when m is not a multiple of n - 1. Using the table and the approaches given in this paper, SSDs can be developed for designs with up to 24 runs and up to 12,190 factors.