A mathematical model for estimating the extent of solute- and water-flux heterogeneity in multiple sample percolation experiments

Multiple sampling is widely used in vadose zone percolation experiments to investigate the extent in which soil structure heterogeneities influence the spatial and temporal distributions of water and solutes. In this note, a simple, robust, mathematical model, based on the beta-statistical distribution, is proposed as a method of quantifying the magnitude of heterogeneity in such experiments. The model relies on fitting two parameters, alpha and zeta to the cumulative elution curves generated in multiple-sample percolation experiments. The model does not require knowledge of the soil structure. A homogeneous or uniform distribution of a solute and/or soil-water is indicated by alpha = zeta = 1, Using these parameters, a heterogeneity index (HI) is defined as root 3 times the ratio of the standard deviation and mean. Uniform or homogeneous flow of water or solutes is indicated by HI = 1 and heterogeneity is indicated by HI > 1. A large value for this index may indicate preferential flow. The heterogeneity index relies only on knowledge of the elution curves generated from multiple sample percolation experiments and is, therefore, easily calculated. The index may also be used to describe and compare the differences in solute and soil-water percolation from different experiments. The use of this index is discussed for several different leaching experiments. (C) 1999 Elsevier Science B.V. All rights reserved.

Thermal and radiation curing of phenylethynyl terminated macromers

The thermal and gamma-irradiation induced curing of two phenylethynyl terminated composite resin systems, DFB/BPF and PETI5A, was investigated. Thermal curing of these matrix resin samples was performed at a temperature of 360 degrees C, gamma irradiation of the samples was conducted at 300 degrees C at a dose rate of 2.2 kGy h(-1). The reaction and subsequent loss of ethynyl groups in the resins for both cure methods was demonstrated by observing the decrease of the 2215 cm(-1) peak in the Raman spectra of the resins. Fully cured resin samples were found to have glass transition temperatures of 244-246 degrees C and 278-280 degrees C for DFB/BPF and PETI5A respectively. Similar relationships between T-g and fractional conversion were observed in both resins. The apparent polymerization rate, R-p, for thermal cure at 360 degrees C, was found to be 4.79 x 10(-2)% s(-1) in PETI5A and 3.22 x 10(-2)% s(-1) in DFB/BPF. Catastrophic degradation under nitrogen was observed to commence near 450 degrees C and 530 degrees C, with 5% weight losses occurring at 455 degrees C and 540 degrees C for DFB/BPF and PETI5A respectively. Gamma radiation induced cure at 300 degrees C was shown to be feasible, with full cure being reached with doses of 40 kGy for DFB/BPF and 100 kGy for PETI5A.

New methods for determining quasi-stationary distributions for Markov chains

We shall be concerned with the problem of determining quasi-stationary distributions for Markovian models directly from their transition rates Q. We shall present simple conditions for a mu-invariant measure m for Q to be mu-invariant for the transition function, so that if m is finite, it can be normalized to produce a quasi-stationary distribution. (C) 2000 Elsevier Science Ltd. All rights reserved.

Topological methods for some boundary value problems

We establish existence of solutions for a finite difference approximation to y = f(x, y, y ') on [0, 1], subject to nonlinear two-point Sturm-Liouville boundary conditions of the form g(i)(y(i),y ' (i)) = 0, i = 0, 1, assuming S satisfies one-sided growth bounds with respect to y '. (C) 2001 Elsevier Science Ltd. All rights reserved.

Application of Petrov-Galerkin methods to transient boundary value problems in chemical engineering: adsorption with steep gradients in bidisperse solids

Petrov-Galerkin methods are known to be versatile techniques for the solution of a wide variety of convection-dispersion transport problems, including those involving steep gradients. but have hitherto received little attention by chemical engineers. We illustrate the technique by means of the well-known problem of simultaneous diffusion and adsorption in a spherical sorbent pellet comprised of spherical, non-overlapping microparticles of uniform size and investigate the uptake dynamics. Solutions to adsorption problems exhibit steep gradients when macropore diffusion controls or micropore diffusion controls, and the application of classical numerical methods to such problems can present difficulties. In this paper, a semi-discrete Petrov-Galerkin finite element method for numerically solving adsorption problems with steep gradients in bidisperse solids is presented. The numerical solution was found to match the analytical solution when the adsorption isotherm is linear and the diffusivities are constant. Computed results for the Langmuir isotherm and non-constant diffusivity in microparticle are numerically evaluated for comparison with results of a fitted-mesh collocation method, which was proposed by Liu and Bhatia (Comput. Chem. Engng. 23 (1999) 933-943). The new method is simple, highly efficient, and well-suited to a variety of adsorption and desorption problems involving steep gradients. (C) 2001 Elsevier Science Ltd. All rights reserved.

Purchasing Power Parities for Industry Comparisons Using Weighted Elteto-Koves-Szulc (EKS) Methods

This study has three main objectives. First, it develops a generalization of the commonly used EKS method to multilateral price comparisons. It is shown that the EKS system can be generalized so that weights can be attached to each of the link comparisons used in the EKS computations. These weights can account for differing levels of reliability of the underlying binary comparisons. Second, various reliability measures and corresponding weighting schemes are presented and their merits discussed. Finally, these new methods are applied to an international data set of manufacturing prices from the ICOP project. Although theoretically superior, it appears that the empirical impact of the weighted EKS method is generally small compared to the unweighted EKS. It is also found that this impact is larger when it is applied at lower levels of aggregation. Finally, the importance of using sector specific PPPs in assessing relative levels of manufacturing productivity is indicated.

Numerical methods for strong solutions of stochastic differential equations: an overview

This paper gives a review of recent progress in the design of numerical methods for computing the trajectories (sample paths) of solutions to stochastic differential equations. We give a brief survey of the area focusing on a number of application areas where approximations to strong solutions are important, with a particular focus on computational biology applications, and give the necessary analytical tools for understanding some of the important concepts associated with stochastic processes. We present the stochastic Taylor series expansion as the fundamental mechanism for constructing effective numerical methods, give general results that relate local and global order of convergence and mention the Magnus expansion as a mechanism for designing methods that preserve the underlying structure of the problem. We also present various classes of explicit and implicit methods for strong solutions, based on the underlying structure of the problem. Finally, we discuss implementation issues relating to maintaining the Brownian path, efficient simulation of stochastic integrals and variable-step-size implementations based on various types of control.

Expression and function of neuronal nicotinic ACh receptors in rat microvascular endothelial cells

The expression and function of nicotinic ACh receptors (nAChRs) in rat coronary microvascular endothelial cells (CMECs) were examined using RT-PCR and whole cell patch-clamp recording methods. RT-PCR revealed expression of mRNA encoding for the subunits alpha(2), alpha(3), alpha(4), alpha(5), alpha(7), beta(2), and beta(4) but not beta(3). Focal application of ACh evoked an inward current in isolated CMECs voltage clamped at negative membrane potentials. The current-voltage relationship of the ACh-induced current exhibited marked inward rectification and a reversal potential (E-rev) close to 0 mV. The cholinergic agonists nicotine, epibatidine, and cytisine activated membrane currents similar to those evoked by ACh. The nicotine-induced current was abolished by the neuronal nAChR antagonist mecamylamine. The direction and magnitude of the shift in E-rev of nicotine-induced current as a function of extracellular Na+ concentration indicate that the nAChR channel is cation selective and follows that predicted by the Goldman-Hodgkin-Katz equation assuming K+/Na+ permeability ratio of 1.11. In fura-2-loaded CMECs, application of ACh, but not of nicotine, elicited a transient increase in intracellular free Ca2+ concentration. Taken together, these results demonstrate that neuronal nAChR activation by cholinergic agonists evokes an inward current in CMECs carried primarily by Na+, which may contribute to the plasma nicotine-induced changes in microvascular permeability and reactivity induced by elevations in plasma nicotine.

The structure and infrastructure of Mexico's science and technology

The structure and infrastructure of the Mexican technical literature was determined. A representative database of technical articles was extracted from the Science Citation Index for the year 2002, with each article containing at least one author with a Mexican address. Many different manual and statistical clustering methods were used to identify the structure of the technical literature (especially the science and technology core competencies). One of the pervasive technical topics identified from the clustering, thin films research, was analyzed further using bibliometrics, in order to identify the infrastructure of this technology. Published by Elsevier Inc.

A Bayesian approach to imposing curvature on distance functions

The estimated parameters of output distance functions frequently violate the monotonicity, quasi-convexity and convexity constraints implied by economic theory, leading to estimated elasticities and shadow prices that are incorrectly signed, and ultimately to perverse conclusions concerning the effects of input and output changes on productivity growth and relative efficiency levels. We show how a Bayesian approach can be used to impose these constraints on the parameters of a translog output distance function. Implementing the approach involves the use of a Gibbs sampler with data augmentation. A Metropolis-Hastings algorithm is also used within the Gibbs to simulate observations from truncated pdfs. Our methods are developed for the case where panel data is available and technical inefficiency effects are assumed to be time-invariant. Two models-a fixed effects model and a random effects model-are developed and applied to panel data on 17 European railways. We observe significant changes in estimated elasticities and shadow price ratios when regularity restrictions are imposed. (c) 2004 Elsevier B.V. All rights reserved.

Estimating variable returns to scale production frontiers with alternative stochastic assumptions

Two stochastic production frontier models are formulated within the generalized production function framework popularized by Zellner and Revankar (Rev. Econ. Stud. 36 (1969) 241) and Zellner and Ryu (J. Appl. Econometrics 13 (1998) 101). This framework is convenient for parsimonious modeling of a production function with returns to scale specified as a function of output. Two alternatives for introducing the stochastic inefficiency term and the stochastic error are considered. In the first the errors are added to an equation of the form h(log y, theta) = log f (x, beta) where y denotes output, x is a vector of inputs and (theta, beta) are parameters. In the second the equation h(log y,theta) = log f(x, beta) is solved for log y to yield a solution of the form log y = g[theta, log f(x, beta)] and the errors are added to this equation. The latter alternative is novel, but it is needed to preserve the usual definition of firm efficiency. The two alternative stochastic assumptions are considered in conjunction with two returns to scale functions, making a total of four models that are considered. A Bayesian framework for estimating all four models is described. The techniques are applied to USDA state-level data on agricultural output and four inputs. Posterior distributions for all parameters, for firm efficiencies and for the efficiency rankings of firms are obtained. The sensitivity of the results to the returns to scale specification and to the stochastic specification is examined. (c) 2004 Elsevier B.V. All rights reserved.

Cholinergic and adrenergic innervation of lingual salt glands of the estuarine crocodile, Crocodylus porosus

Many marine reptiles and birds possess extrarenal salt glands that facilitate the excretion of excess sodium and chloride ions accumulated as a consequence of living in saline environments. Control of the secretory activity of avian salt glands is under neural control, but little information is available on the control of reptilian salt glands. Innervation of the lingual salt glands of the salt water crocodile, Crocodylus porosus, was examined in salt water-acclimated animals using histological methods. Extensive networks of both cholinergic and adrenergic nerve fibres were identified close to salt-secreting lobules and vasculature. The identification of both catecholamine-containing and cholinergic neurons in the salt gland epithelium and close to major blood vessels in the tissue suggests the action of the neurotransmitters on the salt-secreting epithelium itself and the rich vascular network of the lingual salt glands.

Age and growth in olive ridley seaturtles (Lepidochelys olivacea) from the north-central Pacific: a skeletochronological analysis

The olive ridley is the most abundant seaturtle species in the world but little is known of the demography of this species. We used skeletochronological data on humerus diameter growth changes to estimate the age of North Pacific olive ridley seaturtles caught incidentally by pelagic longline fisheries operating near Hawaii and from dead turtles washed ashore on the main Hawaiian Islands. Two age estimation methods [ranking, correction factor (CF)] were used and yielded age estimates ranging from 5 to 38 and 7 to 24 years, respectively. Rank age-estimates are highly correlated (r = 0.93) with straight carapace length (SCL), CF age estimates are not (r = 0.62). We consider the CF age-estimates as biologically more plausible because of the disassociation of age and size. Using the CF age-estimates, we then estimate the median age at sexual maturity to be around 13 years old (mean carapace size c. 60 cm SCL) and found that somatic growth was negligible by 15 years of age. The expected age-specific growth rate function derived using numerical differentiation suggests at least one juvenile growth spurt at about 10–12 years of age when maximum age-specific growth rates, c. 5 cm SCL year−1, are apparent.