Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.

Squeezing-induced complete transparency in two-level systems

The modification of the statistical properties of vacuum fluctuations, via quadrature squeezing, can dramatically reduce the absorptive and dispersive properties of two-level atoms. We show that for some range of parameter values the system exhibits zero absorption accompanied by zero dispersion of the probe field. This complete transparency is attributed to the coherent population oscillations induced by the squeezed vacuum.

Classical and quantum noise in electronic systems

We review the description of noise in electronic circuits in terms of electron transport. The Poisson process is used as a unifying principle. In recent years, much attention has been given to current noise in light-emitting diodes and laser diodes. In these devices, random events associated with electron transport are correlated with photon emission times, thus modifying both the current statistics and the statistics of the emitted light. We give a review of experiments in this area with special emphasis on the ability of such devices to produce subshot-noise currents and light beams. Finally we consider the noise properties of a class of mesoscopic devices based on the quantum tunnelling of an electron into and out of a bound state. We present a simple quantum model of this process which confirms that the current noise in such a device should be subshot-noise.

Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

5-cycle systems of K-v-F with a hole

Recently the problem of the existence of a 5-cycle system of K-v with a hole of size u was completely solved. In this paper we prove necessary and sufficient conditions on v and u for the existence of a 5-cycle system of K-v - F, with a hole of size u.

X-ray properties of the Parkes sample of flat-spectrum radio sources: dust in radio-loud quasars?

We investigate the X-ray properties of the Parkes sample of Bat-spectrum radio sources using data from the ROSAT All-Sky Survey and archival pointed PSPC observations. In total, 163 of the 323 sources are detected. For the remaining 160 sources, 2 sigma upper limits to the X-ray flux are derived. We present power-law photon indices in the 0.1-2.4 keV energy band for 115 sources, which were determined either with a hardness ratio technique or from direct fits to pointed PSPC data if a sufficient number of photons were available. The average photon index is <Gamma > = 1.95(-0.12)(+0.13) for flat-spectrum radio-loud quasars, <Gamma > = 1.70(-0.24)(+0.23) for galaxies, and <Gamma > = 2.40(-0.31)(+0.12) for BL Lac objects. The soft X-ray photon index is correlated with redshift and with radio spectral index in the sense that sources at high redshift and/or with flat (or inverted) radio spectra have flatter X-ray spectra on average. The results are in accord with orientation-dependent unification schemes for radio-loud active galactic nuclei. Webster et al. discovered many sources with unusually red optical continua among the quasars of this sample, and interpreted this result in terms of extinction by dust. Although the X-ray spectra in general do not show excess absorption, we find that low-redshift optically red quasars have significantly lower soft X-ray luminosities on average than objects with blue optical continua. The difference disappears for higher redshifts, as is expected for intrinsic absorption by cold gas associated with the dust. In addition, the scatter in log(f(x)/f(o)) is consistent with the observed optical extinction, contrary to previous claims based on optically or X-ray selected samples. Although alternative explanations for the red optical continua cannot be excluded with the present X-ray data, we note that the observed X-ray properties are consistent with the idea that dust plays an important role in some of the radio-loud quasars with red optical continua.

(m,n)-cycle systems

We describe a method which, in certain circumstances, may be used to prove that the well-known necessary conditions for partitioning the edge set of the complete graph on an odd number of vertices (or the complete graph on an even number of vertices with a 1-factor removed) into cycles of lengths m(1),m(2),...,m(t) are sufficient in the case \{m(1), m(2), ..., m(t)}\=2. The method is used to settle the case where the cycle lengths are 4 and 5. (C) 1998 Elsevier Science B.V. All rights reserved.

Binary vectors for sense and antisense expression of arabidopsis ESTs

Our laboratory is interested in devising methods to identify functions for the vast numbers of arabidopsis genes now available. For this purpose, we have constructed a set of binary vectors that will allow the quick production of transgenic arabidopsis plants containing either sense or antisense copies of EST clones obtained from the PRL2 library. These vectors are based on the pSLJ series containing the bialophos resistance (BAR) gene that confers resistance to the herbicide BASTA. Tn addition, our vectors contain a 35S CaMV promoter-polylinker-nos terminator cassette that allows the direct cloning of arabidopsis ESTs in either antisense (pAOV and pAOV2) or sense (pSOV and pSOV2) orientation. We also describe the construction of two additional vectors conferring BASTA resistance and containing the pBluescript polylinker in both orientations inserted between the 35S CaMV promoter and nos terminator (pKMB and pSMB).

The influence of Lys68 in decapeptide agonists of C5a on C5a receptor binding, activation and selectivity

The potent, conformationally biased C5a agonist peptide YSFKPMPLaR (C5a(65-74), Y65, F67, P69, P71, D-Ala73) was used as a template to gain insight into the nature and importance of lysine at position 68 in the peptide-receptor interaction. A panel of YSFKPMPLaR analogs with systematic substitutions for Lys68 was evaluated for C5a receptor (C5aR) binding affinity and activation in two well-characterized assay systems: human polymorphonuclear leukocytes (PMNs) and human fetal artery. In addition, we determined the activity of these new analogs in transfected rat basophilic leukemia (RBL) cells in which the Glu at position 199 of the C5aR (wtGlu199) was replaced by a Gin (C5aR-Gln199) or a Lys (C5aR-Lys199). Our results indicated that Lys68 in YSFKPMPLaR plays an important role in binding the C5aR expressed on PMNs and RBL cells. Furthermore, the data indicated that Lys68 interacted with Glu199 of the C5aR in PMNs and RBL cells. In human fetal artery, however, Lys68 substitutions had little or no effect on activity, which suggested that the receptor conformation may be different in this tissue. Thus, the interaction between Lys68 of the decapeptide agonist and Glu199 of the C5aR may be cell type-specific and may form the molecular basis for tissue-specific responses to C5a agonists.

Quantum open-systems approach to current noise in resonant tunneling junctions

A quantum Markovian master equation is derived to describe the current noise in resonant tunneling devices. This equation includes both incoherent and coherent quantum tunneling processes. We show how to obtain the population master equation by adiabatic elimination of quantum coherences in the presence of elastic scattering. We calculate the noise spectrum for a double well device and predict subshot noise statistics for strong tunneling between the wells. The method is an alternative to Green's function methods and population master equations for very small coherently coupled quantum dots.

New techniques in nutritional assessment: body composition methods

New techniques in air-displacement plethysmography seem to have overcome many of the previous problems of poor reproducibility and validity. These have made body-density measurements available to a larger range of individuals, including children, elderly and sick patients who often have difficulties in being submerged underwater in hydrodensitometry systems. The BOD POD air-displacement system (BOD POD body composition system; Life Measurement Instruments, Concord, CA, USA) is more precise than hydrodensitometry, is simple and rapid to operate (approximately 1 min measurements) and the results agree closely with those of hydrodensitometry (e.g. +/-3.4% for estimation of body fat). Body line scanners employing the principles of three-dimensional photography are potentially able to measure the surface area and volume of the body and its segments even more rapidly (approximately 10 s), but the validity of the measurements needs to be established. Advances in i.r. spectroscopy and mathematical modelling for calculating the area under the curve have improved precision for measuring enrichment of (H2O)-H-2 in studies of water dilution (CV 0.1-0.9% within the range of 400-1000 mu l/l) in saliva, plasma and urine. The technique is rapid and compares closely with mass spectrometry (bias 1 (SD 2) %). Advances in bedside bioelectrical-impedance techniques are making possible potential measurements of skinfold thicknesses and limb muscle mass electronically. Preliminary results suggest that the electronic method is more reproducible (intra-and inter-individual reproducibility for measuring skinfold thicknesses) and associated with less bias (+ 12%), than anthropometry (+ 40%). In addition to these selected examples, the 'mobility' or transfer of reference methods between centres has made the distinction between reference and bedside or field techniques less distinct than in the past.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.