induced upward migration of the alpine treeline ... - Wiley Online Library

1 downloads 0 Views 1MB Size Report
Warming-induced upward migration of the alpine treeline in the Changbai Mountains, northeast China. Haibo Du1 | Jie Liu1 | Mai-He Li2,3 | Ulf B¨untgen2,4,5,6 ...
Received: 1 August 2017

|

Accepted: 5 October 2017

DOI: 10.1111/gcb.13963

PRIMARY RESEARCH ARTICLE

Warming-induced upward migration of the alpine treeline in the Changbai Mountains, northeast China € ntgen2,4,5,6 | Yue Yang1 | Lei Wang1,7 | Haibo Du1 | Jie Liu1 | Mai-He Li2,3 | Ulf Bu Zhengfang Wu1 | Hong S. He1,8 1

School of Geographical Sciences, Northeast Normal University, Changchun, China

Abstract Treeline responses to environmental changes describe an important phenomenon in

2

Swiss Federal Research Institute WSL, Birmensdorf, Switzerland

global change research. Often conflicting results and generally too short observa-

3

tions are, however, still challenging our understanding of climate-induced treeline

Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang, China 4

Department of Geography, University of Cambridge, Cambridge, UK 5 Department of Geography, Masaryk University, Brno, Czech Republic 6

Global Change Research Institute CAS, Brno, Czech Republic 7 Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun, China 8

School of Natural Resources, University of Missouri, Columbia, MO, USA

dynamics. Here, we use a state-of-the-art dendroecological approach to reconstruct long-term changes in the position of the alpine treeline in relation to air temperature at two sides in the Changbai Mountains in northeast China. Over the past 160 years, the treeline increased by around 80 m, a process that can be divided into three phases of different rates and drives. The first phase was mainly influenced by vegetation recovery after an eruption of the Tianchi volcano in 1702. The slowly upward shift in the second phase was consistent with the slowly increasing temperature. The last phase coincided with rapid warming since 1985, and shows with 33 m per 1°C, the most intense upward shift. The spatial distribution and age structure of trees beyond the current treeline confirm the latest, warming-induced

Correspondence Hong S. He and Zhengfang Wu, School of Geographical Sciences, Northeast Normal University, Changchun, China. Emails: [email protected]; [email protected] Funding information National Key R&D Program of China, Grant/ Award Number: 2016YFA0602301; National Natural Science Foundation of China, Grant/ Award Number: 41601052, 41471085, 41371076; China Postdoctoral Science Foundation, Grant/Award Number: 2016M600225; Fundamental Research Funds for the Central Universities, Grant/ Award Number: 2412016KJ025; Youth and Sports of Czech Republic within the National Sustainability Program I, Grant/ Award Number: LO1415

upward shift. Our results suggest that the alpine treeline will continue to rise, and that the alpine tundra may disappear if temperatures will increase further. This study not only enhances mechanistic understanding of long-term treeline dynamics, but also highlights the effects of rising temperatures on high-elevation vegetation dynamics. KEYWORDS

altitudinal transect, Betula ermanii, Changbai Mountains, climate change, dendroecology, forest growth, treeline dynamics

1 | INTRODUCTION

Hik, 2007; Kammer et al., 2009; Kelly & Goulden, 2008; Kirdyanov gout, Marquet, De Ruffray, & Brisse, 2008). Some et al., 2011; Lenoir, Ge

The upper distribution of alpine ecosystems is highly sensitive to tem-

disparities or even downward shifts from more drought-prone environ-

€rner & Riedl, 2012). Upward shifts of the alpine perature changes (Ko

ments have recently been reported as well (Crimmins, Dobrowski,

treeline in relation to global warming have been observed in many

Greenberg, Abatzoglou, & Mynsberge, 2011; Foster & D’amato, 2015;

mountain systems around the world (Beckage et al., 2008; Danby &

Hartl-Meier, Dittmar, Zang, & Rothe, 2014; Rabasa et al., 2013).

1256

|

© 2017 John Wiley & Sons Ltd

wileyonlinelibrary.com/journal/gcb

Glob Change Biol. 2018;24:1256–1266.

DU

|

ET AL.

1257

The driving factors of treeline dynamics are complex (Malanson

mountain system (Foster & D’amato, 2015) although the treeline

et al., 2007; Tognetti & Palombo, 2013). Warming may promote the

position (i.e. elevation) varies with slope aspect. Moreover, air warm-

upward shift of treelines because productivity and functioning of

ing will release the low temperature constraints on tree flowering,

€ rner & Riedl, such ecotones are generally temperature-limited (Ko

pollination, seed maturation, seed dispersal, seed germination,

2012; Pauli et al., 2012). Microbial activity and plant metabolism are

€ rner & Riedl, 2012; Morin, Auggrowth and survival of seedlings (Ko

both temperature-dependent (Zhu, Cao, Wang, Xiao, & Li, 2012;

spurger, & Chuine, 2007), which will then result in expansion of

Zhu, Xiang, Wang, & Li, 2012). Tree recruitment and growth within

young recruits above the treeline (Truong, Palm, & Felber, 2007).

€ rner & the ecotone are mainly limited by low temperatures (Ko

We thus further hypothesize that tree populations will be more

€ rner & Riedl, 2012; Li, Kr€auchi, & Dobbertin, Paulsen, 2004; Ko

abundant and occur beyond the current treeline with global warm-

2006; Li, Yang, & Kr€auchi, 2003). Recent studies, however, demon-

ing. The third hypothesis reinforces the first two hypotheses if they

strated that species interaction (Liang et al., 2016) and geomorphic

are accepted.

processes (Macias-Fauriaa & Johnsonb, 2013) can slow down

To test these three hypotheses, we analysed the treeline shifts

warming-induced upward shifts, and that life stage (Malis et al.,

and the corresponding temperature changes in the last 160 years on

2016) and land-use legacies (Ameztegui, Coll, Brotons, & Ninot,

two sides of the Changbai Mountains and investigated the distribu-

2016), rather than warming, control treeline shifts. It has also been

tion and age structure of trees above the current treeline (TACT).

suggested that changes in soil nutrients, permafrost depth and snow

We aimed to answer the following questions: (i) How did tempera-

cover, mediated by topography and wind exposure (Frost & Epstein, € 2014; Kullman & Oberg, 2009; Lloyd, 2005; Macias-Fauriaa &

ture change and the treeline shift during the last 160 years? (ii) Did

Johnsonb, 2013; Richardson, 2004; Wilmking et al., 2012), influence

growth of TACT reflect temperature changes on the Changbai

treeline movements.

Mountains?

temperature changes coincide with treeline shifts? (iii) Did the

Previous studies mainly used remote sensing images (Ameztegui et al., 2016), historic photographs (Hagedorn et al., 2014) and treering data (Devi et al., 2008; Elliott, 2011) to identify treeline locations or changes in tree density in relation to global warming. The detected changes were often correlated against changes in tempera-

2 | MATERIALS AND METHODS 2.1 | Study area

ture to derive the causal effects of temperature on treeline dynamics

The Changbai Mountains (41°410 49″ to 42°250 18″N and 127°420 55″

(Gaire, Koirala, Bhuju, & Borgaonkar, 2014), or simply on variation in

to 128°160 48″E) is located in the northeast China at the border to

stem density within the upper ecotone (Liang, Wang, Eckstein, &

North Korea (Cao & Guo, 2000) (Figure 1). Its altitude ranges from

Luo, 2011; Wang et al., 2016). Although remote sensing images have

713 to 2,691 m above sea level (a.s.l.), and the climate is temperate

proven to be efficient at analysing large-scale and consecutive tree-

continental, with annual mean temperatures from 7.3 to 4.9°C and

line variation, they are mainly limited in their temporal scale that is

annual precipitation from 800 to 1,800 mm. The Changbai Moun-

100 years) and high temporal resolution data on treeline position,

2.2 | Historical treeline reconstruction

in tandem with reliable meteorological measurements. In addition,

We established two altitudinal transects on the north side of the

although global temperature has increased over the last decades, the

Changbai Mountains in 2014 (North1 and North2) and two transects

increase was not uniform around the world, with some areas experi-

on the west side (West1 and West2) in 2015 (Figure 1 and Table 1).

encing no increase or even cooling (Harsch et al., 2009). Conse-

North1 was a protected slope and North2 was an exposed slope,

quently, upward shifts of treeline did not occur uniformly (Holtmeier

whereas West1 and West2 were both exposed slopes. We did not

& Broll, 2007) and particularly corresponded to local temperature

sample the protected slope of the west side of the Changbai Moun-

€rner, 2007). Therefore, we hypothesize that the rate of change (Ko

tains because of safety issues. Instead, we used an unmanned aerial

treeline shift is similar for all sides of a mountain because tempera-

vehicle and measured the altitude of the current treeline of the pro-

ture change (warming or cooling) is generally uniform in a single

tected slope (Fig. S6). Selecting exposed and protected slopes

1258

|

DU

ET AL.

F I G U R E 1 Geographical features of the study area at regional and local scales. (a) Satellite image of the study area showing the location of sampling sites (red and black dots, and blue box); (b) location of the study area at the regional scale and (c) observed summer mean temperature (red line) and summer precipitation (blue bars) of the study area during 1959–2015. TACT represents the trees above the current treeline. North1, north side transect 1; North2, north side transect 2; West1, west side transect 1; West2, west side transect 2 [Colour figure can be viewed at wileyonlinelibrary.com]

allowed us to account for the terrain effects on treeline altitudes

model (described later). We found that the age of the current tree-

(Zhang, Zhang, & Mian, 2016; Zhao et al., 2016). These transects

line trees was ~10 years according to this model. We thus defined

were from the current treeline downward to lower altitudes at 10 m

the historic treeline position (altitude) as the location where B. er-

intervals. The plot size (100 m2) was 5 m (altitudinal direc-

manii had established and grown for 10 years. Therefore, the time of

tion) 9 20 m (along contour line). We defined the current treeline

tree establishment plus 10 years was the time when the treeline

position as the canopy cover of ≥0.2 (20%) of trees with an average

reaches a given altitude. Because the average tree’s canopy in the

height of >2 m. For the small trees in the current treeline, tree age

current treeline is about 1 m2, our treeline definition meant that a

was determined by a tree age-DB (diameter at base) regression

treeline plot (100 m2) required a minimum 20 such trees.

DU

|

ET AL.

We sampled increment cores or discs at various altitudes on the two sides of the Change Mountains to determine how tree age

1259

2.3 | Sampling of trees above the current treeline

changes with altitude. We used the locations of the old tree ages to

To analyse the distribution and age structure of B. ermanii trees

delineate the treeline locations of different historic periods. For his-

above the current treeline (TACTs), in 2015 we established a big

toric treelines at lower altitudes, where the current tree height and

B. ermanii plot with 450 9 50 m (altitudinal length 9 width) ranging

canopy cover exceed the treeline criterion very much (Table 1), we

from 2,074 to 2,194 m a.s.l. above the current treeline on the west

assumed that the same treeline criterion was met at some time in the

side of the Changbai Mounts (Figure 1). We counted all B. ermanii

past. Thus, the historic treeline position was determined using the

TACTs in this plot by dividing the plot into nine subplots of

oldest 20 trees in the 100 m2 plot. The corresponding time and tree-

50 9 50 m and then recording B. ermanii in each subplot. In total, we

line reaching that plot location was determined by the ages of the 20

found 598 B. ermanii TACTs. Most TACTs were smaller than 1.3 m in

trees. We randomly sampled at least five such trees to determine

height, and therefore we measured diameter at base (DB) for all

their ages. Tree age changes by altitude were revealed by averaging

TACTs. Because these trees were too small to drill cores, we ran-

the samples in each altitudinal plot. The largest one or two trees,

domly selected 35 nearby B. ermanii, measured DB of each tree, and

however, were abnormally larger than other trees in some plots

took their trunks to laboratory to determine their ages using TSAPwin

(random establishment before large number of trees moved up and

system (LINTAB 6). We built regression models of tree age with DB

formed the treelines), and these trees were excluded from determin-

and tree height for the 35 trees. Compared to the tree age-tree

ing tree ages. In total, 186 trees were sampled in the four transects.

height regression model, the tree age-DB regression model performed

We drilled two tree cores per tree, parallel to contour, at 1.3 m for trees with DBH >5 cm and height >2 m. Tree cores were held in

better (Appendix S1, Fig. S2). Therefore, we estimated the ages of all B. ermanii in this plot using the regression of tree age with DB.

self-made paper tubes in the field and were mounted as wooden stripes in laboratory. We sanded the cores using coarse-to-fine sandpaper. Cores were then measured using the TSAPwin4.7 system (LINTAB 6, with precision of 0.001 mm) and cross-dated using the

2.4 | Deriving weather data of the past five decades for the current treeline areas

COFECHA program (Holmes, 1983). For some cores hitting the pith,

Observed climatic data from 1959 onward (the Tianchi weather sta-

the distance to the centre of the tree was estimated by fitting a cir-

tion, 42°10 N, 128°50 E, 2,623 m near the volcanic lake on the moun-

cular template to the innermost curved ring or by subtracting core

tain top) were provided by the China meteorological data network

length from the radius at core height (Hollingsworth, Walker, Chapin

(http://data.cma.gov.cn/). We performed quality control based on

Iii, & Parsons, 2006). The years of this segment were determined

the method by Alexander et al. (2006) and performed homogeneity

from the regression of tree age with DB (discussed later; also see

test using the RHtestV4 software (Wang & Feng, 2013) (the nearby

Appendix S1). Cores with unrecognized rings were rejected, and ser-

reference stations were also provided by the China meteorological

ies not correlated with the entire dataset were also discarded. We

data network). The data proved to be homogeneous by this test;

then reconstructed treeline changes using the altitudinal changes in

however, the altitude difference between our study sites and Tianchi

tree age. For North1, however, we sampled tree cores along the

station was between 400 and 700 m, and the temperature differ-

ridge upward to the mountain-top rather than from the valley

ence was notable. Therefore, we set up an automatic meteorological

upward to the mountain-top. Many trees in lower altitudes were

station (Weather hawk) at 2,050 m a.s.l. (41°590 22″N, 128°000 11.37″

smaller than in higher altitudes, and therefore data in North1 were

E) in 2015. We compared the temperature recordings between the

not used to analyse the altitudinal changes in B. ermanii’s age. In

two sites and established a regression model (p < .001, R2 = 0.96)

total, 163 cores were available for analysing the altitudinal changes

(Appendix S2, Fig. S3) that was used to construct the temperature

using tree age after cross-dating (Table 1).

data near the current treeline.

T A B L E 1 Characteristics of the four Betula ermanii transects on the Changbai Mountains No. of cores (trees) Transects

Altitude range (m a.s.l.)

No. of plots

a

1892–2115

9

0 (0)

40 (26)

6.4–38.5



4.2–18.8

40–207

North2b

1940–2025

8

71 (37)

16 (11)

13.7–40.4



11.4–17.9

8–188

North1

a

b

DBH (cm)

West1

1991–2062

8

54 (42)

15 (9)

7.2–26.4

West2c

1994–2060

4

38 (19)

16 (9)

8.8–32.6

DB (cm)

2.9–6.3 —

Height (m)

Age (year)

1.4–18.8

10–182

5.2–17.5

26–178

North1, north side transect 1; North2, north side transect 2; West1, west side transect 1; West2, west side transect 2; DBH, diameter at breast height; DB, diameter at base. DB was measured for trees with a height of 200 cm. The DB of 90%

(p < .01) on the Changbai Mountains. This finding further indicates

of trees was