Extreme value distribution for singular measures

In this paper we perform an analytical and numerical study of Extreme Value distributions in discrete dynamical systems that have a singular measure. Using the block maxima approach described in Faranda et al. [2011] we show that, numerically, the Extreme Value distribution for these maps can be associated to the Generalised Extreme Value family where the parameters scale with the information dimension. The numerical analysis are performed on a few low dimensional maps. For the middle third Cantor set and the Sierpinskij triangle obtained using Iterated Function Systems, experimental parameters show a very good agreement with the theoretical values. For strange attractors like Lozi and H\`enon maps a slower convergence to the Generalised Extreme Value distribution is observed. Even in presence of large statistics the observed convergence is slower if compared with the maps which have an absolute continuous invariant measure. Nevertheless and within the uncertainty computed range, the results are in good agreement with the theoretical estimates.

The causal surface and its effect on distribution transparency in a distributed virtual environment

One goal in the development of distributed virtual environments (DVEs) is to create a system such that users are unaware of the distribution-the distribution should be transparent. The paper begins by discussing the general issues in DVEs that might make this possible, and a system that allows some level of distribution transparency is described. The system described suffers from effects of inconsistency, which in turn cause undesirable visual effects. The causal surface is introduced as a solution that removes these visual effects. The paper then introduces two determining factors of distribution transparency relating to user perception and performance. With regard to these factors, two hypotheses are stated relating to the causal surface. A user-trial on forty-five subjects is used to validate the hypotheses. A discussion of the results of the trial concludes that the causal surface solution does significantly improve the distribution transparency in a DVE.

Strengthening the link between climate, hydrological and species distribution modeling to assess the impacts of climate change on freshwater biodiversity

To understand the resilience of aquatic ecosystems to environmental change, it is important to determine how multiple, related environmental factors, such as near-surface air temperature and river flow, will change during the next century. This study develops a novel methodology that combines statistical downscaling and fish species distribution modeling, to enhance the understanding of how global climate changes (modeled by global climate models at coarse-resolution) may affect local riverine fish diversity. The novelty of this work is the downscaling framework developed to provide suitable future projections of fish habitat descriptors, focusing particularly on the hydrology which has been rarely considered in previous studies. The proposed modeling framework was developed and tested in a major European system, the Adour-Garonne river basin (SW France, 116,000 km(2)), which covers distinct hydrological and thermal regions from the Pyrenees to the Atlantic coast. The simulations suggest that, by 2100, the mean annual stream flow is projected to decrease by approximately 15% and temperature to increase by approximately 1.2 °C, on average. As consequence, the majority of cool- and warm-water fish species is projected to expand their geographical range within the basin while the few cold-water species will experience a reduction in their distribution. The limitations and potential benefits of the proposed modeling approach are discussed. Copyright © 2012 Elsevier B.V. All rights reserved.

The role of the international cocoa germplasm database and the international cocoa quarantine centre in information management and distribution of cocoa genetic resources

A range of physiological parameters (canopy light transmission, canopy shape, leaf size, flowering and flushing intensity) were measured from the International Clone Trial, typically over the course of two years. Data were collected from six locations, these being: Brazil, Ecuador, Trinidad, Venezuela, Côte d’Ivoire and Ghana. Canopy shape varied significantly between clones, although it showed little variation between locations. Genotypic variation in leaf size was differentially affected by the growth location; such differences appeared to underlie a genotype by environment interaction in relation to canopy light transmission. Flushing data were recorded at monthly intervals over the course of a year. Within each location, a significant interaction was observed between genotype and time of year, suggesting that some genotypes respond to a greater extent than others to environmental stimuli. A similar interaction was observed for flowering data, where significant correlations were found between flowering intensity and temperature in Brazil and flowering intensity and rainfall in Côte d’Ivoire. The results demonstrate the need for local evaluation of cocoa clones and also suggest that the management practices for particular planting material may need to be fine-tuned to the location in which they are cultivated.

Volume-based size distribution of accumulation and coarse particles (PM0.1–10) from cooking fume during oil heating

Particulate matter generated during the cooking process has been identified as one of the major problems of indoor air quality and indoor environmental health. Reliable assessment of exposure to cooking-generated particles requires accurate information of emission characteristics especially the size distribution. This study characterizes the volume/mass-based size distribution of the fume particles at the oil-heating stage for the typical Chinese-style cooking in a laboratory kitchen. A laser-diffraction size analyzer is applied to measure the volume frequency of fume particles ranged from 0.1 to 10 μm, which contribute to most mass proportion in PM2.5 and PM10. Measurements show that particle emissions have little dependence on the types of vegetable oil used but have a close relationship with the heating temperature. It is found that volume frequency of fume particles in the range of 1.0–4.0 μm accounts for nearly 100% of PM0.1–10 with the mode diameter 2.7 μm, median diameter 2.6 μm, Sauter mean diameter 3.0 μm, DeBroukere mean diameter 3.2 μm, and distribution span 0.48. Such information on emission characteristics obtained in this study can be possibly used to improve the assessment of indoor air quality due to PM0.1–10 in the kitchen and residential flat.

The long-term stability of a possible aqueous ammonium sulfate ocean inside Titan

We model the thermal evolution of a subsurface ocean of aqueous ammonium sulfate inside Titan using a parameterized convection scheme. The cooling and crystallization of such an ocean depends on its heat flux balance, and is governed by the pressure-dependent melting temperatures at the top and bottom of the ocean. Using recent observations and previous experimental data, we present a nominal model which predicts the thickness of the ocean throughout the evolution of Titan; after 4.5 Ga we expect an aqueous ammonium sulfate ocean 56 km thick, overlain by a thick (176 km) heterogeneous crust of methane clathrate, ice I and ammonium sulfate. Underplating of the crust by ice I will give rise to compositional diapirs that are capable of rising through the crust and providing a mechanism for cryovolcanism at the surface. We have conducted a parameter space survey to account for possible variations in the nominal model, and find that for a wide range of plausible conditions, an ocean of aqueous ammonium sulfate can survive to the present day, which is consistent with the recent observations of Titan's spin state from Cassini radar data [Lorenz, R.D., Stiles, B.W., Kirk, R.L., Allison, M.D., del Marmo, P.P., Iess, L., Lunine, J.I., Ostro, S.J., Hensley, S., 2008. Science 319, 1649–1651].

From ensemble forecasts to predictive distribution functions

The translation of an ensemble of model runs into a probability distribution is a common task in model-based prediction. Common methods for such ensemble interpretations proceed as if verification and ensemble were draws from the same underlying distribution, an assumption not viable for most, if any, real world ensembles. An alternative is to consider an ensemble as merely a source of information rather than the possible scenarios of reality. This approach, which looks for maps between ensembles and probabilistic distributions, is investigated and extended. Common methods are revisited, and an improvement to standard kernel dressing, called ‘affine kernel dressing’ (AKD), is introduced. AKD assumes an affine mapping between ensemble and verification, typically not acting on individual ensemble members but on the entire ensemble as a whole, the parameters of this mapping are determined in parallel with the other dressing parameters, including a weight assigned to the unconditioned (climatological) distribution. These amendments to standard kernel dressing, albeit simple, can improve performance significantly and are shown to be appropriate for both overdispersive and underdispersive ensembles, unlike standard kernel dressing which exacerbates over dispersion. Studies are presented using operational numerical weather predictions for two locations and data from the Lorenz63 system, demonstrating both effectiveness given operational constraints and statistical significance given a large sample.

Constraints on the spectral distribution of energy and enstrophy dissipation in forced two-dimensional turbulence

We study two-dimensional (2D) turbulence in a doubly periodic domain driven by a monoscale-like forcing and damped by various dissipation mechanisms of the form νμ(−Δ)μ. By “monoscale-like” we mean that the forcing is applied over a finite range of wavenumbers kmin≤k≤kmax, and that the ratio of enstrophy injection η≥0 to energy injection ε≥0 is bounded by kmin2ε≤η≤kmax2ε. Such a forcing is frequently considered in theoretical and numerical studies of 2D turbulence. It is shown that for μ≥0 the asymptotic behaviour satisfies ∥u∥12≤kmax2∥u∥2, where ∥u∥2 and ∥u∥12 are the energy and enstrophy, respectively. If the condition of monoscale-like forcing holds only in a time-mean sense, then the inequality holds in the time mean. It is also shown that for Navier–Stokes turbulence (μ=1), the time-mean enstrophy dissipation rate is bounded from above by 2ν1kmax2. These results place strong constraints on the spectral distribution of energy and enstrophy and of their dissipation, and thereby on the existence of energy and enstrophy cascades, in such systems. In particular, the classical dual cascade picture is shown to be invalid for forced 2D Navier–Stokes turbulence (μ=1) when it is forced in this manner. Inclusion of Ekman drag (μ=0) along with molecular viscosity permits a dual cascade, but is incompatible with the log-modified −3 power law for the energy spectrum in the enstrophy-cascading inertial range. In order to achieve the latter, it is necessary to invoke an inverse viscosity (μ<0). These constraints on permissible power laws apply for any spectrally localized forcing, not just for monoscale-like forcing.

Inhalation exposure to particulate matter in rooms with underfloor air distribution

This work investigated the personal exposure to indoor particulate matters using the intake fraction metric and provided a possible way to trace the particle inhaled from an indoor particle source. A turbulence model validated by the particle measurements in a room with underfloor air distribution (UFAD) system was used to predict the indoor particle concentrations. Inhalation intake fraction of indoor particles was defined and evaluated in two rooms equipped with the UFAD, i.e., the experimental room and a small office. According to the exposure characteristics and a typical respiratory rate, the intake fraction was determined in two rooms with a continuous and episodic (human cough) source of particles, respectively. The findings showed that the well-mixing assumption of indoor air failed to give an accurate estimation of inhalation exposure and the average concentration at return outlet or within the overall room could not relate well the intake fraction to the amount of particle emitted from an indoor source.

When long-range zero-lag synchronization is feasible in cortical networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (~ 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Distribution and availability of trace elements in municipal solid waste composts

Trace element contamination is one of the main problems linked to the quality of compost, especially when it is produced from urban wastes, which can lead to high levels of some potentially toxic elements such as Cu, Pb or Zn. In this work, the distribution and bioavailability of five elements (Cu, Zn, Pb, Cr and Ni) were studied in five Spanish composts obtained from different feedstocks (municipal solid waste, garden trimmings, sewage sludge and mixed manure). The five composts showed high total concentrations of these elements, which in some cases limited their commercialization due to legal imperatives. First, a physical fractionation of the composts was performed, and the five elements were determined in each size fraction. Their availability was assessed by several methods of extraction (water, CaCl2–DTPA, the PBET extract, the TCLP extract, and sodium pyrophosphate), and their chemical distribution was assessed using the BCR sequential extraction procedure. The results showed that the finer fractions were enriched with the elements studied, and that Cu, Pb and Zn were the most potentially problematic ones, due to both their high total concentrations and availability. The partition into the BCR fractions was different for each element, but the differences between composts were scarce. Pb was evenly distributed among the four fractions defined in the BCR (soluble, oxidizable, reducible and residual); Cu was mainly found in the oxidizable fraction, linked to organic matter, and Zn was mainly associated to the reducible fraction (iron oxides), while Ni and Cr were mainly present almost exclusively in the residual fraction. It was not possible to establish a univocal relation between trace elements availability and their BCR fractionation. Given the differences existing for the availability and distribution of these elements, which not always were related to their total concentrations, we think that legal limits should consider availability, in order to achieve a more realistic assessment of the risks linked to compost use.

Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim  Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location  Europe. Methods  We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000 yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results  Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions  The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Pollen-based reconstructions of Japanese biomes at 0, 6000 and 18,000 14C yr bp

A biomization method, which objectively assigns individual pollen assemblages to biomes ( Prentice et al., 1996 ), was tested using modern pollen data from Japan and applied to fossil pollen data to reconstruct palaeovegetation patterns 6000 and 18,000 14C yr bp Biomization started with the assignment of 135 pollen taxa to plant functional types (PFTs), and nine possible biomes were defined by specific combinations of PFTs. Biomes were correctly assigned to 54% of the 94 modern sites. Incorrect assignments occur near the altitudinal limits of individual biomes, where pollen transport from lower altitudes blurs the local pollen signals or continuous changes in species composition characterizes the range limits of biomes. As a result, the reconstructed changes in the altitudinal limits of biomes at 6000 and 18,000 14C yr bp are likely to be conservative estimates of the actual changes. The biome distribution at 6000 14C yr bp was rather similar to today, suggesting that changes in the bioclimate of Japan have been small since the mid-Holocene. At 18,000 14C yr bp the Japanese lowlands were covered by taiga and cool mixed forests. The southward expansion of these forests and the absence of broadleaved evergreen/warm mixed forests reflect a pronounced year-round cooling.

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.