The curious case of the camelthorn: competition, coexistence and nest-site limitation in a multispecies mutualism

Myrmecophyte plants house ants in domatia in exchange for protection from herbivores. Ant-myrmecophyte mutualisms exhibit two general patterns due to competition between ants for plant occupancy: i) domatia nest-sites are a limiting resource and ii) each individual plant hosts one ant species at a time. However, individual camelthorn trees (Vachellia erioloba) typically host two to four ant species simultaneously, often coexisting in adjacent domatia on the same branch. Such fine-grain spatial coexistence brings into question the conventional wisdom on ant-myrmecophyte mutualisms. Camelthorn ants appear not to be nest-site limited, despite low abundance of suitable domatia, and have random distributions of nest-sites within and across trees. These patterns suggest a lack of competition between ants for domatia and contrast strongly with other ant-myrmecophyte systems. Comparison of this unusual case with others suggests that spatial scale is crucial to coexistence or competitive exclusion involving multiple ant species. Furthermore, coexistence may be facilitated when co-occurring ant species diverge strongly on at least one niche axis. Our conclusions provide recommendations for future ant-myrmecophyte research, particularly in utilising multispecies systems to further our understanding of mutualism biology.

The EarthCARE satellite: the next step forward in global measurements of clouds, aerosols, precipitation, and radiation

The collective representation within global models of aerosol, cloud, precipitation, and their radiative properties remains unsatisfactory. They constitute the largest source of uncertainty in predictions of climatic change and hamper the ability of numerical weather prediction models to forecast high-impact weather events. The joint European Space Agency (ESA)–Japan Aerospace Exploration Agency (JAXA) Earth Clouds, Aerosol and Radiation Explorer (EarthCARE) satellite mission, scheduled for launch in 2018, will help to resolve these weaknesses by providing global profiles of cloud, aerosol, precipitation, and associated radiative properties inferred from a combination of measurements made by its collocated active and passive sensors. EarthCARE will improve our understanding of cloud and aerosol processes by extending the invaluable dataset acquired by the A-Train satellites CloudSat, Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Aqua. Specifically, EarthCARE’s cloud profiling radar, with 7 dB more sensitivity than CloudSat, will detect more thin clouds and its Doppler capability will provide novel information on convection, precipitating ice particle, and raindrop fall speeds. EarthCARE’s 355-nm high-spectral-resolution lidar will measure directly and accurately cloud and aerosol extinction and optical depth. Combining this with backscatter and polarization information should lead to an unprecedented ability to identify aerosol type. The multispectral imager will provide a context for, and the ability to construct, the cloud and aerosol distribution in 3D domains around the narrow 2D retrieved cross section. The consistency of the retrievals will be assessed to within a target of ±10 W m–2 on the (10 km)2 scale by comparing the multiview broadband radiometer observations to the top-of-atmosphere fluxes estimated by 3D radiative transfer models acting on retrieved 3D domains.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

Impacts of climate change scenarios on dissolved oxygen in the River Thames, UK

A water quality model is used to assess the impact of possible climate change on dissolved oxygen (DO) in the Thames. The Thames catchment is densely populated and, typically, many pressures are anthropogenic. However, that same population also relies on the river for potable water supply and as a disposal route for treated wastewater. Thus, future water quality will be highly dependent on future activity. Dynamic and stochastic modelling has been used to assess the likely impacts on DO dynamics along the river system and the probability distributions associated with future variability. The modelling predictions indicate that warmer river temperatures and drought act to reduce dissolved oxygen concentrations in lowland river systems

Extended magnetic reconnection across the dayside magnetopause

The extent of where magnetic reconnection (MR), the dominant process responsible for energy and plasma transport into the magnetosphere, operates across Earth’s dayside magnetopause has previously been only indirectly shown by observations. We report the first direct evidence of X-line structure resulting from the operation of MR at each of two widely separated locations along the tilted, subsolar line of maximum current on Earth’s magnetopause, confirming the operation of MR at two or more sites across the extended region where MR is expected to occur. The evidence results from in-situ observations of the associated ion and electron plasma distributions, present within each magnetic X-line structure, taken by two spacecraft passing through the active MR regions simultaneously.

Enhanced viability of corneal epithelial cells for efficient transport/storage using a structurally-modified calcium alginate hydrogel

Aims: Therapeutic limbal epithelial stem cells could be managed more efficiently if clinically validated batches were transported for ‘on-demand’ use. Materials & methods: In this study, corneal epithelial cell viability in calcium alginate hydrogels was examined under cell culture, ambient and chilled conditions for up to 7 days. Results: Cell viability improved as gel internal pore size increased, and was further enhanced with modification of the gel from a mass to a thin disc. Ambient storage conditions were optimal for supporting cell viability in gel discs. Cell viability in gel discs was significantly enhanced with increases in pore size mediated by hydroxyethyl cellulose. Conclusion: Our novel methodology of controlling alginate gel shape and pore size together provides a more practical and economical alternative to established corneal tissue/cell storage methods.

Inner plasma structure of the low-latitude reconnection layer

We report a clear transition through a reconnection layer at the low-latitude magnetopause which shows a complete traversal across all reconnected field lines during northwestward interplanetary magnetic field (IMF) conditions. The associated plasma populations confirm details of the electron and ion mixing and the time history and acceleration through the current layer. This case has low magnetic shear with a strong guide field and the reconnection layer contains a single density depletion layer on the magnetosheath side which we suggest results from nearly field-aligned magnetosheath flows. Within the reconnection boundary layer, there are two plasma boundaries, close to the inferred separatrices on the magnetosphere and magnetosheath sides (Ssp and Ssh) and two boundaries associated with the Alfvén waves (or Rotational Discontinuities, RDsp and RDsh). The data are consistent with these being launched from the reconnection site and the plasma distributions are well ordered and suggestive of the time elapsed since reconnection of the field lines observed. In each sub-layer between the boundaries the plasma distribution is different and is centered around the current sheet, responsible for magnetosheath acceleration. We show evidence for a velocity dispersion effect in the electron anisotropy that is consistent with the time elapsed since reconnection. In addition, new evidence is presented for the occurrence of partial reflection of magnetosheath electrons at the magnetopause current layer.

Estimating root–soil contact from 3D X-ray microtomographs

Adequate contact with the soil is essential for water and nutrient adsorption by plant roots, but the determination of root–soil contact is a challenging task because it is difficult to visualize roots in situ and quantify their interactions with the soil at the scale of micrometres. A method to determine root–soil contact using X-ray microtomography was developed. Contact areas were determined from 3D volumetric images using segmentation and iso-surface determination tools. The accuracy of the method was tested with physical model systems of contact between two objects (phantoms). Volumes, surface areas and contact areas calculated from the measured phantoms were compared with those estimated from image analysis. The volume was accurate to within 0.3%, the surface area to within 2–4%, and the contact area to within 2.5%. Maize and lupin roots were grown in soil (<2 mm) and vermiculite at matric potentials of −0.03 and −1.6 MPa and in aggregate fractions of 4–2, 2–1, 1–0.5 and < 0.5 mm at a matric potential of −0.03 MPa. The contact of the roots with their growth medium was determined from 3D volumetric images. Macroporosity (>70 µm) of the soil sieved to different aggregate fractions was calculated from binarized data. Root-soil contact was greater in soil than in vermiculite and increased with decreasing aggregate or particle size. The differences in root–soil contact could not be explained solely by the decrease in porosity with decreasing aggregate size but may also result from changes in particle and aggregate packing around the root.

Classification of cortical microcircuits based on micro-electrode-array data from slices of rat barrel cortex

The bewildering complexity of cortical microcircuits at the single cell level gives rise to surprisingly robust emergent activity patterns at the level of laminar and columnar local field potentials (LFPs) in response to targeted local stimuli. Here we report the results of our multivariate data-analytic approach based on simultaneous multi-site recordings using micro-electrode-array chips for investigation of the microcircuitary of rat somatosensory (barrel) cortex. We find high repeatability of stimulus-induced responses, and typical spatial distributions of LFP responses to stimuli in supragranular, granular, and infragranular layers, where the last form a particularly distinct class. Population spikes appear to travel with about 33 cm/s from granular to infragranular layers. Responses within barrel related columns have different profiles than those in neighbouring columns to the left or interchangeably to the right. Variations between slices occur, but can be minimized by strictly obeying controlled experimental protocols. Cluster analysis on normalized recordings indicates specific spatial distributions of time series reflecting the location of sources and sinks independent of the stimulus layer. Although the precise correspondences between single cell activity and LFPs are still far from clear, a sophisticated neuroinformatics approach in combination with multi-site LFP recordings in the standardized slice preparation is suitable for comparing normal conditions to genetically or pharmacologically altered situations based on real cortical microcircuitry.

A Gaussian-mixture ensemble transform filter

We generalize the popular ensemble Kalman filter to an ensemble transform filter, in which the prior distribution can take the form of a Gaussian mixture or a Gaussian kernel density estimator. The design of the filter is based on a continuous formulation of the Bayesian filter analysis step. We call the new filter algorithm the ensemble Gaussian-mixture filter (EGMF). The EGMF is implemented for three simple test problems (Brownian dynamics in one dimension, Langevin dynamics in two dimensions and the three-dimensional Lorenz-63 model). It is demonstrated that the EGMF is capable of tracking systems with non-Gaussian uni- and multimodal ensemble distributions. Copyright © 2011 Royal Meteorological Society

Explanations of the inconsistencies in survey respondents' forecasts

A comparison of the point forecasts and the probability distributions of inflation and output growth made by individual respondents to the US Survey of Professional Forecasters indicates that the two sets of forecasts are sometimes inconsistent. We evaluate a number of possible explanations, and find that not all forecasters update their histogram forecasts as new information arrives. This is supported by the finding that the point forecasts are more accurate than the histograms in terms of first-moment prediction.

Experimentally driven atomistic model of 1,2 Polybutadiene

We present an efficient method of combining wide angle neutron scattering data with detailed atomistic models, allowing us to perform a quantitative and qualitative mapping of the organisation of the chain conformation in both glass and liquid phases. The structural refinement method presented in this work is based on the exploitation of the intrachain features of the diffraction pattern and its intimate linkage with atomistic models by the use of internal coordinates for bond lengths, valence angles and torsion rotations. Atomic connectivity is defined through these coordinates that are in turn assigned by pre-defined probability distributions, thus allowing for the models in question to be built stochastically. Incremental variation of these coordinates allows for the construction of models that minimise the differences between the observed and calculated structure factors. We present a series of neutron scattering data of 1,2 polybutadiene at the region 120-400K. Analysis of the experimental data yield bond lengths for C-C and C=C of 1.54Å and 1.35Å respectively. Valence angles of the backbone were found to be at 112° and the torsion distributions are characterised by five rotational states, a three-fold trans-skew± for the backbone and gauche± for the vinyl group. Rotational states of the vinyl group were found to be equally populated, indicating a largely atactic chan. The two backbone torsion angles exhibit different behaviour with respect to temperature of their trans population, with one of them adopting an almost all trans sequence. Consequently the resulting configuration leads to a rather persistent chain, something indicated by the value of the characteristic ratio extrapolated from the model. We compare our results with theoretical predictions, computer simulations, RIS models and previously reported experimental results.

Quantifying the importance of galactic cosmic rays in cloud microphysical processes

Galactic Cosmic Rays are one of the major sources of ion production in the troposphere and stratosphere. Recent studies have shown that ions form electrically charged clusters which may grow to become cloud droplets. Aerosol particles charge by the attachment of ions and electrons. The collision efficiency between a particle and a water droplet increases, if the particle is electrically charged, and thus aerosol-cloud interactions can be enhanced. Because these microphysical processes may change radiative properties of cloud and impact Earth's climate it is important to evaluate these processes' quantitative effects. Five different models developed independently have been coupled to investigate this. The first model estimates cloud height from dew point temperature and the temperature profile. The second model simulates the cloud droplet growth from aerosol particles using the cloud parcel concept. In the third model, the scavenging rate of the aerosol particles is calculated using the collision efficiency between charged particles and droplets. The fourth model calculates electric field and charge distribution on water droplets and aerosols within cloud. The fifth model simulates the global electric circuit (GEC), which computes the conductivity and ionic concentration in the atmosphere in altitude range 0–45 km. The first four models are initially coupled to calculate the height of cloud, boundary condition of cloud, followed by growth of droplets, charge distribution calculation on aerosols and cloud droplets and finally scavenging. These models are incorporated with the GEC model. The simulations are verified with experimental data of charged aerosol for various altitudes. Our calculations showed an effect of aerosol charging on the CCN concentration within the cloud, due to charging of aerosols increase the scavenging of particles in the size range 0.1 µm to 1 µm.

Flux transfer events at the dayside magnetopause: transient reconnection or magnetosheath dynamic pressure pulses?

The suggestion is discussed that characteristic particle and field signatures at the dayside magnetopause, termed “flux transfer events” (FTEs), are, in at least some cases, due to transient solar wind and/or magnetosheath dynamic pressure increases, rather than time-dependent magnetic reconnection. It is found that most individual cases of FTEs observed by a single spacecraft can, at least qualitatively, be explained by the pressure pulse model, provided a few rather unsatisfactory features of the predictions are explained in terms of measurement uncertainties. The most notable exceptions to this are some “two-regime” observations made by two satellites simultaneously, one on either side of the magnetopause. However, this configuration has not been frequently achieved for sufficient time, such observations are rare, and the relevant tests are still not conclusive. The strongest evidence that FTEs are produced by magnetic reconnection is the dependence of their occurrence on the north-south component of the interplanetary magnetic field (IMF) or of the magnetosheath field. The pressure pulse model provides an explanation for this dependence (albeit qualitative) in the case of magnetosheath FTEs, but this does not apply to magnetosphere FTEs. The only surveys of magnetosphere FTEs have not employed the simultaneous IMF, but have shown that their occurrence is strongly dependent on the north-south component of the magnetosheath field, as observed earlier/later on the same magnetopause crossing (for inbound/outbound passes, respectively). This paper employs statistics on the variability of the IMF orientation to investigate the effects of IMF changes between the times of the magnetosheath and FTE observations. It is shown that the previously published results are consistent with magnetospheric FTEs being entirely absent when the magnetosheath field is northward: all crossings with magnetosphere FTEs and a northward field can be attributed to the field changing sense while the satellite was within the magnetosphere (but close enough to the magnetopause to detect an FTE). Allowance for the IMF variability also makes the occurrence frequency of magnetosphere FTEs during southward magnetosheath fields very similar to that observed for magnetosheath FTEs. Conversely, the probability of attaining the observed occurrence frequencies for the pressure pulse model is 10−14. In addition, it is argued that some magnetosheath FTEs should, for the pressure pulse model, have been observed for northward IMF: the probability that the number is as low as actually observed is estimated to be 10−10. It is concluded that although the pressure model can be invoked to qualitatively explain a large number of individual FTE observations, the observed occurrence statistics are in gross disagreement with this model.

Variability of the interplanetary medium at 1 A.U. over 24 years: 1963–1986

A survey is presented of hourly averages of observations of the interplanetary medium, made by satellites close to the Earth (i.e. at l a.u.) in the years 1963-1986. This survey therefore covers two complete solar cycles (numbers 20 and 21). The distributions and solar-cycle variations of IMF field strength, B, and its northward component (in GSM coordinates), B(z), and of the solar-wind density, n, speed, v, and dynamic pressure, P, are discussed. Because of their importance to the terrestrial magnetosphere/ionosphere, particular attention is given to B(z) and P. The solar-cycle variation in the magnitude and variability of B(z) previously reported for cycle 20, is also found for cycle 21. However, the solar-wind data show a number of differences between cycles 20 and 21. The average dynamic pressure is found to show a solar-cycle variation and a systematic increase over the period of the survey. The minimum of dynamic pressure at sunspot maximum is mainly due to reduced solar-wind densities in cycle 20, but lower solar-wind speed in cycle 21 is a more significant factor. The distribution of the duration of periods of stable polarity of the IMF B(z) component shows that the magnetosphere could achieve steady state for only a small fraction of the time and there is some evidence for a solar-cycle variation in this fraction. It is also found that the polarity changes in the IMF B(z) fall into two classes: one with an associated change in solar-wind dynamic pressure, the other without such a change. However, in only 20% of cases does the dynamic pressure change exceed 50%.