A statistical approach to evaluating distance metrics and analog assignments for pollen records

The modern analog technique typically uses a distance metric to determine the dissimilarity between fossil and modern biological assemblages. Despite this quantitative approach, interpretation of distance metrics is usually qualitative and rules for selection of analogs tend to be ad hoc. We present a statistical tool, the receiver operating characteristic (ROC) curve, which provides a framework for identifying analogs from distance metrics. If modern assemblages are placed into groups (e.g., biomes), this method can (1) evaluate the ability of different distance metrics to distinguish among groups, (2) objectively identify thresholds of the distance metric for determining analogs, and (3) compute a likelihood ratio and a Bayesian probability that a modern group is an analog for an unknown (fossil) assemblage. Applied to a set of 1689 modern pollen assemblages from eastern North America classified into eight biomes, ROC analysis confirmed that the squared-chord distance (SCD) outperforms most other distance metrics. The optimal threshold increased when more dissimilar biomes were compared. The probability of an analog vs no-analog result (a likelihood ratio) increased sharply when SCD decreased below the optimal threshold, indicating a nonlinear relationship between SCD and the probability of analog. Probabilities of analog computed for a postglacial pollen record at Tannersville Bog (Pennsylvania, USA) identified transitions between biomes and periods of no analog.     

Resolving mechanisms of toxicity while pursuing ecotoxicological relevance?

In this age of modern biology, aquatic toxicological research has pursued mechanisms of action of toxicants. This has provided potential tools for ecotoxicologic investigations. However, problems of biocomplexity and issues at higher levels of biological organization remain a challenge. In the 1980s and 1990s and continuing to a lesser extent today, organisms residing in highly contaminated field sites or exposed in the laboratory to calibrated concentrations of individual compounds were carefully analyzed for their responses to priority pollutants. Correlation of biochemical and structural analyses in cultured cells and tissues, as well as the in vivo exposures led to the production and application of biomarkers of exposure and effect and to our awareness of genotoxicity and its chronic manifestations, such as neoplasms, in wild fishes. To gain acceptance of these findings in the greater environmental toxicology community, “validation of the model” versus other, better-established often rodent models, was necessary and became a major focus. Resultant biomarkers were applied to heavily contaminated and reference field sites as part of effects assessment and with investigations following large-scale disasters such as oil spills or industrial accidents.

Over the past 15 years, in the laboratory, small aquarium fish models such as medaka (Oryzias latipes), zebrafish (Danio rerio), platyfish (Xiphophorus species), fathead minnow (Pimephales promelas), and sheepshead minnow (Cyprinodon variegatus) were increasingly used establishing mechanisms of toxicants. Today, the same organisms provide reliable information at higher levels of biological organization relevant to ecotoxicology. We review studies resolving mechanisms of toxicity and discuss ways to address biocomplexity, mixtures of contaminants, and the need to relate individual level responses to populations and communities.     

Searches for very high energy gamma rays from blazars with CANGAROO-III telescope in 2005-2009

We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data, even though they are non-simultaneous, are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size.     

A Refined Model for Solid Particle Rock Erosion

A procedure for the estimation of distribution parameters of a Weibull distribution model K₁ = f(K_Ic^12/4/σ_C^23/4) for solid particle erosion, as recently suggested in Rock Mech Rock Eng, doi: 10.1007/s00603-014-0658-x, 2014, is derived. The procedure is based on examinations of elastic–plastically responding rocks (rhyolite, granite) and plastically responding rocks (limestone, schist). The types of response are quantified through SEM inspections of eroded surfaces. Quantitative numbers for the distribution parameter K₁ are calculated for 30 rock materials, which cover a wide range of mechanical properties. The ranking according to the parameter K₁ is related to qualitative rock classification schemes. A modified proposal for the erosion of schist due to solid particle impingement at normal incidence is introduced.     

The oxygen isotopic composition of the fraction less than 2 microns (clay fraction) in Recent sediments from western Europe

The oxygen isotopic composition of the clay fraction in Recent argillaceous sediments in western Europe shows regional differences, the δ¹⁸_SMOW values varying between + 15.8 and + 21.7‰. This phenomenon is used for sediment-transport studies. In the Rhine-Meuse estuary a landward transport of marine sediments is observed. The clay fraction of the sediments in the Wadden Sea is only for a minor part derived from the rivers Rhine and Meuse.     

Faunal relationships between the New Zealand plateau and the New Zealand sector of Antarctica based on echinoderm distribution

The problem of New Zealand‐Antarctic faunal relationships is discussed on the basis of echinoderm distribution. The limitations of the data upon which past zoogeographical speculations have been based are pointed out.

Macquarie Island (occupying an intermediate geographical position between the New Zealand Plateau and the Antarctic) shows definite relationships with New Zealand, and the submarine Macquarie Ridge may have provided a connecting migration route. However, only four species of echinoderms are shared between the Ross Sea‐Balleny Islands area (the New Zealand sector of the Antarctic Region) and the New Zealand Pliateau‐Macquarie Island area (the New Zealand Region). Although the as yet unsampled part of the Macquarie‐Balleny Ridge may reveal other faunal similarities, the present systematic sampling has made possible a sounder understanding of the zoogeographical affinities of the two regions.     

Bizarre geochemistry of komatiites from the Crixás greenstone belt,Brazil

Komatiite lava flows in the Crixás greenstone belt, Goiás, Brazil, have textures and volcanic structures typical of Archean komatiites, but are geochemically most unusual. The flows are porphyritic and massive, or layered with spinifex upper parts and olivine cumulate lower parts. MgO contents range from 18 to 40%. In such lavas, only olivine (and minor chromite) can have crystallized, but neither major nor trace elements fall on olivine control lines. In MgO variation diagrams, CaO and Sr fall on lines with slopes steeper than olivine control lines; SiO₂, FeO, Na₂O, K₂O and Y show little systematic variation; Zr shows a large variation that does not correlate with MgO; and Al₂O₃ decreases markedly with decreasing MgO. The aberrant behaviour is highlighted by the REE (rare earth elements) in spinifex and olivine cumulate layers from three flows: in the spinifex layers, chondrite-normalized REE patterns are hump-shaped with maxima at Nd or Sm ((La/Sm)_N=0.6, (Gd/Yb)_N=1.6–2.3), whereas cumulate zones in the same flows have steadily sloping patterns, with LREE enriched relative to HREE ((La/Sm)_N=1.3, (Gd/Yb)_N=1.4). Neither normal magmatic processes acting during emplacement of the komatiites, nor thermal erosion and wall-rock assimilation can explain these effects, and we speculate that elements commonly thought of as “immobile” (e.g. Al, Zr, REE) migrated during hydrothermal alteration or metamorphism. A Pb-Pb whole rock isochron gave an age of 2,728±140 Ma and selected Sm-Nd analyses an apparent isochron age of 2,825±98 Ma (ɛ_Nd≈0). The Pb-Pb age is believed to be the approximate time of emplacement. Interpretation of the Sm-Nd data is complicated by the evidence of mobility of REE.     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.