Determinants of species richness in the Park Grass experiment

The Park Grass Experiment at Rothamsted in southeast England was started in 1856, making it the longest-running experiment in plant ecology anywhere in the world. Experimental inputs include a range of fertilizers (nitrogen, phosphorus, potassium, and organic manures) applied annually, with lime applied occasionally, and these have led to an increase in biomass and, where nitrogen was applied in the form of ammonium sulfate, to substantial decreases in soil pH. The number of species per plot varies from three to 44 per 200 m2, affording a unique opportunity to study the determinants of plant species richness and to estimate the effect sizes attributable to different factors. The response of species richness to biomass depends on the amount and type of nitrogen applied; richness declined monotonically with increasing biomass on plots receiving no nitrogen or receiving nitrogen in the form of sodium nitrate, but there was no relationship between species richness and biomass on plots acidified by ammonium sulfate application. The response to lime also depended on the type of nitrogen applied; there was no relationship between lime treatment and species richness, except in plots receiving nitrogen in the form of ammonium sulfate, where species richness increased sharply with increasing soil pH. The addition of phosphorus reduced species richness, and application of potassium along with phosphorus reduced species richness further, but the biggest negative effects were when nitrogen and phosphorus were applied together. The analysis demonstrates how multiple factors contribute to the observed diversity patterns and how environmental regulation of species pools can operate at the same spatial and temporal scale as biomass effects

Scaling Regimes, Crossovers, and Lattice Corrections in 2D Heisenberg Antiferromagnets

We study scaling behavior in 2D, S=1/2 and S=1 Heisenberg antiferromagnets using the data on full q-dependences of the equal time structure factor and the static susceptibility, calculated through high temperature expansions. We also carry out comparisons with a model of two coupled S=1/2 planes with the interlayer coupling tuned to the T=0 critical point. We separately determine the spin-wave velocity c and mass $m=c/\xi$, in addition to the correlation length, $\xi$, and find that c is temperature dependent; only for T\alt JS, it approaches its known T=0 value $c_0$. Despite this temperature dependent spin-wave velocity, full q- and $\omega$-dependences of the dynamical susceptibility $\chi(\bf q,\omega)$ agree with the universal scaling functions computable for the $\sigma$-model, for temperatures upto $T_0 \sim 0.6c_0/a$. Detailed comparisons show that below $T_0$ the S=1 model is in the renormalized classical (RC) regime, the two plane model is in the quantum critical (QC) regime, and the S=1/2 model exhibits a RC-QC crossover, centered at T=0.55J. In particular, for the S=1/2 model above this crossover and for the two-plane model at all T, the spin-wave mass is in excellent agreement with the universal QC prediction, $m\simeq 1.04\,T$. In contrast, for the S=1/2 model below the RC-QC crossover, and for the S=1 model at all T, the behavior agrees with the known RC expression. For all models nonuniversal behavior occurs above $T\sim 0.6c_0/a$. Our results strongly support the conjecture of Chubukov and Sachdev that the S=1/2 model is close to the T=0 critical point to exhibit QC behavior.Comment: 13 pages, REVTeX with attached PostScript (see file for addl info

Meshfree finite differences for vector Poisson and pressure Poisson equations with electric boundary conditions

We demonstrate how meshfree finite difference methods can be applied to solve vector Poisson problems with electric boundary conditions. In these, the tangential velocity and the incompressibility of the vector field are prescribed at the boundary. Even on irregular domains with only convex corners, canonical nodal-based finite elements may converge to the wrong solution due to a version of the Babuska paradox. In turn, straightforward meshfree finite differences converge to the true solution, and even high-order accuracy can be achieved in a simple fashion. The methodology is then extended to a specific pressure Poisson equation reformulation of the Navier-Stokes equations that possesses the same type of boundary conditions. The resulting numerical approach is second order accurate and allows for a simple switching between an explicit and implicit treatment of the viscosity terms.Comment: 19 pages, 7 figure

Incipient Separation in Shock Wave Boundary Layer Interactions as Induced by Sharp Fin

The incipient separation induced by the shock wave turbulent boundary layer interaction at the sharp fin is the subject of present study. Existing theories for the prediction of incipient separation, such as those put forward by McCabe (1966) and Dou and Deng (1992), can have thus far only predicting the direction of surface streamline and tend to over-predict the incipient separation condition based on the Stanbrook's criterion. In this paper, the incipient separation is firstly predicted with Dou and Deng (1992)'s theory and then compared with Lu and Settles (1990)' experimental data. The physical mechanism of the incipient separation as induced by the shock wave/turbulent boundary layer interactions at sharp fin is explained via the surface flow pattern analysis. Furthermore, the reason for the observed discrepancy between the predicted and experimental incipient separation conditions is clarified. It is found that when the wall limiting streamlines behind the shock wave becomes\ aligning with one ray from the virtual origin as the strength of shock wave increases, the incipient separation line is formed at which the wall limiting streamline becomes perpendicular to the local pressure gradient. The formation of this incipient separation line is the beginning of the separation process. The effects of Reynolds number and the Mach number on incipient separation are also discussed. Finally, a correlation for the correction of the incipient separation angle as predicted by the theory is also given.Comment: 34 pages; 9 figure

Characterisation of haematological profiles and whole blood relative gene expression levels in Holstein-Friesian and Jersey bull calves undergoing gradual weaning

peer-reviewedHaematological profiles indicate the health status of an animal and can be used to identify sub-clinical stress responses. The objectives of the study were to examine (i) the effect of breed and plane of nutrition, on haematological profiles of artificially reared Holstein-Friesian and Jersey bull calves in response to gradual weaning, and (ii) the effect of breed on immune response genes in bovine whole blood using real-time quantitative PCR. Holstein-Friesian and Jersey bull calves were group housed indoors and individually fed using an automatic feeder. They were allocated to a high, medium or low plane of nutrition, based on milk replacer (MR) and concentrate. The nutrition treatments were calculated using National Research Council guidelines in order to achieve a high, medium or low growth rate for each respective breed. During the weaning phase MR was gradually reduced over a 14-day (d) period (d −13 to d 0). Calves were blood sampled on d −14, −6, −3, 0, 1, 3, 8 and 14 relative to weaning (d 0) for subsequent haematological analysis. On d −14, 1 and 8, a subset of eight Holstein-Friesian calves randomly selected from the medium nutrition treatment and eight Jersey calves randomly selected from the high nutrition treatment, were blood sampled for gene expression profiling, targeting biomarkers of weaning stress. These two treatment groups were chosen to examine the effect of breed on expression of the genes of interest, as energy intake and animal performance were similar. There was no effect of breed×plane of nutrition interaction nor effect of plane of nutrition on any variable measured (P>0.05). Gradual weaning produced differential biological responses in the two breeds evidenced by breed×time interactions for lymphocyte, monocyte and red blood cell number, plasma haemoglobin and haptoglobin concentrations (P<0.05). The typical stress response consisting of neutrophilia and lymphopaenia was not observed for any treatment. An immune response to gradual weaning was observed as the relative gene expression level of the pro-apoptotic gene, Fas, increased on d 1 relative to d −14 (P<0.05). Relative gene expression levels were greater in Jersey calves compared with Holstein-Friesian for the pro-inflammatory cytokine CXCL8 (P=0.05) and the glucocorticoid receptor, GRα (P<0.05). The increased levels of these transcripts suggest that Jersey calves may have a more sensitive immune system compared with Holstein-Friesian

Covid-19, Ebola, and HIV - Leveraging Lessons to Maximize Impact

As researchers who have long experience with HIV/AIDS prevention, vaccines, and therapies, some of whom also have experience with Ebola, we believe it is critical to build the response to the Covid-19 pandemic on lessons from the HIV pandemic and recent Ebola outbreaks

Crossover and scaling in a nearly antiferromagnetic Fermi liquid in two dimensions

We consider two-dimensional Fermi liquids in the vicinity of a quantum transition to a phase with commensurate, antiferromagnetic long-range order. Depending upon the Fermi surface topology, mean-field spin-density-wave theory predicts two different types of such transitions, with mean-field dynamic critical exponents $z=1$ (when the Fermi surface does not cross the magnetic zone boundary, type $A$) and $z=2$ (when the Fermi surface crosses the magnetic zone boundary, type $B$). The type $A$ system only displays $z=1$ behavior at all energies and its scaling properties are similar (though not identical) to those of an insulating Heisenberg antiferromagnet. Under suitable conditions precisely stated in this paper, the type $B$ system displays a crossover from relaxational behavior at low energies to type $A$ behavior at high energies. A scaling hypothesis is proposed to describe this crossover: we postulate a universal scaling function which determines the entire, temperature-, wavevector-, and frequency-dependent, dynamic, staggered spin susceptibility in terms of 4 measurable, $T=0$, parameters (determining the distance, energy, and order parameter scales, plus one crossover parameter). The scaling function contains the full scaling behavior in all regimes for both type $A$ and $B$ systems. The crossover behavior of the uniform susceptibility and the specific heat is somewhat more complicated and is also discussed. Explicit computation of the crossover functions is carried out in a large $N$ expansion on a mean-field model. Some new results for the critical properties on the ordered side of the transition are also obtained in a spin-density wave formalism. The possible relevance of our results to the doped cuprate compounds is briefly discussed.Comment: 20 pages, REVTeX, 6 figures (uuencoded compressed PostScript file for figures is appended

Far-Ultraviolet and Far-Infrared Bivariate Luminosity Function of Galaxies: Complex Relation between Stellar and Dust Emission

Far-ultraviolet (FUV) and far-infrared (FIR) luminosity functions (LFs) of galaxies show a strong evolution from $z = 0$ to $z = 1$, but the FIR LF evolves much stronger than the FUV one. The FUV is dominantly radiated from newly formed short-lived OB stars, while the FIR is emitted by dust grains heated by the FUV radiation field. It is known that dust is always associated with star formation activity. Thus, both FUV and FIR are tightly related to the star formation in galaxies, but in a very complicated manner. In order to disentangle the relation between FUV and FIR emissions, we estimate the UV-IR bivariate LF (BLF) of galaxies with {\sl GALEX} and {\sl AKARI} All-Sky Survey datasets. Recently we invented a new mathematical method to construct the BLF with given marginals and prescribed correlation coefficient. This method makes use of a tool from mathematical statistics, so called "copula". The copula enables us to construct a bivariate distribution function from given marginal distributions with prescribed correlation and/or dependence structure. With this new formulation and FUV and FIR univariate LFs, we analyze various FUV and FIR data with {\sl GALEX}, {\sl Spitzer}, and {\sl AKARI} to estimate the UV-IR BLF. The obtained BLFs naturally explain the nonlinear complicated relation between FUV and FIR emission from star-forming galaxies. Though the faint-end of the BLF was not well constrained for high-$z$ samples, the estimated linear correlation coefficient $\rho$ was found to be very high, and is remarkably stable with redshifts (from 0.95 at $z = 0$ to 0.85 at $z = 1.0$). This implies the evolution of the UV-IR BLF is mainly due to the different evolution of the univariate LFs, and may not be controlled by the dependence structure.Comment: 10 pages, 7 figures, Earth, Planets and Space, in pres

The ODD protocol for describing agent-based and other simulation models: A second update to improve clarity, replication, and structural realism

© 2020, University of Surrey. All rights reserved. The Overview, Design concepts and Details (ODD) protocol for describing Individual-and Agent-Based Models (ABMs) is now widely accepted and used to document such models in journal articles. As a standardized document for providing a consistent, logical and readable account of the structure and dynamics of ABMs, some research groups also find it useful as a workflow for model design. Even so, there are still limitations to ODD that obstruct its more widespread adoption. Such limitations are discussed and addressed in this paper: the limited availability of guidance on how to use ODD; the length of ODD documents; limitations of ODD for highly complex models; lack of sufficient details of many ODDs to enable reimplementation without access to the model code; and the lack of provision for sections in the document structure covering model design ratio-nale, the model’s underlying narrative, and the means by which the model’s fitness for purpose is evaluated. We document the steps we have taken to provide better guidance on: structuring complex ODDs and an ODD summary for inclusion in a journal article (with full details in supplementary material; Table 1); using ODD to point readers to relevant sections of the model code; update the document structure to include sections on model rationale and evaluation. We also further advocate the need for standard descriptions of simulation experiments and argue that ODD can in principle be used for any type of simulation model. Thereby ODD would provide a lingua franca for simulation modelling