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.

Diurnal cycle of fossil and nonfossil carbon using radiocarbon analyses during CalNex

Radiocarbon (14C) analysis is a unique tool to distinguish fossil/nonfossil sources of carbonaceous aerosols. We present 14C measurements of organic carbon (OC) and total carbon (TC) on highly time resolved filters (3–4 h, typically 12 h or longer have been reported) from 7 days collected during California Research at the Nexus of Air Quality and Climate Change (CalNex) 2010 in Pasadena. Average nonfossil contributions of 58% ± 15% and 51% ± 15% were found for OC and TC, respectively. Results indicate that nonfossil carbon is a major constituent of the background aerosol, evidenced by its nearly constant concentration (2–3 μgC m−3). Cooking is estimated to contribute at least 25% to nonfossil OC, underlining the importance of urban nonfossil OC sources. In contrast, fossil OC concentrations have prominent and consistent diurnal profiles, with significant afternoon enhancements (~3 μgC m−3), following the arrival of the western Los Angeles (LA) basin plume with the sea breeze. A corresponding increase in semivolatile oxygenated OC and organic vehicular emission markers and their photochemical reaction products occurs. This suggests that the increasing OC is mostly from fresh anthropogenic secondary OC (SOC) from mainly fossil precursors formed in the western LA basin plume. We note that in several European cities where the diesel passenger car fraction is higher, SOC is 20% less fossil, despite 2–3 times higher elemental carbon concentrations, suggesting that SOC formation from gasoline emissions most likely dominates over diesel in the LA basin. This would have significant implications for our understanding of the on-road vehicle contribution to ambient aerosols and merits further study.

Fossil vs. non-fossil sources of fine carbonaceous aerosols in four Chinese cities during the extreme winter haze episode of 2013

During winter 2013, extremely high concentrations (i.e., 4–20 times higher than the World Health Organization guideline) of PM2.5 (particulate matter with an aerodynamic diameter < 2.5 μm) mass concentrations (24 h samples) were found in four major cities in China including Xi'an, Beijing, Shanghai and Guangzhou. Statistical analysis of a combined data set from elemental carbon (EC), organic carbon (OC), 14C and biomass-burning marker measurements using Latin hypercube sampling allowed a quantitative source apportionment of carbonaceous aerosols. Based on 14C measurements of EC fractions (six samples each city), we found that fossil emissions from coal combustion and vehicle exhaust dominated EC with a mean contribution of 75 ± 8% across all sites. The remaining 25 ± 8% was exclusively attributed to biomass combustion, consistent with the measurements of biomass-burning markers such as anhydrosugars (levoglucosan and mannosan) and water-soluble potassium (K+). With a combination of the levoglucosan-to-mannosan and levoglucosan-to-K+ ratios, the major source of biomass burning in winter in China is suggested to be combustion of crop residues. The contribution of fossil sources to OC was highest in Beijing (58 ± 5%) and decreased from Shanghai (49 ± 2%) to Xi'an (38 ± 3%) and Guangzhou (35 ± 7%). Generally, a larger fraction of fossil OC was from secondary origins than primary sources for all sites. Non-fossil sources accounted on average for 55 ± 10 and 48 ± 9% of OC and total carbon (TC), respectively, which suggests that non-fossil emissions were very important contributors of urban carbonaceous aerosols in China. The primary biomass-burning emissions accounted for 40 ± 8, 48 ± 18, 53 ± 4 and 65 ± 26% of non-fossil OC for Xi'an, Beijing, Shanghai and Guangzhou, respectively. Other non-fossil sources excluding primary biomass burning were mainly attributed to formation of secondary organic carbon (SOC) from non-fossil precursors such as biomass-burning emissions. For each site, we also compared samples from moderately to heavily polluted days according to particulate matter mass. Despite a significant increase of the absolute mass concentrations of primary emissions from both fossil and non-fossil sources during the heavily polluted events, their relative contribution to TC was even decreased, whereas the portion of SOC was consistently increased at all sites. This observation indicates that SOC was an important fraction in the increment of carbonaceous aerosols during the haze episode in China.

Linking megathrust earthquakes to brittle deformation in a fossil accretionary complex

Seismological data from recent subduction earthquakes suggest that megathrust earthquakes induce transient stress changes in the upper plate that shift accretionary wedges into an unstable state. These stress changes have, however, never been linked to geological structures preserved in fossil accretionary complexes. The importance of coseismically induced wedge failure has therefore remained largely elusive. Here we show that brittle faulting and vein formation in the palaeo-accretionary complex of the European Alps record stress changes generated by subduction-related earthquakes. Early veins formed at shallow levels by bedding-parallel shear during coseismic compression of the outer wedge. In contrast, subsequent vein formation occurred by normal faulting and extensional fracturing at deeper levels in response to coseismic extension of the inner wedge. Our study demonstrates how mineral veins can be used to reveal the dynamics of outer and inner wedges, which respond in opposite ways to megathrust earthquakes by compressional and extensional faulting, respectively.

Wet deposition of fossil and non-fossil derived particulate carbon: insights from radiocarbon measurement

Radiocarbon (14C) measurements of both organic carbon (OC) and elemental carbon (EC) allow a more detailed source apportionment, leading to a full and unambiguous distinction and quantification of the contributions from non-fossil and fossil sources. A thermal-optical method with a commercial OC/EC analyzer to isolate water-insoluble OC (WIOC) and EC for their subsequent 14C measurement was applied for the first time to filtered precipitation samples collected at a costal site in Portugal and at a continental site in Switzerland. Our results show that WIOC in precipitation is dominated by non-fossil sources such as biogenic and biomass-burning emissions regardless of rain origins and seasons, whereas EC sources are shared by fossil-fuel combustion and biomass burning. In addition, monthly variation of WIOC in Switzerland was characterized by higher abundance in warm than in cold seasons, highlighting the importance of biogenic emissions to particulate carbon in rainwater. Samples with high particulate carbon concentrations in Portugal were found to be associated with increased biogenic input. Despite the importance of non-fossil sources, fossil emissions account for approximately 20% of particulate carbon in wet deposition for our study, which is in line with fossil contribution in bulk rainwater dissolved organic carbon as well as aerosol WIOC and EC estimated by the 14C approach from other studies.

Pinus nigra (European black pine) as the dominant species of the last glacial pinewoods in south-western to central Iberia: a morphological study of modern and fossil pollen

Aim Our aim was to discriminate different species of Pinus via pollen analysis in order to assess the responses of particular pine species to orbital and millennial-scale climate changes, particularly during the last glacial period. Location Modern pollen grains were collected from current pine populations along transects from the Pyrenees to southern Iberia and the Balearic Islands. Fossil pine pollen was recovered from the south-western Iberian margin core MD95-2042. Methods We measured a set of morphological traits of modern pollen from the Iberian pine species Pinus nigra, P. sylvestris, P. halepensis, P. pinea and P. pinaster and of fossil pine pollen from selected samples of the last glacial period and the early to mid-Holocene. Classification and regression tree (CART) analysis was used to establish a model from the modern dataset that discriminates pollen from the different pine species and allows identification of fossil pine pollen at the species level. Results The CART model was effective in separating pollen of P. nigra and P. sylvestris from that of the Mediterranean pine group (P. halepensis, P. pinea and P. pinaster). The pollen of Pinus nigra diverged from that of P. sylvestris by having a more flattened corpus. Predictions using this model suggested that fossil pine pollen is mainly from P. nigra in all the samples analysed. Pinus sylvestris was more abundant in samples from Greenland stadials than Heinrich stadials, whereas Mediterranean pines increased in samples from Greenland interstadials and during the early to mid-Holocene. Main conclusions Morphological parameters can be successfully used to increase the taxonomic resolution of fossil pine pollen at the species level for the highland pines (P. nigra and P. sylvestris) and at the group of species level for the Mediterranean pines. Our study indicates that P. nigra was the dominant component of the last glacial south-western/central Iberian pinewoods, although the species composition of these woodlands varied in response to abrupt climate changes.

Younger Dryas and Allerød summer temperatures at Gerzensee (Switzerland) inferred from fossil pollen and cladoceran assemblages

Linear- and unimodal-based inference models for mean summer temperatures (partial least squares, weighted averaging, and weighted averaging partial least squares models) were applied to a high-resolution pollen and cladoceran stratigraphy from Gerzensee, Switzerland. The time-window of investigation included the Allerød, the Younger Dryas, and the Preboreal. Characteristic major and minor oscillations in the oxygen-isotope stratigraphy, such as the Gerzensee oscillation, the onset and end of the Younger Dryas stadial, and the Preboreal oscillation, were identified by isotope analysis of bulk-sediment carbonates of the same core and were used as independent indicators for hemispheric or global scale climatic change. In general, the pollen-inferred mean summer temperature reconstruction using all three inference models follows the oxygen-isotope curve more closely than the cladoceran curve. The cladoceran-inferred reconstruction suggests generally warmer summers than the pollen-based reconstructions, which may be an effect of terrestrial vegetation not being in equilibrium with climate due to migrational lags during the Late Glacial and early Holocene. Allerød summer temperatures range between 11 and 12°C based on pollen, whereas the cladoceran-inferred temperatures lie between 11 and 13°C. Pollen and cladocera-inferred reconstructions both suggest a drop to 9–10°C at the beginning of the Younger Dryas. Although the Allerød–Younger Dryas transition lasted 150–160 years in the oxygen-isotope stratigraphy, the pollen-inferred cooling took 180–190 years and the cladoceran-inferred cooling lasted 250–260 years. The pollen-inferred summer temperature rise to 11.5–12°C at the transition from the Younger Dryas to the Preboreal preceded the oxygen-isotope signal by several decades, whereas the cladoceran-inferred warming lagged. Major discrepancies between the pollen- and cladoceran-inference models are observed for the Preboreal, where the cladoceran-inference model suggests mean summer temperatures of up to 14–15°C. Both pollen- and cladoceran-inferred reconstructions suggest a cooling that may be related to the Gerzensee oscillation, but there is no evidence for a cooling synchronous with the Preboreal oscillation as recorded in the oxygen-isotope record. For the Gerzensee oscillation the inferred cooling was ca. 1 and 0.5°C based on pollen and cladocera, respectively, which lies well within the inherent prediction errors of the inference models.

The palaeolimnology of Soppensee (Central Switzerlabnd), as evidenced by diatom, pollen, and fossil-pigment analyses

The development of Soppensee (Central Switzerland, 596 m a.s.l.) has been reconstructed using algal remains such as diatoms, chlorophytes and fossil pigments, as well as the pollen and spores of macrophytes. Sediment accumulation in Soppensee began at the end of the last glacial period, approximately 15,000 yrs ago. During the Oldest Dryas biozone (> 12,700 radiocarbon yrs B.P.) the lake had low primary productivity. After reforestation with birch and later pine, around 12,700 B.P., phases of summer anoxia occurred in the lake. These anoxic conditions were most probably caused by additional carbon input from the catchment, as well as longer phases of stratification due to reduced wind exposure caused by the sheltering effect of increased tree cover. From the Younger Dryas biozone (10,800 to 10,000 radiocarbon yrs B.P.) onwards, Soppensee became meromictic for several millennia.

The recent eutrophication of Baldeggersee (Switzerland) as assessed by fossil diatom assemblages

Diatom analyses with an annual resolution were carried out on varves of the hypertrophic Baldeggersee (Central Swiss Plateau) for the timespan ad 1885 to 1993. They reveal seven major changes in the dominant planktonic diatoms. As a result of progressive nutrient enrichment, Baldeggersee changed in the 1910s from a Cyclotella to a Tabellaria fenestrata dominated assemblage, and eventually in the 1950s to a Stephanodiscus parvus dominated diatom assemblage. The timing and direction of diatom-assemblage changes in the varved sediment compare well with sedimentological and limnological observations. Partitioning of the variance in the diatom data revealed that TP is a stronger explanatory variable than temperature for these changes. A diatom-inferred total phosphorus (TP) reconstruction indicates three major steps in eutrophication, occurring at 1909, the mid-1950s and the mid-1970s. Comparison with TP measurements in the water column demonstrates that the diatom-TP inference model used is able to hindcast past TP concentrations reliably. The major steps in eutrophication led to decreases in diatom diversity and also resulted in a progressive increase of calcite grain-size. The lake restoration programme established since 1982 shows no direct impact on the composition of the diatom assemblages. However, the decrease in phosphorus loads since the mid-1970s is reflected in the diatom assemblages and in decreasing diatom-inferred TP concentrations.

Fossil and non-fossil source contributions to atmospheric carbonaceous aerosols during extreme spring grassland fires in Eastern Europe

In early spring the Baltic region is frequently affected by high-pollution events due to biomass burning in that area. Here we present a comprehensive study to investigate the impact of biomass/grass burning (BB) on the evolution and composition of aerosol in Preila, Lithuania, during springtime open fires. Non-refractory submicron particulate matter (NR-PM1) was measured by an Aerodyne aerosol chemical speciation monitor (ACSM) and a source apportionment with the multilinear engine (ME-2) running the positive matrix factorization (PMF) model was applied to the organic aerosol fraction to investigate the impact of biomass/grass burning. Satellite observations over regions of biomass burning activity supported the results and identification of air mass transport to the area of investigation. Sharp increases in biomass burning tracers, such as levoglucosan up to 683 ngm-3 and black carbon (BC) up to 17 μgm-3 were observed during this period. A further separation between fossil and non-fossil primary and secondary contributions was obtained by coupling ACSM PMF results and radiocarbon (14C) measurements of the elemental (EC) and organic (OC) carbon fractions. Non-fossil organic carbon (OCnf/ was the dominant fraction of PM1, with the primary (POCnf/ and secondary (SOCnf/ fractions contributing 26–44% and 13–23% to the total carbon (TC), respectively. 5–8% of the TC had a primary fossil origin (POCf/, whereas the contribution of fossil secondary organic carbon (SOCf/ was 4–13 %. Nonfossil EC (ECnf/ and fossil EC (ECf/ ranged from 13–24 and 7–13 %, respectively. Isotope ratios of stable carbon and nitrogen isotopes were used to distinguish aerosol particles associated with solid and liquid fossil fuel burning.

Fossil and Nonfossil Sources of Organic and Elemental Carbon Aerosols in the Outflow from Northeast China

Source quantification of carbonaceous aerosols in the Chinese outflow regions still remains uncertain despite their high mass concentrations. Here, we unambiguously quantified fossil and nonfossil contributions to elemental carbon (EC) and organic carbon (OC) of total suspended particles (TSP) from a regional receptor site in the outflow of Northeast China using radiocarbon measurement. OC and EC concentrations were lower in summer, representing mainly marine air, than in other seasons, when air masses mostly traveled over continental regions in Mongolia and northeast China. The annual-mean contribution from fossil-fuel combustion to EC was 76 ± 11% (0.1−1.3 μg m−3). The remaining 24 ± 11% (0.03−0.42 μg m−3) was attributed to biomass burning, with slightly higher contribution in the cold period (∼31%) compared to the warm period (∼21%) because of enhanced emissions from regional biomass combustion sources in China. OC was generally dominated by nonfossil sources, with an annual average of 66 ± 11% (0.5−2.8 μg m−3), approximately half of which was apportioned to primary biomass burning sources (34 ± 6%). In winter, OC almost equally originated from primary OC (POC) emissions and secondary OC (SOC) formation from fossil fuel and biomass-burning sources. In contrast, summertime OC was dominated by primary biogenic emissions as well as secondary production from biogenic and biomass-burning sources, but fossil-derived SOC was the smallest contributor. Distinction of POC and SOC was performed using primary POC-to-EC emission ratios separated for fossil and nonfossil emissions.