Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

A modified discrete random pore model allowing for different initial surface reactivity

We investigate here a modification of the discrete random pore model [Bhatia SK, Vartak BJ, Carbon 1996;34:1383], by including an additional rate constant which takes into account the different reactivity of the initial pore surface having attached functional groups and hydrogens, relative to the subsequently exposed surface. It is observed that the relative initial reactivity has a significant effect on the conversion and structural evolution, underscoring the importance of initial surface chemistry. The model is tested against experimental data on chemically controlled char oxidation and steam gasification at various temperatures. It is seen that the variations of the reaction rate and surface area with conversion are better represented by the present approach than earlier random pore models. The results clearly indicate the improvement of model predictions in the low conversion region, where the effect of the initially attached functional groups and hydrogens is more significant, particularly for char oxidation. It is also seen that, for the data examined, the initial surface chemistry is less important for steam gasification as compared to the oxidation reaction. Further development of the approach must also incorporate the dynamics of surface complexation, which is not considered here.

Random numbers from astronomical imaging

This article describes a method to turn astronomical imaging into a random number generator by using the positions of incident cosmic rays and hot pixels to generate bit streams. We subject the resultant bit streams to a battery of standard benchmark statistical tests for randomness and show that these bit streams are statistically the same as a perfect random bit stream. Strategies for improving and building upon this method are outlined.

On the applicability of the Aris' equivalent sphere in multicomponent reacting systems

The validity of the concept of equivalent sphere introduced by Aris in 1957 to multicomponent reacting systems is investigated in this paper. A network of C6 hydrocarbon reforming reaction and a fixed bed reactor are taken as the model reaction network and the reactor configuration, respectively.

Isolation of genotype- and chromosome-specific DNA markers in a biotype of Colletotrichum gloeosporioides using random amplified polymorphic DNA analysis

Genetic markers that distinguish fungal genotypes are important tools for genetic analysis of heterokaryosis and parasexual recombination in fungi. Random amplified polymorphic DNA (RAPD) markers that distinguish two races of biotype B of Colletotrichum gloeosporioides infecting the legume Stylosanthes guianensis were sought. Eighty-five arbitrary oligonucleotide primers were used to generate 895 RAPD bands but only two bands were found to be specifically amplified from DNA of the race 3 isolate. These two RAPD bands were used as DNA probes and hybridised only to DNA of the race 3 isolate. Both RAPD bands hybridised to a dispensable 1.2 Mb chromosome of the race 3 isolate. No other genotype-specific chromosomes or DNA sequences were identified in either the race 2 or race 3 isolates. The RAPD markers hybridised to a 2 Mb chromosome in all races of the genetically distinct biotype A pathogen which infects other species of Stylosanthes as well as S. guianensis. The experiments indicate that RAPD analysis is a potentially useful tool for obtaining genotype-and chromosome-specific DNA probes in closely related isolates of one biotype of this fungal pathogen.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Availability of two forms of dissolved nitrogen to the coral Pocillopora damicornis and its symbiotic zooxanthellae

The relative contribution of dissolved nitrogen (ammonium and dissolved free amino acids DFAAs) to the nitrogen budget of the reef-building coral Pocillopora damicornis was assessed for colonies growing on control and ammonium-enriched reefs at One Tree Island (southern Great Barrier Reef) during the ENCORE (Enrichment of Nutrient on Coral Reef; 1993 to 1996) project. P. damicornis acquired ammonium at rates of between 5.1 and 91.8 nmol N cm(-2) h(-1) which were not affected by nutrient treatment except in the case of one morph. In this case, uptake rates decreased from 80.5 to 42.8 nmol cm(-2) h(-1) (P < 0.05) on exposure to elevated ammonium over 12 mo. The presence or absence of light during measurement did not influence the uptake of ammonium ions. Nitrogen budgets revealed that the uptake of ammonium from concentrations of 0.11 to 0.13 mu M could completely satisfy the demand of growing P. damicornis for new nitrogen. P. damicornis also took up DFAAs at rates ranging from 4.9 to 9.8 nmol N cm(-2) h(-1). These rates were higher in the dark than in the light (9.0 vs 5.1 nmol m(-2) h(-1), P < 0.001). Uptake rates were highest for the amino acids serine, arginine and alanine, and lowest for tyrosine. DFAA concentrations within the ENCORE microatolls that received ammonium were undetectable, whereas they ranged up to 100 nM within the control microatolls. The contribution of DFAAs to the nitrogen budget of P. damicornis constituted only a small fraction of the nitrogen potentially contributed by ammonium under field conditions. Even at the highest field concentrations measured during this study, DFAAs could contribute only similar or equal to 11.3% of the nitrogen demand of P. damicornis. This contribution, however, may be an important source of nitrogen when other sources such as ammonium are scarce or during periods when high concentrations of DFAAs become sporadically available (e.g. cell breakage during fish-grazing).

The N-15 natural abundance (delta N-15) of ecosystem samples reflects measures of water availability

We assembled a globally-derived data set for site-averaged foliar delta(15)N, the delta(15)N of whole surface mineral soil and corresponding site factors (mean annual rainfall and temperature, latitude, altitude and soil pH). The delta(15)N of whole soil was related to all of the site variables (including foliar delta(15)N) except altitude and, when regressed on latitude and rainfall, provided the best model of these data, accounting for 49% of the variation in whole soil delta(15)N. As single linear regressions, site-averaged foliar delta(15)N was more strongly related to rainfall than was whole soil delta(15)N. A smaller data set showed similar, negative correlations between whole soil delta(15)N, site-averaged foliar delta(15)N and soil moisture variations during a single growing season. The negative correlation between water availability (measured here by rainfall and temperature) and soil or plant delta(15)N fails at the landscape scale, where wet spots are delta(15)N-enriched relative to their drier surroundings. Here we present global and seasonal data, postulate a proximate mechanism for the overall relationship between water availability and ecosystem delta(15)N and, newly, a mechanism accounting for the highly delta(15)N-depleted values found in the foliage and soils of many wet/cold ecosystems. These hypotheses are complemented by documentation of the present gaps in knowledge, suggesting lines of research which will provide new insights into terrestrial N-cycling. Our conclusions are consistent with those of Austin and Vitousek (1998) that foliar (and soil) delta(15)N appear to be related to the residence time of whole ecosystem N.

Life-cycle preferences over consumption and health: when is cost-effectiveness analysis equivalent to cost-benefit analysis?

This paper studies life-cycle preferences over consumption and health status. We show that cost-effectiveness analysis is consistent with cost-benefit analysis if the Lifetime utility function is additive over time, multiplicative in the utility of consumption and the utility of health status, and if the utility of consumption is constant over time. We derive the conditions under which the lifetime utility function takes this form, both under expected utility theory and under rank-dependent utility theory, which is currently the most important nonexpected utility theory. If cost-effectiveness analysis is consistent with cost-benefit analysis, it is possible to derive tractable expressions for the willingness to pay for quality-adjusted life-years (QALYs). The willingness to pay for QALYs depends on wealth, remaining life expectancy, health status, and the possibilities for intertemporal substitution of consumption. (C) 1999 Elsevier Science B.V. All rights reserved.

Effects of mate availability on female longevity, fecundity and egg development of Homichloda barkeri (Jacoby) (Coleoptera : Chrysomelidae)

This study examined the effect of the number of maternal matings on egg production, embryonic development and female longevity in Homichloda barkeri (Jacoby) (Cole optera: Chrysomelidae). Single-mated females lived longer than multiple-mated and unmated females. The number of eggs produced per day and the proportion of eggs that developed between single- and multiple-mated females was not affected by mating frequency.

Neuromuscular-skeletal constraints upon the dynamics of unimanual and bimanual coordination

In the first of three experiments, 11 participants generated pronation and supination movements of the forearm, in time with an auditory metronome. The metronome frequency was increased in eight steps (0.25 Hz) from a base frequency of 1.75 Hz. On alternating trials, participants were required to coordinate either maximum pronation or maximum supination with each beat of the metronome. In each block of trials, the axis of rotation was either coincident with the long axis of the forearm, above this axis, or below this axis. The stability of the pronate-on-the-beat pattern, as indexed by the number of pattern changes, and the time of onset of pattern change, was greatest when the axis of rotation of the movement was below the long axis of the forearm. In contrast, the stability of the supinate-on-the-beat pattern was greatest when the axis of rotation of the movement was above the long axis of the forearm. In a second experiment, we examined how changes in the position of the axis of rotation alter the activation patterns of muscles that contribute to pronation and supination of the forearm. Variations in the relative dominance of the pronation and supination phases of the movement cycle across conditions were accounted for primarily by changes in the activation profile of flexor carpi radialis (FCR) and extensor carpi radialis longus (ECR). In the Final experiment we examined how these constraints impact upon the stability of bimanual coordination. Thirty-two participants were assigned at random to one of four conditions, each of which combined an axis of rotation configuration (bottom or top) for each limb. The participants generated both inphase (both limbs pronating simultaneously, and supinating simultaneously) and antiphase (left limb pronating and right limb supinating simultaneously, and vice versa) patterns of coordination. When the position of the axis of rotation was equivalent for the left and the right limb, transitions from antiphase to inphase patterns of coordination were Frequently observed. In marked contrast, when the position of the axis of rotation for the left and right limb was contradistinct, transitions From inphase to antiphase patterns of coordination occurred. The results demonstrated that when movements are performed in an appropriate mechanical context, inphase patterns of coordination are less stable than antiphase patterns.

The fine structure of (v, 3) directed triple systems: v equivalent to 2 (mod 3)

The fine structure of a directed triple system of index lambda is the vector (c(1), c(2),...,C-lambda), where c(i) is the number of directed triples appearing precisely i times in the system. We determine necessary and sufficient conditions for a vector to be the fine structure of a directed triple system of index 3 for upsilon = 2 (mod 3).

Content and context of monocular regions determine perceived depth in random dot, unpaired background and phantom stereograms

Perceived depth was measured for three-types of stereograms with the colour/texture of half-occluded (monocular) regions either similar to or dissimilar to that of binocular regions or background. In a two-panel random dot stereogram the monocular region was filled with texture either similar or different to the far panel or left blank. In unpaired background stereograms the monocular region either matched the background or was different in colour or texture and in phantom stereograms the monocular region matched the partially occluded object or was a different colour or texture. In all three cases depth was considerably impaired when the monocular texture did not match either the background or the more distant surface. The content and context of monocular regions as well as their position are important in determining their role as occlusion cues and thus in three-dimensional layout. We compare coincidence and accidental view accounts of these effects. (C) 2002 Elsevier Science Ltd. All rights reserved.