Evolutionary Ecology of Mountain Birch in Subarctic Stress Gradients: Interplay of Biotic and Abiotic Factors in Plant-Plant Interactions and Evolutionary Processes

In nature, variation for example in herbivory, wind exposure, moisture and pollution impact often creates variation in physiological stress and plant productivity. This variation is seldom clear-cut, but rather results in clines of decreasing growth and productivity towards the high-stress end. These clines of unidirectionally changing stress are generally known as ‘stress gradients’. Through its effect on plant performance, stress has the capacity to fundamentally alter the ecological relationships between individuals, and through variation in survival and reproduction it also causes evolutionary change, i.e. local adaptations to stress and eventually speciation. In certain conditions local adaptations to environmental stress have been documented in a matter of just a few generations. In plant-plant interactions, intensities of both negative interactions (competition) and positive ones (facilitation) are expected to vary along stress gradients. The stress-gradient hypothesis (SGH) suggests that net facilitation will be strongest in conditions of high biotic and abiotic stress, while a more recent ‘humpback’ model predicts strongest net facilitation at intermediate levels of stress. Plant interactions on stress gradients, however, are affected by a multitude of confounding factors, making studies of facilitation-related theories challenging. Among these factors are plant ontogeny, spatial scale, and local adaptation to stress. The last of these has very rarely been included in facilitation studies, despite the potential co-occurrence of local adaptations and changes in net facilitation in stress gradients. Current theory would predict both competitive effects and facilitative responses to be weakest in populations locally adapted to withstand high abiotic stress. This thesis is based on six experiments, conducted both in greenhouses and in the field in Russia, Norway and Finland, with mountain birch (Betula pubescens subsp. czerepanovii) as the model species. The aims were to study potential local adaptations in multiple stress gradients (both natural and anthropogenic), changes in plant-plant interactions under conditions of varying stress (as predicted by SGH), potential mechanisms behind intraspecific facilitation, and factors confounding plant-plant facilitation, such as spatiotemporal, ontogenetic, and genetic differences. I found rapid evolutionary adaptations (occurring within a time-span of 60 to 70 years) towards heavy-metal resistance around two copper-nickel smelters, a phenomenon that has resulted in a trade-off of decreased performance in pristine conditions. Heavy-metal-adapted individuals had lowered nickel uptake, indicating a possible mechanism behind the detected resistance. Seedlings adapted to heavy-metal toxicity were not co-resistant to others forms of abiotic stress, but showed co-resistance to biotic stress by being consumed to a lesser extent by insect herbivores. Conversely, populations from conditions of high natural stress (wind, drought etc.) showed no local adaptations, despite much longer evolutionary time scales. Due to decreasing emissions, I was unable to test SGH in the pollution gradients. In natural stress gradients, however, plant performance was in accordance with SGH, with the strongest host-seedling facilitation found at the high-stress sites in two different stress gradients. Factors confounding this pattern included (1) plant size / ontogenetic status, with seedling-seedling interactions being competition dominated and host-seedling interactions potentially switching towards competition with seedling growth, and (2) spatial distance, with competition dominating at very short planting distances, and facilitation being strongest at a distance of circa ¼ benefactor height. I found no evidence for changes in facilitation with respect to the evolutionary histories of plant populations. Despite the support for SGH, it may be that the ‘humpback’ model is more relevant when the main stressor is resource-related, while what I studied were the effects of ‘non-resource’ stressors (i.e. heavy-metal pollution and wind). The results have potential practical applications: the utilisation of locally adapted seedlings and plant facilitation may increase the success of future restoration efforts in industrial barrens as well as in other wind-exposed sites. The findings also have implications with regard to the effects of global change in subarctic environments: the documented potential by mountain birch for rapid evolutionary change, together with the general lack of evolutionary ‘dead ends’, due to not (over)specialising to current natural conditions, increase the chances of this crucial forest-forming tree persisting even under the anticipated climate change.

Seasonal trends in the stable isotopic composition of snow and meltwater runoff in a subarctic catchment at Okstindan, Norway

Seasonal variations in the stable isotopic composition of snow and meltwater were investigated in a sub-arctic, mountainous, but non-glacial, catchment at Okstindan in northern Norway based on analyses of delta(18)O and deltaD. Samples were collected during four field periods (August 1998; April 1999; June 1999 and August 1999) at three sites lying on an altitudinal transect (740-970 m a.s.l.). Snowpack data display an increase in the mean values of delta(18)O (increasing from a mean value of - 13.51 to - 11.49% between April and August), as well as a decrease in variability through the melt period. Comparison with a regional meteoric water line indicates that the slope of the delta(18)O - deltaD line for the snowpacks decreases over the same period, dropping from 7.49 to approximately 6.2. This change points to the role of evaporation in snowpack ablation and is confirmed by the vertical profile of deuterium excess. Snowpack seepage data, although limited, also suggest reduced values of deltaD, as might be associated with local evaporation during meltwater generation. In general, meltwaters were depleted in delta(18)O relative to the source snowpack at the peak of the melt (June), but later in the year (August) the difference between the two was not statistically significant. The diurnal pattern of isotopic composition indicates that the most depleted meltwaters coincide with the peak in temperature and, hence, meltwater production.

Tracking the origins of Yakutian horses and the genetic basis for their fast adaptation to subarctic environments

Yakutia, Sakha Republic, in the Siberian Far East, represents one of the coldest places on Earth, with winter record temperatures dropping below -70 °C. Nevertheless, Yakutian horses survive all year round in the open air due to striking phenotypic adaptations, including compact body conformations, extremely hairy winter coats, and acute seasonal differences in metabolic activities. The evolutionary origins of Yakutian horses and the genetic basis of their adaptations remain, however, contentious. Here, we present the complete genomes of nine present-day Yakutian horses and two ancient specimens dating from the early 19th century and ∼5,200 y ago. By comparing these genomes with the genomes of two Late Pleistocene, 27 domesticated, and three wild Przewalski's horses, we find that contemporary Yakutian horses do not descend from the native horses that populated the region until the mid-Holocene, but were most likely introduced following the migration of the Yakut people a few centuries ago. Thus, they represent one of the fastest cases of adaptation to the extreme temperatures of the Arctic. We find cis-regulatory mutations to have contributed more than nonsynonymous changes to their adaptation, likely due to the comparatively limited standing variation within gene bodies at the time the population was founded. Genes involved in hair development, body size, and metabolic and hormone signaling pathways represent an essential part of the Yakutian horse adaptive genetic toolkit. Finally, we find evidence for convergent evolution with native human populations and woolly mammoths, suggesting that only a few evolutionary strategies are compatible with survival in extremely cold environments.

(Table 1) Dissolved organic carbon, water temperature and conductivity in three subarctic ponds near Kilpisjärvi, North Finland

Daphnia was collected from five subarctic ponds which differed greatly in their DOC contents and, consequently, their underwater light (UV) climates. Irrespective of which Daphnia species was present, and contrary to expectations, the ponds with the lowest DOC concentrations (highest UV radiation levels) contained Daphnia with the highest eicosapentaenoic acid (EPA) concentrations. In addition, EPA concentrations in these Daphnia generally decreased in concert with seasonally increasing DOC concentrations. Daphnia from three of the ponds was also tested for its tolerance to solar ultraviolet radiation (UVR) with respect to survival. Daphnia pulex from the clear water pond showed, by far, the best UV-tolerance, followed by D. longispina from the moderately humic and D. longispina from the very humic pond. In addition, we measured sublethal parameters related to UV-damage such as the degree to which the gut of Daphnia appeared green (as a measure of their ability to digest algae), and whether their guts appeared damaged. We developed a simple, noninvasive scoring system to quantify the proportion of the gut in which digestive processes were presumably active. This method allowed repeated measurement of the same animals over the course of the experiment. We demonstrated, for the first time, that sublethal damage of the gut precedes mortality caused by exposure to UVR. In a parallel set of experiments we fed UV-exposed and non-exposed algae to UV-exposed and non-exposed daphnids. UVR pretreatment of algae enhanced the negative effects of exposure to natural solar UV-irradiation in Daphnia. These UV-related effects were generally not specific to the species of Daphnia.

Dust flux records from the Subarctic North Pacific

We present a new record of eolian dust flux to the western Subarctic North Pacific (SNP) covering the past 27000 years based on a core from the Detroit Seamount. Comparing the SNP dust record to the NGRIP ice core record shows significant differences in the amplitude of dust changes to the two regions during the last deglaciation, while the timing of abrupt changes is synchronous. If dust deposition in the SNP faithfully records its mobilization in East Asian source regions, then the difference in the relative amplitude must reflect climate-related changes in atmospheric dust transport to Greenland. Based on the synchronicity in the timing of dust changes in the SNP and Greenland, we tie abrupt deglacial transitions in the 230Th-normalized 4He flux record to corresponding transitions in the well-dated NGRIP dust flux record to provide a new chronostratigraphic technique for marine sediments from the SNP. Results from this technique are complemented by radiocarbon dating, which allows us to independently constrain radiocarbon paleoreservoir ages. We find paleoreservoir ages of 745 ± 140 yr at 11653 yr BP, 680 ± 228 yr at 14630 yr BP and 790 ± 498 yr at 23290 yr BP. Our reconstructed paleoreservoir ages are consistent with modern surface water reservoir ages in the western SNP. Good temporal synchronicity between eolian dust records from the Subantarctic Atlantic and equatorial Pacific and the ice core record from Antarctica supports the reliability of the proposed dust tuning method to be used more widely in other global ocean regions.

Radiocarbon dates, accumulation rates, d13C and July temperature of profiles obtained from Selwyn and Misaw Lake, subarctic Canada

The rapid warming of arctic regions during recent decades has been recorded by instrumental monitoring, but the natural climate variability in the past is still sparsely reconstructed across many areas. We have reconstructed past climate changes in subarctic west-central Canada. Stable carbon and oxygen isotope ratios (d13C, d18O) were derived from a single Sphagnum fuscum plant component; alpha-cellulose isolated from stems. Periods of warmer and cooler conditions identified in this region, described in terms of a "Mediaeval Climatic Anomaly" and "Little Ice Age" were registered in the temperature reconstruction based on the d13C record. Some conclusions could be drawn about wet/dry shifts during the same time interval from the d18O record, humification indices and the macrofossil analysis. The results were compared with other proxy data from the vicinity of the study area. The amplitude of the temperature change was similar to that in chironomid based reconstructions, showing c. 6.5 ±2.3 °C variability in July temperatures during the past 6.2 ka.

Radiocarbon ages and stable isotope measurements on sediment cores of the subarctic Pacific and its marginal seas

Under modern conditions only North Pacific Intermediate Water is formed in the northwest Pacific Ocean. This situation might have changed in the past. Recent studies with general circulation models indicate a switch to deep-water formation in the northwest Pacific during Heinrich Stadial 1 (17.5-15.0 ka) of the last glacial termination. Reconstructions of past ventilation changes based on paleoceanographic proxy records are still insufficient to test whether a deglacial mode of deep-water formation in the North Pacific Ocean existed. Here we present deglacial ventilation records based on radiocarbon-derived ventilation ages in combination with epibenthic stable carbon isotopes from the northwest Pacific including the Okhotsk Sea and Bering Sea, the two potential source regions for past North Pacific ventilation changes. Evidence for most rigorous ventilation of the intermediate-depth North Pacific occurred during Heinrich Stadial 1 and the Younger Dryas, simultaneous to significant reductions in Atlantic Meridional Overturning Circulation. Concurrent changes in d13C and ventilation ages point to the Okhotsk Sea as driver of millennial-scale changes in North Pacific Intermediate Water ventilation during the last deglaciation. Our records additionally indicate that changes in the d13C intermediate-water (700-1750 m water depth) signature and radiocarbon-derived ventilation ages are in antiphase to those of the deep North Pacific Ocean (>2100 m water depth) during the last glacial termination. Thus, intermediate- and deep-water masses of the northwest Pacific have a differing ventilation history during the last deglaciation.

Establishment, expansion and regeneration of aspen (Populus tremula) in the Torneträsk area, subarctic Sweden

In subarctic Sweden, recent decadal colonization and expansion of aspen (Populus tremula L.) were recorded. Over the past 100 years, aspen became c. 16 times more abundant, mainly as a result of increased sexual regeneration. Moreover, aspen now reach tree-size (>2 m) at the alpine treeline, an ecotone that has been dominated by mountain birch (Betula pubescens ssp. czerepanovii) for at least the past 4000 years. We found that sexual regeneration in aspen probably occurred seven times or more within the last century. Whereas sexual regeneration occurred during moist years following a year with an exceptionally high June-July temperature, asexual regeneration was favored by warm and dry summers. Disturbance to the birch forest by cyclic moth population outbreaks was critical in aspen establishment in the subalpine area. At the treeline, aspen colonization was less determined by these moth outbreaks, and was mainly restricted by summer temperature. If summer warming persists, aspen spread may continue in subarctic Sweden, particularly at the treeline. However, changing disturbance regimes, future herbivore population dynamics and the responses of aspen's competitors birch and pine to a changing climate may result in different outcomes.

Stable isotope record of foraminifera from the subarctic Pacific

Benthic (Uvigerina spp., Cibicidoides spp., Gyroidinoides spp.) and planktonic (N. pachyderma sinistral, G. bulloides) stable isotope records from three core sites in the central Gulf of Alaska are used to infer mixed-layer and deepwater properties of the late glacial Subarctic Pacific. Glacial-interglacial amplitudes of the planktonic delta18O records are 1.1-1.3 per mil, less than half the amplitude observed at core sites at similar latitudes in the North Atlantic; these data imply that a strong, negative deltaw anomaly existed in the glacial Subarctic mixed layer during the summer, which points to a much stronger low-salinity anomaly than exists today. If true, the upper water column in the North Pacific would have been statically more stable than today, thus suppressing convection even more efficiently. This scenario is further supported by vertical (i.e., planktic versus benthic) delta18O and delta13C gradients of >1 per mil, which suggest that a thermohaline link between Pacific deep waters and the Subarctic Pacific mixed layer did not exist during the late glacial. Epibenthic delta13C in the Subarctic Pacific is more negative than at tropical-subtropical Pacific sites but similar to that recorded at Southern Ocean sites, suggesting ventilation of the deep central Pacific from mid-latitude sources, e.g., from the Sea of Japan and Sea of Okhotsk. Still, convection to intermediate depths could have occurred in the Subarctic during the winter months when heat loss to the atmosphere, sea ice formation, and wind-driven upwelling of saline deep waters would have been most intense. This would be beyond the grasp of our planktonic records which only document mixed-layer temperature-salinity fields extant during the warmer seasons. Also we do not have benthic isotope records from true intermediate water depths of the Subarctic Pacific.

Seawater carbonate chemistry and phytoplankton and eubacterial community compositions in the northwest subarctic Pacific

On-deck CO2-Fe-manipulated incubation experiments were conducted using surface seawater collected from the Western Subarctic Gyre of the NW Pacific in the summer of 2008 to elucidate the impacts of ocean acidification and Fe enrichment on the abundance and community composition of phytoplankton and eubacteria in the study area. During the incubation, excluding the initial period, the mean partial pressures of CO2 in non-Fe-added bottles were 230, 419, 843, and 1124 µatm, whereas those in Fe-added treatments were 152, 394, 791, and 1008 µatm. Changes in the abundance and community composition of phytoplankton were estimated using HPLC pigment signatures with the program CHEMTAX and flow cytometry. A DGGE fingerprint technique targeting 16S rRNA gene fragments was also used to estimate changes in eubacterial phylotypes during incubation. The Fe addition induced diatom blooms, and subsequently stimulated the growth of heterotrophic bacteria such as Roseobacter, Phaeobacter, and Alteromonas in the post-bloom phase. In both the Fe-limited and Fe-replete treatments, concentrations of 19'-hexanoyloxyfucoxanthin, a haptophyte marker, and the cell abundance of coccolithophores decreased at higher CO2 levels (750 and 1000 ppm), whereas diatoms exhibited little response to the changes in CO2 availability. The abundances of Synechococcus and small eukaryotic phytoplankton (<10 µm) increased at the higher CO2 levels. DGGE band positions revealed that Methylobacterium of Alphaproteobacteria occurred solely at lower CO2 levels (180 and 380 ppm) during the post-bloom phase. These results suggest that increases in CO2 level could affect not only the community composition of phytoplankton but also that of eubacteria. As these microorganisms play critical roles in the biological carbon pump and microbial loop, our results indicate that the progression of ocean acidification can alter the biogeochemical processes in the study area.

Radiocarbon dating of peat profiles from the subarctic ecoclimatic region, west-central Canada

Long-term vegetation succession and permafrost dynamics in subarctic peat plateaus of west-central Canada have been studied through detailed plant macrofossil analysis and extensive AMS radiocarbon dating of two peat profiles. Peatland inception at these sites occurred around 5800-5100 yr BP (6600-5900 cal. BP) as a result of paludification of upland forests. At the northern peat plateau site, located in the continuous permafrost zone, palaeobotanical evidence suggests that permafrost was already present under the forested upland prior to peatland development. Paludification was initiated by permafrost collapse, but re-aggradation of permafrost occurred soon after peatland inception. At the southern site, located in the discontinuous permafrost zone, the aggradation of permafrost occurred soon after peatland inception. In the peat plateaus, permafrost conditions have remained very stable until present. Sphagnum fuscum-dominated stages have alternated with more xerophytic communities characterized by ericaceous shrubs. Local peat fires have occurred, but most of these did not cause degradation of the permafrost. Starting from 2800-1100 yr BP (2900-1000 cal. BP) consistently dry surface conditions have prevailed, possibly related to continued frost heave or nearby polygon crack formation.

Impacts of elevated CO2 on particulate and dissolved organic matter production: microcosm experiments using iron-deficient plankton communities in open subarctic waters

Response of phytoplankton to increasing CO2 in seawater in terms of physiology and ecology is key to predicting changes in marine ecosystems. However, responses of natural plankton communities especially in the open ocean to higher CO2 levels have not been fully examined. We conducted CO2 manipulation experiments in the Bering Sea and the central subarctic Pacific, known as high nutrient and low chlorophyll regions, in summer 2007 to investigate the response of organic matter production in iron-deficient plankton communities to CO2 increases. During the 14-day incubations of surface waters with natural plankton assemblages in microcosms under multiple pCO2 levels, the dynamics of particulate organic carbon (POC) and nitrogen (PN), and dissolved organic carbon (DOC) and phosphorus (DOP) were examined with the plankton community compositions. In the Bering site, net production of POC, PN, and DOP relative to net chlorophyll-a production decreased with increasing pCO2. While net produced POC:PN did not show any CO2-related variations, net produced DOC:DOP increased with increasing pCO2. On the other hand, no apparent trends for these parameters were observed in the Pacific site. The contrasting results observed were probably due to the different plankton community compositions between the two sites, with plankton biomass dominated by large-sized diatoms in the Bering Sea versus ultra-eukaryotes in the Pacific Ocean. We conclude that the quantity and quality of the production of particulate and dissolved organic matter may be altered under future elevated CO2 environments in some iron-deficient ecosystems, while the impacts may be negligible in some systems.