Size dependence of biomass spectra and abundance spectra: the optimal distributions

Based on the hypothesis of self-optimization, we derive four models of biomass spectra and abundance spectra in communities with size-dependent metabolic rates. In Models 1 and 2, the maximum diversity of population abundance in different size classes subject to the constraints of constant mean body mass and constant mean respiration rate is assumed to be the strategy for ecosystems to organize their size structure. In Models 3 and 4, the organizing strategy is defined as the maximum diversity of biomass in different size classes without constraints on mean body mass and subject to the constant mean specific respiration rate of all individuals, i.e. the average specific respiration rate over all individuals of a community or group, which characterizes the mean rate of energy consumption in a community. Models 1 and 2 generate peaked distributions of biomass spectral density whereas Model 3 generates a fiat distribution. In Model 4, the distributions of biomass spectral density and of abundance spectral density depend on the Lagrangian multipler (lambda (2)). When lambda (2) tends to zero or equals zero, the distributions of biomass spectral density and of abundance spectral density correspond to those from Model 3. When lambda (2) has a large negative value, the biomass spectrum is similar to the empirical fiat biomass spectrum organized in logarithmic size intervals. When lambda (2) > 0, the biomass spectral density increases with body mass and the distribution of abundance spectral density is an unimodal curve. (C) 2001 Elsevier Science B.V. All rights reserved.

Oscillations by minimal bacterial suicide circuits reveal hidden facets of host-circuit physiology.

Synthetic biology seeks to enable programmed control of cellular behavior though engineered biological systems. These systems typically consist of synthetic circuits that function inside, and interact with, complex host cells possessing pre-existing metabolic and regulatory networks. Nevertheless, while designing systems, a simple well-defined interface between the synthetic gene circuit and the host is frequently assumed. We describe the generation of robust but unexpected oscillations in the densities of bacterium Escherichia coli populations by simple synthetic suicide circuits containing quorum components and a lysis gene. Contrary to design expectations, oscillations required neither the quorum sensing genes (luxR and luxI) nor known regulatory elements in the P(luxI) promoter. Instead, oscillations were likely due to density-dependent plasmid amplification that established a population-level negative feedback. A mathematical model based on this mechanism captures the key characteristics of oscillations, and model predictions regarding perturbations to plasmid amplification were experimentally validated. Our results underscore the importance of plasmid copy number and potential impact of "hidden interactions" on the behavior of engineered gene circuits - a major challenge for standardizing biological parts. As synthetic biology grows as a discipline, increasing value may be derived from tools that enable the assessment of parts in their final context.

Fitting host-parasitoid models with CV² > 1 using hierarchical generalized linear models

The powerful general Pacala-Hassell host-parasitoid model for a patchy environment, which allows host density–dependent heterogeneity (HDD) to be distinguished from between-patch, host density–independent heterogeneity (HDI), is reformulated within the class of the generalized linear model (GLM) family. This improves accessibility through the provision of general software within well–known statistical systems, and allows a rich variety of models to be formulated. Covariates such as age class, host density and abiotic factors may be included easily. For the case where there is no HDI, the formulation is a simple GLM. When there is HDI in addition to HDD, the formulation is a hierarchical generalized linear model. Two forms of HDI model are considered, both with between-patch variability: one has binomial variation within patches and one has extra-binomial, overdispersed variation within patches. Examples are given demonstrating parameter estimation with standard errors, and hypothesis testing. For one example given, the extra-binomial component of the HDI heterogeneity in parasitism is itself shown to be strongly density dependent.

A unified 3D model of the VAR process

A 3D time-dependent model of the VAR process has been developed using CFD techniques. The model solves the coupled field equations for fluid flow, heat transfer (including phase change) and electromagnetic field, for both the electrode and the ingot. The motion of the electic arc 'preferred spot' can be specified based on observations. Correlations are sought between the local gap height, resulting from instantaneous liquid pool surface shape and electrode tip shape, and the arc motion. The detailed behaviour of the melting film on the electrode tip is studies using a spectral free surface technique, which allows investigation of the drops' detachment and drip shorts.

Models of open populations with space-limited recruitment: extension of theory and application to the barnacle Chthamalus montagui.

1. Barnacles are a good model organism for the study of open populations with space-limited recruitment. These models are applicable to other species with open supply of new individuals and resource limitation. The inclusion of space in models leads to reductions in recruitment with increasing density, and thus predictions of population size and stability are possible. 2. Despite the potential generality of a demographic theory for open space-limited populations, the models currently have a narrow empirical base. In this study, a model for an open population with space-limited recruitment was extended to include size-specific survival and promotions to any size class. The assumptions of this model were tested using data from a pan-European study of the barnacle Chthamalus montagui Southward. Two models were constructed: a 6-month model and a periodic annual model. Predicted equilibria and their stabilities were compared between shores. 3. Tests of model assumptions supported the extension of the theory to include promotions to any size class. Mortality was found to be size-specific and density independent. Studied populations were open, with recruitment proportional to free space. 4. The 6-month model showed a significant interaction between time and location for equilibrium free space. This may have been due to contrasts in the timing of structuring processes (i.e. creating and filling space) between Mediterranean and Atlantic systems. Integration of the 6-month models into a periodic annual model removed the differences in equilibrium-free space between locations. 5. Model predictions show a remarkable similarity between shores at a European scale. Populations were persistent and all solutions were stable. This reflects the apparent absence of density-dependent mortality and a high adult survivorship in C. montagui. As populations are intrinsically stable, observations of fluctuations in density are directly attributable to variations in the environmental forcing of recruitment or mortality

Space and phase resolved plasma parameters in an industrial dual-frequency capacitively coupled radio-frequency discharge

The dynamics of high energetic electrons (>= 11.7 eV) in a modified industrial confined dual-frequency capacitively coupled RF discharge (Exelan, Lam Research Inc.), operated at 1.937 MHz and 27.118 MHz, is investigated by means of phase resolved optical emission spectroscopy. Operating in a He-O-2. plasma with small rare gas admixtures the emission is measured, with one-dimensional spatial resolution along the discharge axis. Both the low and high frequency RF cycle are resolved. The diagnostic is based on time dependent measurements of the population densities of specifically chosen excited rare gas states. A time dependent model, based on rate equations, describes the dynamics of the population densities of these levels. Based on this model and the comparison of the excitation of various rare gas states, with different excitation thresholds, time and space resolved electron temperature, propagation velocity and qualitative electron density as well as electron energy distribution functions are determined. This information leads to a better understanding of the dual-frequency sheath dynamics and shows, that separate control of ion energy and electron density is limited.

Cyanopolyynes in hot cores: modelling G305.2+0.2

We present results from a time-dependent gas-phase chemical model of a hot core based on the physical conditions of G305.2+0.2. While the cyanopolyyne HC3N has been observed in hot cores, the longer chained species, HC5N, HC7N and HC9N, have not been considered as the typical hot-core species. We present results which show that these species can be formed under hot core conditions. We discuss the important chemical reactions in this process and, in particular, show that their abundances are linked to the parent species acetylene which is evaporated from icy grain mantles. The cyanopolyynes show promise as ‘chemical clocks’ which may aid future observations in determining the age of hot core sources. The abundance of the larger cyanopolyynes increases and decreases over relatively short time-scales, ~10^2.5 yr. We present results from a non-local thermodynamic equilibrium statistical equilibrium excitation model as a series of density, temperature and column density dependent contour plots which show both the line intensities and several line ratios. These aid in the interpretation of spectral-line data, even when there is limited line information available. In particular, non-detections of HC5N and HC7N in Walsh et al. are analysed and discussed.

Range expansion in an invasive small mammal: influence of life-history and habitat quality.

Invasive species pose a major threat to biodiversity but provide an opportunity to describe the processes that lead to changes in a species’ range. The bank vole (Myodes glareolus) is an invasive rodent that was introduced to Ireland in the early twentieth century. Given its continuing range expansion, the substantial empirical data on its spread thus far, and the absence of any eradication program, the bank vole in Ireland represents a unique model system for studying the mechanisms influencing the rate of range expansion in invasive small mammals. We described the invasion using a reaction–diffusion model informed by empirical data on life history traits and demographic parameters. We subsequently modelled the processes involved in its range expansion using a rule-based spatially explicit simulation. Habitat suitability interacted with density-dependent parameters to influence dispersal, most notably the density at which local populations started to donate emigrating individuals, the number of dispersing individuals and the direction of dispersal. Whilst local habitat variability influenced the rate of spread, on a larger scale the invasion resembled a simple reaction–diffusion process. Our results suggest a Type 1 range expansion where the rate of expansion is generally constant over time, but with some evidence for a lag period following introduction. We demonstrate that a two-parameter empirical model and a rule-based spatially explicit simulation are sufficient to accurately describe the invasion history of a species that exhibits a complex, density-dependent pattern of dispersal.

Dynamics of the entanglement spectrum in spin chains

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

Home-range size in juveniles of the temperate reef fish, the common triplefin (Forsterygion lapillum)

An organism’s home range dictates the spatial scale on which important processes occur (e.g. competition and predation) and directly affects the relationship between individual fitness and local habitat quality. Many reef fish species have very restricted home ranges after settlement and, here, we quantify home-range size in juveniles of a widespread and abundant reef fish in New Zealand, the common triplefin (Forsterygion lapillum). We conducted visual observations on 49 juveniles (mean size = 35-mm total length) within the Wellington harbour, New Zealand. Home ranges were extremely small, 0.053 m2 ± 0.029 (mean ± s.d.) and were unaffected by adult density, body size or substrate composition. A regression tree indicated that home-range size sharply decreased ~4.5 juveniles m–2 and a linear mixed model confirmed that home-range sizes in high-density areas (>4.5 juveniles m–2) were significantly smaller (34%) than those in low-density areas (after accounting for a significant effect of fish movement on our home-range estimates). Our results suggest that conspecific density may have negative and non-linear effects on home-range size, which could shape the spatial distribution of juveniles within a population, as well as influence individual fitness across local density gradients.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.

Experimentelle und numerische Untersuchungen zur Dispersion von Dichteströmungen in einem stochastischen Modellaquifer

Es ist bekannt, dass die Dichte eines gelösten Stoffes die Richtung und die Stärke seiner Bewegung im Untergrund entscheidend bestimmen kann. Eine Vielzahl von Untersuchungen hat gezeigt, dass die Verteilung der Durchlässigkeiten eines porösen Mediums diese Dichteffekte verstärken oder abmindern kann. Wie sich dieser gekoppelte Effekt auf die Vermischung zweier Fluide auswirkt, wurde in dieser Arbeit untersucht und dabei das experimentelle sowohl mit dem numerischen als auch mit dem analytischen Modell gekoppelt. Die auf der Störungstheorie basierende stochastische Theorie der macrodispersion wurde in dieser Arbeit für den Fall der transversalen Makodispersion. Für den Fall einer stabilen Schichtung wurde in einem Modelltank (10m x 1.2m x 0.1m) der Universität Kassel eine Serie sorgfältig kontrollierter zweidimensionaler Experimente an einem stochastisch heterogenen Modellaquifer durchgeführt. Es wurden Versuchsreihen mit variierenden Konzentrationsdifferenzen (250 ppm bis 100 000 ppm) und Strömungsgeschwindigkeiten (u = 1 m/ d bis 8 m/d) an drei verschieden anisotrop gepackten porösen Medien mit variierender Varianzen und Korrelationen der lognormal verteilten Permeabilitäten durchgeführt. Die stationäre räumliche Konzentrationsausbreitung der sich ausbreitenden Salzwasserfahne wurde anhand der Leitfähigkeit gemessen und aus der Höhendifferenz des 84- und 16-prozentigen relativen Konzentrationsdurchgang die Dispersion berechnet. Parallel dazu wurde ein numerisches Modell mit dem dichteabhängigen Finite-Elemente-Strömungs- und Transport-Programm SUTRA aufgestellt. Mit dem kalibrierten numerischen Modell wurden Prognosen für mögliche Transportszenarien, Sensitivitätsanalysen und stochastische Simulationen nach der Monte-Carlo-Methode durchgeführt. Die Einstellung der Strömungsgeschwindigkeit erfolgte - sowohl im experimentellen als auch im numerischen Modell - über konstante Druckränder an den Ein- und Auslauftanks. Dabei zeigte sich eine starke Sensitivität der räumlichen Konzentrationsausbreitung hinsichtlich lokaler Druckvariationen. Die Untersuchungen ergaben, dass sich die Konzentrationsfahne mit steigendem Abstand von der Einströmkante wellenförmig einem effektiven Wert annähert, aus dem die Makrodispersivität ermittelt werden kann. Dabei zeigten sich sichtbare nichtergodische Effekte, d.h. starke Abweichungen in den zweiten räumlichen Momenten der Konzentrationsverteilung der deterministischen Experimente von den Erwartungswerten aus der stochastischen Theorie. Die transversale Makrodispersivität stieg proportional zur Varianz und Korrelation der lognormalen Permeabilitätsverteilung und umgekehrt proportional zur Strömungsgeschwindigkeit und Dichtedifferenz zweier Fluide. Aus dem von Welty et al. [2003] mittels Störungstheorie entwickelten dichteabhängigen Makrodispersionstensor konnte in dieser Arbeit die stochastische Formel für die transversale Makrodispersion weiter entwickelt und - sowohl experimentell als auch numerisch - verifiziert werden.

Numerical Feasibility Study for Treated Wastewater Recharge as a Tool to Impede Saltwater Intrusion in the Coastal Aquifer of Gaza – Palestine

The ongoing depletion of the coastal aquifer in the Gaza strip due to groundwater overexploitation has led to the process of seawater intrusion, which is continually becoming a serious problem in Gaza, as the seawater has further invaded into many sections along the coastal shoreline. As a first step to get a hold on the problem, the artificial neural network (ANN)-model has been applied as a new approach and an attractive tool to study and predict groundwater levels without applying physically based hydrologic parameters, and also for the purpose to improve the understanding of complex groundwater systems and which is able to show the effects of hydrologic, meteorological and anthropogenic impacts on the groundwater conditions. Prediction of the future behaviour of the seawater intrusion process in the Gaza aquifer is thus of crucial importance to safeguard the already scarce groundwater resources in the region. In this study the coupled three-dimensional groundwater flow and density-dependent solute transport model SEAWAT, as implemented in Visual MODFLOW, is applied to the Gaza coastal aquifer system to simulate the location and the dynamics of the saltwater–freshwater interface in the aquifer in the time period 2000-2010. A very good agreement between simulated and observed TDS salinities with a correlation coefficient of 0.902 and 0.883 for both steady-state and transient calibration is obtained. After successful calibration of the solute transport model, simulation of future management scenarios for the Gaza aquifer have been carried out, in order to get a more comprehensive view of the effects of the artificial recharge planned in the Gaza strip for some time on forestall, or even to remedy, the presently existing adverse aquifer conditions, namely, low groundwater heads and high salinity by the end of the target simulation period, year 2040. To that avail, numerous management scenarios schemes are examined to maintain the ground water system and to control the salinity distributions within the target period 2011-2040. In the first, pessimistic scenario, it is assumed that pumping from the aquifer continues to increase in the near future to meet the rising water demand, and that there is not further recharge to the aquifer than what is provided by natural precipitation. The second, optimistic scenario assumes that treated surficial wastewater can be used as a source of additional artificial recharge to the aquifer which, in principle, should not only lead to an increased sustainable yield of the latter, but could, in the best of all cases, revert even some of the adverse present-day conditions in the aquifer, i.e., seawater intrusion. This scenario has been done with three different cases which differ by the locations and the extensions of the injection-fields for the treated wastewater. The results obtained with the first (do-nothing) scenario indicate that there will be ongoing negative impacts on the aquifer, such as a higher propensity for strong seawater intrusion into the Gaza aquifer. This scenario illustrates that, compared with 2010 situation of the baseline model, at the end of simulation period, year 2040, the amount of saltwater intrusion into the coastal aquifer will be increased by about 35 %, whereas the salinity will be increased by 34 %. In contrast, all three cases of the second (artificial recharge) scenario group can partly revert the present seawater intrusion. From the water budget point of view, compared with the first (do nothing) scenario, for year 2040, the water added to the aquifer by artificial recharge will reduces the amount of water entering the aquifer by seawater intrusion by 81, 77and 72 %, for the three recharge cases, respectively. Meanwhile, the salinity in the Gaza aquifer will be decreased by 15, 32 and 26% for the three cases, respectively.

The counter-propagating Rossby-wave perspective on baroclinic instability. II: Application to the Charney model

The constant-density Charney model describes the simplest unstable basic state with a planetary-vorticity gradient, which is uniform and positive, and baroclinicity that is manifest as a negative contribution to the potential-vorticity (PV) gradient at the ground and positive vertical wind shear. Together, these ingredients satisfy the necessary conditions for baroclinic instability. In Part I it was shown how baroclinic growth on a general zonal basic state can be viewed as the interaction of pairs of ‘counter-propagating Rossby waves’ (CRWs) that can be constructed from a growing normal mode and its decaying complex conjugate. In this paper the normal-mode solutions for the Charney model are studied from the CRW perspective. Clear parallels can be drawn between the most unstable modes of the Charney model and the Eady model, in which the CRWs can be derived independently of the normal modes. However, the dispersion curves for the two models are very different; the Eady model has a short-wave cut-off, while the Charney model is unstable at short wavelengths. Beyond its maximum growth rate the Charney model has a neutral point at finite wavelength (r=1). Thereafter follows a succession of unstable branches, each with weaker growth than the last, separated by neutral points at integer r—the so-called ‘Green branches’. A separate branch of westward-propagating neutral modes also originates from each neutral point. By approximating the lower CRW as a Rossby edge wave and the upper CRW structure as a single PV peak with a spread proportional to the Rossby scale height, the main features of the ‘Charney branch’ (0model where the positive PV gradient exists only at the tropopause

Significance of cross correlations in the stress relaxation of polymer melts

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent "colloidal" relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. (c) 2007 American Institute of Physics