Skip to main content

Treeline dynamics in response to climate change in the Min Mountains, southwestern China

Abstract

Background

Abies faxoniana is the dominant plant species of the forest ecosystem on the eastern edge of Qinghai-Tibet Plateau, where the treeline is strongly defined by climate. The tree-ring chronologies and age structure of Abies faxoniana were developed in the treeline ecotones on the northwestern and southeastern aspects of the Min Mountains in the Wanglang Nature Reserve to examine the treeline dynamics of recent decades in response to climate change.

Results

On the northwestern aspect, correlation analysis showed that the radial growth was significantly and positively correlated with precipitation in current January and monthly mean temperature in current April, but significantly and negatively correlated with monthly mean temperature in previous August. On the southeastern aspect, the radial growth was significantly negatively correlated with monthly mean temperature in previous July and August.

Conclusions

The different responses of radial growth to climatic variability on both the aspects might be mainly due to the micro-environmental conditions. The recruitment benefited from the warm temperature in current April, July and September on the northwestern aspect. The responses of radial growth and recruitment to climatic variability were similar on the northwestern slope. Recruitment was greatly restricted by competition with dense bamboos on the southeastern aspect.

Background

Upper treeline ecotones within alpine and arctic ecosystems are mainly controlled by climate (e.g. Lavoie and Payette, 1994; Kullman, 1999; Grace et al., 2002), and are considered as sensitive proxy biomonitors for revealing the impact of climate variability on the distribution of high-elevation mountain forests (e.g. Camarero and Gutiérrez, 2004; Danby and Hik, 2007; Elliott, 2011). Many studies have shown that climatic variables limit the tree’s radial growth and recruitment at altitudinal treelines (Cullen et al., 2001; Takahashi et al., 2003; Wilmking et al., 2004; Elliott and Kipfmueller, 2011). The close relationships between climatic factors and tree radial growth are widely used to reconstruct past patterns of climate changes (e.g. Fritts, 1976; Cullen et al., 2001; Liu et al., 2009). The frequency of seedling recruitment at treeline ecotones can also reflect climatic change well (e.g. Kullman, 1993; Daniels and Veblen, 2004).

Various studies have revealed that treeline position is sensitive to temperature changes and climate warming has caused an increase in treeline elevation over time (e.g. Brubaker, 1986; Lloyd and Fastie, 2003; Danby and Hik, 2007; Leonelli et al., 2011) and is likely to cause further increases in treeline in the future (Munier et al., 2010). Increasing temperatures can provide a possible mechanism for abrupt increases in recruitment (Hessl and Baker, 1997; Elliott and Kipfmueller, 2011). Tree growth and survival at some upper timberlines are strongly limited by the low-temperature among the main factors controlling the treeline altitude (Tranquillini, 1979; Stevens and Fox, 1991; Körner, 2003; Holtmeier, 2009). In some altitudinal regions, radial growth of trees is driven by summer temperatures (LaMarche and Fritts, 1971; Eckstein and Aniol, 1981; Schweingruber et al., 1991; Bradley and Jones, 1993). There have been strong links between increased recruitment and warmer temperatures during the growing season and the cool seasons at treeline ecotones (Elliott and Baker, 2004; Danby and Hik, 2007; Holtmeier and Broll, 2007; Harsch et al., 2009; Kullman and Öberg, 2009).

As well as temperature, precipitation can also influence the treeline dynamics greatly. Global warming could exacerbate possible water limitations (Andersen et al., 2009) and cause heat-induced moisture stress without a concomitant increase in precipitation at altitudinal treeline (Weisberg and Baker, 1995; Daniels and Veblen, 2004). It has shown that moisture stress could limit seedling recruitment (Hessl and Baker, 1997; Lloyd and Graumlich, 1997) and tree growth (Jacoby and D’Arrigo, 1995; Barber et al., 2000; Lloyd and Fastie, 2002) at some upper treelines. Trees may lose the ability to grow continuously in the warmer temperature condition if insufficient soil water leads to drought stress, therefore precipitation may show a high positive correlation with tree ring-width (Bunn et al., 2005). Inversely, less precipitation in late spring and early summer may favor tree establishment by prolonging the growing season (Elliott and Kipfmueller, 2011).

It has been suggested that dendroclimatic data alone cannot determine the causes of changes in the structure of ecosystems and populations (Moiseev, 2002). In fact, tree radial growth and seedling recruitment are interrelated and must be considered together in order to gain an accurate understanding of treeline dynamics. Studies have uncovered treeline dynamics in relation to climatic change at the population and community level by studying dendroecological techniques coupled with stand age structures, climatic factors and ecological attributes (Ruffner and Abrams, 1998; Daniels and Veblen, 2004; Bunn et al., 2005; Wang et al., 2006; Jump et al., 2007; Elliott and Kipfmueller, 2011). At some altitudinal treelines, the climate conditions that facilitate the radial growth are similar to those that are conductive to recruitment (e.g. Szeicz and Macdonald, 1995; Camarero and Gutiérrez 1999; Gervais and MacDonald, 2000; Jump et al., 2007; Dang et al., 2009). Yet, the two processes of recruitment and growth may respond differently to climatic factors in some other treelines (e.g. Earle, 1993; Daniels and Veblen, 2004; Wang et al., 2006). The sensitivity of treelines to climate change varies with local and regional topographical conditions and thus differs as to its extent, intensity and the process of change (Holtmeier and Broll, 2005).

In mountainous areas, slope aspect has been considered as an important role for exploring the variability of upper treelines to climate change (e.g. Danby and Hik, 2007; Dang et al., 2009; Elliott and Kipfmueller, 2010; Elliott and Kipfmueller, 2011). For instance, soil moisture conditions on different slopes may exert notable differences in the spatiotemporal patterns of tree regeneration at upper treelines (e.g. Daniels and Veblen, 2004; Elliott and Kipfmueller, 2011). Treeline elevation and stand density may increase differently between slope aspects due to the differential presence of permafrost (e.g. Danby and Hik, 2007). The environmental factors mediated by slope aspect should be considered when assessing possible treeline response to climate change (Elliott and Kipfmueller, 2010).

To date, the manner by which climate variability affects radial growth and seedling recruitment of many upper treeline species in different geographic locations is not completely understood (Wang et al., 2006; Dang et al., 2009). Qinghai-Tibet Plateau (QTP) is considered as one of the most sensitive areas to global climate change in China (Hou et al., 2008). Some dendrochronological studies have been conducted in the northeastern QTP (Liu et al., 2006; Li et al., 2008; Fang et al., 2009), but few of these have explored the main climatic factors of temperature and precipitation to determine how each (alone and in combination) influence the subalpine treeline on the eastern edge of QTP. Based on recent data, the temperature in QTP has increased significantly over the past 50 years (Ding et al., 2009). Thus, the study presented herein was designed to determine how radial growth and recruitment of Abies faxoniana responded to the variability of temperature and precipitation on both northwestern and southeastern aspects in the Min Mountains on the eastern edge of QTP. We predicted that: climate warming enhanced both the radial growth and the seedling recruitment of A. faxoniana in the latest decades.

Methods

Study area

This study was conducted in the Wanglang National Nature Reserve (32°49′-33°02′N, 103°55′-104°10′E) in the Min Mountains, on the eastern edge of QTP in western Sichuan Province, southwestern China (Figure 1). Elevations in the reserve range from 2300 to 4980 m. The terrain is steep and deeply dissected, with sharp environmental gradients.

Figure 1
figure 1

Locations of sampling sites (▲) at the treeline ecotones on the northwestern (NW) and southeastern (SE) of the Min Mountains in the Wanglang Nature Reserve, southwestern China.

The Wanglang National Nature Reserve belongs to Danba-Songpan semi-humid climate, which is characterized by dry, cold winters and wet, cool summers. The mean annual temperature is 2.3°C and the mean annual total precipitation is about 1100 mm (Taylor et al., 2006). The average July temperature is 12.7°C and the average January temperature is −6.1°C, with a recorded extreme high temperature of 26.2°C and extreme low temperature of −17.8°C. The annual rainfall time may last more than 195 days and is concentrated in May, June and July. Between 1951 and 2009, the annual mean temperature increased at an average rate of 0.0137°C /yr (R2 = 0.1631, P = 0.0015) (Figure 2a). However, insignificant descending trend for the annual precipitation was found during this period (Figure 2b).

Figure 2
figure 2

Variations in annual mean temperature (a), annual precipitation (b) during the period 1951–2009 at the study area, calculated with data from Songpan meteorological station and using MTCLIM program.

The subalpine fir, A. faxoniana, is the dominant coniferous tree species in the study area, having the most widely distributed range and largest volume. Most of the reserve from an elevation of ~2700 m to the upper treeline is covered by this species. The conifer forests below 2700 m were clear-felled in the 1950s, and some separate trees on or near the valley bottoms above this altitude were also cut until the reserve was established in 1964 (Taylor et al., 2006). However, the uppermost trees, often situated in fairly inaccessible locales, have remained generally undisturbed by humans or livestock. The treeline forests have been generating without disturbance, with no evidences of fires or previous logging.

Site selection and field sampling

From June to August of 2010, two plots (each was 100 m × 20 m, 0.2 ha) were established at the subalpine treeline ecotones on the northwestern (32°59′27.3″ N, 104°01′44.0″ E, 3297 m a.s.l.) and southeastern (32°54′09.5″ N, 104°02′51.6″ E, 3,225 m a.s.l.) aspects of the Min Mountains in the Wanglang Nature Reserve (Figure 1). The main reason for choosing contrasting slope aspects is that the main mountain ridges all run northeast to southwest in this area. Both plots on northwestern and southeastern aspects can represent the treeline ecotones exactly here. The longer side of each plot was parallel to the isoline. The plots were selected based on the criterion that they should represent the fir forest structure at subalpine treeline.

In the plot, A. faxoniana individuals were divided into three height classes (Wang et al., 2006): trees (> 2 m), saplings (0.5-2 m) and seedlings (< 0.5 m). The d.b.h. (diameter at breast height) of each fir tree in the plot was measured. At least one core of each tree was extracted at breast height (1.3 m above the ground) in the direction parallel to the contour line, using an increment borer. One additional core was extracted from the opposite side of some trees. These trees included the ones which were selected for developing chronologies and the ones with broken or rotten increment cores. For all the saplings and seedlings in the plot, the number of branch whorls and bud scars on the main stem was counted and recorded (Daniels and Veblen, 2004). In total, 208 increment cores from 119 trees were sampled and 615 saplings and seedlings were measured at both the subalpine treelines. In addition, 68 saplings and seedlings with normal growth were randomly uprooted after recording their height and the number of whorls and scars, and cross-sectional disks were cut from their base stems to determine their accurate ages. The time required for seedlings to reach the d.b.h. was estimated through the age-height regression (Dang et al., 2009). Shade-tolerant species can vary greatly in growth rate, for example, individuals in the seedling bank of some Abies species typically grow slowly and can persist for a very long period (Antos et al., 2005). Consequently, fir seedlings and saplings with severe release or suppression events (Nowacki and Abrams, 1997) were excluded from developing age-height regression (Dang et al., 2010).

Climate data

Climate data (1951–2009) were obtained from the nearest meteorological station, Songpan (32°39′N, 103°34′E, 2,580 m a.s.l.), approximately 55 km southwest of the sampling sites. Monthly mean temperature and monthly total precipitation of the sampling sites were simulated by Mountain Climate Simulator (MTCLIM, version 4.3; School of Forestry, University of Montana, Bozeman) based on daily maximum and minimum temperature and precipitation of base station (Thornton et al., 1997; Dang et al., 2009).

Chronology development

In the laboratory, all the increment cores were mounted in slotted wooden boards, air-dried and sanded with successively finer sandpaper to produce a polished transverse surface for visual cross-dating (Stokes and Smiley, 1996). The annual ring-widths were measured to the nearest 0.01 mm, and the rings were counted using a LINTAB II measuring system. False rings, missing rings or measurement mistakes were all identified by cross-dating the patterns of wide and narrow rings among trees with the software package TSAP-Win (Rinn, 2003). The quality of cross-dating was controlled using the software COFECHA (Holmes, 1983). Cores which were too short, fragmented or rotted were discarded in order to improve the common signals in tree ring-width sequences. In total, 73 cores from 38 trees and 59 cores from 41 trees were used to build the chronologies for the northwestern and southeastern treelines, respectively.

The ring-width chronologies were standardized using the program ARSTAN (Cook, 1985). Either negative exponential curves or straight lines were used to remove the effects of tree age (Fritts, 1976). If both of the curves were failed, a cubic smoothing spline with a 50% frequency-response cut-off of 80 years was applied. In addition, a univariate autoregressive model was used to remove the time series effect. The residual indices from autoregressive modeling of the detrended series were averaged by year with a robust mean calculation (Cook, 1985). The results were residual tree-ring chronologies, representing the common signal for the site. Descriptive statistics that were calculated for assessing the quality of the chronology included mean sensitivity (MS), standard deviation (SD), signal to noise ratio (SNR) and expressed population signal (EPS).

Dendroclimatic analysis

In this study, correlation coefficients between tree-ring chronologies with monthly mean temperature and monthly precipitation derived from the MTCLM program were used to identify the relationships between radial growth variability of A. faxoniana and climate variables (Tardif et al., 2001). Because climatic conditions in previous growing season often influence the radial growth in the following year (Fritts, 1976), temperature and precipitation beginning in previous June until current September (from 1951 to 2009) were used to analyze the relationships between annual radial growth and climate variables. The analyses were produced using the DENDROCLIM 2002 software (Biondi and Waikul, 2004).

Age structure and recruitment analysis

For the individual fir trees, the ages were determined using the cross-dated years from cores and the time to reach the coring height. If the cores passed close to the pith or missed inner parts, the number of rings missing from the pith was estimated by the geometric method (Duncan, 1989). Based on the relationship between age (y, year) and height (x, cm) of the fir seedlings/saplings (Figure 3):

y = 0.180 x + 4.615 R 2 = 0.728 , P < 0.001 , n = 68
(1)

the fir trees required 29 years to reach breast height (1.3 m above the ground). The ages of the seedlings and saplings were determined by the number of whorls and scars on the main stem (Dang et al.,2009). The ageing methodology may underestimate age by several years (Lv and Zhang, 2011). From the statistics of differences between the discs’ ages and the number of whorls/scars of the randomly uprooted seedlings and saplings, the results showed that only 20% of the samples were accurate. Of the samples, 15% were being underestimated one year, 20% were being underestimated two years and 45% were being underestimated three or more years. So the age frequency data for seedlings and saplings were smoothed using the following function (a method similar to the one of Daniels and Veblen, 2004):

x t = 0.20 f t + 0.15 f t + 1 + 0.20 f t + 2 + 0.45 f t + 3
(2)

where x is the smoothed age frequency, f is the original age frequency, and t is the year of seedling recruitment.

Figure 3
figure 3

Relationship between ages and height of A. faxoniana seedlings and saplings at the treeline ecotone in the Wanglang Nature Reserve.

Rates of recruitment vs mortality cannot be differentiated from a static age structure alone. So we detrended the time series of annually recruited seedlings and saplings to account for mortality by the best-fit theoretical distributions, the exponential and power functions (Hett and Loucks, 1976). The difference between the theoretical age frequency and observed age frequency provided a time series of residuals (Szeicz and Macdonald, 1995). Then the “seedling residuals” were used to assess climate-recruitment relationships. The influence of climate change on recruitment was tested using Pearson’s correlation coefficients between recruitment and climate data derived from MTCLIM for the 1951–2009 periods (Daniels and Veblen, 2004; Jump et al., 2007). The climate parameters were the same as those described for the dendroclimatic analysis of the fir trees. The time period of current year January-December was used in this analysis (Jump et al., 2007).

Results

Chronologies and descriptive statistics

Following standard procedures, residual chronologies for A. faxoniana in the treelines on both the northwestern and southeastern aspects and the statistics are presented in Figure 4 and Table 1. Both chronologies for the two aspects showed similar patterns with characteristic narrow rings in the same pointer years such as 1893, 1935–1937, 1967, 1982 and 1989 (Figure 4). All the descriptive statistics are higher in the chronology for the northwestern aspect than in that for the southeastern aspect (Table 1). Although the MS and SD are not very high, the common interval analysis indicates high SNR and EPS for the two chronologies, respectively (Table 1).

Figure 4
figure 4

Radial growth patterns of A. faxoniana at the treeline ecotones on the northwestern (a) and southeastern (b) aspects in the Wanglang Nature Reserve. Thin line and thick line represent the annual series and 10-year smoothing spline of each ring-width chronology, respectively. Dash line represents the number of cores included in each residual chronology.

Table 1 Descriptive parameters of residual chronologies of A. faxoniana at the treeline ecotones on the northwestern (NW) and southeastern (SE) aspects of the Min Mountains in the Wanglang National Nature Reserve

Radial growth trends with climate

The main climatic influence on the radial growth of A. faxoniana appears to be temperatures in the previous summer, especially in previous August, which were significantly negatively correlated with radial growth on both aspects (Figure 5). On the northwestern aspect, the radial growth of A. faxoniana was significantly negatively correlated with monthly mean temperature for the previous August, and was significantly positively correlated with monthly mean temperature for the current April (Figure 5a); and precipitation in the current January and September were significantly correlated positively and negatively with the current radial growth, respectively (Figure 5a). On the southeastern aspect, temperatures in the previous July and August were significantly negatively correlated with the radial growth (Figure 5b); and precipitation showed no significant correlation with ring-width indices (Figure 5b).

Figure 5
figure 5

Correlations between climate variables from 1951 to 2009 and chronologies for A. faxoniana at the treelines on the northwestern (a) and southeastern (b) aspects in the Wanglang Nature Reserve from June of the previous year to September of the current year. The gray bars represent monthly mean temperature and the black bars represent monthly total precipitation. * P < 0.05.

Distribution of age structure

The age structures of A. faxoniana (10-year intervals) on both northwestern and southeastern aspects are shown in Figure 6. On the northwestern aspect, the 21-30-year-old age class (fir trees that established in the 1980s) accounted for the largest age class (25%). The majority of fir trees (90%) were found during the period from the 1970s to the 2000s and their ages ranged from 1 to 40 years (Figure 6a).

Figure 6
figure 6

Age frequency distribution patterns (10-year intervals) of A. faxoniana at the treelines on the northwestern (a) and southeastern (b) aspects in the Wanglang Nature Reserve.

On the southeastern aspect, the 1-10-year-old age class (fir trees that established in the 2000s) formed the largest age class (11%). The majority of fir trees (78%) were successfully recruited into the treeline ecotone from 1750 s to 1900s (Figure 6b). However, there was no successful recruitment in successive decades from 1910s to 1970s.

Recent recruitment and climate

The short-time recruitment records in recent decades did not match with the available meteorological records (1951–2009) in the treeline on the southeastern aspect (Figure 6b), so we only assessed the climate-recruitment relationships on the northwestern aspect.

During the last 60 years, significantly positive correlations existed between recruitment and monthly mean temperatures for the current April, July and September. However, monthly precipitation had no significant correlations with the recruitment (Figure 7).

Figure 7
figure 7

Correlations between climate variables and recruitment residuals for A. faxoniana on the northwestern aspect from January to December (1951–2009). Legends as for Figure 5. #P < 0.01 and * P < 0.05.

Discussion

Correlations between radial growth and climate

Generally, the responses of radial growth to temperatures in previous summer were similar at the treeline ecotones on both the aspects. Significant negative correlations with previous July/August mean temperatures likely reflected the negative carry-over effects (Figure 5). Negative correlations between tree growth and previous-year July temperatures have been reported at high-elevation forests on southern (Liang et al., 2010) and southeastern (Lv and Zhang, 2011) QTP. High temperatures in the previous growing season (July/August) may reduce storage of assimilates for growth of next year (Liang et al., 2010). Moreover, fir radial growth was significantly positively correlated with monthly mean temperature in current April on the northwestern aspect (Figure 5a). Li et al. (2010) have also shown that radial growth of A. faxoniana was positively correlated with the temperature in current March at the treeline site in Wolong National Natural Reserve (30°53′N, 102°59′E), which lies west of our study area. Warmer temperatures early in the growing season benefit A. faxoniana growth by inducing early snowmelt and increasing available soil water (Wang et al., 2006; Dang et al., 2009), which causes early initiation of cambial activity and increases photosynthates (Splechtna et al., 2000; Case and Peterson, 2005).

The correlations between the precipitation and radial growth were different on both aspects (Figure 5). There was no significant correlation between precipitation and radial growth on the southeastern aspect (Figure 5b). On the northwestern aspect, the precipitation of the current January facilitated the radial growth (Figure 5a). This might represent relationships with snowpack and the subsequent effects on soil moisture (D’Arrigo et al. 2001. Significant negative correlation was found between radial growth and precipitation in the current September on the northwestern aspect, indicating that excessively sufficient precipitation will restrict the radial growth. Frequent rainfall in the end of the growing season of A. faxoniana might reduce solar radiation and effective photosynthesis, thereby, shorten the growing season (Wang et al., 2007; Jiang et al., 2010). The relationship between radial growth and precipitation on both the aspects might indicate that the precipitation was not the main limiting factor affecting radial growth during the main growing season.

The differential responses of radial growth to climate factors on both sites might be mainly due to the difference in micro-environmental conditions between the contrasting slope aspects. The studied area is dominated by the southeast monsoon (Pu et al., 2008; Yao et al., 2010). The southeastern aspect is wetter and warmer than the northwestern aspect, and also with stronger solar radiation. The sensitivity and response of treelines to climate variability may vary both in local and regional topographical conditions for the environmental variation (Holtmeier and Broll, 2005; Elliott and Kipfmueller, 2011). In our study, on the northwestern aspect with relatively weak solar radiation, the radial growth showed significantly negative correlation with precipitation in the end of the growing season (Figure 5a). The precipitation in January only favored the radial growth significantly on the northwestern aspect (Figure 5a) might because the difference of soil moisture on both the aspects. January snow might improve the soil moisture on the relatively drought northwestern slope in spring, therefore favored fir growth at the beginning of growing season.

Recruitment and climate correlations on the northwestern aspect

The age structure of a stand can provide a fairly accurate picture of temporal variations in the recruitment rate (Kullman, 1991) with the dynamics of climate change (Payette and Filion, 1985), because tree recruitment is more sensitive than tree mortality to climate variability (Camarero and Gutiérrez, 2004). Most of the fir trees in the treeline on the northwestern aspect were younger than 40 years, but only about 10% were successfully recruited into the stand before 1960 (Figure 6a), which showed that fir recruitment had a sporadic mode from 1670 to 1960. In the Min Mountains, the fir trees possessed good soil seed banks at treeline ecotones (Fang, 2006), which suggests that seed production has not been an important influencing factor on the age-structure distribution (Wang et al., 2006; Dang et al., 2009). Cheng et al. (2005) have revealed that the high mortality of A. faxoniana at the treeline of the Min Mountains on the eastern edge of QTP might be controlled mainly by temperature, wind, snowpack depth and winter drought.

Recruitment was mainly affected by the temperatures in the spring, summer and early autumn seasons, and it was found to have significantly positive correlations with mean temperatures in April, July and September (Figure 7). Temperatures in spring facilitating fir seedling establishment were also found in the treeline ecotones of the Shennongjia Mountains (Dang et al., 2009). High temperatures in spring are very important for germination (Camarero and Gutiérrez, 1999) and might simulate tree recruitment from soil seed banks (Dang et al., 2009). The significant positive correlations between summer temperatures and recruitment were also reported from conifer seedlings in alpine treeline ecotone of the Snowy Range in Wyoming USA (Germino and Smith, 1999) and from Abies spectabilis forest in the alpine timberline of the Mt. Everest in southern QTP, China (Lv and Zhang, 2011). Higher summer temperatures will strengthen the photosynthesis rate, which could encourage both growth and nonstructural carbon storage for fir seedlings to survive harsh winter climate (Camarero and Gutiérrez, 1999).

Recruitment on the southeastern aspect

There was a large gap in recruitment during the middle of the 20th century on the southeastern aspect (Figure 6). Competition from the dense bamboos might be mainly responsible for the rare recruitment (Dang et al., 2009). The bamboo cover in the southeastern plot (42.25%) was much higher than that in the northwestern plot (19.53%) (Zhao et al., 2012). Some studies have demonstrated that bamboos with a relatively high cover seem sufficient to impede tree establishment in subalpine forests (Takahashi, 1997; Holz and Veblen, 2006). In our study, the microenvironment on the southeastern aspect might be more suitable for bamboos’ cloning growth. The restriction of dense bamboos might exceed the facilitation of warm climate to the recruitment in recent decades.

Radial growth versus recruitment

Over the most recent decades, the radial growth of A. faxoniana had no obvious increasing trends with small fluctuations on both the aspects (Figure 4). However, the recruitment has increased sharply, especially on the northwestern aspect (Figure 6a). Obvious increases of tree recruitment with the recent climate change were also found in other treeline ecotones, such as European Alps (Leonelli et al., 2011), Swedish Scandes (Kullman, 2005) and U.S. Rocky Mountain (Elliott and Kipfmueller, 2011). Our result was accordant with the treeline dynamics about A. spectabilis on the southern QTP (Lv and Zhang, 2011), but was different from those obtained in similar studies carried out in the central Tianshan Mountains (Wang et al., 2006) and the Shennongjia Mountains (Dang et al., 2009). Wang et al. (2006) and Dang et al. (2009) have found more radial growth but less recruitment with Picea schrenkiana and Abies fargesii in response to warm climate over several of the most recent decades. The different results suggest that the response of treeline ecotones to climate change varies with both local site conditions and the individual species (Luckman and Kavanagh, 1998).

At many upper treeline ecotones the climatic conditions that facilitate seedling recruitment are frequently similar to those conducive to radial growth of trees (Szeicz and Macdonald, 1995; Camarero and Gutiérrez, 1999; Jump et al., 2007; Dang et al., 2009). In our study the climatic conditions that facilitated the A. faxoniana seedling recruitment were also similar to those that enhanced the radial growth of fir trees. For example, high temperatures in current April enhanced the recruitment of A. faxoniana and facilitated the radial growth (Figures 5 and 7); more January precipitation in current year strengthened both seedling recruitment and radial growth (Figures 5 and 7).

Conclusions

In summary, we investigated the climatic response of radial growth and recruitment of A. faxoniana in the treeline ecotones on the eastern edge of QTP. For the geographical novelty, QTP is attracting more and more attention on conducting similar studies (Liang et al., 2010; Lv and Zhang, 2011). In the future, continuous studies should be conducted, such as carrying out a second census within the same plots or setting up plots covering larger elevational gradients, in order to understand the dynamics of treeline position with the climatic variability on the eastern edge of QTP. Moreover, such studies will undoubtedly advance our knowledge of the other treelines on the eastern edge of QTP and reveal common and species-specific responses to climate change that may be exploited in conservation and protection efforts.

References

  • Andersen T, Carstensen J, Hernandez-Garcia E, Duarte CM: Ecological thresholds and regime shifts: approaches to identification. Trends Ecol Evol 2009, 24: 49–57. 10.1016/j.tree.2008.07.014

    Article  PubMed  Google Scholar 

  • Antos JA, Guest HJ, Parish R: The tree seedling bank in an ancient montane forest: stress tolerators in a productive habitat. J Ecol 2005, 93: 536–543. 10.1111/j.1365-2745.2005.00968.x

    Article  Google Scholar 

  • Barber VA, Juday GP, Finney BP: Reduced growth of Alaskan white spruce in the twentieth century from temperature-induced drought stress. Nature 2000, 405: 668–673. 10.1038/35015049

    Article  CAS  PubMed  Google Scholar 

  • Biondi F, Waikul K: DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree-ring chronologies. Comput. Geosci. (UK) 2004, 30: 303–311. 10.1016/j.cageo.2003.11.004

    Article  Google Scholar 

  • Bradley RS, Jones PD: ‘Little Ice Age’ summer temperature variations: Their nature and relevance to recent global warming trends. Holocene 1993, 3: 367–376. 10.1177/095968369300300409

    Article  Google Scholar 

  • Brubaker LB: Responses of tree populations to climatic change. Plant Ecol 1986, 67: 119–130. 10.1007/BF00037362

    Article  Google Scholar 

  • Bunn AG, Waggoner LA, Graumlich LJ: Topographic mediation of growth in high elevation foxtail pine ( Pinus balfouriana Grev. et Balf.) forests in the Sierra Nevada, USA. Global Ecol Biogeogr. 2005, 14: 103–114. 10.1111/j.1466-822X.2005.00141.x

    Article  Google Scholar 

  • Camarero JJ, Gutiérrez E: Structure and recent recruitment at alpine forest-pasture ecotones in the Spanish central Pyrenees. Ecoscience 1999, 6: 451–464.

    Google Scholar 

  • Camarero JJ, Gutiérrez E: Pace and pattern of recent treeline dynamics: Response of ecotones to climatic variability in the Spanish Pyrenees. Clim Change 2004, 63: 181–200.

    Article  Google Scholar 

  • Case MJ, Peterson DL: Fine-scale variability in growth-climate relationships of Douglas-fir, North Cascade Range. Washington. Can. J. For. Res. 2005, 35: 2743–2755. 10.1139/x05-191

    Article  Google Scholar 

  • Cheng W, Wu N, Luo P: Survival analysis of Abies faxoniana populations near timberline on the upper Minjiang River. Acta Phytoecol. Sin. 2005, 29: 349–353.

    Google Scholar 

  • Cook ER: A time series analysis approach to tree ring standardization (Dendrochronology, Forestry, Dendroclimatology, Autoregressive Process). Doctor’s thesis. Tucson, Arizona: The University of Arizona; 1985.

    Google Scholar 

  • Cullen LE, Palmer JG, Duncan RP, Stewart GH: Climate change and tree-ring relationships of Nothofagus menziesii tree-line forests. Can J For Res 2001, 31: 1981–1991.

    Article  Google Scholar 

  • D’Arrigo RD, Schuster WSF, Lawrence DM, Cook ED, Wiljanen M: Climate-growth relationships of eastern hemlock and chestnut oak from Black Rock Forest in the highlands of southeastern New York. Tree-Ring Res. 2001, 57: 131–191.

    Google Scholar 

  • Danby RK, Hik DS: Variability, contingency and rapid change in recent subarctic alpine tree line dynamics. J Ecol 2007, 95: 352–363. 10.1111/j.1365-2745.2006.01200.x

    Article  Google Scholar 

  • Dang HS, Zhang KR, Zhang YJ, Tan SD, Jiang MX, Zhang QF: Tree-line dynamics in relation to climate variability in the Shennongjia Mountains, central China. Can J For Res 2009, 39: 1848–1858. 10.1139/X09-106

    Article  Google Scholar 

  • Dang HS, Zhang YJ, Zhang KR, Jiang MX, Zhang QF: Age structure and regeneration of subalpine fir ( Abies fargesii ) forests across an altitudinal range in the Qinling Mountains. China. Forest Ecol. Manag. 2010, 259: 547–554. 10.1016/j.foreco.2009.11.011

    Article  Google Scholar 

  • Daniels LD, Veblen TT: Spatiotemporal influences of climate on altitudinal treeline in northern Patagonia. Ecology 2004, 85: 1284–1296. 10.1890/03-0092

    Article  Google Scholar 

  • Ding YH, Lin ED, He JK: Climate change research over China-Science, Impact, Adaptation and Strategy Policy. Beijing: China Environmental Science Press; 2009:455.

    Google Scholar 

  • Duncan RP: An evaluation of errors in tree age estimates based on increment cores in kahikatea (Dacrycarpus dacrydioides). New Zealand Nat. Sci. 1989, 16: 31–37.

    Google Scholar 

  • Earle CJ: Forest dynamics in a forest-tundra ecotone, Medicine Bow Mountains, Wyoming. Doctor’s thesis. Seattle: University of Washington; 1993.

    Google Scholar 

  • Eckstein D, Aniol RW: Dendroclimatological reconstruction of the summer temperatures for an alpine region. Mitt Forstl Bundesvers, Wien. 1981, 142: 391–398.

    Google Scholar 

  • Elliott GP: Influences of 20th-century warming at the upper treeline contingent on local-scale interactions: Evidence from a latitudinal gradient in the Rocky Mountains. USA. Global Ecol. Biogeogr. 2011, 20: 46–57. 10.1111/j.1466-8238.2010.00588.x

    Article  Google Scholar 

  • Elliott GP, Baker WL: Quaking aspen (Populus tremuloides Michx.) at treeline: A century of change in the San Juan Mountains, Colorado. USA. J. Biogeogr. 2004, 31: 733–745. 10.1111/j.1365-2699.2004.01064.x

    Article  Google Scholar 

  • Elliott GP, Kipfmueller KF: Multi-scale influences of slope aspect and spatial pattern on ecotonal dynamics at upper treeline in the Southern Rocky Mountains. USA. Arct. Antarct. Alp. Res. 2010, 42: 45–56. 10.1657/1938-4246-42.1.45

    Article  Google Scholar 

  • Elliott GP, Kipfmueller KF: Multiscale influences of climate on upper treeline dynamics in the southern rocky mountains, USA: evidence of intraregional variability and bioclimatic thresholds in response to twentieth-century warming. Ann Assoc Am Geogr 2011, 101: 1181–1203. 10.1080/00045608.2011.584288

    Article  Google Scholar 

  • Fang JQ: Soil seed banks and seedling recruitments of the dominant tree species near the alpine treeline ecotone on the upper reaches of Minjiang River, China. Master’s thesis. Chengdu: Chengdu Institute of Biology, Chinese Academy of Science; 2006.

    Google Scholar 

  • Fang KY, Gou XH, Chen FH, Peng JF, D’Arrigo R, Wright W, Li MH: Response of regional tree-line forests to climate change: evidence from the northeastern Tibetan Plateau. Trees-Struct. Funct. 2009, 23: 1321–1329. 10.1007/s00468-009-0373-5

    Article  Google Scholar 

  • Fritts HC: Tree rings and climate. New York: Academic Press; 1976:567.

    Google Scholar 

  • Germino MJ, Smith WK: Sky exposure, crown architecture, and low-temperature photoinhibition in conifer seedlings at alpine treeline. Plant Cell Environ 1999, 22: 407–415. 10.1046/j.1365-3040.1999.00426.x

    Article  Google Scholar 

  • Gervais BR, MacDonald GM: A 403-year record of July temperatures and treeline dynamics of Pinus sylvestris from the Kola Peninsula, northwest Russia. Arct Antarct Alp Res 2000, 32: 295–302. 10.2307/1552528

    Article  Google Scholar 

  • Grace J, Berninger F, Nagy L: Impacts of climate change on the treeline. Ann. Bot.-London 2002, 90: 537–544. 10.1093/aob/mcf222

    Article  CAS  Google Scholar 

  • Harsch MA, Hulme PE, McGlone MS, Duncan RP: Are treelines advancing? A global meta-analysis of treeline response to climate warming. Ecol Lett 2009, 12: 1040–1049. 10.1111/j.1461-0248.2009.01355.x

    Article  PubMed  Google Scholar 

  • Hessl AE, Baker WL: Spruce and fir regeneration and climate in the forest-tundra ecotone of Rocky Mountain National Park, Colorado, U.S.A. Arct Antarct Alp Res 1997, 29: 173–183. 10.2307/1552044

    Article  Google Scholar 

  • Hett JM, Loucks OL: Age structure models of balsam fir and eastern hemlock. J Ecol 1976, 64: 1029–1044. 10.2307/2258822

    Article  Google Scholar 

  • Holmes RL: Computer-assisted quality control in tree-ring dating and measurement. Tree-ring Bull. 1983, 43: 69–78.

    Google Scholar 

  • Holtmeier FK: Mountain timberlines: Ecology, patchiness and dynamics. New York: Springer; 2009:384.

    Book  Google Scholar 

  • Holtmeier FK, Broll G: Sensitivity and response of northern hemisphere altitudinal and polar treelines to environmental change at landscape and local scales. Global Ecol. Biogeogr. 2005, 14: 395–410. 10.1111/j.1466-822X.2005.00168.x

    Article  Google Scholar 

  • Holtmeier FK, Broll G: Treeline advance-driving processes and adverse factors. Landscape Online 2007, 1: 1–33.

    Article  Google Scholar 

  • Holz A, Veblen TT: Tree regeneration responses to Chusquea montana bamboo die-off in a subalpine Nothofagus forest in the southern Andes. J Veg Sci 2006, 17: 19–28.

    Google Scholar 

  • Hou Y, Wang CY, Zhang YB, Zou CJ: Effects of elevated CO 2 concentration and temperature on root structure of subalpine Betula albosinensis seedlings in western Sichuan, southwestern China. Journal of Beijing Forestry University 2008, 30: 29–33.

    Google Scholar 

  • Jacoby GC, D’Arrigo RD: Tree ring width and density evidence of climatic and potential forest change in Alaska. Global Biogeochem. Cy. 1995, 9: 227–234. 10.1029/95GB00321

    Article  CAS  Google Scholar 

  • Jiang SH, Ren LL, Yong B, Liu L, Yang HF: Reconstruction of precipitation series for Ma’an Mountain area from July to September using tree rings. Journal of Hohai University (Natural Sciences) 2010, 38: 1–35.

    Google Scholar 

  • Jump AS, Hunt JM, Penuelas J: Climate relationships of growth and recruitment across the altitudinal range of Fagus sylvatica in the Montseny Mountains, northeast Spain. Ecoscience 2007, 14: 507–518. 10.2980/1195-6860(2007)14[507:CROGAE]2.0.CO;2

    Article  Google Scholar 

  • Körner C: Alpine plant life. Berlin: Springer; 2003:344.

    Book  Google Scholar 

  • Kullman L: Structural change in a subalpine birth woodland in north Sweden during the past century. J Biogeogr 1991, 18: 53–62. 10.2307/2845244

    Article  Google Scholar 

  • Kullman L: Tree limit dynamics of betula pubescens ssp. tortuosa in relation to climate variability: Evidence from central Sweden. J Veg Sci 1993, 4: 765–772. 10.2307/3235613

    Article  Google Scholar 

  • Kullman L: Early holocene tree growth at a high elevation site in the northernmost scandes of Sweden (Lapland): A palaeobiogeographical case study based on megafossil evidence. Geogr. Ann. A 1999, 81: 63–74. 10.1111/j.0435-3676.1999.00049.x

    Article  Google Scholar 

  • Kullman L: Pine (Pinus sylvestris) treeline dynamics during the past millennium - a population study in west-central Sweden. Ann Bot Fenn 2005, 42: 95–106.

    Google Scholar 

  • Kullman L, Öberg L: Post-little Ice Age tree line rise and climate warming in the Swedish scandes: a landscape ecological perspective. J Ecol 2009, 97: 415–429. 10.1111/j.1365-2745.2009.01488.x

    Article  Google Scholar 

  • LaMarche VC, Fritts HC: Tree rings, glacial advance, and climate in the Alps. Z. Gletscherk. Glazialgeol. 1971, 7: 125–131.

    Google Scholar 

  • Lavoie C, Payette S: Recent fluctuations of the lichen-spruce forest limit in subarctic Quebec. J Ecol 1994, 82: 725–734. 10.2307/2261438

    Article  Google Scholar 

  • Leonelli G, Pelfini M, di Cella UM, Garavaglia V: Climate warming and the recent treeline shift in the European Alps: the role of geomorphological factors in high-altitude sites. Ambio 2011, 40: 264–273. 10.1007/s13280-010-0096-2

    Article  PubMed Central  PubMed  Google Scholar 

  • Li JB, Cook ER, D’Arrigo R, Chen FH, Gou XH, Peng JF, Huang JG: Common tree growth anomalies over the northeastern Tibetan Plateau during the last six centuries: implications for regional moisture change. Global Change Biol. 2008, 14: 2096–2107. 10.1111/j.1365-2486.2008.01603.x

    Article  Google Scholar 

  • Li ZS, Liu GH, Zhang QB, Hu CJ, Huo SZ, Liu XL, He F: Tree ring reconstruction of summer temperature variations over the past 159 years in Wolong National Natural Reserve, western Sichuan, China. Acta Phytoecologica Sinica 2010, 34: 628–641.

    Google Scholar 

  • Liang EY, Wang YF, Xu Y, Liu B, Shao XM: Growth variation in Abies georgei var. smithii along altitudinal gradients in the Sygena mountains, southeastern Tibetan Plateau. Trees 2010, 24: 363–373. 10.1007/s00468-009-0406-0

    Article  Google Scholar 

  • Liu LS, Shao XM, Liang EY: Climate signals from tree ring chronologies of the upper and lower treelines in the Dulan region of the Northeastern Qinghai-Tibetan Plateau. Journal of Integrative Plant Biol. 2006, 48: 278–285. 10.1111/j.1744-7909.2006.00158.x

    Article  Google Scholar 

  • Liu WH, Gou XH, Yang MX, Zhang Y, Fang KY, Yang T, Jin LY: Drought reconstruction in the Qilian Mountains over the last two centuries and its implications for large-scale moisture patterns. Adv. Atmos. Sci. 2009, 26: 621–629. 10.1007/s00376-009-9028-0

    Article  Google Scholar 

  • Lloyd AH, Fastie CL: Spatial and temporal variability in the growth and climate response of treeline trees in Alaska. Clim Change 2002, 52: 481–509. 10.1023/A:1014278819094

    Article  Google Scholar 

  • Lloyd AH, Fastie CL: Recent changes in treeline forest distribution and structure in interior Alaska. Ecoscience 2003, 10: 176–185.

    Google Scholar 

  • Lloyd AH, Graumlich LJ: Holocene dynamics of treeline forests in the Sierra Nevada. Ecology 1997, 78: 1199–1210. 10.1890/0012-9658(1997)078[1199:HDOTFI]2.0.CO;2

    Article  Google Scholar 

  • Luckman BH, Kavanagh TA: Documenting the effects of recent climate change at treeline in the Canadian Rockies. In The impacts of climate variability on forests. Lecture notes in Earth Sciences. Edited by: Beniston M, Innes JL. Berlin: Springer; 1998:121–144.

    Chapter  Google Scholar 

  • Lv LX, Zhang QB: Asynchronous recruitment history of Abies spectabilis along an altitudinal gradient in the Mt. Everest region. J. Plant Ecol. 2011, 5: 147–156.

    Article  Google Scholar 

  • Moiseev PA: Effects of climatic changes on radial increment and age structure formation in high-mountain larch forests of the Kuznetsk Ala Tau. Russ J Ecol 2002, 33: 7–13. 10.1023/A:1013659319452

    Article  Google Scholar 

  • Munier A, Hermanutz L, Jacobs JD, Lewis K: The interacting effects of temperature, ground disturbance, and herbivory on seedling recruitment: implications for treeline advance with climate warming. Plant Ecol 2010, 210: 19–30. 10.1007/s11258-010-9724-y

    Article  Google Scholar 

  • Nowacki GJ, Abrams MD: Radial-growth averaging criteria for reconstructing disturbance histories from presettlement-origin oaks. Ecol Monogr 1997, 67: 225–249.

    Google Scholar 

  • Payette S, Filion L: White spruce expansion at the treeline and recent climate change. Can J For Res 1985, 15: 241–251. 10.1139/x85-042

    Article  Google Scholar 

  • Pu YL, Liu SQ, Zhang SR, Long GF, Lu CT: Slope-directive variation of mountain soil basic attributes in the northern Hengduan Mountains regions. Journal of Soil and Water Conservation 2008, 22: 112–117.

    Google Scholar 

  • Rinn F: TSAPWin: Time series analysis and presentation for dendrochronology and related applications, version 0.55 user reference. Heidelberg, Germany: RINNTECH; 2003.

    Google Scholar 

  • Ruffner CM, Abrams MD: Relating land-use history and climate to the dendroecology of a 326-year-old Quercus prinus talus slope forest. Can J For Res 1998, 28: 347–358. 10.1139/x97-220

    Article  Google Scholar 

  • Schweingruber FH, Briffa KR, Jones PD: Yearly maps of summer temperatures in Western Europe from A.D. 1750 to 1975 and western North America from 1600 to 1982: Results of a radiodensitometrical study on tree rings. Vegetatio 1991, 92: 5–71.

    Google Scholar 

  • Splechtna BE, Dobrys J, Klinka K: Tree-ring characteristics of subalpine fir ( Abies lasiocarpa (Hook.) Nutt.) in relation to elevation and climatic fluctuations. Ann For Sci 2000, 57: 89–100. 10.1051/forest:2000105

    Article  Google Scholar 

  • Stevens GC, Fox JF: The causes of treeline. Annu Rev Ecol Syst 1991, 22: 177–191. 10.1146/annurev.es.22.110191.001141

    Article  Google Scholar 

  • Stokes MA, Smiley TL: An introduction to tree-ring dating. Tucson, Arizona: University of Arizona Press; 1996:77.

    Google Scholar 

  • Szeicz JM, Macdonald GM: Recent white spruce dynamics at the subarctic alpine treeline of north-western Canada. J Ecol 1995, 83: 873–885. 10.2307/2261424

    Article  Google Scholar 

  • Takahashi K: Regeneration and coexistence of two subalpine conifer species in relation to dwarf bamboo in the understorey. J Veg Sci 1997, 8: 529–536. 10.2307/3237203

    Article  Google Scholar 

  • Takahashi K, Azuma H, Yasue K: Effects of climate on the radial growth of tree species in the upper and lower distribution limits of an altitudinal ecotone on Mount Norikura, central Japan. Ecol Res 2003, 18: 549–558. 10.1046/j.1440-1703.2003.00577.x

    Article  Google Scholar 

  • Tardif J, Conciatori F, Bergeron Y: Comparative analysis of the climatic response of seven boreal tree species from northwestern Quebec. Canada. Tree-Ring Res. 2001, 57: 169–181.

    Google Scholar 

  • Taylor AH, Jiang SW, Zhao LJ, Liang CP, Miao CJ, Huang JY: Regeneration patterns and tree species coexistence in old-growth Abies-Picea forests in southwestern China. Forest Ecol. Manag. 2006, 223: 303–317. 10.1016/j.foreco.2005.11.010

    Article  Google Scholar 

  • Thornton PE, Running SW, White MA: Generating surfaces of daily meteorological variables over large regions of complex terrain. J Hydrol 1997, 190: 214–251. 10.1016/S0022-1694(96)03128-9

    Article  Google Scholar 

  • Tranquillini W: Physiological ecology of the alpine treeline. New York: Springer; 1979.

    Book  Google Scholar 

  • Wang T, Zhang QB, Ma KP: Treeline dynamics in relation to climatic variability in the central Tianshan Mountains, northwestern China. Global Ecol. Biogeogr. 2006, 15: 406–415. 10.1111/j.1466-822X.2006.00233.x

    Article  Google Scholar 

  • Wang YX, Lu ZH, Su HX, Sang WG: Response analysis between tree-ring widths and climatic factors along different regions in Tianshan Mountains northwestern China. Journal of China University of Mining & Technology 2007, 36: 251–256.

    Google Scholar 

  • Weisberg PJ, Baker WL: Spatial variation in tree regeneration in the forest-tundra ecotone, Rocky Mountain National Park. Colorado. Can. J. For. Res. 1995, 25: 1326–1339. 10.1139/x95-145

    Article  Google Scholar 

  • Wilmking M, Juday GP, Barber VA, Zald H: Recent climate warming forces contrasting growth responses of white spruce at treeline in Alaska through temperature thresholds. Global Change Biol. 2004, 10: 1–13. 10.1111/j.1365-2486.2003.00715.x

    Article  Google Scholar 

  • Yao YH, Zhang BP, Han F, Pang Y: Spatial pattern and exposure effect of altitudinal belts in the Hengduan Mountains. J. Mount. Sci. 2010, 28: 11–20.

    Google Scholar 

  • Zhao ZJ, Eamus D, Yu Q, Li Y, Yang HW, Li JQ: Climate constraints on growth and recruitment patterns of Abies faxoniana over altitudinal gradients in the Wanglang Natural Reserve, eastern Tibetan Plateau. Aust J Bot 2012, 60: 602–614. 10.1071/BT12051

    Article  Google Scholar 

Download references

Acknowledgements

This study was supported by grants from the Special Research Program for Public-Welfare Forestry of China (No. 200804001). We thank the Wanglang National Nature Reserve for permission and assistance in the field sampling. We are grateful to Yue Li, Hong-Wei Yang, Yu Ruan and Shen Zhao who helped in fieldwork, and to Prof. Xiu-Hai Zhao and Dr. Lu-Shuang Gao for guidance in the LINTAB system.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jun-Qing Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZJZ carried out the indoor work and tree-ring analysis, ZJZ, GZS, LYT, DWK, MJW and WK worked together in the field, ZJZ, GZS and JQL conceived and designed the experiments, and drafted the manuscript, WXG, MJBZ and QY revised the manuscript; All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhao, ZJ., Shen, GZ., Tan, LY. et al. Treeline dynamics in response to climate change in the Min Mountains, southwestern China. Bot Stud 54, 15 (2013). https://doi.org/10.1186/1999-3110-54-15

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1999-3110-54-15

Keywords