Skip to main content

Holocene variations of wildfire occurrence as a guide for sustainable management of the northeastern Canadian boreal forest

Abstract

Background

Cumulative impacts of wildfires and forest harvesting can cause shifts from closed-crown forest to open woodland in boreal ecosystems. To lower the probability of occurrence of such catastrophic regime shifts, forest logging must decrease when fire frequency increases, so that the combined disturbance rate does not exceed the Holocene maximum. Knowing how climate warming will affect fire regimes is thus crucial to sustainably manage the forest. This study aimed to provide a guide to determine sustainable forest harvesting levels, by reconstructing the Holocene fire history at the northern limit of commercial forestry in Quebec using charcoal particles preserved in lake sediments.

Methods

Sediment cores were sampled from four lakes located close to the northern limit of commercial forestry in Quebec. The cores were sliced into consecutive 0.5 cm thick subsamples from which 1 cm3 was extracted to count and measure charcoal particles larger than 150 microns. Age-depth models were obtained for each core based on accelerator mass spectroscopy (AMS) radiocarbon dates. Holocene fire histories were reconstructed by combining charcoal counts and age-depth models to obtain charcoal accumulation rates and, after statistical treatment, long-term trends in fire occurrence (expressed as number of fires per 1000 years).

Results

Fire occurrence varied between the four studied sites, but fires generally occurred more often during warm and dry periods of the Holocene, especially during the Holocene Thermal Maximum (7000–3500 cal. BP), when fire occurrence was twice as high as at present.

Conclusions

The current fire regime in the study area is still within the natural range of variability observed over the Holocene. However, climatic conditions comparable to the Holocene Thermal Maximum could be reached within the next few decades, thus substantially reducing the amount of wood available to the forest industry.

Background

The northern part of the boreal biome contains some of the largest remaining tracts of intact forest in the world (Bradshaw et al. 2009). However, forestry operations have expanded in a south-north fashion over the last few decades (Bouchard and Pothier 2011; Kivinen et al. 2012), and the cumulative impacts of natural disturbances and forest harvesting may affect the ecological integrity of boreal ecosystems. For example, closed-crown spruce-moss forests have shifted to open woodlands in eastern Canada under the combined action of wildfire and clearcut logging (Dussart and Payette 2002; Girard et al. 2008). To lower the probability of occurrence of such catastrophic ecosystem shifts, harvesting rate must decrease when fire frequency increases (Bergeron et al. 2006; Raulier et al. 2013), so that the total disturbance rate (harvesting + wildfire) does not exceed the maximum disturbance rate recorded during the Holocene.

Climate models predict a major rise in mean annual temperature across boreal regions, but only a slight – if any – precipitation increase (IPCC 2014). This could lead to dryer conditions, as higher evapotranspiration would offset the precipitation increase (Hély et al. 2010). Hence, forest fires will likely be more frequent and cover larger areas by the end of the 21st century (Goldammer 2013), meaning that the margin of manoeuvre for forestry operations will decrease.

In Quebec (eastern Canada), government authorities have established a northern limit to commercial forestry operations in the boreal forest (Jobidon et al. 2015). North of this limit, forest productivity is deemed too low, and fire frequency too high, to support a profitable forest industry. In a comparative study of Holocene fire regimes along a transect crossing the limit of commercial forestry, Oris et al. (2014b) noted that all regions were characterized by considerably higher fire occurrences during the Holocene Thermal Maximum (HTM; ca. 7000–3000 cal. years BP), followed by a sharp decrease in fire occurrence at the beginning of the Neoglacial period (ca. 3000 cal. years BP) in forests located close to and north of the limit, compared to a more gradual decrease south of the limit. They thus concluded that fire occurrence was more sensitive to climate change close to and north of the limit than further south. These regional patterns, however, hide important between-site variability (see Figure S2 in Oris et al. (2014b)), and additional studies are needed to verify if the same regional fire history can be reconstructed from other sites.

We used high-resolution charcoal analysis to reconstruct the Holocene variability of fire occurrence from 4 small lakes located near the northern limit of commercial forestry in Quebec, ca. 100–200 km east of the lakes studied by Oris et al. (2014b). As wildfire dynamics are mainly controlled by regional climatic patterns in boreal forests (Flannigan and Wotton 2001; Ali et al. 2009; Boulanger et al. 2013), we expected to find a regional fire history similar to that reported by Oris et al. (2014b), despite inter-site variability. Having a clear understanding of the Holocene regional variability in fire occurrence is mandatory to adjust forest harvesting levels so that the total disturbance rate does not exceed the Holocene maximum. We show that fire occurrence at the northern limit of commercial forestry in Quebec was higher during warm and dry periods of the Holocene. According to predictions, climate warming over the next few decades will increase fire occurrence close to the maximum values recorded during the Holocene, thus substantially reducing the margin of manoeuvre of the forest industry.

Methods

Study area

The study area is located within the western black spruce-moss bioclimatic subdomain (Saucier et al. 1998), in northern Quebec, 50–100 km north of Chibougamau (Fig. 1). Mean annual temperature is –0.4 °C and mean annual precipitation is 961.4 mm (Environment Canada 2011). The regional climate history displayed a general progression from a cool early Holocene (<7000 cal. BP) to a warm and dry middle Holocene (Holocene Thermal Maximum: 7000–3500 cal. BP) followed by a cool and wet Neoglacial period (>3500 cal. BP), briefly interrupted (around 1000 cal. BP) by a Medieval Warm Period (Viau and Gajewski 2009). Forests of this region have been dominated by black spruce (Picea mariana (Mill.) BSP) throughout the Holocene (Garralla and Gajewski 1992; Payette 1993).

Fig. 1
figure 1

Location of the 4 sampled lakes (triangles) at the northern limit of commercial forestry (thick black line). The dark gray area in the inset represents the western spruce-moss bioclimatic subdomain

Sampling

Four lakes were selected close to the northern limit of commercial forestry (Fig. 1). The lakes are within 100 km of each other and have similar characteristics: size (<4 ha), absence of inlet or outlet, elevation (376–432 m a.s.l.), water depth (500–571 cm), and local vegetation (dominated by black spruce). Lake sediment cores were sampled in March 2009 at the deepest point of each lake. The water-sediment interface was sampled using a Kajak-Brinkhurst (KB) gravity corer (Glew 1991), and a Livingstone corer was used to sample deeper sediments (Wright et al. 1984). The cores were sliced into consecutive 0.5 cm thick subsamples to maximize temporal resolution. Subsamples were stored in plastic bags in the refrigerator until analysis.

Charcoal quantification

A volume of 1 cm3 was taken from each subsample and soaked in a 3 % sodium hexametaphosphate – (NaPO3)6 – dispersing solution for 48 h before sieving through a 150 microns mesh. The material remaining on the sieve was washed with a 10 % sodium hypochlorite (NaOCl) solution to whiten organic matter particles and facilitate distinction of charcoal particles. Charcoal particles larger than 150 microns were counted and measured (area) under a dissecting microscope (20 ×) equipped with a digital camera coupled to an image analysis software (Winseedle, Regent Instruments Inc., Canada). Charcoal particles larger than 150 microns are rarely carried more than 1 km from the fire and are usually assumed to represent local fires (Wein et al. 1987), although long-distance transport (>30 km) of such particles has been reported (Oris et al. 2014a). We nevertheless elected to work with charcoal particles larger than 150 microns to ease comparisons with Oris et al. (2014b).

Age-depth models

Radiocarbon dates were obtained by accelerator mass spectroscopy (AMS) performed on samples from different levels of the sedimentary profiles (Table 1). Macroremains were rare (except for Aurélie lake), and most dates were obtained from bulk sediments (gyttja). All dates were calibrated using the CALIB 6.0.1 program (http://calib.qub.ac.uk/calib/). Age-depth models were produced using MCAgeDepth software version 0.1 (Higuera et al. 2009; available for free at http://code.google.com/p/mcagedepth/) (Fig. 2).

Table 1 Radiocarbon dates at different depths in the 4 sediment profiles
Fig. 2
figure 2

Age-depth models (cubic spline interpolation) for lakes Aurélie (A), Nans (B), Richard (C), and Twin (D). Radiocarbon dates are shown (black dots) with 95 % confidence intervals (whiskers). One of the dates from the lake Nans profile was discarded (grey dot) as it was too young and was likely contaminated by more recent material during sampling

Fire history reconstruction

Charcoal accumulation rates (CHAR) expressed in mm2∙cm−2∙yr−1 were calculated by multiplying charcoal concentration (mm2∙cm−3) by sedimentation rate (cm∙yr−1) for each subsample. To avoid possible influence of sedimentation variations within and between sites, a constant sedimentation rate (20 years, the average of the 4 lakes) was used to compare lakes and to compute the regional fire history (all lakes combined). Version 1.1. of the CharAnalysis software (Higuera et al. 2009; available for free at http://sites.google.com/site/charanalysis/) was used to remove background noise from the CHAR series with a LOWESS function. The remaining peak component was composed of two Gaussian subpopulations: (1) noise inherent to statistical analysis and part of the background noise that would not have been removed previously, and (2) fire peaks, each one representing the occurrence of one or more local fire events (Gavin et al. 2006). The noise subpopulation was removed from the fire peak subpopulation by applying a threshold corresponding to the 99th percentile of the noise distribution. General trends in fire occurrence were expressed as number of fires per 1000 years, smoothed with a 100 years moving window using K1D software version 1.2 (Gavin 2010). The composite record illustrating regional fire occurrence was obtained by averaging data from the 4 sites, with a 95 % bootstrap confidence interval.

Results

As expected, Holocene trends in fire occurrence (number of fires per 1000 years) varied between lakes (Fig. 3). At Nans Lake, fire occurrence was higher between 5000 and 3800 cal. BP and around 1000 cal. BP. At Richard Lake, fire occurrence was higher between 4600 and 3600 cal. BP and between 1400 and 800 cal. BP. Most fire events at Aurélie Lake occurred before 3400 cal. BP, and the highest occurrence was recorded between 6800 and 4000 cal. BP. At Twin Lake, fire occurrence was higher before 3200 cal. BP.

Fig. 3
figure 3

Fire occurrence (number of fire events per 1000 years) for lakes Aurélie (thick full line), Nans (thin full line), Richard (dashed line), and Twin (dotted line)

The composite curve of regional fire occurrence shows that fire occurrences were generally higher during the middle Holocene (ca. 7000–3500 cal. BP) than during the late Holocene (3500–0 cal. BP), where a brief increase of fire occurrence was nevertheless recorded around 1000 cal. BP (Fig. 4).

Fig. 4
figure 4

Regional average of fire occurrence (number of fire events per 1000 years) combining the 4 sampled lakes. The 95 % confidence interval is shown in grey

Discussion

The long-term fire history at the northern limit of commercial forestry in Quebec fluctuated during the last 7000 years under the influence of regional climate. Fire occurrence was higher between 7000 and 3500 cal. BP, especially between 5500 and 3500 ca. years B.P. This period corresponds to the Holocene Thermal Maximum, previously shown to have been warmer and drier in northern Quebec (Garralla and Gajewski 1992; Viau and Gajewski 2009) and in various locations across North America (Bartlein et al. 1998; Viau et al. 2006). A sharp decrease in fire occurrence occurred in the study area at the beginning of the Neoglacial period (3500–0 cal. BP), especially at lakes Aurélie and Twin, likely in response to a cooler and wetter climate (Garralla and Gajewski 1992; Viau and Gajewski 2009). This general trend is similar to that found further west by Oris et al. (2014b), also at the northern limit of commercial forestry. A brief period of increased fire occurrence was noted in the study area around 1000 cal. BP, coinciding with the Medieval Warm Period (Gajewski 1988; Hunt 2006). The imprint of the Medieval Warm Period was not recorded at all sites, a specificity also noted further south and west in the spruce-moss forest (Ali et al. 2009; Oris et al. 2014b). This could be due to the short duration of this climatic period (a few centuries) making it less likely to affect fire regimes. Alternatively, local factors, such as topography, landscape connectivity, vegetation flammability, or weather could have been more important than regional climate to explain fire occurrence during that period (Gavin et al. 2006; Hu et al. 2006; Long et al. 2007; Ali et al. 2009).

Fire occurrence at the northern limit of commercial forestry in Quebec was variable between sites, but was generally higher during warm and dry periods of the Holocene, especially during the Holocene Thermal Maximum (7000–3500 cal. years BP), when fire occurrence was twice as high as at present. Current fire occurrence is still within the range of natural variability, probably because the amount of precipitation is still high enough to balance the effects of global warming. Fire return intervals (FRI) calculated from the composite chronology of fire occurrences presented here (ca. 167–333 years, based on 3–6 fires per millennium) are slightly higher than – but overlap with – those obtained from spruce-moss forests further south (Cyr et al. 2009: 111–267 years; Ali et al. 2009: 90–230 years) and west (Oris et al. 2014b; 196–312 years). The present-time FRI in our study area (ca. 333 years) is twice as long as that presented by Le Goff et al. (2007) for a particularly dry area located further west (132–153 years). It is however well within the range presented by Mansuy et al. (2010) for a wider study area encompassing both mesic and xeric sites (90–715 years). This is in line with reconstructions from northern Ontario (Canada) that showed significant differences between FRI values from xeric and mesic sites (Senici et al. 2013).

Predicted high temperatures for the middle of the 21st century could be close to those recorded during the Holocene Thermal Maximum (Plummer et al. 2006). If increased temperature is not accompanied by a marked increase in precipitation, the mid-Holocene scenario could be repeated and fire occurrence could increase at the northern limit of commercial forestry in Quebec. The extent to which forest managers will be able to replace wildfire by harvesting activities will be more limited (Bergeron et al. 2006, 2010; Raulier et al. 2013). Indeed, important fire suppression efforts have been deployed in the study area at times when fire activity was at its lowest levels of the past 7000 years.

Conclusions

We showed that fire occurrence was higher during warm and dry periods of the Holocene at the northern limit of commercial forestry in Quebec. Climate predictions for the middle of the 21st century point towards conditions similar to those that prevailed during the Holocene Thermal Maximum, when maximum fire occurrence was recorded. Because sustainable forest management implies that total disturbance (harvesting + wildfire) should not exceed the Holocene maximum, a substantial reduction of the amount of wood available to the forestry industry is to be expected over the next few decades.

References

  • Ali AA, Carcaillet C, Bergeron Y (2009) Long-term fire frequency variability in the eastern Canadian boreal forest: the influences of climate vs. local factors. Global Change Biol 15:1230–1241

    Article  Google Scholar 

  • Bartlein PJ, Anderson KH, Anderson PM, Edwards ME, Mock CJ, Thompson RS, Webb RS, Webb T, Whitlock C (1998) Paleoclimate simulations for North America over the past 21,000 years: features of the simulated climate and comparisons with paleoenvironmental data. Quaternary Sci Rev 17:549–585

    Article  Google Scholar 

  • Bergeron Y, Cyr D, Drever CR, Flannigan M, Gauthier S, Kneeshaw D, Lauzon È, Leduc A, Le Goff H, Lesieur D, Logan K (2006) Past, current, and future fire frequencies in Quebec’s commercial forests: implications for the cumulative effects of harvesting and fire on age-class structure and natural disturbance-based management. Can J Forest Res 36:2737–2744

    Article  Google Scholar 

  • Bergeron Y, Cyr D, Girardin MP, Carcaillet C (2010) Will climate change drive 21st century burn rates in Canadian boreal forest outside of its natural variability: collating global climate model experiments with sedimentary charcoal data. Int J Wildland Fire 19:1127–1139

    Article  Google Scholar 

  • Bouchard M, Pothier D (2011) Long-term influence of fire and harvesting on boreal forest age structure and forest composition in Eastern Québec. Forest Ecol Manag 261:811–820

    Article  Google Scholar 

  • Boulanger Y, Gauthier S, Gray DR, Le Goff H, Lefort P, Morissette J (2013) Fire regime zonation under current and future climate over Eastern Canada. Ecol Appl 23:904–923

    Article  PubMed  Google Scholar 

  • Bradshaw CJA, Warkentin IG, Sodhi NS (2009) Urgent preservation of boreal carbon stocks and biodiversity. Trends Ecol Evol 24:541–548

    Article  PubMed  Google Scholar 

  • Cyr D, Gauthier S, Bergeron Y, Carcaillet C (2009) Forest management is driving the eastern North American boreal forest outside its natural range of variability. Front Ecol Environ 7:519–524

    Article  Google Scholar 

  • Dussart E, Payette S (2002) Ecological impact of clearcutting on black spruce-moss forests in Southern Québec. Ecoscience 9:533–543

    Google Scholar 

  • Environment Canada (2011) Normales climatiques au Canada 1971-2000., National Climate Archive. http://www.climate.weatheroffice.ec.gc.ca/climate_normals/

    Google Scholar 

  • Flannigan MD, Wotton BM (2001) Climate, weather, and area burned. In: Johnson EA, Miyanishi K (eds) Forest fires: behavior and ecological effects. Academic Press, New York

    Google Scholar 

  • Gajewski K (1988) Late Holocene climate changes in eastern North America estimated from pollen data. Quaternary Res 29:255–262

    Article  Google Scholar 

  • Garralla S, Gajewski K (1992) Holocene vegetation history of the boreal forest near Chibougamau, Central Quebec. Can J Bot 70:1364–1368

    Article  Google Scholar 

  • Gavin DG (2010) K1D: multivariate Ripley’s K-function for one-dimensional data. Version 1.2. User guide. Department of Geography, University of Oregon, Eugene, OR, USA

    Google Scholar 

  • Gavin DG, Sheng HF, Lertzman K, Corbett P (2006) Weak climatic control of stand-scale fire history during the late Holocene. Ecology 87:1722–1732

    Article  PubMed  Google Scholar 

  • Girard F, Payette S, Gagnon R (2008) Rapid expansion of lichen woodlands within the closed-crown boreal forest zone over the last 50 years caused by stand disturbances in Eastern Canada. J Biogeogr 35:529–537

    Article  Google Scholar 

  • Glew JR (1991) Miniature gravity corer for recovering short sediment cores. J Paleolimnol 5:285–287

    Article  Google Scholar 

  • Goldammer J (2013) Vegetation fires and global change: Challenges for concerted international action. A white paper directed to the United Nations and International Organizations. A publication of the Global Fire Monitoring Center (GFMC). Kessel Publishing House, Remagen-Oberwinter, Germany

    Google Scholar 

  • Hély C, Girardin MP, Ali AA, Carcaillet C, Brewer S, Bergeron Y (2010) Eastern boreal North American wildfire risk of the past 7000 years: a model-data comparison. Geophys Res Lett 37:L14709

    Google Scholar 

  • Higuera PE, Brubaker LB, Anderson PM, Hu FS, Brown TA (2009) Vegetation mediated the impacts of postglacial climatic change on fire regimes in the south-central Brooks Range, Alaska. Ecol Monogr 79:201–219

    Article  Google Scholar 

  • Hu FS, Brubaker LB, Gavin DG, Higuera PE, Lynch JA, Rupp TS, Tinner W (2006) How climate and vegetation influence the fire regime of the Alaskan Boreal biome: the Holocene perspective. Mitig Adapt Strateg Glob Change 11:829–846

    Article  Google Scholar 

  • Hunt BG (2006) The medieval warm period, the little ice age and simulated climatic variability. Clim Dynam 27:677–694

    Article  Google Scholar 

  • IPCC (2014) Climate change 2014: Synthesis report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. IPCC, Geneva, Switzerland

    Google Scholar 

  • Jobidon R, Bergeron Y, Robitaille A, Raulier F, Gauthier S, Imbeau L, Saucier J-P, Boudreault C (2015) A biophysical approach to delineate a northern limit to commercial forestry: the case of Quebec’s boreal forest. Can J Forest Res 45:515–528

    Article  Google Scholar 

  • Kivinen S, Berg A, Moen J, Östlund L, Olofsson J (2012) Forest fragmentation and landscape transformation in a reindeer husbandry area in Sweden. Environ Manag 49:295–304

    Article  Google Scholar 

  • Le Goff H, Flannigan M, Bergeron Y, Girardin MP (2007) Historical fire regime shifts related to climate teleconnections in the Waswanipi area, central Quebec, Canada. Int J Wildland Fire 16:607–618

    Article  Google Scholar 

  • Long C, Whitlock C, Bartlein PJ (2007) Holocene-scale variations in vegetation and fire history of the Oregon Coast Range, USA. The Holocene 17:917–926

    Article  Google Scholar 

  • Mansuy N, Gauthier S, Robitaille A, Bergeron Y (2010) The effects of surficial deposit–drainage combinations on spatial variations of fire cycles in the boreal forest of eastern Canada. Int J Wildland Fire 19:1083–1098

    Article  Google Scholar 

  • Oris F, Ali AA, Asselin H, Paradis L, Bergeron Y, Finsinger W (2014a) Charcoal dispersion and deposition in boreal lakes from 3 years of monitoring: Differences between local and regional fires. Geophys Res Lett 41:6743–6752

    Article  CAS  Google Scholar 

  • Oris F, Asselin H, Finsinger W, Hély C, Blarquez O, Ferland M-E, Bergeron Y, Ali AA (2014b) Long-term fire history in northern Quebec: implications for the northern limit of commercial forests. J Appl Ecol 51:675–683

    Article  Google Scholar 

  • Payette S (1993) The range limit of boreal tree species in Québec-Labrador: an ecological and palaeoecological interpretation. Rev Palaeobot Palyno 79:7–30

    Article  Google Scholar 

  • Plummer DA, Caya D, Frigon A, Côté H, Giguère M, Paquin D, Biner S, Harvey R, De Elia R (2006) Climate and climate change over North America as simulated by the Canadian RCM. J Climate 19:3112–3132

    Article  Google Scholar 

  • Raulier F, Le Goff H, Gauthier S, Rapanoela R, Bergeron Y (2013) Introducing two indicators for fire risk consideration in the management of boreal forests. Ecol Indic 24:451–461

    Article  Google Scholar 

  • Saucier J-P, Bergeron J-F, Grondin P, Robitaille A (1998) Les régions écologiques du Québec méridional (troisième version). Ministère des Ressources naturelles du Québec. Ministère des Ressources naturelles du Québec, Québec, Canada

    Google Scholar 

  • Senici D, Lucas A, Chen HYH, Bergeron Y, Larouche A, Brossier B, Blarquez O, Ali AA (2013) Multi-millennial fire frequency and tree abundance differ between xeric and mesic boreal forests in central Canada. J Ecol 101:356–367

    Article  Google Scholar 

  • Viau AE, Gajewski K (2009) Reconstructing millennial-scale, regional paleoclimates of boreal Canada during the Holocene. J Climate 22:316–330

    Article  Google Scholar 

  • Viau AE, Gajewski K, Sawada MC, Fines P (2006) Millennial-scale temperature variations in North America during the Holocene. J Geophys Res 111:D09102

    Google Scholar 

  • Wein RW, Burzynski MP, Sreenivasa BA, Tolonen K (1987) Bog profile evidence of fire and vegetation dynamics since 3000 years BP in the Acadian forest. Can J Bot 65:1180–1186

    Article  Google Scholar 

  • Wright HE Jr, Mann DH, Glaser PH (1984) Piston corers for peat and lake sediments. Ecology 65:657–659

    Article  Google Scholar 

Download references

Acknowledgements

We thank Loic Bircker, Laurent Bremond, Aurélie Genries and Raynald Julien for their assistance during fieldwork. Financial support was provided by the Natural Sciences and Engineering Research Council of Canada and by the CNRS (Paleo2, INSU).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hugo Asselin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AE-G carried out the field and laboratory work, performed the analyses and drafted the manuscript. HA, SG, YB and AAA conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

El-Guellab, A., Asselin, H., Gauthier, S. et al. Holocene variations of wildfire occurrence as a guide for sustainable management of the northeastern Canadian boreal forest. For. Ecosyst. 2, 15 (2015). https://doi.org/10.1186/s40663-015-0039-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40663-015-0039-2

Keywords