Mortality and climatic/environmental correlates of Widdringtonia cedargensis trees in the Cederberg, South Africa

Background: Conifer populations appear disproportionately threatened by global change. Most examples are, however, drawn from the northern hemisphere and long-term rates of population decline are not well documented as historical data are often lacking. We use a large and long-term (1931-2013) repeat photography dataset together with environmental data and fire records to account for the decline of the critically endangered Widdringtonia cedarbergensis. Eighty-seven historical and repeat photo-pairs were analysed to establish 20th century changes in W. cedarbergensis demography. A generalized linear mixed-effects model was fitted to determine the relative importance of environmental factors and fire-return interval on mortality for the species. Results: From an initial total of 1313 live trees in historical photographs, 74% had died and only 44 (3.4%) had recruited in the repeat photographs, leaving 387 live individuals. Juveniles (mature adults) had decreased (increased) from 27% (73%) to 8% (92%) over the intervening period. Our model demonstrates that mortality is related to greater fire frequency, higher temperatures, lower elevations, less rocky habitats and aspect (i.e. east-facing slopes had the least mortality). Conclusions: Our results show that W. cedarbergensis populations have declined significantly over the recorded period, with a pronounced decline in the last 30 years. Individuals that established in open habitats at lower, hotter elevations and experienced a greater fire frequency appear to be more vulnerable to mortality than individuals growing within protected, rocky environments at higher, cooler locations with less frequent fires. Climate models predict increasing temperatures for our study area (and likely increases in wildfires). If these predictions are realised, further declines in the species can be expected. Urgent management interventions, including seedling out-planting in fire-protected high elevation sites, reducing fire frequency in higher elevation populations, and assisted migration, should be considered.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Pacific Ocean

Total organic carbon (TOC) was analyzed on four transects along 140°W in 1992 using a high temperature combustion/discrete injection (HTC/DI) analyzer. For two of the transects, the analyses were conducted on-board ship. Mixed-layer concentrations of organic carbon varied from about 80 µM C at either end of the transect (12°N and 12°S) to about 60 µM C at the equator. Total organic carbon concentrations decreased rapidly below the mixed-layer to about 38-40 µM C at 1000 m across the transect. Little variation was observed below this depth; deep water concentrations below 2000 m were virtually monotonic at about 36 µM C. Repeat measurements made on subsequent cruises consistently found the same concentrations at 1000 m or deeper, but substantial variations were observed in the mixed-layer and the upper water column above 400 m depth. Linear mixing models of total organic carbon versus sigmaT exhibited zones of organic carbon formation and consumption. TOC was found to be inversely correlated with apparent oxygen utilization (AOU) in the region between the mixed-layer and the oxygen minimum. In the mixed-layer, TOC concentrations varied seasonally. Part of the variations in TOC at the equator was driven by changes in the upwelling rate in response to variations in physical forcing related to an El Niño and to the passage of tropical instability waves. TOC export fluxes, calculated from simple box models, averaged 8±4 mmol C/m**2/day at the equator and also varied seasonally. These export fluxes account for 50-75% of the total carbon deficit and are consistent with other estimates and model predictions.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Whale sightings, group sizes and krill biomass in West Greenland in 2005

This study combined data on fin whale Balaenoptera physalus, humpback whale Megaptera novaeangliae, minke whale B. acutorostrata, and sei whale B. borealis sightings from large-scale visual aerial and ship-based surveys (248 and 157 sightings, respectively) with synoptic acoustic sampling of krill Meganyctiphanes norvegica and Thysanoessa sp. abundance in September 2005 in West Greenland to examine the relationships between whales and their prey. Krill densities were obtained by converting relationships of volume backscattering strengths at multiple frequencies to a numerical density using an estimate of krill target strength. Krill data were vertically integrated in 25 m depth bins between 0 and 300 m to obtain water column biomass (g/m**2) and translated to density surfaces using ordinary kriging. Standard regression models (Generalized Additive Modeling, GAM, and Generalized Linear Modeling, GLM) were developed to identify important explanatory variables relating the presence, absence, and density of large whales to the physical and biological environment and different survey platforms. Large baleen whales were concentrated in 3 focal areas: (1) the northern edge of Lille Hellefiske bank between 65 and 67°N, (2) north of Paamiut at 63°N, and (3) in South Greenland between 60 and 61° N. There was a bimodal pattern of mean krill density between depths, with one peak between 50 and 75 m (mean 0.75 g/m**2, SD 2.74) and another between 225 and 275 m (mean 1.2 to 1.3 g/m**2, SD 23 to 19). Water column krill biomass was 3 times higher in South Greenland than at any other site along the coast. Total depth-integrated krill biomass was 1.3 x 10**9 (CV 0.11). Models indicated the most important parameter in predicting large baleen whale presence was integrated krill abundance, although this relationship was only significant for sightings obtained on the ship survey. This suggests that a high degree of spatio-temporal synchrony in observations is necessary for quantifying predator-prey relationships. Krill biomass was most predictive of whale presence at depths >150 m, suggesting a threshold depth below which it is energetically optimal for baleen whales to forage on krill in West Greenland.

(Table 2) Carbon and nitrogen isotopic ratios in blood sampled from arctic foxes (Vulpes lagopus) on Bylot Island between 2003-2008

Inter-individual variation in diet within generalist animal populations is thought to be a widespread phenomenon but its potential causes are poorly known. Inter-individual variation can be amplified by the availability and use of allochthonous resources, i.e., resources coming from spatially distinct ecosystems. Using a wild population of arctic fox as a study model, we tested hypotheses that could explain variation in both population and individual isotopic niches, used here as proxy for the trophic niche. The arctic fox is an opportunistic forager, dwelling in terrestrial and marine environments characterized by strong spatial (arctic-nesting birds) and temporal (cyclic lemmings) fluctuations in resource abundance. First, we tested the hypothesis that generalist foraging habits, in association with temporal variation in prey accessibility, should induce temporal changes in isotopic niche width and diet. Second, we investigated whether within-population variation in the isotopic niche could be explained by individual characteristics (sex and breeding status) and environmental factors (spatiotemporal variation in prey availability). We addressed these questions using isotopic analysis and Bayesian mixing models in conjunction with linear mixed-effects models. We found that: i) arctic fox populations can simultaneously undergo short-term (i.e., within a few months) reduction in both isotopic niche width and inter-individual variability in isotopic ratios, ii) individual isotopic ratios were higher and more representative of a marine-based diet for non-breeding than breeding foxes early in spring, and iii) lemming population cycles did not appear to directly influence the diet of individual foxes after taking their breeding status into account. However, lemming abundance was correlated to proportion of breeding foxes, and could thus indirectly affect the diet at the population scale.

Measurements of gastropods shell sizes from 23 fossil and extant long-lived lakes

Aim: To investigate shell size variation among gastropod faunas of fossil and recent long-lived European lakes and discuss potential underlying processes. Location: 23 long-lived lakes of the Miocene to Recent of Europe. Methods: Based on a dataset of 1412 species of both fossil and extant lacustrine gastropods, we assessed differences in shell size in terms of characteristics of the faunas (species richness, degree of endemism, differences in family composition) and the lakes (surface area, latitude and longitude of lake centroid, distance to closest neighbouring lake) using multiple and linear regression models. Because of a strong species-area relationship, we used resampling to determine whether any observed correlation is driven by that relationship. Results: The regression models indicated size range expansion rather than unidirectional increase or decrease as the dominant pattern of size evolution. The multiple regression models for size range and maximum and minimum size were statistically significant, while the model with mean size was not. Individual contributions and linear regressions indicated species richness and lake surface area as best predictors for size changes. Resampling analysis revealed no significant effects of species richness on the observed patterns. The correlations are comparable across families of different size classes, suggesting a general pattern. Main conclusions: Among the chosen variables, species richness and lake surface area are the most robust predictors of shell size in long-lived lake gastropods. Although the most outstanding and attractive examples for size evolution in lacustrine gastropods derive from lakes with extensive durations, shell size appears to be independent of the duration of the lake as well as longevity of a species. The analogue of long-lived lakes as 'evolutionary islands' does not hold for developments of shell size because different sets of parameters predict size changes.

Phycoerythrin concentration during POLARSTERN cruises ANT-XXIV/4, ANT-XXV/1, and ANT-XXVI/4

Phycobiliproteins are a family of water-soluble pigment proteins that play an important role as accessory or antenna pigments and absorb in the green part of the light spectrum poorly used by chlorophyll a. The phycoerythrins (PEs) are one of four types of phycobiliproteins that are generally distinguished based on their absorption properties. As PEs are water soluble, they are generally not captured with conventional pigment analysis. Here we present a statistical model based on in situ measurements of three transatlantic cruises which allows us to derive relative PE concentration from standardized hyperspectral underwater radiance measurements (Lu). The model relies on Empirical Orthogonal Function (EOF) analysis of Lu spectra and, subsequently, a Generalized Linear Model with measured PE concentrations as the response variable and EOF loadings as predictor variables. The method is used to predict relative PE concentrations throughout the water column and to calculate integrated PE estimates based on those profiles.