Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial–interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia – fossil records, species distribution models and phylogeographic surveys – in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Conservation of threatened habitat types under future climate change – Lessons from plant-distribution models and current extinction trends in southern Germany

A higher risk of future range losses as a result of climate change is expected to be one of the main drivers of extinction trends in vascular plants occurring in habitat types of high conservation value. Nevertheless, the impact of the climate changes of the last 60 years on the current distribution and extinction patterns of plants is still largely unclear. We applied species distribution models to study the impact of environmental variables (climate, soil conditions, land cover, topography), on the current distribution of 18 vascular plant species characteristic of three threatened habitat types in southern Germany: (i) xero-thermophilous vegetation, (ii) mesophilous mountain grasslands (mountain hay meadows and matgrass communities), and (iii) wetland habitats (bogs, fens, and wet meadows). Climate and soil variables were the most important variables affecting plant distributions at a spatial level of 10 × 10 km. Extinction trends in our study area revealed that plant species which occur in wetland habitats faced higher extinction risks than those in xero-thermophilous vegetation, with the risk for species in mesophilous mountain grasslands being intermediary. For three plant species characteristic either of mesophilous mountain grasslands or wetland habitats we showed exemplarily that extinctions from 1950 to the present day have occurred at the edge of the species’ current climatic niche, indicating that climate change has likely been the main driver of extinction. This is largely consistent with current extinction trends reported in other studies. Our study indicates that the analysis of past extinctions is an appropriate means to assess the impact of climate change on species and that vulnerability to climate change is both species- and habitat-specific.

Counterbalancing effects of competition for resources and facilitation against grazing in alpine snowbed communities

Alpine snowbeds are habitats where the major limiting factors for plant growth are herbivory and a small time window for growth due to late snowmelt. Despite these limitations, snowbed vegetation usually forms a dense carpet of palatable plants due to favourable abiotic conditions for plant growth within the short growing season. These environmental characteristics make snowbeds particularly interesting to study the interplay of facilitation and competition. We hypothesised an interplay between resource competition and facilitation against herbivory. Further, we investigated whether these predicted neighbour effects were species-specific and/or dependent on ontogeny, and whether the balance of positive and negative plant–plant interactions shifted along a snowmelt gradient. We determined the neighbour effects by means of neighbour removal experiments along the snowmelt gradient, and linear mixed model analyses. The results showed that the effects of neighbour removal were weak but generally consistent among species and snowmelt dates, and depended on whether biomass production or survival was considered. Higher total biomass and increased fruiting in removal plots indicated that plants competed for nutrients, water, and light, thereby supporting the hypothesis of prevailing competition for resources in snowbeds. However, the presence of neighbours reduced herbivory and thereby also facilitated survival. For plant growth the facilitative effects against herbivores in snowbeds counterbalanced competition for resources, leading to a weak negative net effect. Overall the neighbour effects were not species-specific and did not change with snowmelt date. Our finding of counterbalancing effects of competition and facilitation within a plant community is of special theoretical value for species distribution models and can explain the success of models that give primary importance to abiotic factors and tend to overlook interrelations between biotic and abiotic effects on plants.

Correspondence establishment in statistical modeling of shapes with arbitrary topology

Correspondence establishment is a key step in statistical shape model building. There are several automated methods for solving this problem in 3D, but they usually can only handle objects with simple topology, like that of a sphere or a disc. We propose an extension to correspondence establishment over a population based on the optimization of the minimal description length function, allowing considering objects with arbitrary topology. Instead of using a fixed structure of kernel placement on a sphere for the systematic manipulation of point landmark positions, we rely on an adaptive, hierarchical organization of surface patches. This hierarchy can be built on surfaces of arbitrary topology and the resulting patches are used as a basis for a consistent, multi-scale modification of the surfaces' parameterization, based on point distribution models. The feasibility of the approach is demonstrated on synthetic models with different topologies.

Isolation and characterization of nine polymorphic microsatellite markers for the deep-sea shrimp Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea)

Background: The shrimp Nematocarcinus lanceopes Bate, 1888 is found in the deep sea around Antarctica and sub-Antarctic islands. Previous studies on mitochondrial data and species distribution models provided evidence for a homogenous circum-Antarctic population of N. lanceopes. However, to analyze the fine-scale population genetic structure and to examine influences of abiotic environmental conditions on population composition and genetic diversity, a set of fast evolving nuclear microsatellite markers is required. Findings: We report the isolation and characterization of nine polymorphic microsatellite markers from the Antarctic deep-sea shrimp species Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea). Microsatellite markers were screened in 55 individuals from different locations around the Antarctic continent. All markers were polymorphic with 9 to 25 alleles per locus. The observed heterozygosity ranged from 0.545 to 0.927 and the expected heterozygosity from 0.549 to 0.934. Conclusions: The reported markers provide a novel tool to study genetic structure and diversity in Nematocarcinus lanceopes populations in the Southern Ocean and monitor effects of ongoing climate change in the region on the populations inhabiting these.

A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics

Mountain vegetation is strongly affected by temperature and is expected to shift upwards with climate change. Dynamic vegetation models are often used to assess the impact of climate on vegetation and model output can be compared with paleobotanical data as a reality check. Recent paleoecological studies have revealed regional variation in the upward shift of timberlines in the Northern and Central European Alps in response to rapid warming at the Younger Dryas/Preboreal transition ca. 11700years ago, probably caused by a climatic gradient across the Alps. This contrasts with previous studies that successfully simulated the early Holocene afforestation in the (warmer) Central Alps with a chironomid-inferred temperature reconstruction from the (colder) Northern Alps. We use LandClim, a dynamic landscape vegetation model to simulate mountain forests under different temperature, soil and precipitation scenarios around Iffigsee (2065m a.s.l.) a lake in the Northwestern Swiss Alps, and compare the model output with the paleobotanical records. The model clearly overestimates the upward shift of timberline in a climate scenario that applies chironomid-inferred July-temperature anomalies to all months. However, forest establishment at 9800 cal. BP at Iffigsee is successfully simulated with lower moisture availability and monthly temperatures corrected for stronger seasonality during the early Holocene. The model-data comparison reveals a contraction in the realized niche of Abies alba due to the prominent role of anthropogenic disturbance after ca. 5000 cal. BP, which has important implications for species distribution models (SDMs) that rely on equilibrium with climate and niche stability. Under future climate projections, LandClim indicates a rapid upward shift of mountain vegetation belts by ca. 500m and treeline positions of ca. 2500m a.s.l. by the end of this century. Resulting biodiversity losses in the alpine vegetation belt might be mitigated with low-impact pastoralism to preserve species-rich alpine meadows.

Automatic multi-resolution shape modeling of multi-organ structures.

Point Distribution Models (PDM) are among the most popular shape description techniques and their usefulness has been demonstrated in a wide variety of medical imaging applications. However, to adequately characterize the underlying modeled population it is essential to have a representative number of training samples, which is not always possible. This problem is especially relevant as the complexity of the modeled structure increases, being the modeling of ensembles of multiple 3D organs one of the most challenging cases. In this paper, we introduce a new GEneralized Multi-resolution PDM (GEM-PDM) in the context of multi-organ analysis able to efficiently characterize the different inter-object relations, as well as the particular locality of each object separately. Importantly, unlike previous approaches, the configuration of the algorithm is automated thanks to a new agglomerative landmark clustering method proposed here, which equally allows us to identify smaller anatomically significant regions within organs. The significant advantage of the GEM-PDM method over two previous approaches (PDM and hierarchical PDM) in terms of shape modeling accuracy and robustness to noise, has been successfully verified for two different databases of sets of multiple organs: six subcortical brain structures, and seven abdominal organs. Finally, we propose the integration of the new shape modeling framework into an active shape-model-based segmentation algorithm. The resulting algorithm, named GEMA, provides a better overall performance than the two classical approaches tested, ASM, and hierarchical ASM, when applied to the segmentation of 3D brain MRI.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

pH-dependent distribution of chlorin e6 derivatives across phospholipid bilayers probed by NMR spectroscopy

The pH-dependent membrane adsorption and distribution of three chlorin derivatives, chlorin e6 (CE), rhodin G7 (RG), and monoaspartyl-chlorin e6 (MACE), in the physiological pH range (pH 6-8) were probed by NMR spectroscopy. Unilamellar vesicles consisting of dioleoyl-phosphatidyl-choline (DOPC) were used as membrane models. The chlorin derivatives were characterized with respect to their aggregation behavior, the pK(a) values of individual carboxylate groups, the extent of membrane adsorption, and their flip-flop rates across the bilayer membrane for pH 6-8. External membrane adsorption was found to be lower for RG than for CE and MACE. Both electrostatic interactions and the extent of aggregation seemed to be the main determinants of membrane adsorption. Rate constants for chlorin transfer across the membrane were found to correlate strongly with the pH of the surrounding medium, in particular, for CE and RG. In acidic solution, CE and RG transfer across the membrane was strongly accelerated, and in basic solution, all compounds were retained, mostly in the outer monolayer. In contrast, MACE flip-flop across the membrane remained very low even at pH 6. The protonation of ionizable groups is suggested to be a major determinant of chlorin transfer rates across the bilayer. pK(a) values of CE and RG were found to be between 6 and 8, and two of the carboxylate groups in MACE had pK(a) values below 6. For CE and RG, the kinetic profiles at acidic pH indicated that the initial fast membrane distribution was followed by secondary steps that are discussed in this article.

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.