Turbulence in Small Sub-Tropical Estuary with Semi-Diurnal Tides

In natural estuaries, contaminant transport is driven by the turbulent momentum mixing. The predictions of scalar dispersion can rarely be predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted at high frequency and continuously for up to 50 hours per investigation in a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such small estuarine systems with shallow water depths (less than 0.5 m at low tides), and a thorough post-processing technique was applied. The estuarine flow is always a fluctuating process. The bulk flow parameters fluctuated with periods comparable to tidal cycles and other large-scale processes. But turbulence properties depended upon the instantaneous local flow properties. They were little affected by the flow history, but their structure and temporal variability were influenced by a variety of mechanisms. This resulted in behaviour which deviated from that for equilibrium turbulent boundary layer induced by velocity shear only. A striking feature of the data sets is the large fluctuations in all turbulence characteristics during the tidal cycle. This feature was rarely documented, but an important difference between the data sets used in this study from earlier reported measurements is that the present data were collected continuously at high frequency during relatively long periods. The findings bring new lights in the fluctuating nature of momentum exchange coefficients and integral time and length scales. These turbulent properties should not be assumed constant.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Common multiples of complete graphs

A graph H is said to divide a graph G if there exists a set S of subgraphs of G, all isomorphic to H, such that the edge set of G is partitioned by the edge sets of the subgraphs in S. Thus, a graph G is a common multiple of two graphs if each of the two graphs divides G.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

Effects of viscous dissipation and boundary conditions on forced convection in a channel occupied by a saturated porous medium

Forced convection with viscous dissipation in a parallel plate channel filled by a saturated porous medium is investigated numerically. Three different viscous dissipation models are examined. Two different sets of wall conditions are considered: isothermal and isoflux. Analytical expressions are also presented for the asymptotic temperature profile and the asymptotic Nusselt number. With isothermal walls, the Brinkman number significantly influences the developing Nusselt number but not the asymptotic one. At constant wall heat flux, both the developing and the asymptotic Nusselt numbers are affected by the value of the Brinkman number. The Nusselt number is sensitive to the porous medium shape factor under all conditions considered.

Review of recent results on excursion set models

Many images consist of two or more 'phases', where a phase is a collection of homogeneous zones. For example, the phases may represent the presence of different sulphides in an ore sample. Frequently, these phases exhibit very little structure, though all connected components of a given phase may be similar in some sense. As a consequence, random set models are commonly used to model such images. The Boolean model and models derived from the Boolean model are often chosen. An alternative approach to modelling such images is to use the excursion sets of random fields to model each phase. In this paper, the properties of excursion sets will be firstly discussed in terms of modelling binary images. Ways of extending these models to multi-phase images will then be explored. A desirable feature of any model is to be able to fit it to data reasonably well. Different methods for fitting random set models based on excursion sets will be presented and some of the difficulties with these methods will be discussed.

Innovation as a Meta-Attribute for Graduate Engineers

This paper reviews the attitudes, skills and knowledge that engineering innovators should possess. It critically analyses and compares sets of graduate attributes from the USA, Australia and Malaysia in terms of which of these relate to the ability to innovate. Innovation can be described as an integrative, meta attribute that overarches most of the other graduate attributes. Due to the “graduate attribute paradox”, it is shown how meeting the stated attributes of graduates by industry does not necessarily satisfy the requirements of industry. It is argued that the culture of the engineering school is an important influence on fostering innovation in engineers.

Effects of viscous dissipation and boundary conditions on forced convection in a channel occupied by a saturated porous medium

Forced convection with viscous dissipation in a parallel plate channel filled by a saturated porous medium is investigated numerically. Three different viscous dissipation models are examined. Two different sets of wall conditions are considered: isothermal and isoflux. Analytical expressions are also presented for the asymptotic temperature profile and the asymptotic Nusselt number. With isothermal walls, the Brinkman number significantly influences the developing Nusselt number but not the asymptotic one. At constant wall heat flux, both the developing and the asymptotic Nusselt numbers are affected by the value of the Brinkman number. The Nusselt number is sensitive to the porous medium shape factor under all conditions considered.

Alkyl halides, super hydrogen production and the pathogenesis of pneumatosis cystoides coli

Background and aims-The colons of patients with pneumatosis cystoides coli produce excessive H-2. Exposure to alkyl halides could explain this. Six consecutive patients who had pneumatosis cystoides coli while taking chloral hydrate (1-5+ g/day) are reported. Patients 2 and 3 were investigated after they had ceased chloral hydrate treatment. One produced methane, the other did not. (Pneumatosis cystoides coli patients are non-methanogenic according to the literature.) Both had overnight fasting breath H-2 of less than 10 ppm. A literature review disclosed just one patient who was using chloral at the time of diagnosed pneumatosis cystoides coli, but an epidemic of the disease in workers exposed to trichloroethylene. Methods-(i) In vitro experiments with human faeces: chloral or closely related alkyl halides were added to anaerobic faecal cultures derived from four methane-producing and three non-methanogenic human subjects. H-2 and CH4 gases were measured. (ii) In vivo animal experiment: chloral hydrate was added to drinking water of four Wistar rats, and faecal HI compared with control rats. Results-Alkyl halides increased H-2 up to 900 times in methanogenic and 10 times in non-methanogenic faecal cultures. The K-i of chloral was 0.2 mM. Methanogenesis was inhibited in concert with the increase in net H-2. In the rat experiment, chloral hydrate increased H-2 10 times, but did not cause pneumatosis. Conclusions-Chloral and trichloroethylene are alkyl halides chemically similar to chloroform, a potent inhibitor of H-2 consumption by methanogens and acetogens. These bacteria are the most important H-2-consuming species in the colon. It is postulated that exposure to these alkyl halides increases net H-2 production, which sets the scene for counterperfusion supersaturation and the formation of gas cysts. In recent times, very low prescribing rates for chloral have caused primary pneumatosis cystoides to become extremely rare. As with primary pneumatosis, secondary pneumatosis cystoides, which occurs if there is small bowel bacterial overgrowth distal to a proximally located gut obstruction, is predicted by counterperfusion supersaturation. Inherent unsaturation due to metabolism of O-2 is a safety factor, which could explain why gas bubbles do not form more often in tissue with high H-2 tension.

Adsorption characteristics of carbonized bark for phenol and pentachlorophenol

The potential of using carbonized slash pine bark as a substitute for activated carbon was examined in this study. The bark was carbonized by slow heating in nitrogen for 6.5 h to 672 degrees C. The BET-N-2 surface area, average micropore and mesopore diameter, and micropore volume were 332 m(2) g(-1) 21.7 Angstrom, and 0.125 cm(3) g(-1), respectively. The adsorption capacities for phenol and pentachlorophenol (PCP) at pH 2 and pH 8 were evaluated. The Langmuir equation provided a slightly better fit than the Freundlich equation to two sets of phenol data. The calculated Freundlich constants, K = 0.41 - 0.58 mmol/g/(mmol dm(-3))(1/n) and 1/n = 0.30 - 0.41, were lower and higher, respectively, than literature values for activated carbons. The adsorption capacity of the carbonized bark was much lower for PCP than for phenol. The protonated and anionic PCP isotherms were Type II or III, respectively, in the Brunauer classification. The BET equation provided the best fit to protonated PCP isotherm data. The anionic PCP data were fitted to both the BET model and an equation used in the literature to represent phosphate adsorption on activated carbons. Nonlinear regression of the data for both phenol and PCP adsorption with the Freundlich, Langmuir and BET equations generally gave more accurate parameters, compared with the use of linearized equations to obtain the parameters. (C) 1998 SCI.

Evolutionary distinctiveness and status of the endangered Lake Eacham rainbowfish (Melanotaenia eachamensis)

The Lake Eacham rainbowfish (Melanotaenia eachamensis) was declared extinct in the wild in the late 1980s after it disappeared from its only known locality, an isolated crater lake in northeast Queensland. Doubts have been raised about whether this taxon is distinct from surrounding populations of the eastern rainbowfish (Melanotaenia splendida splendida). We examined the evolutionary distinctiveness of M. eachamensis, obtained from captive stocks, relative to M. s. splendida through analysis of variation in mtDNA sequences, nuclear microsatellites, and morphometric characters Captive M. eachamensis had mtDNAs that were highly divergent from those in most populations of M. s. splendida. A broader geographic survey using RFLPs revealed some populations initially identified as M. s. splendida, that carried eachamensis mtDNA, whereas some others had mixtures of eachamensis and splendida mtDNA. The presence of eachamensis-like mtDNA in these populations could in principle be due to (1) sorting of ancestral polymorphisms, (2) introgression of M. eachamensis mtDNA into M. s. splendida, or (3) incorrect species boundaries, such that some populations currently assigned to M. s. splendida are M. eachamensis or are mixtures of the two species. These alternatives hypotheses were evaluated through comparisons of four nuclear microsatellite loci and morphometrics and meristics. In analyses of both data sets, populations of M. s. splendida with eachamensis mtDNA were more similar to captive M. eachamensis than to M. s. splendida with splendida mtDNA, supporting hypothesis 3. These results are significant for the management of M. eachamensis in several respects. First the combined molecular and morphological evidence indicates that M. eachamensis is a distinct species and a discrete evolutionarily significant unit worthy of conservation effort. Second it appears that the species boundary between M. eachamensis and M. s. splendida has been misdiagnosed such that there are extant populations on the Atherton Tableland as well as areas where both forms coexist. Accordingly we suggest that M. eachamensis be listed as vulnerable, rather than critical (or extinct in the wild). Third, the discovery of extant but genetically divergent populations of M. eachamensis on the Atherton Tableland broadens the options for future reintroductions to Lake Eacham.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - II. Testing simulation capabilities of wheat growth

To simulate cropping systems, crop models must not only give reliable predictions of yield across a wide range of environmental conditions, they must also quantify water and nutrient use well, so that the status of the soil at maturity is a good representation of the starting conditions for the next cropping sequence. To assess the suitability for this task a range of crop models, currently used in Australia, were tested. The models differed in their design objectives, complexity and structure and were (i) tested on diverse, independent data sets from a wide range of environments and (ii) model components were further evaluated with one detailed data set from a semi-arid environment. All models were coded into the cropping systems shell APSIM, which provides a common soil water and nitrogen balance. Crop development was input, thus differences between simulations were caused entirely by difference in simulating crop growth. Under nitrogen non-limiting conditions between 73 and 85% of the observed kernel yield variation across environments was explained by the models. This ranged from 51 to 77% under varying nitrogen supply. Water and nitrogen effects on leaf area index were predicted poorly by all models resulting in erroneous predictions of dry matter accumulation and water use. When measured light interception was used as input, most models improved in their prediction of dry matter and yield. This test highlighted a range of compensating errors in all modelling approaches. Time course and final amount of water extraction was simulated well by two models, while others left up to 25% of potentially available soil water in the profile. Kernel nitrogen percentage was predicted poorly by all models due to its sensitivity to small dry matter changes. Yield and dry matter could be estimated adequately for a range of environmental conditions using the general concepts of radiation use efficiency and transpiration efficiency. However, leaf area and kernel nitrogen dynamics need to be improved to achieve better estimates of water and nitrogen use if such models are to be use to evaluate cropping systems. (C) 1998 Elsevier Science B.V.

The critical power function is dependent on the duration of the predictive exercise tests chosen

The linear relationship between work accomplished (W-lim) and time to exhaustion (t(lim)) can be described by the equation: W-lim = a + CP.t(lim). Critical power (CP) is the slope of this line and is thought to represent a maximum rate of ATP synthesis without exhaustion, presumably an inherent characteristic of the aerobic energy system. The present investigation determined whether the choice of predictive tests would elicit significant differences in the estimated CP. Ten female physical education students completed, in random order and on consecutive days, five art-out predictive tests at preselected constant-power outputs. Predictive tests were performed on an electrically-braked cycle ergometer and power loadings were individually chosen so as to induce fatigue within approximately 1-10 mins. CP was derived by fitting the linear W-lim-t(lim) regression and calculated three ways: 1) using the first, third and fifth W-lim-t(lim) coordinates (I-135), 2) using coordinates from the three highest power outputs (I-123; mean t(lim) = 68-193 s) and 3) using coordinates from the lowest power outputs (I-345; mean t(lim) = 193-485 s). Repeated measures ANOVA revealed that CPI123 (201.0 +/- 37.9W) > CPI135 (176.1 +/- 27.6W) > CPI345 (164.0 +/- 22.8W) (P < 0.05). When the three sets of data were used to fit the hyperbolic Power-t(lim) regression, statistically significant differences between each CP were also found (P < 0.05). The shorter the predictive trials, the greater the slope of the W-lim-t(lim) regression; possibly because of the greater influence of 'aerobic inertia' on these trials. This may explain why CP has failed to represent a maximal, sustainable work rate. The present findings suggest that if CP is to represent the highest power output that an individual can maintain for a very long time without fatigue then CP should be calculated over a range of predictive tests in which the influence of aerobic inertia is minimised.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Velocimetry and thermometry of supersonic flow around a cylindrical body

The supersonic flow around a cylindrical body has been studied using two optical techniques. For both sets of measurements, the cylinder was mounted from the side of the tunnel, allowing investigation of the bow shock region as well as in the wake. A new technique, laser-enhanced ionization flow tagging, was used for streamwise velocity determinations behind the body. From these measurements, it was found that the downstream velocity outside the wake was (1.90 +/- 0.06) km/s, whereas inside the wake the velocity was about 0-500 m/s in the upstream direction. Planar laser induced fluorescence of nitric oxide was employed for temperature determinations. It was established that the freestream temperature was (2120 +/- 100) K, decreasing to around (1550 +/- 400) K in the wake.