Thermal effects in the evolution of initially layered mantle material

A simplified model for anisotropic mantle convection based on a novel class of rheologies, originally developed for folding instabilities in multilayered rock (MUHLHAUS et al., 2002), is extended ¨ through the introduction of a thermal anisotropy dependent on the local layering. To examine the effect of the thermal anisotropy on the evolution of mantle material, a parallel implementation of this model was undertaken using the Escript modelling toolkit and the Finley finite-element computational kernel (DAVIES et al., 2004). For the cases studied, there appears too little if any effect. For comparative purposes, the effects of anisotropic shear viscosity and the introduced thermal anisotropy are also presented. These results contribute to the characterization of viscous anisotropic mantle convection subject to variation in thermal conductivities and shear viscosities.

Evolution of Wolbachia pipientis transmission dynamics in insects

Wolbachia pipientis is an intracellular bacterial parasite of arthropods that enhances its transmission by manipulating host reproduction, most commonly by inducing cytoplasmic incompatibility. The discovery of isolates with modified cytoplasmic incompatibility phenotypes and others with novel virulence properties is an indication of the potential breadth of evolutionary strategies employed by Wolbachia.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

Evolution of a mate recognition system after hybridization between two Drosophila species

I investigated the genetic relationship between male and female components of the mate recognition system and how this relationship influenced the subsequent evolution of the two traits, in a series of replicate populations of interspecific hybrids. Thirty populations of hybrids between Drosophila serrata and Drosophila birchii were established and maintained for 24 generations. At the fifth generation after hybridization, the mating success of hybrid individuals with the D. serrata parent was determined. The genetic correlation between male and female components of the male recognition system, as a consequence of pleiotropy or tight physical linkage, was found to be significant but low (r = 0.388). This result suggested that pleiotropy may play only a minor role in the evolution of mate recognition in this system. At the twenty-fourth generation after hybridization, the mating success of the hybrids was again determined. The evolution of male and female components was investigated by analyzing the direction of evolution of each hybrid line with respect to its initial position in relation to the genetic regression. Male and female components appeared to converge on a single equilibrium point, rather than evolving along trajectories with slope equal to the genetic regression, toward a line of equilibria.

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.

Genetic structure and evolution of species in the mangrove genus Avicennia (Avicenniaceae) in the Indo-West Pacific

Allozyme variation in species of the mangrove genus Avicennia was screened in 25 populations collected from 22 locations in the Indo-West Pacific and eastern North America using 11 loci. Several fixed gene differences supported the specific status of Avicennia alba, A. integra, A. marina, and A. rumphiana from the Indo-West Pacific, and A. germinans from the Atlantic-East Pacific. The three varieties of A. marina, var. marina, var. eucalyptifolia, and var. australasica, had higher genetic similarities (Nei's I) and no fixed gene differences, confirming their conspecific status. Strong genetic structuring was observed in A. marina, with sharp changes in gene frequencies at the geographical margins of varietal distributions. The occurrence of alleles found otherwise in only one variety, in only immediately adjacent populations of another variety, provided evidence of introgession between varieties. The varieties appear to have diverged recently in the Pleistocene and are apparently not of ancient Cretaceous origin, as suggested earlier. Despite evidence of high degrees of outcrossing, gene flow among populations was relatively low (N(e)m less than or equal to 1-2), except where populations were geographically continuous, questioning assumptions that these widespread mangrove species achieve high levels of long-distance dispersal.

Crystal structure of human T cell leukemia virus type 1 gp21 ectodomain crystallized as a maltose-binding protein chimera reveals structural evolution of retroviral transmembrane proteins

Retroviral entry into cells depends on envelope glycoproteins, whereby receptor binding to the surface-exposed subunit triggers membrane fusion by the transmembrane protein (TM) subunit. We determined the crystal structure at 2.5-Angstrom resolution of the ectodomain of gp21, the TM from human T cell leukemia virus type 1. The gp21 fragment was crystallized as a maltose-binding protein chimera, and the maltose-binding protein domain was used to solve the initial phases by the method of molecular replacement. The structure of gp21 comprises an N-terminal trimeric coiled coil, an adjacent disulfide-bonded loop that stabilizes a chain reversal, and a C-terminal sequence structurally distinct from HIV type 1/simian immunodeficiency virus gp41 that packs against the coil in an extended antiparallel fashion. Comparison of the gp21 structure with the structures of other retroviral TMs contrasts the conserved nature of the coiled coil-forming region and adjacent disulfide-bonded loop with the variable nature of the C-terminal ectodomain segment. The structure points to these features having evolved to enable the dual roles of retroviral TMs: conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface. The structure of gp21 implies that the N-terminal fusion peptide is in close proximity to the C-terminal transmembrane domain and likely represents a postfusion conformation.

Some new examples for nonuniqueness of the evolution problem of harmonic maps

We find some new examples to show nonuniquence for the heat flow of harmonic maps where weak solutions satisfy the same monotonicity property.

Rapid evolution towards equal sex ratios in a system with heterogamety

The equal sex ratios found in many species with heterogametic sex determination may be a consequence of selection for equality or the result of the Mendelian segregation of the two sex chromosomes. A lack of genetic variation in sex ratio in species with heterogamety has been the major obstacle in distinguishing between these two hypotheses. We overcome this obstacle by generating hybrids between two species of Drosophila. The resulting hybrid lines had biased sex ratios, allowing us to observe the evolution of sex ratio in replicate populations. Sex ratio converged towards 1:1 after 16 generations of natural selection. These changes in sex ratio were not due to differences in viability between the sexes and the loci underlying the variation in sex ratio were not sex-linked. Equal sex ratios may therefore be the result of natural selection as Fisher predicted.