Severtsov Institute of Ecology and Evolution

, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

▪ Abstract We explore empirical and theoretical evidence for the functional significance of plant-litter diversity and the extraordinary high diversity of decomposer organisms in the process of litter decomposition and the consequences for biogeochemical cycles. Potential mechanisms for the frequently observed litter-diversity effects on mass loss and nitrogen dynamics include fungi-driven nutrient transfer among litter species, inhibition or stimulation of microorganisms by specific litter compounds, and positive feedback of soil fauna due to greater habitat and food diversity. Theory predicts positive effects of microbial diversity that result from functional niche complementarity, but the few existing experiments provide conflicting results. Microbial succession with shifting enzymatic capabilities enhances decomposition, whereas antagonistic interactions among fungi that compete for similar resources slow litter decay. Soil-fauna diversity manipulations indicate that the number of trophic levels, species identity, and the presence of keystone species have a strong impact on decomposition, whereas the importance of diversity within functional groups is not clear at present. In conclusion, litter species and decomposer diversity can significantly influence carbon and nutrient turnover rates; however, no general or predictable pattern has emerged. Proposed mechanisms for diversity effects need confirmation and a link to functional traits for a comprehensive understanding of how biodiversity interacts with decomposition processes and the consequences of ongoing biodiversity loss for ecosystem functioning.

Assessing Biodiversity Declines Understanding human impact on biodiversity depends on sound quantitative projection. Pereira et al. (p. 1496 , published online 26 October) review quantitative scenarios that have been developed for four main areas of concern: species extinctions, species abundances and community structure, habitat loss and degradation, and shifts in the distribution of species and biomes. Declines in biodiversity are projected for the whole of the 21st century in all scenarios, but with a wide range of variation. Hoffmann et al. (p. 1503 , published online 26 October) draw on the results of five decades' worth of data collection, managed by the International Union for Conservation of Nature Species Survival Commission. A comprehensive synthesis of the conservation status of the world's vertebrates, based on an analysis of 25,780 species (approximately half of total vertebrate diversity), is presented: Approximately 20% of all vertebrate species are at risk of extinction in the wild, and 11% of threatened birds and 17% of threatened mammals have moved closer to extinction over time. Despite these trends, overall declines would have been significantly worse in the absence of conservation actions.

Recent case studies showing substantial declines of insect abundances have raised alarm, but how widespread such patterns are remains unclear. We compiled data from 166 long-term surveys of insect assemblages across 1676 sites to investigate trends in insect abundances over time. Overall, we found considerable variation in trends even among adjacent sites but an average decline of terrestrial insect abundance by ~9% per decade and an increase of freshwater insect abundance by ~11% per decade. Both patterns were largely driven by strong trends in North America and some European regions. We found some associations with potential drivers (e.g., land-use drivers), and trends in protected areas tended to be weaker. Our findings provide a more nuanced view of spatiotemporal patterns of insect abundance trends than previously suggested.

We use eddy covariance measurements of net ecosystem productivity (NEP) from 21 FLUXNET sites (153 site-years of data) to investigate relationships between phenology and productivity (in terms of both NEP and gross ecosystem photosynthesis, GEP) in temperate and boreal forests. Results are used to evaluate the plausibility of four different conceptual models. Phenological indicators were derived from the eddy covariance time series, and from remote sensing and models. We examine spatial patterns (across sites) and temporal patterns (across years); an important conclusion is that it is likely that neither of these accurately represents how productivity will respond to future phenological shifts resulting from ongoing climate change. In spring and autumn, increased GEP resulting from an 'extra' day tends to be offset by concurrent, but smaller, increases in ecosystem respiration, and thus the effect on NEP is still positive. Spring productivity anomalies appear to have carry-over effects that translate to productivity anomalies in the following autumn, but it is not clear that these result directly from phenological anomalies. Finally, the productivity of evergreen needleleaf forests is less sensitive to phenology than is productivity of deciduous broadleaf forests. This has implications for how climate change may drive shifts in competition within mixed-species stands.

Land-use transitions can enhance the livelihoods of smallholder farmers but potential economic-ecological trade-offs remain poorly understood. Here, we present an interdisciplinary study of the environmental, social and economic consequences of land-use transitions in a tropical smallholder landscape on Sumatra, Indonesia. We find widespread biodiversity-profit trade-offs resulting from land-use transitions from forest and agroforestry systems to rubber and oil palm monocultures, for 26,894 aboveground and belowground species and whole-ecosystem multidiversity. Despite variation between ecosystem functions, profit gains come at the expense of ecosystem multifunctionality, indicating far-reaching ecosystem deterioration. We identify landscape compositions that can mitigate trade-offs under optimal land-use allocation but also show that intensive monocultures always lead to higher profits. These findings suggest that, to reduce losses in biodiversity and ecosystem functioning, changes in economic incentive structures through well-designed policies are urgently needed.

The widespread extinctions of large mammals at the end of the Pleistocene epoch have often been attributed to the depredations of humans; here we present genetic evidence that questions this assumption. We used ancient DNA and Bayesian techniques to reconstruct a detailed genetic history of bison throughout the late Pleistocene and Holocene epochs. Our analyses depict a large diverse population living throughout Beringia until around 37,000 years before the present, when the population's genetic diversity began to decline dramatically. The timing of this decline correlates with environmental changes associated with the onset of the last glacial cycle, whereas archaeological evidence does not support the presence of large populations of humans in Eastern Beringia until more than 15,000 years later.

Abstract Climate and litter quality are primary drivers of terrestrial decomposition and, based on evidence from multisite experiments at regional and global scales, are universally factored into global decomposition models. In contrast, soil animals are considered key regulators of decomposition at local scales but their role at larger scales is unresolved. Soil animals are consequently excluded from global models of organic mineralization processes. Incomplete assessment of the roles of soil animals stems from the difficulties of manipulating invertebrate animals experimentally across large geographic gradients. This is compounded by deficient or inconsistent taxonomy. We report a global decomposition experiment to assess the importance of soil animals in C mineralization, in which a common grass litter substrate was exposed to natural decomposition in either control or reduced animal treatments across 30 sites distributed from 43°S to 68°N on six continents. Animals in the mesofaunal size range were recovered from the litter by Tullgren extraction and identified to common specifications, mostly at the ordinal level. The design of the trials enabled faunal contribution to be evaluated against abiotic parameters between sites. Soil animals increase decomposition rates in temperate and wet tropical climates, but have neutral effects where temperature or moisture constrain biological activity. Our findings highlight that faunal influences on decomposition are dependent on prevailing climatic conditions. We conclude that (1) inclusion of soil animals will improve the predictive capabilities of region‐ or biome‐scale decomposition models, (2) soil animal influences on decomposition are important at the regional scale when attempting to predict global change scenarios, and (3) the statistical relationship between decomposition rates and climate, at the global scale, is robust against changes in soil faunal abundance and diversity.

Soil organisms, including earthworms, are a key component of terrestrial ecosystems. However, little is known about their diversity, their distribution, and the threats affecting them. We compiled a global dataset of sampled earthworm communities from 6928 sites in 57 countries as a basis for predicting patterns in earthworm diversity, abundance, and biomass. We found that local species richness and abundance typically peaked at higher latitudes, displaying patterns opposite to those observed in aboveground organisms. However, high species dissimilarity across tropical locations may cause diversity across the entirety of the tropics to be higher than elsewhere. Climate variables were found to be more important in shaping earthworm communities than soil properties or habitat cover. These findings suggest that climate change may have serious implications for earthworm communities and for the functions they provide.

Abstract. In 2009, the International Soil Moisture Network (ISMN) was initiated as a community effort, funded by the European Space Agency, to serve as a centralised data hosting facility for globally available in situ soil moisture measurements (Dorigo et al., 2011b, a). The ISMN brings together in situ soil moisture measurements collected and freely shared by a multitude of organisations, harmonises them in terms of units and sampling rates, applies advanced quality control, and stores them in a database. Users can freely retrieve the data from this database through an online web portal (https://ismn.earth/en/, last access: 28 October 2021). Meanwhile, the ISMN has evolved into the primary in situ soil moisture reference database worldwide, as evidenced by more than 3000 active users and over 1000 scientific publications referencing the data sets provided by the network. As of July 2021, the ISMN now contains the data of 71 networks and 2842 stations located all over the globe, with a time period spanning from 1952 to the present. The number of networks and stations covered by the ISMN is still growing, and approximately 70 % of the data sets contained in the database continue to be updated on a regular or irregular basis. The main scope of this paper is to inform readers about the evolution of the ISMN over the past decade, including a description of network and data set updates and quality control procedures. A comprehensive review of the existing literature making use of ISMN data is also provided in order to identify current limitations in functionality and data usage and to shape priorities for the next decade of operations of this unique community-based data repository.

Terrestrial gross primary productivity (GPP) varies greatly over time and space. A better understanding of this variability is necessary for more accurate predictions of the future climate-carbon cycle feedback. Recent studies have suggested that variability in GPP is driven by a broad range of biotic and abiotic factors operating mainly through changes in vegetation phenology and physiological processes. However, it is still unclear how plant phenology and physiology can be integrated to explain the spatiotemporal variability of terrestrial GPP. Based on analyses of eddy-covariance and satellite-derived data, we decomposed annual terrestrial GPP into the length of the CO2 uptake period (CUP) and the seasonal maximal capacity of CO2 uptake (GPPmax). The product of CUP and GPPmax explained >90% of the temporal GPP variability in most areas of North America during 2000-2010 and the spatial GPP variation among globally distributed eddy flux tower sites. It also explained GPP response to the European heatwave in 2003 (r(2) = 0.90) and GPP recovery after a fire disturbance in South Dakota (r(2) = 0.88). Additional analysis of the eddy-covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPPmax than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPPmax and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.

Abstract Solar‐induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite‐observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory‐2 (OCO‐2) provides the first opportunity to examine the SIF–GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO‐2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO‐2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm ( R 2 = 0.72, p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes ( R 2 = 0.57–0.79, p < 0.0001) except evergreen broadleaf forests ( R 2 = 0.16, p < 0.05) at the daily timescale. A higher slope was found for C 4 grasslands and croplands than for C 3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome‐specific SIF–GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO‐2 SIF generally had a better performance for predicting GPP than satellite‐derived vegetation indices and a light use efficiency model. The universal SIF–GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO‐2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

A precise knowledge of the spatial distribution of taxa is essential for decision-making processes in land management and biodiversity conservation, both for present and under future global change scenarios. This is a key base for several scientific disciplines (e.g. macro-ecology, biogeography, evolutionary biology, spatial planning, or environmental impact assessment) that rely on species distribution maps. An atlas summarizing the distribution of European amphibians and reptiles with 50 × 50 km resolution maps based on ca. 85 000 grid records was published by the Societas Europaea Herpetologica (SEH) in 1997. Since then, more detailed species distribution maps covering large parts of Europe became available, while taxonomic progress has led to a plethora of taxonomic changes including new species descriptions. To account for these progresses, we compiled information from different data sources: published in books and websites, ongoing national atlases, personal data kindly provided to the SEH, the 1997 European Atlas, and the Global Biodiversity Information Facility (GBIF). Databases were homogenised, deleting all information except species names and coordinates, projected to the same coordinate system (WGS84) and transformed into a 50 × 50 km grid. The newly compiled database comprises more than 384 000 grid and locality records distributed across 40 countries. We calculated species richness maps as well as maps of Corrected Weighted Endemism and defined species distribution types (i.e. groups of species with similar distribution patterns) by hierarchical cluster analysis using Jaccard’s index as association measure. Our analysis serves as a preliminary step towards an interactive, dynamic and online distributed database system (NA2RE system) of the current spatial distribution of European amphibians and reptiles. The NA2RE system will serve as well to monitor potential temporal changes in their distributions. Grid maps of all species are made available along with this paper as a tool for decision-making and conservation-related studies and actions. We also identify taxonomic and geographic gaps of knowledge that need to be filled, and we highlight the need to add temporal and altitudinal data for all records, to allow tracking potential species distribution changes as well as detailed modelling of the impacts of land use and climate change on European amphibians and reptiles.

Soil organisms drive major ecosystem functions by mineralising carbon and releasing nutrients during decomposition processes, which supports plant growth, aboveground biodiversity and, ultimately, human nutrition. Soil ecologists often operate with functional groups to infer the effects of individual taxa on ecosystem functions and services. Simultaneous assessment of the functional roles of multiple taxa is possible using food-web reconstructions, but our knowledge of the feeding habits of many taxa is insufficient and often based on limited evidence. Over the last two decades, molecular, biochemical and isotopic tools have improved our understanding of the feeding habits of various soil organisms, yet this knowledge is still to be synthesised into a common functional framework. Here, we provide a comprehensive review of the feeding habits of consumers in soil, including protists, micro-, meso- and macrofauna (invertebrates), and soil-associated vertebrates. We have integrated existing functional group classifications with findings gained with novel methods and compiled an overarching classification across taxa focusing on key universal traits such as food resource preferences, body masses, microhabitat specialisation, protection and hunting mechanisms. Our summary highlights various strands of evidence that many functional groups commonly used in soil ecology and food-web models are feeding on multiple types of food resources. In many cases, omnivory is observed down to the species level of taxonomic resolution, challenging realism of traditional soil food-web models based on distinct resource-based energy channels. Novel methods, such as stable isotope, fatty acid and DNA gut content analyses, have revealed previously hidden facets of trophic relationships of soil consumers, such as food assimilation, multichannel feeding across trophic levels, hidden trophic niche differentiation and the importance of alternative food/prey, as well as energy transfers across ecosystem compartments. Wider adoption of such tools and the development of open interoperable platforms that assemble morphological, ecological and trophic data as traits of soil taxa will enable the refinement and expansion of the multifunctional classification of consumers in soil. The compiled multifunctional classification of soil-associated consumers will serve as a reference for ecologists working with biodiversity changes and biodiversity-ecosystem functioning relationships, making soil food-web research more accessible and reproducible.

Summary Based on review and original data, this synthesis investigates carbon pools and fluxes of Siberian and European forests (600 and 300 million ha, respectively). We examine the productivity of ecosystems, expressed as positive rate when the amount of carbon in the ecosystem increases, while (following micrometeorological convention) downward fluxes from the atmosphere to the vegetation (NEE = Net Ecosystem Exchange) are expressed as negative numbers. Productivity parameters are Net Primary Productivity (NPP=whole plant growth), Net Ecosystem Productivity (NEP = CO 2 assimilation minus ecosystem respiration), and Net Biome Productivity (NBP = NEP minus carbon losses through disturbances bypassing respiration, e.g. by fire and logging). Based on chronosequence studies and national forestry statistics we estimate a low average NPP for boreal forests in Siberia: 123 gC m –2 y –1 . This contrasts with a similar calculation for Europe which suggests a much higher average NPP of 460 gC m –2 y –1 for the forests there. Despite a smaller area, European forests have a higher total NPP than Siberia (1.2–1.6 vs. 0.6–0.9 × 10 15 gC region –1 y –1 ). This arises as a consequence of differences in growing season length, climate and nutrition. For a chronosequence of Pinus sylvestris stands studied in central Siberia during summer, NEE was most negative in a 67‐y old stand regenerating after fire (– 192 mmol m –2 d –1 ) which is close to NEE in a cultivated forest of Germany (– 210 mmol m –2 d –1 ). Considerable net ecosystem CO 2 ‐uptake was also measured in Siberia in 200‐ and 215‐y old stands (NEE:174 and – 63 mmol m –2 d –1 ) while NEP of 7‐ and 13‐y old logging areas were close to the ecosystem compensation point. Two Siberian bogs and a bog in European Russia were also significant carbon sinks (– 102 to – 104 mmol m –2 d –1 ). Integrated over a growing season (June to September) we measured a total growing season NEE of – 14 mol m –2 summer –1 (– 168 gC m –2 summer –1 ) in a 200‐y Siberian pine stand and – 5 mol m –2 summer –1 (– 60 gC m –2 summer –1 ) in Siberian and European Russian bogs. By contrast, over the same period, a spruce forest in European Russia was a carbon source to the atmosphere of (NEE: + 7 mol m –2 summer –1 = + 84 gC m –2 summer –1 ). Two years after a windthrow in European Russia, with all trees being uplifted and few successional species, lost 16 mol C m –2 to the atmosphere over a 3‐month in summer, compared to the cumulative NEE over a growing season in a German forest of – 15.5 mol m –2 summer –1 (– 186 gC m –2 summer –1 ; European flux network annual averaged – 205 gC m –2 y –1 ). Differences in CO 2 ‐exchange rates coincided with differences in the Bowen ratio, with logging areas partitioning most incoming radiation into sensible heat whereas bogs partitioned most into evaporation (latent heat). Effects of these different surface energy exchanges on local climate (convective storms and fires) and comparisons with the Canadian BOREAS experiment are discussed. Following a classification of disturbances and their effects on ecosystem carbon balances, fire and logging are discussed as the main processes causing carbon losses that bypass heterotrophic respiration in Siberia. Following two approaches, NBP was estimated to be only about 13–16 mmol m –2 y –1 for Siberia. It may reach 67 mmol m –2 y –1 in North America, and about 140–400 mmol m –2 y –1 in Scandinavia. We conclude that fire speeds up the carbon cycle, but that it results also in long‐term carbon sequestration by charcoal formation. For at least 14 years after logging, regrowth forests remain net sources of CO 2 to the atmosphere. This has important implications regarding the effects of Siberian forest management on atmospheric concentrations. For many years after logging has taken place, regrowth forests remain weaker sinks for atmospheric CO 2 than are nearby old‐growth forests.

The origin and evolution of the domestic dog remains a controversial question for the scientific community, with basic aspects such as the place and date of origin, and the number of times dogs were domesticated, open to dispute. Using whole genome sequences from a total of 58 canids (12 gray wolves, 27 primitive dogs from Asia and Africa, and a collection of 19 diverse breeds from across the world), we find that dogs from southern East Asia have significantly higher genetic diversity compared to other populations, and are the most basal group relating to gray wolves, indicating an ancient origin of domestic dogs in southern East Asia 33 000 years ago. Around 15 000 years ago, a subset of ancestral dogs started migrating to the Middle East, Africa and Europe, arriving in Europe at about 10 000 years ago. One of the out of Asia lineages also migrated back to the east, creating a series of admixed populations with the endemic Asian lineages in northern China before migrating to the New World. For the first time, our study unravels an extraordinary journey that the domestic dog has traveled on earth.

The family of proteins homologous to the green fluorescent protein (GFP) from Aequorea victoria exhibits striking diversity of features, including several different types of autocatalytically synthesized chromophores. Here we report 11 new members of the family, among which there are 3 red-emitters possessing unusual features, and discuss the similarity relationships within the family in structural, spectroscopic, and evolutionary terms. Phylogenetic analysis has shown that GFP-like proteins from representatives of subclass Zoantharia fall into at least four distinct clades, each clade containing proteins of more than one emission color. This topology suggests multiple recent events of color conversion. Combining this result with previous mutagenesis and structural data, we propose that (i) different chromophore structures are alternative products synthesized within a similar autocatalytic environment, and (ii) the phylogenetic pattern and color diversity in reef Anthozoa is a result of a balance between selection for GFP-like proteins of particular colors and mutation pressure driving the color conversions.

Abstract Pollen and plant macrofossil data from northern Eurasia were used to reconstruct the vegetation of the last glacial maximum (LGM: 18,000 ± 2000 14 C yr bp ) using an objective quantitative method for interpreting pollen data in terms of the biomes they represent ( Prentice et al. , 1996 ). The results confirm previous qualitative vegetation reconstructions at the LGM but provide a more comprehensive analysis of the data. Tundra dominated a large area of northern Eurasia (north of 57°N) to the west, south and east of the Scandinavian ice sheet at the LGM. Steppe‐like vegetation was reconstructed in the latitudinal band from western Ukraine, where temperate deciduous forests grow today, to western Siberia, where taiga and cold deciduous forests grow today. The reconstruction shows that steppe graded into tundra in Siberia, which is not the case today. Taiga grew on the northern coast of the Sea of Azov, about 1500 km south of its present limit in European Russia. In contrast, taiga was reconstructed only slightly south of its southern limit today in south‐western Siberia. Broadleaved trees were confined to small refuges, e.g. on the eastern coast of the Black Sea, where cool mixed forest was reconstructed from the LGM data. Cool conifer forests in western Georgia were reconstructed as growing more than 1000 m lower than they grow today. The few scattered sites with LGM data from the Tien‐Shan Mountains and from northern Mongolia yielded biome reconstructions of steppe and taiga, which are the biomes growing there today.

Hantaviruses are rodent-borne, emerging viruses that cause life-threatening human diseases in Eurasia and the Americas. We detected hantavirus genome sequences in an African wood mouse (Hylomyscus simus) captured in Sangassou, Guinea. Sequence and phylogenetic analyses of the genetic material demonstrate a novel hantavirus species, which we propose to name "Sangassou virus."

Mid-to-high latitude forests play an important role in the terrestrial carbon cycle, but the representation of photosynthesis in boreal forests by current modelling and observational methods is still challenging. In particular, the applicability of existing satellite-based proxies of greenness to indicate photosynthetic activity is hindered by small annual changes in green biomass of the often evergreen tree population and by the confounding effects of background materials such as snow. As an alternative, satellite measurements of sun-induced chlorophyll fluorescence (SIF) can be used as a direct proxy of photosynthetic activity. In this study, the start and end of the photosynthetically active season of the main boreal forests are analysed using spaceborne SIF measurements retrieved from the GOME-2 instrument and compared to that of green biomass, proxied by vegetation indices including the Enhanced Vegetation Index (EVI) derived from MODIS data. We find that photosynthesis and greenness show a similar seasonality in deciduous forests. In high-latitude evergreen needleleaf forests, however, the length of the photosynthetically active period indicated by SIF is up to 6 weeks longer than the green biomass changing period proxied by EVI, with SIF showing a start-of-season of approximately 1 month earlier than EVI. On average, the photosynthetic spring recovery as signalled by SIF occurs as soon as air temperatures exceed the freezing point (2-3 °C) and when the snow on the ground has not yet completely melted. These findings are supported by model data of gross primary production and a number of other studies which evaluated in situ observations of CO2 fluxes, meteorology and the physiological state of the needles. Our results demonstrate the sensitivity of space-based SIF measurements to light-use efficiency of boreal forests and their potential for an unbiased detection of photosynthetic activity even under the challenging conditions interposed by evergreen boreal ecosystems.