QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

The structure of the optimal income tax in the quasi-linear model

Existing numerical characterizations of the optimal income tax have been based on a limited number of model specifications. As a result, they do not reveal which properties are general. We determine the optimal tax in the quasi-linear model under weaker assumptions than have previously been used; in particular, we remove the assumption of a lower bound on the utility of zero consumption and the need to permit negative labor incomes. A Monte Carlo analysis is then conducted in which economies are selected at random and the optimal tax function constructed. The results show that in a significant proportion of economies the marginal tax rate rises at low skills and falls at high. The average tax rate is equally likely to rise or fall with skill at low skill levels, rises in the majority of cases in the centre of the skill range, and falls at high skills. These results are consistent across all the specifications we test. We then extend the analysis to show that these results also hold for Cobb-Douglas utility.

Time evolution of the Wigner function in discrete quantum phase space for a soluble quasi-spin model

The discrete phase space approach to quantum mechanics of degrees of freedom without classical counterparts is applied to the many-fermions/quasi-spin Lipkin model. The Wi:ner function is written for some chosen states associated to discrete angle and angular momentum variables, and the rime evolution is numerically calculated using the discrete von Neumnnn-Liouville equation. Direct evidences in the lime evolution of the Wigner function are extracted that identify a tunnelling effect. A connection with a SU(2)-based semiclassical continuous approach to the Lipkin model is also presented.

Quantum spin tunneling in a soluble quasi-spin model: an angle-based potential description

We propose an approach which allows one to construct and use a potential function written in terms of an angle variable to describe interacting spin systems. We show how this can be implemented in the Lipkin-Meshkov-Glick, here considered a paradigmatic spin model. It is shown how some features of the energy gap can be interpreted in terms of a spin tunneling. A discrete Wigner function is constructed for a symmetric combination of two states of the model and its time evolution is obtained. The physical information extracted from that function reinforces our description of phase oscillations in a potential. (c) 2004 Elsevier B.V. All rights reserved.

Langevin dynamics, large deviations and instantons for the quasi-geostrophic model and two-dimensional Euler equations

We investigate a class of simple models for Langevin dynamics of turbulent flows, including the one-layer quasi-geostrophic equation and the two-dimensional Euler equations. Starting from a path integral representation of the transition probability, we compute the most probable fluctuation paths from one attractor to any state within its basin of attraction. We prove that such fluctuation paths are the time reversed trajectories of the relaxation paths for a corresponding dual dynamics, which are also within the framework of quasi-geostrophic Langevin dynamics. Cases with or without detailed balance are studied. We discuss a specific example for which the stationary measure displays either a second order (continuous) or a first order (discontinuous) phase transition and a tricritical point. In situations where a first order phase transition is observed, the dynamics are bistable. Then, the transition paths between two coexisting attractors are instantons (fluctuation paths from an attractor to a saddle), which are related to the relaxation paths of the corresponding dual dynamics. For this example, we show how one can analytically determine the instantons and compute the transition probabilities for rare transitions between two attractors.

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Effective carrying capacity and analytical solution of a particular case of the Richards-like two-species population dynamics model

We consider a generalized two-species population dynamic model and analytically solve it for the amensalism and commensalism ecological interactions. These two-species models can be simplified to a one-species model with a time dependent extrinsic growth factor. With a one-species model with an effective carrying capacity one is able to retrieve the steady state solutions of the previous one-species model. The equivalence obtained between the effective carrying capacity and the extrinsic growth factor is complete only for a particular case, the Gompertz model. Here we unveil important aspects of sigmoid growth curves, which are relevant to growth processes and population dynamics. (C) 2011 Elsevier B.V. All rights reserved.

Overcoming the rare species modelling paradox: test of a novel hierarchical framework with an Iberian endemic plant

Rare species have restricted geographic ranges, habitat specialization, and/or small population sizes. Datasets on rare species distribution usually have few observations, limited spatial accuracy and lack of valid absences; conversely they provide comprehensive views of species distributions allowing to realistically capture most of their realized environmental niche. Rare species are the most in need of predictive distribution modelling but also the most difficult to model. We refer to this contrast as the "rare species modelling paradox" and propose as a solution developing modelling approaches that deal with a sufficiently large set of predictors, ensuring that statistical models aren't overfitted. Our novel approach fulfils this condition by fitting a large number of bivariate models and averaging them with a weighted ensemble approach. We further propose that this ensemble forecasting is conducted within a hierarchic multi-scale framework. We present two ensemble models for a test species, one at regional and one at local scale, each based on the combination of 630 models. In both cases, we obtained excellent spatial projections, unusual when modelling rare species. Model results highlight, from a statistically sound approach, the effects of multiple drivers in a same modelling framework and at two distinct scales. From this added information, regional models can support accurate forecasts of range dynamics under climate change scenarios, whereas local models allow the assessment of isolated or synergistic impacts of changes in multiple predictors. This novel framework provides a baseline for adaptive conservation, management and monitoring of rare species at distinct spatial and temporal scales.

Unification of multi-species vertebrate anatomy ontologies for comparative biology in Uberon.

BACKGROUND: Elucidating disease and developmental dysfunction requires understanding variation in phenotype. Single-species model organism anatomy ontologies (ssAOs) have been established to represent this variation. Multi-species anatomy ontologies (msAOs; vertebrate skeletal, vertebrate homologous, teleost, amphibian AOs) have been developed to represent 'natural' phenotypic variation across species. Our aim has been to integrate ssAOs and msAOs for various purposes, including establishing links between phenotypic variation and candidate genes. RESULTS: Previously, msAOs contained a mixture of unique and overlapping content. This hampered integration and coordination due to the need to maintain cross-references or inter-ontology equivalence axioms to the ssAOs, or to perform large-scale obsolescence and modular import. Here we present the unification of anatomy ontologies into Uberon, a single ontology resource that enables interoperability among disparate data and research groups. As a consequence, independent development of TAO, VSAO, AAO, and vHOG has been discontinued. CONCLUSIONS: The newly broadened Uberon ontology is a unified cross-taxon resource for metazoans (animals) that has been substantially expanded to include a broad diversity of vertebrate anatomical structures, permitting reasoning across anatomical variation in extinct and extant taxa. Uberon is a core resource that supports single- and cross-species queries for candidate genes using annotations for phenotypes from the systematics, biodiversity, medical, and model organism communities, while also providing entities for logical definitions in the Cell and Gene Ontologies. THE ONTOLOGY RELEASE FILES ASSOCIATED WITH THE ONTOLOGY MERGE DESCRIBED IN THIS MANUSCRIPT ARE AVAILABLE AT: http://purl.obolibrary.org/obo/uberon/releases/2013-02-21/ CURRENT ONTOLOGY RELEASE FILES ARE AVAILABLE ALWAYS AVAILABLE AT: http://purl.obolibrary.org/obo/uberon/releases/

Testing the limits of quasi-geostrophic theory: application to observed laboratory flows outside the quasi-geostrophic regime

We compare laboratory observations of equilibrated baroclinic waves in the rotating two-layer annulus, with numerical simulations from a quasi-geostrophic model. The laboratory experiments lie well outside the quasi-geostrophic regime: the Rossby number reaches unity; the depth-to-width aspect ratio is large; and the fluid contains ageostrophic inertia–gravity waves. Despite being formally inapplicable, the quasi-geostrophic model captures the laboratory flows reasonably well. The model displays several systematic biases, which are consequences of its treatment of boundary layers and neglect of interfacial surface tension and which may be explained without invoking the dynamical effects of the moderate Rossby number, large aspect ratio or inertia–gravity waves. We conclude that quasi-geostrophic theory appears to continue to apply well outside its formal bounds.

Nonlinear stability of continuously stratified quasi-geostrophic flow

Nonlinear stability theorems analogous to Arnol'd's second stability theorem are established for continuously stratified quasi-geostrophic flow with general nonlinear boundary conditions in a vertically and horizontally confined domain. Both the standard quasi-geostrophic model and the modified quasi-geostrophic model (incorporating effects of hydrostatic compressibility) are treated. The results establish explicit upper bounds on the disturbance energy, the disturbance potential enstrophy, and the disturbance available potential energy on the horizontal boundaries, in terms of the initial disturbance fields. Nonlinear stability in the sense of Liapunov is also established.

The coevolution of two phytoplankton species on a single resource: allelopathy as a pseudo-mixotrophy

Without the top-down effects and the external/physical forcing, a stable coexistence of two phytoplankton species under a single resource is impossible — a result well known from the principle of competitive exclusion. Here I demonstrate by analysis of a mathematical model that such a stable coexistence in a homogeneous media without any external factor would be possible, at least theoretically, provided (i) one of the two species is toxin producing thereby has an allelopathic effect on the other, and (ii) the allelopathic effect exceeds a critical level. The threshold level of allelopathy required for the coexistence has been derived analytically in terms of the parameters associated with the resource competition and the nutrient recycling. That the extra mortality of a competitor driven by allelopathy of a toxic species gives a positive feed back to the algal growth process through the recycling is explained. And that this positive feed back plays a pivotal role in reducing competition pressures and helping species succession in the two-species model is demonstrated. Based on these specific coexistence results, I introduce and explain theoretically the allelopathic effect of a toxic species as a ‘pseudo-mixotrophy’—a mechanism of ‘if you cannot beat them or eat them, just kill them by chemical weapons’. The impact of this mechanism of species succession by pseudo-mixotrophy in the form of alleopathy is discussed in the context of current understanding on straight mixotrophy and resource-species relationship among phytoplankton species.

Applying a new methodology to verify transmission line model performance - the equivalent impedance test

The objective of this paper is to present a methodology to analyze a transmission line model used in electromagnetic transitory simulators, called equivalent impedance test. Initially the definition of equivalent impedance reference test is shown. Soon after this methodology is applied to a transmission line model, the Quasi-Modes model. The studies were accomplished in a hypothetical non-transposed three-phase transmission fine of 440 kV. The line length is 500 km, and it was modeled through cascades of pi-circuits (with 50 pi's circuits, each with 10 km length).

Impaired immune evasion in HIV through intracellular delays and multiple infection of cells

With its high mutation rate, HIV is capable of escape from recognition, suppression and/or killing by CD8(+) cytotoxic T lymphocytes (CTLs). The rate at which escape variants replace each other can give insights into the selective pressure imposed by single CTL clones. We investigate the effects of specific characteristics of the HIV life cycle on the dynamics of immune escape. First, it has been found that cells in HIV-infected patients can carry multiple copies of proviruses. To investigate how this process affects the emergence of immune escape, we develop a mathematical model of HIV dynamics with multiple infections of cells. Increasing the frequency of multiple-infected cells delays the appearance of immune escape variants, slows down the rate at which they replace the wild-type variant and can even prevent escape variants from taking over the quasi-species. Second, we study the effect of the intracellular eclipse phase on the rate of escape and show that escape rates are expected to be slower than previously anticipated. In summary, slow escape rates do not necessarily imply inefficient CTL-mediated killing of HIV-infected cells, but are at least partly a result of the specific characteristics of the viral life cycle.

Two-temperature model of the coronal irradiated pellets

A two electron-temperature, quasi-steady model of the corona of a laser-ablated pellet is considered. Ablation pressure, critical radius and mass flow rate are determined. Results are close to those obtained with heat flux saturation well below the free-streaming limit.