Climate response of tree-ring width in Larix sibirica growing in the drought-stressed forest-steppe ecotone of northern Mongolia

Abstract

• Introduction

Siberian larch (Larix sibirica) has its southern distribution limit in Mongolia in an area of rapidly rising temperatures. Direct effects of climate on tree-ring formation due to drought stress or indirect effects via the control of insect herbivore populations are little studied.

• Methods

The hypotheses were tested that stem increment of Siberian larch is reduced by (1) drought during the growing season and (2) high snow cover, as the latter is thought to protect hibernating herbivores, including gypsy moth (Lymantria dispar).

• Results

Tree-ring width increases with decreasing summer temperature, increasing precipitation during the growing season and decreasing winter precipitation.

• Conclusions

The susceptibility of stem wood formation to drought during the growing season suggests that a future climate warming will decrease productivity of Siberian larch, thus affecting its existence within the forest-steppe ecotone of Mongolia. Narrow tree rings in years following winters with low snowfall support the hypothesis that winter precipitation exerts an indirect effect on the growth of Siberian larch by controlling the survival rates of gypsy moth eggs.

1 Introduction

The southern distribution limit of the vast Siberian taiga in the northern Mongolian forest-steppe ecotone is of special interest for the study of the impact of global warming on vegetation for four reasons: (1) late twentieth century warming in central and semiarid northeastern Asia significantly exceeded the global average (Batima et al. 2005); (2) northern Mongolia is a marginal habitat for forests already under the current climate, as tree growth is principally limited by drought (Gunin et al. 1999; Dulamsuren et al. 2009a); (3) the position of the forest-steppe borderline already migrated through Mongolia along with changing climate throughout the Holocene (Dinesman et al. 1989); and (4) the livelihood of Mongolia’s rural population strongly depends on natural resources and is, in turn, influenced by climate-driven changes in vegetation patterns or productivity (Havstad et al. 2008). Even the main tree species in Mongolia’s forest-steppe ecotone, Siberian larch (Larix sibirica Ledeb.), covering 80% of the country’s forested area (Tsogtbaatar 2004), regularly suffers from drought stress (Dulamsuren et al. 2009a), shows reduced stem growth (Dulamsuren et al. 2010a), and has difficulties to regenerate in that area (Dulamsuren et al. 2008, 2010b; Hauck et al. 2008).

To assess the potential impact of global warming on the performance of Siberian larch forests in Mongolia, it is necessary to know which climate parameters exert an influence on tree growth. We, therefore, analyzed the impact of the seasonal and interannual variation of temperature and precipitation on tree-ring width in Siberian larch in the western Khentey Mountains, an area where weather recording started in the early 1960s. Since summer drought was shown to result in drought stress in Siberian larch of the western Khentey Mountains (Dulamsuren et al. 2009a) and larch trees show reduced stem growth since the 1950s (Dulamsuren et al. 2010a), we tested the hypothesis that high summer temperatures and low summer precipitation are correlated with narrow tree rings.

In addition to climate parameters which may exert a direct physiological effect on water relations, photosynthesis, and wood formation, indirect effects of climate on tree growth should be considered. Such indirect effects include the impact of herbivores or pathogens. Gypsy moth (Lymantria dispar L.) is a major pest in the forest-steppe ecotone of Mongolia, which has the potential to reduce tree-ring widths in Siberian larch (Ghent and Onken 2004; Hauck et al. 2008; Dulamsuren et al. 2010c). Since the eggs of gypsy moth are susceptible to severe frost (Waggoner 1985), we also tested the hypothesis that winter temperatures are inversely correlated with tree-ring width.

Increasing attention has recently been paid to the fact that climate–growth relationships can differ between the individual trees within one and the same stand (Wilmking et al. 2004). Several reasons may cause such divergent behavior, including genetic differences within the tree population or the small-scale variation of site parameters. Carrer and Urbinati (2004) additionally showed the existence of age-dependent variation in the response of tree-ring width on climate. We therefore distinguished three different age classes and between trees growing directly at the forest edge to the steppe or at the forest interior to test the hypotheses that the position along the forest-steppe ecotone and tree age influence the correlation of tree-ring width with climate parameters.

2 Materials and methods

2.1 Study area

Field work was carried out in the western Khentey Mountains, northern Mongolia, Mandal county (Mandal somon) near Khonin Nuga Research Station (49°04′48″ N, 107°17′15″ E), located 130 km north of Ulan Bator in the valley of the River Eroo, a tributary of the Orkhon River. Elevation of the study area ranges from 900 m in the river valleys up to 1,600 m asl on the mountain tops. Geologically, the Khentey mainly consists of Proterozoic and Paleozoic rocks, especially of granite. Permafrost is found under forests on north-facing slopes. Study sites were located near Khonin Nuga Research Station on northwest-facing slopes of the Sharlan Valley (49°6′ N, 107°19′ E; in the following called site A) at an elevation of 1,010–1,072 m, Mt. Bayantogol (49°5′ N, 107°17′ E; site B) at 1,020–1,120 m, and Mt. Baziin Am (49°3′ N, 107°15′ E; site C) at 990–1,090 m.

The vegetation of the study area was currently surveyed by Dulamsuren et al. (2005a). Light taiga forests with L. sibirica dominate in the lower montane belt and partly in floodplains. Mountain forests of the lower montane belt are limited to the better water-supplied north-facing slopes, whereas south-facing slopes are covered with different types of montane meadow steppe, mountain steppe, and small savanna-like Ulmus pumila stands (Dulamsuren et al. 2005b, 2009b).

2.2 Climate data

The climate of the Khentey Mountains is characterized by the Asiatic anticyclone in winter, which typically has its center southwest of Lake Baikal and causes dry and cold winters. In summer, warm air masses from the south flow into northern Mongolia resulting in the formation of cyclones when they meet the cold air from Siberia (Dulamsuren and Hauck 2008). The nearest permanent weather station to the study area with a continuous dataset is located 87 km NW of Khonin Nuga Research Station (Station “Eroo”) at an elevation of c. 900 m (49°48′ N, 106°42′ E). From this weather station, data of air temperature and precipitation are available since 1961 (Fig. 1) and were used as explanatory variables. Mean temperatures from 1961 to 2005 amounted to −2.1 ± 0.2°C for the whole year, −27.4 ± 0.5°C for January, and 18.7 ± 0.2°C for July. Variation of mean temperature is higher in winter (mean January temperatures, −35.3 to −17.0°C) than in summer (July, 16.4 to 22.3°C). Mean annual precipitation was 277 ± 11 mm. A marked precipitation peak in summer (Fig. 1) with 78 ± 6 mm in July (but only 3 ± 1 mm in January) is typical of the area (Dulamsuren and Hauck 2008).

Fig. 1

Climatic diagram for Eroo (49°48′ N, 106°42′ E, 900 m asl). Mean monthly precipitation and temperatures were calculated for the period 1961 to 2005

2.3 Wood core sampling and processing

Wood cores were sampled from the trees of the three sites A, B, and C in September 2007. Six plots of 20 m × 20 m were selected per site. They were located in minimum distances of 30 m from the forest line to the steppe and of 100 m from one another. In these plots, wood cores of 272 trees of L. sibirica were collected (Table 1) with an increment borer of 5 mm of the inner diameter. In addition to the trees growing inside the forest, single trees located immediately at the forest line were sampled each in one oblong plot per site. The length of these plots varied between 50 and 100 m. The number of these additionally sampled trees amounted to 93 trees (Table 1). At 1 m above the ground, the borer was driven into the wood parallel to the contour lines of the mountain slopes to avoid compression wood. Annual tree-ring width was measured with a precision of 10 μm on a movable object table (Lintab 5, Rinntech, Heidelberg, Germany), the movements of which are electronically transmitted to a computer system equipped with TSAP (Time Series Analysis and Presentation)-Win software (Rinntech). The light earlywood and the dark latewood were visually differentiated and separately measured.

Table 1 Descriptive statistics of the sampled L. sibirica trees, including the number of sample trees (N), mean tree-ring width, values of the expressed population signal (EPS), mean sensitivity and first-order auto correlation (arithmetic means ± standard deviation)

2.4 Evaluation of the tree-ring data

Evaluation of tree-ring data was, like data recording, conducted with TSAP-Win software. Tree-ring series were controlled for missing rings and false rings during crossdating, as especially the failure of tree rings is not rare in the semiarid environment of Mongolia at the drought limit of forests. Crossdating was based on the use of two parameters: coefficient of agreement (“Gleichläufigkeit” [GL]; Eckstein and Bauch 1969) and (standard) t values. The GL and t values measure the similarity between tree-ring series in the high- and low-frequency domain, respectively. Trees of the same site were pooled by calculating mean values of annual increment. Tree-ring series used for the calculation of means had GL values >65% (P ≤ 0.05) and t values >3 (>6 in c. 90% of the samples). Trend lines were calculated using moving 5-year averages, since we aimed at keeping as much of the interannual (high-frequency) variation as possible, because this is crucial for any climate–response analysis.

Standardized tree-ring series were used in correlation analyses with climate parameters. Standardization is principally achieved by dividing the observed tree-ring width (w _i ) through the expected annual increment. Finding the correct function for the age-related growth trend (for estimating the expected growth rate) is easier in semiarid environments, including our study area, than in regions with ample water supply, as the stand density in semiarid woodlands is relatively low and, thus, the individual trees are less influenced by changes in the stand structure than in dense and moist forests (Cook 1985). Therefore, the same type of standardization could be applied for all tree-ring series. The annual tree-ring index (z _i ) of year i was calculated with the equation ${{\hbox{z}}_{\rm{i}}} = 100 \times \left( {{{\hbox{w}}_{\rm{i}}} - { }{{\hbox{m}}_{\rm{i}}}} \right)/\sum {{{\hbox{w}}_{\rm{i}}}}$ with m _i being the 5-year moving average of year i. Three age classes were distinguished in the analyses including “old trees” with the oldest tree ring (at the sampling height of 1 m) older than 90 years, “middle-aged trees” (oldest tree ring being between 50 and 90 years old), and “young trees” (oldest tree ring younger than 50 years old).

Tree-ring series were representative of the studied stand as shown by the calculation of the expressed population signal (EPS) using sums of squares of within-core variation and error sums of squares calculated in a two-way analysis of variance (ANOVA; Wigley et al. 1984). The EPS clearly exceeded the 0.85 threshold (Wigley et al. 1984), except for the few middle-aged trees available at site C (Table 1). Therefore, in contrast to all other tree-ring series, the reliability of these particular data from site C is limited.

Year-to-year variability of tree-ring width was expressed as mean sensitivity, as was the influence of the increment in the previous year on the increment of the current year as first-order autocorrelation coefficient (Fritts 1976). Mean sensitivity and autocorrelation coefficients were calculated over the whole lifespan of each sample tree. Afterwards, mean values of these parameters were calculated over age classes and sites. These means were tested for statistical significance using Duncan’s multiple range test, after testing for normal distribution with the Shapiro–Wilk test. These analyses and the abovementioned ANOVAs were calculated with SAS 6.04 software (SAS Institute Inc., Cary, NC, USA).

2.5 Analysis of climate–growth relationships

Relationships between climate parameters and tree-ring width were analyzed both with a multivariate approach using multiple regression analysis and a bivariate approach calculating Pearson’s correlations coefficients. Bivariate correlations were calculated in addition to multiple regression analysis to show the positive or negative character of the influence of the individual climate parameters on tree-ring width. As the repeated analysis of the same dataset by multiple and bivariate regression analyses carries the inherent danger of increasing type I errors, results of the bivariate regression analyses are interpreted with caution if the relevant determinants do not occur in the multivariate models. Independent climate parameters in the multivariate analysis included quarterly and annual mean temperatures and mean precipitation of the current and the previous year. The bivariate analysis was restricted to quarterly and monthly climate values; months without any significant correlation are not displayed in the relevant figures. Tree-ring parameters (dependent variables) included total tree-ring width as well as the widths of earlywood and latewood. We refrained from dividing the data in two intervals for calibration and verification (Fritts et al. 1971; Fritts 1976), because the available series of climate data covers only 45 years. The renouncement of a verification period and the choice of relatively simple correlation statistics were also preferred by Bräuning (1999), who had to deal with similarly short time series of climate data in Tibet. Trees which formed the oldest tree ring after 1961 were omitted from the analysis. Regression analyses were computed with SAS 6.04 software.

3 Results

3.1 Descriptive statistics

The tree-ring chronology from the western Khentey Mountains ranges over little more than 200 years (Fig. 2). In many cases, mean sensitivity is significantly higher in trees from the forest edge than in trees from the forest interior (Table 1). There is also a trend for lower first-order autocorrelation at the edge than in the interior of the forest, which is rarely significant, though. Mean sensitivity decreased with increasing autocorrelation in old interior trees from sites A (r = –0.45, P = 0.02) and C (r = –0.73, P = 0.004) and in middle-aged interior trees from sites A (r = –0.33, P = 0.05) and B (r = –0.44, P < 0.001). At the forest edge or in young trees, such relationship was lacking.

Fig. 2

Tree-ring chronology of Siberian larch calculated from all old trees (cambial age >90 years) growing in the forest interior in the western Khentey Mountains (sites A, B, and C)

3.2 Climate response

High summer temperatures and low precipitation during the growing season are correlated with the occurrence of narrow tree rings. Analyzing the effect of climate on tree-ring width separately for trees of different age groups, sites, distances to the forest edge and for width of the total tree ring, earlywood and latewood, summer temperature occurs in 36 (71%) out of the 51 most significant regression models with three climate parameters and in 47 (92%) regression models with four climate parameters (Table S1, available online at www.afs-journal.org). In a monthwise analysis, significant correlations with summer temperature are found with June, July, and/or August temperature depending on the plot (Fig. 3, S1, S2), whereas significant correlations are lacking in the hottest, but simultaneously most rain-laden month, July (Fig. 1). Correlations with summer temperature occur more regularly if calculations are based on quarterly values (from June to August; Fig. S3 to S5). Precipitation during the growing season is correlated with tree-ring width in 45 (88%) of the most significant regression models with three independent variables and in 48 (94%) of the most significant models with four independent variables (Table S1). Mostly, these correlations refer to the precipitation at the end of the growing season in the year prior to tree-ring formation in September and October (Fig. 4, S6, S7) at a point in time when monthly precipitation is limited to 30 mm (September) or <10 mm (October; Fig. 1). Summer precipitation (June to August) is only included in 10% of the regression models with three variables and in 14% of the models with four independent variables (Table S1). These correlations refer both to June and July and to the year of and prior to tree-ring formation (Fig. 4, S6, S7; Table S1). Correlation of tree-ring width with spring precipitation is negligible (Table S1).

Fig. 3

Response of tree-ring width to monthly mean temperatures of the year of and the year prior to tree-ring formation at site A (Pearson’s product-moment correlation coefficients; *P ≤ 0.05). a, c, e Forest interior. b, d, f Forest edge. a, b Total tree-ring width. c, d Earlywood. e, f Latewood. Months without any significant correlation are omitted

Fig. 4

Response of tree-ring width to monthly precipitation of the year of and the year prior to tree-ring formation at site A (Pearson’s product-moment correlation coefficients; *P ≤ 0.05). a, c, e Forest interior. b, d, f Forest edge. a, b Total tree-ring width. c, d Earlywood. e, f Latewood. Months without any significant correlation are omitted

Low snowfall in the winter before the growing season when the tree-ring is formed, or more rarely even in the winter before that winter, is correlated with wide tree rings (Fig. 4, S6, S7). Nearly all most significant models with three (49; 96%) or four (47; 92%) independent climate parameters include winter precipitation (Table S1), which never exceeds 5 mm per month (Fig. 1). Winter temperature is not relevant to tree-ring formation (Fig. 3, S1, S2; Table S1).

No consistent pattern for an influence of tree age, the site or the plot position in the forest-steppe ecotone (forest edge vs. interior) on the correlations between tree-ring width and climate parameters is observed.

4 Discussion

The results support the hypothesis that drought periods during the growing season limit stem wood formation in Siberian larch. The negative influence of low precipitation and high summer temperatures match with low shoot water potentials, close to the point of zero turgor, measured during drought periods in trees of Siberian larch in the study area (Dulamsuren et al. 2009a). Furthermore, the finding agrees with the fact that Siberian larch has its drought limit in Mongolia (Gunin et al. 1999; Miehe et al. 2007). It is noteworthy that the minimum shoot water potentials in the study of Dulamsuren et al. (2009a) decreased throughout the growing season with a minimum in early fall. Similarly, Li et al. (2007b) found steadily decreasing δ¹³C values in needles of Siberian larch from June to late September in the Kherlen valley, northeastern Mongolia. The lack of recovery of the plant water status in fall despite significantly lower temperatures than in summer is attributable to the low fall precipitation (Fig. 1). Li et al. (2007a) showed by the analysis of δ¹⁸O and δD (deuterium) signatures in precipitation, soil, and shoot water that forests of Siberian larch in northern Mongolia primarily depend on surface water from the current precipitation. Carbon assimilated at the end of the growing season is apparently used for wood formation in the following growing season, explaining the increase of tree-ring width with increasing precipitation in September and even October of the year prior to tree-ring formation. Siberian larch is photosynthetically active in the study area until late September or early October (Dugarjav 2006; Li et al. 2007b; Dulamsuren et al. 2009a), though there is considerable interannual variation due to the highly variable climate in northern Mongolia (Nandintsetseg et al. 2007; Dulamsuren and Hauck 2008).

Decreasing tree-ring width with increasing winter precipitation supports the hypothesis that the snow cover protects herbivores of Siberian larch hibernating on the ground from frost damage and thereby reduces the annual wood formation in the subsequent growing seasons. Gypsy moth is the most common and most detrimental herbivore of Siberian larch in northern Mongolia (Hauck et al. 2008). Heavily infested larch trees were shown to show strongly reduced tree-ring widths for at least 2 years after the gypsy moth invasion (Dulamsuren et al. 2010c). Gypsy moth preferentially oviposits on rock outcrops on south-facing steppe slopes near forest edges (Hauck et al. 2008). Eggs of gypsy moth were shown to die at steady temperatures below −26°C (Waggoner 1985) and therefore benefit from insulation by snow (Nealis et al. 1999; Andresen et al. 2001). Freezing rapidly kills the eggs below the −26°C threshold, as all gypsy moth eggs used in a laboratory experiment died within 46 h at −26°C, but already within 20 min at −27.8°C (Waggoner 1985). Periods with very low winter temperatures are characteristic of the Mongolian forest-steppe ecotone. Dulamsuren and Hauck (2008) measured minimum air temperatures (at 1.5 m above the ground) of −43°C on south-facing slopes of the western Khentey Mountains in immediate proximity to egg masses of gypsy moth. Periods with daily minima below −26°C occur for weeks, whereas intervals where even the daily maximum temperature does not exceed this value are limited to single or a few consecutive days. Mongolia’s low winter temperatures and the extremely low snowfall (total precipitation from October to April in the study area is <30 mm; Fig. 1) rule out that the high significance of winter precipitation for tree-ring width in the regression models is directly connected to carbon assimilation or water relations in L. sibirica. Snow cover plays no role at reducing frost damage in near-surface roots (Korotaev 1994), as otherwise tree-ring width should increase, not decrease, with increasing winter precipitation. Limitation of root activity in spring by persistent snow packs as postulated for Larix gmelinii stands in central Siberia (Kujansuu et al. 2007) is not to be expected in Mongolia given the low amount of winter precipitation.

Trends for higher mean sensitivity and lower autocorrelation in tree-ring series from the forest edge than the forest interior show that trees at the edge to the steppe are exposed to more strongly varying climatic conditions than interior trees. The regression analyses, however, suggest that wood formation in Siberian larch is influenced by drought during the growing season and snow cover in winter largely irrespective of their distance from the forest line to the steppe, their age or the study site. Both direct (drought years) and indirect (climate-influenced herbivore invasions) climate effects can exert an enduring effect on carbon allocation in larch trees. This is suggested by the numerous correlations of tree-ring width found with climate parameters from the year prior tree-ring formation (Kujansuu et al. 2007).

5 Conclusion

Our data suggest that both drought periods during the growing season and winter snow pack influence wood formation in Siberian larch in the northern Mongolian forest-steppe ecotone across a wide range of habitats differing in microclimate (as assessed by the comparative study of trees from the forest interior and the edge to the steppe) and age classes. Aridity has steadily increased during the late twentieth and the early twenty-first centuries in the western Khentey Mountains (Dulamsuren et al. 2010a). The results of the regression analyses confirm that growth declines observed in Siberian larch in the western Khentey since the mid-twentieth century are attributable to increasing drought during the growing season as already established by Dulamsuren et al. (2010a). The significance of gypsy moth invasions will probably also increase due to climate warming, which may, among others, cause higher survival rates in winter (Logan et al. 2003). Thus, it is probable that global warming will increasingly deteriorate the growth conditions of Siberian larch in the northern Mongolian forest-steppe ecotone both by causing drought stress and by triggering more frequent gypsy moth outbreaks.