An historical perspective on the development of the thermodynamic equation of seawater - 2010

Oceanography is concerned with understanding the mechanisms controlling the movement of seawater and its contents. A fundamental tool in this process is the characterization of the thermophysical properties of seawater as functions of measured temperature and electrical conductivity, the latter used as a proxy for the concentration of dissolved matter in seawater. For many years a collection of algorithms denoted the Equation of State 1980 (EOS-80) has been the internationally accepted standard for calculating such properties. However, modern measurement technology now allows routine observations of temperature and electrical conductivity to be made to at least one order of magnitude more accurately than the uncertainty in this standard. Recently, a new standard has been developed, the Thermodynamical Equation of Seawater 2010 (TEOS-10). This new standard is thermodynamically consistent, valid over a wider range of temperature and salinity, and includes a mechanism to account for composition variations in seawater. Here we review the scientific development of this standard, and describe the literature involved in its development, which includes many of the articles in this special issue.

Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity

Background: Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host. Results: Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each pollinator species was most common in a different geographic region, emphasising the importance of wide geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in apparently similar or identical resource niches. Conclusion: Our study offers a striking example of extreme deviation from reciprocal partner specificity over the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host’s range. We show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the true structure of cryptic biodiversity underpinning intimate ecological interactions.

Conditioning of the weak-constraint variational data assimilation problem for numerical weather prediction

4-Dimensional Variational Data Assimilation (4DVAR) assimilates observations through the minimisation of a least-squares objective function, which is constrained by the model flow. We refer to 4DVAR as strong-constraint 4DVAR (sc4DVAR) in this thesis as it assumes the model is perfect. Relaxing this assumption gives rise to weak-constraint 4DVAR (wc4DVAR), leading to a different minimisation problem with more degrees of freedom. We consider two wc4DVAR formulations in this thesis, the model error formulation and state estimation formulation. The 4DVAR objective function is traditionally solved using gradient-based iterative methods. The principle method used in Numerical Weather Prediction today is the Gauss-Newton approach. This method introduces a linearised `inner-loop' objective function, which upon convergence, updates the solution of the non-linear `outer-loop' objective function. This requires many evaluations of the objective function and its gradient, which emphasises the importance of the Hessian. The eigenvalues and eigenvectors of the Hessian provide insight into the degree of convexity of the objective function, while also indicating the difficulty one may encounter while iterative solving 4DVAR. The condition number of the Hessian is an appropriate measure for the sensitivity of the problem to input data. The condition number can also indicate the rate of convergence and solution accuracy of the minimisation algorithm. This thesis investigates the sensitivity of the solution process minimising both wc4DVAR objective functions to the internal assimilation parameters composing the problem. We gain insight into these sensitivities by bounding the condition number of the Hessians of both objective functions. We also precondition the model error objective function and show improved convergence. We show that both formulations' sensitivities are related to error variance balance, assimilation window length and correlation length-scales using the bounds. We further demonstrate this through numerical experiments on the condition number and data assimilation experiments using linear and non-linear chaotic toy models.

Generalized convective quasi-equilibrium principle

A generalization of Arakawa and Schubert's convective quasi-equilibrium principle is presented for a closure formulation of mass-flux convection parameterization. The original principle is based on the budget of the cloud work function. This principle is generalized by considering the budget for a vertical integral of an arbitrary convection-related quantity. The closure formulation includes Arakawa and Schubert's quasi-equilibrium, as well as both CAPE and moisture closures as special cases. The formulation also includes new possibilities for considering vertical integrals that are dependent on convective-scale variables, such as the moisture within convection. The generalized convective quasi-equilibrium is defined by a balance between large-scale forcing and convective response for a given vertically-integrated quantity. The latter takes the form of a convolution of a kernel matrix and a mass-flux spectrum, as in the original convective quasi-equilibrium. The kernel reduces to a scalar when either a bulk formulation is adopted, or only large-scale variables are considered within the vertical integral. Various physical implications of the generalized closure are discussed. These include the possibility that precipitation might be considered as a potentially-significant contribution to the large-scale forcing. Two dicta are proposed as guiding physical principles for the specifying a suitable vertically-integrated quantity.

Conditioning and preconditioning of the optimal state estimation problem

Optimal state estimation is a method that requires minimising a weighted, nonlinear, least-squares objective function in order to obtain the best estimate of the current state of a dynamical system. Often the minimisation is non-trivial due to the large scale of the problem, the relative sparsity of the observations and the nonlinearity of the objective function. To simplify the problem the solution is often found via a sequence of linearised objective functions. The condition number of the Hessian of the linearised problem is an important indicator of the convergence rate of the minimisation and the expected accuracy of the solution. In the standard formulation the convergence is slow, indicating an ill-conditioned objective function. A transformation to different variables is often used to ameliorate the conditioning of the Hessian by changing, or preconditioning, the Hessian. There is only sparse information in the literature for describing the causes of ill-conditioning of the optimal state estimation problem and explaining the effect of preconditioning on the condition number. This paper derives descriptive theoretical bounds on the condition number of both the unpreconditioned and preconditioned system in order to better understand the conditioning of the problem. We use these bounds to explain why the standard objective function is often ill-conditioned and why a standard preconditioning reduces the condition number. We also use the bounds on the preconditioned Hessian to understand the main factors that affect the conditioning of the system. We illustrate the results with simple numerical experiments.

Proanthocyanidins inhibit Ascaris suum glutathione-S-transferase activity and increase susceptibility of larvae to levamisole and ivermectin in vitro

Proanthocyanidins (PAC) are a class of plant secondary metabolites commonly found in the diet that have shown potential to control gastrointestinal nematode infections. The anti-parasitic mechanism(s) of PAC remain obscure, however the protein-binding properties of PAC suggest that disturbance of key enzyme functions may be a potential mode of action. Glutathione-S-transferases (GSTs) are essential for parasite detoxification and have been investigated as drug and vaccine targets. Here, we show that purified PAC strongly inhibit the activity of both recombinant and native GSTs from the parasitic nematode Ascaris suum. As GSTs are involved in detoxifying xenobiotic substances within the parasite, we hypothesised that this inhibition may render parasites hyper-susceptible to anthelmintic drugs. Migration inhibition assays with A. suum larvae demonstrated that the potency of levamisole (LEV) and ivermectin (IVM) were significantly increased in the presence of PAC purified from pine bark (4.6-fold and 3.2-fold reduction in IC50 value for LEV and IVM, respectively). Synergy analysis revealed that the relationship between PAC and LEV appeared to be synergistic in nature, suggesting a specific enhancement of LEV activity, whilst the relationship between PAC and IVM was additive rather than synergistic, suggesting independent actions. Our results demonstrate that these common dietary compounds may increase the efficacy of synthetic anthelmintic drugs in vitro, and also suggest one possible mechanism for their well-known anti-parasitic activity.

The reactions of specific neuron types to intestinal ischemia in the guinea pig enteric nervous system

Damage following ischemia and reperfusion (I/R) is common in the intestine and can be caused during abdominal surgery, in several disease states and following intestinal transplantation. Most studies have concentrated on damage to the mucosa, although published evidence also points to effects on neurons. Moreover, alterations of neuronally controlled functions of the intestine persist after I/R. The present study was designed to investigate the time course of damage to neurons and the selectivity of the effect of I/R damage for specific types of enteric neurons. A branch of the superior mesenteric artery supplying the distal ileum of anesthetised guinea pigs was occluded for 1 h and the animals were allowed to recover for 2 h to 4 weeks before tissue was taken for the immunohistochemical localization of markers of specific neuron types in tissues from sham and I/R animals. The dendrites of neurons with nitric oxide synthase (NOS) immunoreactivity, which are inhibitory motor neurons and interneurons, were distorted and swollen by 24 h after I/R and remained enlarged up to 28 days. The total neuron profile areas (cell body plus dendrites) increased by 25%, but the sizes of cell bodies did not change significantly. Neurons of type II morphology (intrinsic primary afferent neurons), revealed by NeuN immunoreactivity, were transiently reduced in cell size, at 24 h and 7 days. These neurons also showed signs of minor cell surface blebbing. Calretinin neurons, many of which are excitatory motor neurons, were unaffected. Thus, this study revealed a selective damage to NOS neurons that was observed at 24 h and persisted up to 4 weeks, without a significant change in the relative numbers of NOS neurons.

Impaired dendritic cell differentiation and maturation in the absence of C3

Human monocytes can be differentiated into immature dendritic cells (DCs) in the presence of serum and cytokines. One of the main functions of immature DCs is to capture and process antigens. Following maturation, they differentiate into antigen presenting cells. The role of complement in the differentiation process from monocytes towards immature DCs remains elusive. Here we demonstrate that complement 3 (C3) has a regulatory impact on the expression of specific DC surface molecules and DC-derived cytokine production during DC differentiation. We isolated human adherent peripheral blood mononuclear cells, which were cultured in the presence of GM-CSF plus IL-4 in medium supplemented with normal human serum or C3 deficient serum. The lack of C3 during DC differentiation negatively impacted the expression of C-type lectin receptor DC-SIGN, the antigen presenting molecules HLA-DR and CD1a, and the costimulatory molecules CD80 and CD86. Further, the spontaneous production of IL-6 and IL-12 was reduced in the absence of C3. Moreover, the maturation of immature DCs in response to LPS challenge was impaired in the absence of C3 as evidenced by reduced MHC-II, co-stimulatory molecule expression as well as modulated IL-12 and TNF-alpha production. Collectively, our results provide evidence for a novel role of C3 as a critical cofactor in human DC differentiation and maturation. (C) 2007 Elsevier Ltd. All rights reserved.

Immune response to cytolethal distending toxin of Aggregatibacter actinomycetemcomitans in periodontitis patients

Background and Objective: Cytolethal distending toxin (CDT) is a genotoxin produced by Aggregatibacter actinomycetemcomitans. In spite of its association with pathogenesis, little is known about the humoral immune response against the CDT. This study aimed to test whether subgingival colonization and humoral response to A. actinomycetemcomitans would lead to a response against CDT. Material and Methods: Sera from periodontally healthy, localized and generalized aggressive periodontitis and chronic periodontitis subjects (n = 80) were assessed for immunoglobulin G titers to A. actinomycetemcomitans serotypes a/b/c and to each CDT subunit (CdtA, CdtB and CdtC) by ELISA. A. actinomycetemcomitans subgingival levels and neutralization of CDT activity were also analyzed. Results: Sera from 75.0% localized and 81.8% generalized aggressive periodontitis patients reacted to A. actinomycetemcomitans. A response to serotype b was detected in localized (66.7%) and generalized aggressive periodontitis (54.5%). Reactivity to A. actinomycetemcomitans correlated with subgingival colonization (R = 0.75, p < 0.05). There was no correlation between A. actinomycetemcomitans colonization or response to serotypes and the immunoglobulin G response to CDT subunits. Titers of immunoglobulin G to CdtA and CdtB did not differ among groups; however, sera of all generalized aggressive periodontitis patients reacted to CdtC. Neutralization of CDT was not correlated with levels of antibodies to CDT subunits. Conclusion: Response to CdtA and CdtB did not correlate with the periodontal status of the subject in the context of an A. actinomycetemcomitans infection. However, a response to CdtC was found in sera of generalized but not of localized aggressive periodontitis subjects. Differences in response to CdtC between generalized and localized aggressive periodontitis subjects indicate that CDT could be expressed differently by the infecting strains. Alternatively, the antibody response to CdtC could require the colonization of multiple sites.

Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

C(0) and bi-Lipschitz K-equivalence of mappings

In this paper we investigate the classification of mappings up to K-equivalence. We give several results of this type. We study semialgebraic deformations up to semialgebraic C(0) K-equivalence and bi-Lipschitz K-equivalence. We give an algebraic criterion for bi-Lipschitz K-triviality in terms of semi-integral closure (Theorem 3.5). We also give a new proof of a result of Nishimura: we show that two germs of smooth mappings f, g : R(n) -> R(n), finitely determined with respect to K-equivalence are C(0)-K-equivalent if and only if they have the same degree in absolute value.

Dynamics of a class of ODEs more general than almost periodic

A temporally global solution, if it exists, of a nonautonomous ordinary differential equation need not be periodic, almost periodic or almost automorphic when the forcing term is periodic, almost periodic or almost automorphic, respectively. An alternative class of functions extending periodic and almost periodic functions which has the property that a bounded temporally global solution solution of a nonautonomous ordinary differential equation belongs to this class when the forcing term does is introduced here. Specifically, the class of functions consists of uniformly continuous functions, defined on the real line and taking values in a Banach space, which have pre-compact ranges. Besides periodic and almost periodic functions, this class also includes many nonrecurrent functions. Assuming a hyperbolic structure for the unperturbed linear equation and certain properties for the linear and nonlinear parts, the existence of a special bounded entire solution, as well the existence of stable and unstable manifolds of this solution are established. Moreover, it is shown that this solution and these manifolds inherit the temporal behaviour of the vector field equation. In the stable case it is shown that this special solution is the pullback attractor of the system. A class of infinite dimensional examples involving a linear operator consisting of a time independent part which generates a C(0)-semigroup plus a small time dependent part is presented and applied to systems of coupled heat and beam equations. (C) 2010 Elsevier Ltd. All rights reserved.

Transformed generalized linear models

The estimation of data transformation is very useful to yield response variables satisfying closely a normal linear model, Generalized linear models enable the fitting of models to a wide range of data types. These models are based on exponential dispersion models. We propose a new class of transformed generalized linear models to extend the Box and Cox models and the generalized linear models. We use the generalized linear model framework to fit these models and discuss maximum likelihood estimation and inference. We give a simple formula to estimate the parameter that index the transformation of the response variable for a subclass of models. We also give a simple formula to estimate the rth moment of the original dependent variable. We explore the possibility of using these models to time series data to extend the generalized autoregressive moving average models discussed by Benjamin er al. [Generalized autoregressive moving average models. J. Amer. Statist. Assoc. 98, 214-223]. The usefulness of these models is illustrated in a Simulation study and in applications to three real data sets. (C) 2009 Elsevier B.V. All rights reserved.

Power series generalized nonlinear models

We introduce in this paper a new class of discrete generalized nonlinear models to extend the binomial, Poisson and negative binomial models to cope with count data. This class of models includes some important models such as log-nonlinear models, logit, probit and negative binomial nonlinear models, generalized Poisson and generalized negative binomial regression models, among other models, which enables the fitting of a wide range of models to count data. We derive an iterative process for fitting these models by maximum likelihood and discuss inference on the parameters. The usefulness of the new class of models is illustrated with an application to a real data set. (C) 2008 Elsevier B.V. All rights reserved.

Statistical analysis of proficiency testing results under elliptical distributions

The use of inter-laboratory test comparisons to determine the performance of individual laboratories for specific tests (or for calibration) [ISO/IEC Guide 43-1, 1997. Proficiency testing by interlaboratory comparisons - Part 1: Development and operation of proficiency testing schemes] is called Proficiency Testing (PT). In this paper we propose the use of the generalized likelihood ratio test to compare the performance of the group of laboratories for specific tests relative to the assigned value and illustrate the procedure considering an actual data from the PT program in the area of volume. The proposed test extends the test criteria in use allowing to test for the consistency of the group of laboratories. Moreover, the class of elliptical distributions are considered for the obtained measurements. (C) 2008 Elsevier B.V. All rights reserved.