A Dual Stable Isotope Approach Unravels Common Climate Signals and Species-Specific Responses to Environmental Change Stored in Multi-Century Tree-Ring Series from the Tibetan Plateau

Tree-rings are recorders of environmental signals and are therefore often used to reconstruct past environmental conditions. In this paper, we present four annually resolved, multi-centennial tree-ring isotope series from the southeastern Tibetan plateau. The investigation site, where juniper and spruce trees jointly occur, is one of the highest known tree-stands in the world. Tree ring cellulose oxygen (δ18O) and carbon (δ13C) isotopes were analyzed for a common period of 1685–2007 AD to investigate climate–isotope relationships. Therefore, various climate parameters from a local meteorological station and from the CRU 4.02 dataset were used. Tree-ring δ18O of both species revealed highly significant sensitivities with a high degree of coherence to hydroclimate variables during the growing season. The obtained δ18O–climate relationships can even be retained using a species mean. In contrast, the individual δ13C series indicated a weaker and non-uniform response to the tested variables. Underlying species-specific responses and adaptations to the long-term trend in atmospheric CO2 bias even after a trend correction identified dominant environmental factors triggering the tree-ring δ13C at our site. However, analysis of individual intrinsic water-use efficiency in juniper and spruce trees indicated a species-specific adaptation strategy to climate change.


Introduction
Often-quoted in a broader context as the "Third Pole" [1], an important part of the "Asian Water Towers" [2], or the "Roof of the World" [3], the Tibetan plateau (TP) exerts a key role for the climatic, hydrological and biological conditions of its highlands, mountain ranges as well as for the adjoining Asian lowland areas.Beyond the fact of its substantial role within the global climate system as the key driver of the Asian monsoonal systems [4], the current changes within the regional cryosphere, hydrosphere, biosphere and associated responses of environmental systems on the TP are subjects of intensive research.Especially, the responses, threats and adaptations of different (socio-)ecosystems to the substantial rise in temperature during the recent decades and associated changes of water availability underline the necessity of a precise quantification of these changes [5][6][7][8][9].Contrastingly, for High Asia in general and for the Tibetan plateau in particular, only a sparse climate station network with relatively short time-series of high-resolution meteorological data is available.Therefore, existing datasets, which are mainly available at the earliest since the 1950s, hamper a substantial quantification and classification of the recent climate change within a retrospective context and underline the needs for paleoclimatic datasets [4,[10][11][12].
Besides studies focusing on tree-ring width (TRW), more and more studies display the advantages of using stable isotopes in tree-rings for ecological as well as paleoclimate research on the Tibetan plateau.During the last decade, the analysis and use of tree-ring δ 18 O has become one of the main proxies used for paleoclimate research in High Asia.Recent works on the Tibetan Plateau, in particular, have demonstrably contributed to a better understanding of regional climate and environmental history [10][11][12][13][14][15][16][17].In comparison, studies on the δ 13 C variations in tree-rings are still underrepresented due to more complex isotope fractionation effects affecting the δ 13 C in tree-rings.In addition, multi-century δ 13 C tree-ring time series usually have to undergo a trend correction procedure primarily due to the rising atmospheric CO 2 concentrations (CO 2 ATM ) since the beginning of the industrialization (the "Suess-effect").Especially the respective correction factor that has to be considered due to plant-internal reactions addressing rising CO 2 ATM has considerable impact on the adjustment of carbon isotope time series [18,19] and may strongly modify trends of δ 13 C tree-ring data after correction.
Recently, progress has been made in studies on analyzing both stable carbon and oxygen isotopes, pointing to the paleoclimatic potential inherent to such an approach [20][21][22][23].In addition, studies underline the potential and benefit within the analysis of multiple stable isotopes across species and across sites [24,25].The reported advantage is mainly based on the possibility of a common assessment and interpretation of stored climate signals with coupled ecophysiological reactions across tree species.However, up to now, the mentioned studies focus either on the analysis of two tree species, but only for over a rather short time period of a few decades, or the establishment of longer (mutli-century to millennial) time series for two isotopes, but only investigating a single tree species.The main reason for this desideratum is presumably caused by the comparably high laboratory and time expenditure inherent to tree-ring isotope studies, although many sites on the TP possibly bear the potential of additionally analyzing several tree species.For a profound retrospective assessment, classification and categorization of the recently ongoing environmental changes and tree-specific adaptations, longer isotope time-series across species are, however, urgently required.
When considering existing studies on tree-ring δ 13 C and δ 18 O from the TP, their spatial representativity and probable inter-isotope dependencies sometimes remain difficult.Due to the lack of other existing, closely situated tree-ring isotope series especially for carbon isotope series, it remains often challenging to corroborate a wider regional representativity.Moreover, the derived, detected and stored environmental signal in the tree-ring proxy is reported to reflect a more site-and species-related physiological response to the ongoing environmental change [26,27].This includes in particular observed adaptations in different tree species with regard to changes in their intrinsic water-use efficiency [25,28].The sometimes heterogenous picture with regard to the (reconstructed) environmental signal stored in δ 13 C and δ 18 O series from the TP, the different levels in signal strength and differing spatial representativities of the time series still contain challenges for actual and future research.
By analyzing multi-centennial variations of different stable isotopes (δ 13 C and δ 18 O) in two different tree species from a common site on the southeastern part of the TP, our study evaluated the potential, benefit, and limitations inherent in such an approach by: (i) deciphering and explaining common and contrasting long-term isotopic variations of δ 13 C and δ 18 O between two tree species from high-altitude, climate-sensitive tree-stands; (ii) assessing the added value for an analysis of δ 13 C and δ 18 O from two local tree species; and (iii) testing if common/differing signals of the proxy series from one site can be used to better assess, decipher, and explain diverging responses of tree species to recent environmental change.
The main aim of this study was therefore a calibration study on the potential of multiple isotope parameters derived from tree-rings from a subalpine tree-stand in High Asia.

Study Site
The study site is situated around 20 km southwest to the municipality of Qamdo on the upper reaches of the Mekong River on the southeastern Tibetan plateau (TP) (Figure 1).The  The main aim of this study was therefore a calibration study on the potential of multiple isotope parameters derived from tree-rings from a subalpine tree-stand in High Asia.

Study Site
The study site is situated around 20 km southwest to the municipality of Qamdo on the upper reaches of the Mekong River on the southeastern Tibetan plateau (TP) (Figure 1).The local climate as recorded by the nearby climate station in the municipality of Qamdo (31°09´ N/97°10´ E; available period 1954-2012) is characterized by a mean annual temperature of 8.8 °C and an annual precipitation sum of 482 mm, indicating a sub-humid continental climate (Figure 1).The main summer season lasts from June to August and represents the period with highest temperatures (TJJA = 16.7 °C).Summer months with the maximum monthly temperatures are July (TJuly = 17.4 °C) and August (TAug = 16.6 °C), respectively.The summer months July and August have the highest amounts of precipitation (PrecipJJA = 296 mm; PrecipJuly = 107 mm, PrecipAug = 101 mm).In contrast, winter months are cold and show little to no rainfall/snow (TDec = 0 °C, PrecipDec = 2mm).Overall, local climate conditions display a typical summer monsoon climate with warm and humid summers contrasted by cold and dry winter conditions [4].The investigated tree-stand is located at the local subalpine vegetation belt between 4400 m and 4600 m asl. on a convex slope with slightly changing southerly expositions (Figure 2).Developed from debris material, predominant soils are well-drained and shallow without influence of stagnating groundwater [12,29,30].The upper timberline on southerly exposed slopes is formed by open, steppeforest-like formations of Tibetan juniper (Juniperus tibetica Kom.), being one of the highest known tree-stands of the world at 4600 m asl.[31].On average, juniper tree individuals reach heights of up The investigated tree-stand is located at the local subalpine vegetation belt between 4400 m and 4600 m asl. on a convex slope with slightly changing southerly expositions (Figure 2).Developed from debris material, predominant soils are well-drained and shallow without influence of stagnating groundwater [12,29,30].The upper timberline on southerly exposed slopes is formed by open, steppe-forest-like formations of Tibetan juniper (Juniperus tibetica Kom.), being one of the highest known tree-stands of the world at 4600 m asl.[31].On average, juniper tree individuals reach heights of up to 12 m, with some individuals exceeding 20 m.As being typical for steppe forests, the juniper tree stand is characterized by a canopy cover of 10%.On eastern and southeastern expositions, denser stands of 20-25 m high Balfour spruce trees (Picea balfouriana) with a canopy cover of maximum 25% are dominant, reaching their maximum upper distribution limit at 4500 m asl.In contrast to the juniper stands, spruce trees are in parts covered with lichens of the species Alectoriella ssp., pointing to more humid microclimatic conditions presumably caused by lower radiation and higher humidity (cf. Figure 2).to 12 m, with some individuals exceeding 20 m.As being typical for steppe forests, the juniper tree stand is characterized by a canopy cover of 10%.On eastern and southeastern expositions, denser stands of 20-25 m high Balfour spruce trees (Picea balfouriana) with a canopy cover of maximum 25% are dominant, reaching their maximum upper distribution limit at 4500 m asl.In contrast to the juniper stands, spruce trees are in parts covered with lichens of the species Alectoriella ssp., pointing to more humid microclimatic conditions presumably caused by lower radiation and higher humidity (cf. Figure 2).

Sampling, Dating and Sample Preparation
During several field campaigns, undisturbed and healthy mature tree-individuals of Picea and Juniperus species were sampled.We collected discs from dead trees and increment cores from living trees using an increment borer (Suunto).In the laboratory, samples were surfaced and tree ring widths were measured with an accuracy of 0.01 mm using a Lintab 6 system (Rinntech, Heidelberg/Germany). Subsequently, samples were cross-dated and validated against an existing regional master chronology [29,30].For the analysis of the ratios of stable oxygen and carbon isotopes from tree-ring cellulose (thereafter expressed as δ 18 OTRC and δ 13 CTRC), at least 5-7 cores from both juniper and spruce individuals were selected.The obtained minimum sample depth of five trees was within the generally accepted threshold of analyzed trees per species per sample site for isotope studies [32] and was consistently available throughout the common period 1685-2007 AD covered by both species (Table 1).

Sampling, Dating and Sample Preparation
During several field campaigns, undisturbed and healthy mature tree-individuals of Picea and Juniperus species were sampled.We collected discs from dead trees and increment cores from living trees using an increment borer (Suunto).In the laboratory, samples were surfaced and tree ring widths were measured with an accuracy of 0.01 mm using a Lintab 6 system (Rinntech, Heidelberg/Germany). Subsequently, samples were cross-dated and validated against an existing regional master chronology [29,30].For the analysis of the ratios of stable oxygen and carbon isotopes from tree-ring cellulose (thereafter expressed as δ 18 O TRC and δ 13 C TRC ), at least 5-7 cores from both juniper and spruce individuals were selected.The obtained minimum sample depth of five trees was within the generally accepted threshold of analyzed trees per species per sample site for isotope studies [32] and was consistently available throughout the common period 1685-2007 AD covered by both species (Table 1).All samples chosen for isotope analysis met the following criteria: (i) similar tree ages for juniper and spruce trees; (ii) no missing rings or occurrence of problematic plateaus with extremely depressed growth and narrow rings to obtain sufficient α-cellulose for the analysis of both δ 13 C TRC and δ 18 O TRC ; and (iii) the analytical exclusion of the innermost 50 years of each sample to avoid a non-climatic juvenile trend of the most recent δ 13 C TRC values [33].Especially Criterion (ii) was difficult to meet within multi-century old junipers known for lobate growth and in parts extremely narrow ring widths [11,12].To obtain sufficient material for the analysis of both isotopes from such samples, it was therefore inevitable to pool material of corresponding calendar years prior to the subsequent laboratory steps, as reported and successfully applied within existing studies [11][12][13][14]17,32,33]. To avoid a mass-related over-representation of a single tree individual within our pooled samples, cut wood samples of every tree-ring were weighed individually prior to cellulose extraction.Subsequently, same amounts of each calendar year were combined for the final pooled sample year.In regard to the aim of this study of an inter tree-species comparison with the precondition of a comparable analytical setting and subsequent sound statistical analyses, sample pooling was also performed for the spruce samples.

Evaluation of δ 18 O TRC and δ 13 C TRC
In a first analytical step, dated growth rings of both species were cut under a binocular using a razor blade.Calendar year-wise pooled wholewood samples of juniper and spruce were processed to α-cellulose following the method in [34].After the extraction process, the obtained cellulose samples were homogenized using an ultrasonic sound system and finally freeze-dried in a lyophilization unit [35].Cellulose samples of each species of 200 µg were weighed using a micro balance and individually packed in tin (for carbon isotope analyses) and silver capsules (for oxygen isotope analyses), respectively.Measurements of the 12 C/ 13 C ratio (expressed as δ 13 C) and 16 O/ 18 O ratios (expressed as δ 18 O) of spruce and juniper samples were determined using a continuous flow DELTA V Advantage Isotope Ratio Mass Spectrometer (IRMS, Thermo Fisher Scientific Inc.) (Waltham, MA, USA).To avoid a hygroscopically caused exchange with ambient air, all samples for the assessment of oxygen analyses were stored prior to the IRMS analysis for at least 48 hours in a vacuum oven.
For the analysis of δ 18 O TRC isotopes, samples were pyrolyzed to CO in a HEKATech pyrolysis reactor coupled to the IRMS at 1080 • C. For the measurements of the δ 13 C variations in tree-ring cellulose, samples were combusted at 1020 • C under an excess of oxygen to CO 2 .All measurements were calibrated by periodically interposing internal and international laboratory standards resulting in an overall analytical precision for the δ 18 O TRC measurements of <0.2‰ and for the δ 13 C TRC measurements of <0.1‰, respectively.
As a result, we were able to set up for both investigated stable isotopes juniper chronologies with periods reaching from 489 to 2011 AD.At this point, it has to be noted that the Juniperus δ 18 O TRC time series was already used for a climate reconstruction [12].However, by using new material, we were able to update the series until 2011 AD.Thereto, pooled samples of the new individuals were analyzed with an overlap of 20 years to the existing series.Due to their lower maximal biological age, the obtained spruce isotope chronologies cover the period from 1685 to 2007 AD.Hence, the common and further considered period covered by both species lasts from 1685 to 2007 AD.

Examination of Isotope-Climate Relationships
As being a crucial step for the calibration of tree-ring proxy data in remote high elevation areas, we conducted a thorough survey of available and representative climate datasets.For a first examination of climate-proxy relationships, instrumental data from the nearby climate station in Qamdo (available since 1954) were used.The dataset includes a wide variety of climate parameters including precipitation, temperature (T min , T max , and T mean ), wind speed (ws), vapor pressure deficit (vp), sunshine hours (sh), and relative humidity (rH).Especially in high mountain areas such as the TP with steep topoclimatic gradients within short distances, precipitation data from climate stations located in valley bottoms often show a limited representativity for study sites located at the upper tree line.Therefore we additionally used gridded means of climate variables from the four nearest 0.5 • grid points of the CRU TS 4.02 dataset (Figure 1; https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4. 02/cruts.1811131722.v4.02/) [36].Therein, provided climate parameters are precipitation, temperature (T max , T mean , and T min ), diurnal temperature range (dtr), wet day frequencies (wet), cloud cover (cld), frost day frequency (frs), and potential evapotranspiration (pet).When using gridded datasets from this area, however, it must always to be taken into account that the early periods are interpolated from a rather low number of available climate stations.Therefore, we truncated the first five decades and only used the CRU data for the common period with stationary data for 1954-2007 (n = 53) where highly significant relationships between both datasets appear (p < 0.01).
For the quantification of the statistical relationships, we computed Pearson's correlation coefficients (r) for all available variables in both datasets.The resultant correlations are visualized in a correlogram [37] where correlation coefficients are colored according to their positive/negative value including the respective level of significance.To determine the influence of the applied CO 2 trend correction on the correlations with climate, time series of all isotopes and tree species were filtered using a multiple local regression filter with windows set to 1/10 of the time series length.Wavelet analysis was employed to explore the existence and significance of high-and low-frequency dynamics and the threshold for the high-pass filter (R package WaveletComp).Running correlations for all isotopes and tree species were computed using a 50-year moving window.

δ 13 C TRC Detrending and iWUE Calculation
As reported by many other studies, δ 13 C TRC time-series are often characterized by a significant decreasing trend.This decline is mainly caused by the steadily increasing combustion of fossil fuels since the beginning of the industrialization after ca.1850 AD [18,33].As a result, the absolute amount of the atmospheric CO 2 concentration has risen steadily since 1800 from 270 ppm up to presently 410 ppm ( [38], https://www.esrl.noaa.gov/gmd/ccgg/trends/data.html).Through primarily released lighter 12 C isotopes, fossil fuel combustion has led to a coupled shift in the atmospheric 12 CO 2 / 13 CO 2 relationship resulting in a depletion of the atmospheric δ 13 CO 2 values.As a result, plants have and do still respond with plant physiological responses to the increased atmospheric CO 2 which in turn are biasing the environmental signals stored in tree-ring cellulose [19,32].To minimize this anthropogenic bias and to extract the environmental signal incorporated within the δ 13 C TRC , a trend-correction of the δ 13 C TRC raw values has to be performed prior to the assessment of climate-δ 13 C TRC -relationships.Within our study, we used the rather conservative non-linear δ 13 C-trend correction for CO 2 ATM , as introduced in [18].
To initially test a possible individual response of our two studied tree species towards changing CO 2 concentrations, we calculated for the period 1900−2007 their respective intrinsic water-use efficiencies (iWUE).As reported in [39], the iWUE is defined as the ratio of net photosynthetic assimilation rate (A) to stomatal conductance (g s ) for water vapor.It is calculated as: where C i and C a are the CO 2 concentrations in the intercellular air space of the leaves and in the atmosphere, respectively.C i can be obtained by calculating the carbon isotope discrimination (∆) for C3 plants, as introduced in [40].The ∆ is expressed therein as where δ 13 C a is the stable isotope value of the ambient air (in ‰) and δ 13 C p the stable isotope value of plant cellulose (in ‰).Values for δ 13 C a (including C a ) for the period 1900−2007 were derived from the datasets in [18,28]; the respective δ 13 C p denote the measured δ 13 C TRC values from our study.
Since ∆ is also a function of the difference between the intercellular and the atmospheric CO 2 , it can be additionally calculated as: where a reflects the isotopic discrimination during diffusion of CO 2 from the atmosphere into the intercellular space (~4.4‰.), and b denotes the isotopic discrimination caused by the discrimination against CO 2 in the chloroplasts by the enzyme RuBisCO (~27‰.).2).The mean values of the oxygen isotope series differ slightly within a range of about 0.6‰ (δ 18 O TRC Juniper = 23.43‰,δ 18 O TRC Spruce = 22.98‰).Since both tree-species are situated within a short elevational range, this is assumed to be caused by slightly different plant-internal fractionations during cellulose synthesis and/or-on a microclimatic scale-a slightly differing moisture availability during the vegetation period.Probably caused by a higher proportion of earlywood [11,12], the juniper chronology is characterized by a higher autocorrelation (AC1 = 0.61) than the spruce time-series (AC1 = 0.26).It is apparent that the inter-annual variability of the individual time-series is considerably high (Figure 3).In particular, during the recent decades, the number of positive and negative excursions in both time series exceeding a deviation of more than 1.5 SD increased compared to their long-term mean.The highly significant (p < 0.01) correlation coefficient of r = 0.60 and Gleichläufigkeit (sign test) of 70% between the isotope chronologies underline their strong common signal (Table 2).The similarity between the series is also shown by highly consistent variations on lower frequencies (11-year and 21-year Fast Fourier Transform (FFT) filter).While regarding the respective low frequency domains of the two series, the statistical relationships are even higher and stable within their respective levels of significance with p < 0.01 (11-year FFT filtered: r = 0.73, GLK = 91%; 21-year FFT filtered: r = 0.75, GLK 78%).This leads to the assumption that both species are probably controlled by the same environmental influences and therefore show a very similar temporal variability and environmental sensitivity.In contrast to the δ 18 O TRC -variations, the courses of the δ 13 C TRC time series highly differ between species.In general, both time series show an offset of around 2‰ in their individual mean values (mean δ 13 C TRC juniper = −18.83‰,mean δ 13 C TRC Spruce = −20.92‰).Mainly caused by the steadily rising influence of 13 C-depleted atmospheric CO 2 , both series are highly influenced by a markedly negative trend since the end of the 19th century.After 1900 AD, this anthropogenic trend becomes more apparent for the juniper trees, resulting in an absolute minimum of the δ 13 C TRC juniper series in 2010 at −21.59‰ (Figure 4).For the spruce trees, a corresponding trend can be observed, resulting in the lowest overall value of −24.21‰ in 2005.As mentioned, this non-climatic trend has to be corrected before evaluating the respective climate-proxy relationship.Figure 4 displays the results of the performed δ 13 C trend correction on the two isotope chronologies.By regarding the common period (1685-2007) of both CO 2 -corrected time series, it is noticeable that both tree species exhibit an overall negative trend (Figure 4), which is more prominent for spruce.During the pre-industrial period, long-term linear trends in both species display a similar, only slightly negative trend.After 1850, these trends are further enhanced with different timing and magnitude between the species.Compared to the juniper trees, the linear long-term trend of spruce δ 13 C TRC indicate a more pronounced decrease after 1850 AD.While regarding the 11-year FFT filter of both species (Figure 4), this can be explained by an apparent positive trend in the juniper δ 13 C TRC until the early 20th century.Thereafter, the junipers also show a comparable negative trend to the absolute minimum of the time series in 2007.Presumably, these differing trends point to species-specific physiological strategies against the rising atmospheric CO 2 content.Compared to the δ 18 OTRC time-series, the correlation between the species-specific corrected δ 13 CTRC chronologies is lower (r = 0.44, GLK = 52%; Table 2).Although the statistical interrelation increases in the low frequency domains (11-year FFT filtered: r = 0.64, GLK = 64%; 21-year FFT filtered: r = 0.69, GLK = 65%), this is assumable mainly caused by-despite the applied CO2-correction-the still apparent negative trend since 1800 AD.The high influence and species-specific bias of this trend is additionally displayed in Table 2 by calculations with the uncorrected δ 13 CTRC raw values.

Climate-δ 18 OTRC Relationships
The climate-isotope relationships for the two investigated tree species based on calculations with the climate station at Qamdo are displayed in the Figure 5A,B.According to the high similarity between the δ 18 OTRC-series of both species, it was expected that the resulting climate correlations are highly synchronous.To check the potential benefit and additional information inherent in a combination of tree species, we tested the applicability of a mean isotope series.To evaluate if the climate signal is dominant over species-specific plant-physiological growth characteristics, juniper and spruce δ 18 OTRC chronologies were averaged.Additionally, the first principal components (PC) of the two tree species for both isotopes were computed, with prior centering and scaling applied, as δ 13 CTRC values differ significantly in their mean.Since for two centered and standardized variables, the correlation matrix always has the same eigenvectors, the correlation between its first PC and the mean is 1.Therefore, we used the calculated species δ 18 OTRC mean and performed the same correlation analyses as for the individual species (Figure 5C).
It is apparent that hydroclimatic parameters including precipitation, vapor pressure deficit and relative humidity (rH) have the highest negative influence on the δ 18 O variations in spruce and juniper tree-rings.Values of juniper δ 18 OTRC reveal highly significant correlations with different monthly and seasonal means during the vegetation period (lasting from May to September) of different hydroclimatic variables (Figure 5A).Within those, correlations with relative humidity Compared to the δ 18 O TRC time-series, the correlation between the species-specific corrected δ 13 C TRC chronologies is lower (r = 0.44, GLK = 52%; Table 2).Although the statistical interrelation increases in the low frequency domains (11-year FFT filtered: r = 0.64, GLK = 64%; 21-year FFT filtered: r = 0.69, GLK = 65%), this is assumable mainly caused by-despite the applied CO 2 -correction-the still apparent negative trend since 1800 AD.The high influence and species-specific bias of this trend is additionally displayed in Table 2 by calculations with the uncorrected δ 13 C TRC raw values.

Climate-δ 18 O TRC Relationships
The climate-isotope relationships for the two investigated tree species based on calculations with the climate station at Qamdo are displayed in the Figure 5A,B.According to the high similarity between the δ 18 O TRC -series of both species, it was expected that the resulting climate correlations are highly synchronous.To check the potential benefit and additional information inherent in a combination of tree species, we tested the applicability of a mean isotope series.To evaluate if the climate signal is dominant over species-specific plant-physiological growth characteristics, juniper and spruce δ 18 O TRC chronologies were averaged.Additionally, the first principal components (PC) of the two tree species for both isotopes were computed, with prior centering and scaling applied, as δ 13 C TRC values differ significantly in their mean.Since for two centered and standardized variables, the correlation matrix always has the same eigenvectors, the correlation between its first PC and the mean is 1.Therefore, we used the calculated species δ 18 O TRC mean and performed the same correlation analyses as for the individual species (Figure 5C).
It is apparent that hydroclimatic parameters including precipitation, vapor pressure deficit and relative humidity (rH) have the highest negative influence on the δ 18 O variations in spruce and juniper tree-rings.Values of juniper δ 18 O TRC reveal highly significant correlations with different monthly and seasonal means during the vegetation period (lasting from May to September) of different hydroclimatic variables (Figure 5A).Within those, correlations with relative humidity averaged over the complete vegetation period have the strongest (negative) impact on δ 18 O TRC (r MJJAS = −0.67,p < 0.01.Considering the high autocorrelation of the juniper δ 18 O TRC series, the implementation of the previous year's vegetation period even slightly increases correlations with rH (r MJJAS + pyMJJAS = −0.68,p < 0.01).Precipitation (r Precip MJJAS = −0.46,p> 0.01), nd vapor pressure (r MJJAS = −0.59,p < 0.01) are also significantly correlated to δ 18 O TRC during the summer season.Positive significant correlations occur with maximum temperature during May to September (r MJJAS = 0.45, p < 0.01).Results for single months display similar sensitivities, although showing continuously lower correlation strength.Nearly identical climate-proxy relationships are apparent for spruce and the Qamdo climate station (Figure 5B).The corresponding correlogram reveals nearly the same pattern as for juniper, although the individual absolute values of the correlations are mostly higher.Highly significant relationships occur with relative humidity and vapor pressure deficit during the summer months (rH: r MJJA = −0.74,p < 0.01, vp: r MJJA = −0.62,p < 0.01), precipitation during the summer monsoon season (r JJA = −0.62,p < 0.01), and maximum temperatures in August (r Aug = 0.52, p < 0.01).It is apparent that the hydroclimatic conditions during the vegetation period are decisive for the δ 18 O TRC variations in the tree-ring.
The combined species' mean chronology showed largely identical results as displayed for the individual species' correlations (Figure 5C).Interestingly, statistical relationships and significance levels mainly represent mainly the levels of the results from the spruce trees (e.g., rH: r MJJA = −0.74;vp: r MJJA = −0.62,p < 0.01), or are even increase in contrast to the species-specific results (rH: r JJA = −0.75,p < 0.01).The climate-isotope relationships are kept stable within such a combination of inter-species δ 18 O TRC series, clearly indicating the strong influence of the same external/climate factor(s) on both species.This result underlines the use of an inter-species approach for using δ 18 O TRC to reconstruct past climate variability and enhancing the representativity of tree-sites.
Proxy-climate calculations with climate variables from the CRU dataset also reflect the proposed general high similarities of responses among species (Figure S1).Interestingly, differences occur concerning the respective parameter and season the investigated species react most sensitive.Whereas the juniper δ 18 O TRC is highly sensitive to climate parameters related to drought or moisture stress during the summer monsoon period (pet: r MJJA = −0.54,p < 0.01, T Max : r MJJA = −0.51p < 0.01, precipitation r MJJA = 0.53, p < 0.01), spruce δ 18 O TRC is more related to precipitation variations (r MJJA = −0.71,p < 0.01).Especially for the junipers, the location on a southerly exposed slope results in a more drought and insolation related sensitivity.In general, the correlations for junipers are more consistent than those of the spruce trees (Figure S1A,B).Presumably, the CRU dataset: (i) better reflects climate conditions at the upper tree-line; and (ii) is more representative for juniper than for spruce.
The assessment of the CRU-δ 18 O TRC relationship with a species' mean exceeds the already strong results revealed with local instrumental climate data.It is obvious that, by using a species' mean δ 18 O TRC value, the species-specific correlations coefficients can be further increased (Figure S1C).For example, the correlation with potential evapotranspiration during the summer season (May to August) increases up to an r-value of r MJJA = 0.60 (p < 0.01), while species-specific results are noticeably lower (pet juniper: r MJJA = 0.50; pet spruce: r MJJA = 0.54).Similar effects can be obtained for maximum temperature during the same season, with a species-mean r-value of r MJJA = 0.52 in comparison to the δ 18 O TRC of juniper (r MJJA = 0.50, p < 0.01) or spruce (r MJJA = 0.47, p < 0.01).Apparently, by averaging the δ 18 O TRC values of individual species, biases within the δ 18 O TRC variations related to species-specific responses are dampened in favor to a common site signal.

Climate-δ 13 C TRC Relationships
As expectable from the rather low inter-series correlation between the δ 13 C TRC time series, the two conifer species show only minor similarities in their respective climate-proxy relationships (Figure 6).For junipers, the correlations with climate station data reveal mostly weak statistical relationships and therefore a surprisingly weak storage of environmental signals in the δ 13 C TRC (Figure 6A).Highly significant positive correlations with climate parameters and juniper δ 13 C TRC -variations can only be obtained for different temperature variables during spring season (Tmin; r MAM = 0.36, p < 0.01), and the hydrological year (mean Temperature tmp: r HydY = 0.35, p < 0.01).It is apparent that all displayed correlations stay on much lower significance levels than the δ 18 O TRC -juniper relationships.Figure 6B reveals the correlations between δ 13 C TRC values of spruce trees and the climate station in Qamdo.Interestingly, the results are in parts opposite to the correlations for the juniper trees.Especially the change of the sign in the correlation, e.g., for hydroclimate variables such as relative humidity and precipitation, is remarkable.This justifies and confirms that, in the case of δ 13 C TRC , the use of a cross-species mean for evaluating proxy-isotope relationships is not appropriate.Highest positive correlations were obtained for wind speed throughout the year, indicating the strong influence of wind on transpiration and gas exchange on leaf level.Significantly negative relationships are apparent for vapor pressure during the hydrological year (November to October) and spring (r HydY = −0.54,p < 0.01; r MAM = −0.58,p < 0.01), and for relative humidity during the hydrological year and the vegetation period (r HydY = −0.53,p < 0.01; r MJJA = −0.50,p < 0.01).It can be concluded that only the spruce δ 13 C TRC time series at our study site fulfills the preconditions to reconstruct past environmental history and variability.
As for the δ 18 O TRC , we tested climate-δ 13 C TRC relationships with the CRU TS 4.02 dataset for the period 1901-2007 (Figure S2).Interestingly, correlation coefficients as well as the respective significance levels for juniper are more stable throughout the year than for spruce.This agrees with the results obtained by the analyses with δ 18 O TRC (cf.Section 3.3), however, with considerably lower correlation coefficients.For juniper δ 13 C TRC , highest correlations are manifested with different temperature variables and potential evaporation (Figure S2A).Highly significant positive correlations of spruce δ 13 C TRC occur with diurnal temperature range (dtr), and negative correlations with minimum temperatures (tmn) (Figure S1B).

Implications of Species-Specific Proxy-Climate Responses of Multi-Centennial Isotope Series
Within the subalpine setting and design of our study, across-species δ 18 O TRC variations show a very uniform picture in their individual relationships to the tested environmental parameters.The derived climate-proxy correlations clearly indicate that, despite minor differences in exposure and altitude, the same climate and environmental signals are stored in the individual δ 18 O TRC variations.The investigated species are regardless of the used climate dataset (climate station or CRU 4.02 dataset) most sensitive to humidity related climate parameters during the vegetation period (May to August), clearly underlining a common external forcing.Among those, relative humidity, precipitation, and maximum temperature have the highest impact.These relationships can be directly linked to the influence of the Asian summer monsoon system, which leads via high amounts of precipitation during the vegetation period season to a relative depletion of δ 18 O in precipitation (=amount effect) [41].As a result, negative relationships of δ 18 O TRC with hydroclimate parameters such as rH and precipitation occur and are in accordance with plant-physiological models [42,43].Our findings are corroborated through similar results from numerous studies from the TP and the surrounding areas [10][11][12][13][14][15][16][17]27].This underlines the straightforward and reasonable fingerprints of summer monsoon signals in δ 18 O TRC time-series from the summer monsoon-influenced area.
Minor tree-specific differences are apparent in the absolute values of the derived correlation coefficients, where spruce trees in general show higher values.However, the displayed correlation patters are very similar in terms of seasonality and signs.The higher sensitivities of spruce trees point to a higher general climate sensitivity compared to the long-living junipers.Since it is not possible to carry out separate analyses of earlywood (EW) and latewood (LW) due to narrow annual rings in junipers with only few rows of latewood, it remains open whether varying climate sensitivities of EW and LW are responsible for the lower signal strength in juniper.An indication of such different sensitivities was found for instance in other studies and species, e.g., from southwestern and southeastern China [23,25].It has to be noted, however, that these studies are based on significantly shorter time periods of only few years to decades, and use species with significantly wider annual rings compared to ours (Abies and Picea).
A striking result of our study is that both tree species indicate over the full common period of approximately 400 years a very high conformity in different frequency ranges.Hereby, a very high agreement can be found for single positive and negative extreme years as well as for lower frequency ranges, underlining the uniform response to a common forcing.This is confirmed by the climate-proxy relationship of the cross-species mean isotope chronology compared to the individual species chronologies.When analyzing the climate-isotope response of the average value, species-specific correlations are retained highly significant in both, their response and temporal linkage (individual months and seasons).Consequently, the spatial representativity of both oxygen time series should be expectedly similar.In fact, the analysis of spatial correlations with the three different δ 18 O TRC datasets using the KNMI climate explorer (https://climexp.knmi.nl)allows this relationship to be further quantified (Figure 7).The discussed higher sensitivity of spruce trees compared to junipers also results in consistently higher spatial correlations with different humidity and drought parameters.These are insofar as from a theoretical perspective of isotope fractionation and fixation reasonable, since they display an often reported negative (MJJA precipitation amd MJJA scPDSI) or positive relationship (increased evaporation during the summer season via T Max ) mainly caused by a precipitation amount or temperature effect during the summer monsoon season.Interestingly, Figure 7 displays that the high spatial correlations remain stable with the calculation of a species' mean and thus validates a general high common sensitivity to hydroclimate variables and large spatial representativity of multi-century oxygen isotope time-series across species.This uniform behavior across species further underlines the suitability and applicability of annually resolved, tree-ring based stable oxygen isotope time series to serve as an independent dating control for other archives (e.g., ice cores), as suggested in [44].It is striking that, in contrast to the investigated oxygen isotope series, the derived δ 13 CTRC series reveal a low common signal across species combined with a missing uniform response to the tested environmental parameters.This confirms results of previous studies on δ 13 CTRC from China, where no common cross-regional proxy signal has been found [25,44,[45][46][47][48]. Instead, these studies suggest that the respective δ 13 CTRC series rather display more site-related ecophysiological responses on a changing environment.This can be underlined by the variety of interpretations according to which the respective δ 13 CTRC series serves as a proxy for changes in the plant-internal water-use strategy [25,28], as a proxy for temperature during the vegetation period [45], and for summer temperatures [46].Only few studies indicate a dependency of δ 13 CTRC variations on moisture conditions (including precipitation and rH) [47,48].This is probably connected to a generally drier setting at the investigation site and therefore a higher drought limitation at these sites compared to our location.
Although we can infer for both studied species at our study site similar macro-and mesoclimatic conditions and soil characteristics, differing physiological factors and processes seem to overprint or bias the detection of dominant environmental signals in our tree-ring δ 13 C.As mentioned in [33, [38][39][40], 13 C variations in tree-rings are on the one hand controlled by external environmental factors that regulate the uptake of atmospheric CO2 via the stomata of the leaves.Stomatal conductivity is in turn mainly determined by factors such as rH, temperature, moisture and wind.On the other hand, the It is striking that, in contrast to the investigated oxygen isotope series, the derived δ 13 C TRC series reveal a low common signal across species combined with a missing uniform response to the tested environmental parameters.This confirms results of previous studies on δ 13 C TRC from China, where no common cross-regional proxy signal has been found [25,[44][45][46][47][48].Instead, these studies suggest that the respective δ 13 C TRC series rather display more site-related ecophysiological responses on a changing environment.This can be underlined by the variety of interpretations according to which the respective δ 13 C TRC series serves as a proxy for changes in the plant-internal water-use strategy [25,28], as a proxy for temperature during the vegetation period [45], and for summer temperatures [46].Only few studies indicate a dependency of δ 13 C TRC variations on moisture conditions (including precipitation and rH) [47,48].This is probably connected to a generally drier setting at the investigation site and therefore a higher drought limitation at these sites compared to our location.
Although we can infer for both studied species at our study site similar macro-and mesoclimatic conditions and soil characteristics, differing physiological factors and processes seem to overprint or bias the detection of dominant environmental signals in our tree-ring δ 13 C.As mentioned in [33, [38][39][40], 13 C variations in tree-rings are on the one hand controlled by external environmental factors that regulate the uptake of atmospheric CO 2 via the stomata of the leaves.Stomatal conductivity is in turn mainly determined by factors such as rH, temperature, moisture and wind.On the other hand, the carbon isotopes of CO 2 are subject to a fractionation at the enzyme RUBISCO dependent on the respective CO 2 concentration within the leaves.This, e.g., results in a discrimination in favor to the lighter 12 C and against the heavier 13 C isotopes during photosynthesis.Therefore, sensitivities and responses of the individual δ 13 C TRC to humidity parameters like rH and maximum temperature would be expectable.However, this is not straightforward displayed in our results.Instead, results of climate-isotope correlations with relative humidity indicate an indistinct, species-specific response during the vegetation period (positive for junipers, negative for spruce).
The displayed negative results for correlations between δ 13 C TRC of spruce and rH are in so far expectable and reasonable since a higher rH during the vegetation period linked to more open stomata leads to a depletion in δ 13 C TRC .This relationship stays stable for different seasonal means and single months during the vegetation period.Therefore.a seasonally higher moisture availability during the summer monsoon season apparently does not seem to be decisive for the regulations of the stomata.Presumably, moisture is not the main limiting factor for spruce during the summer monsoon season, although monsoonally triggered higher temperatures and humidity are then predominant.A higher moisture availability during the summer season leads to a further increase in the discrimination against 13 C.For juniper, the explanation of the positive relationship between δ 13 C TRC and rH remains at this point challenging.As mentioned, conditions with higher humidity during the summer monsoon season would expectably be associated to more favorable conditions resulting in wider opened stomata.This will in turn result in lower δ 13 C TRC values.Thus, a negative correlation between δ 13 C TRC and the relative humidity would be expected as consequently displayed by the spruce δ 13 C TRC .Since correlations with a more radiation related parameter such as sunshine hours for the junipers do not display any significant result, possible short-range differences in humidity between both tree stands can be excluded as the reason for this.At this stage, explanations for such differences between species cannot be evaluated without knowledge of possibly underlying and possibly biasing plant-internal factors.
It remains striking that the climate δ 13 C TRC relationships are compared to the results derived from the δ 18 O TRC significantly lower and somehow inconsistent in terms of a uniform (site and species) signal and representativity.In addition, the underlying long-term trend especially apparent for the spruce trees clearly hampers the decoding of climate signals in our δ 13 C TRC series.To evaluate the influence of low-frequency variations on the different species and isotopes, high-pass filtered versions of all time-series were correlated against the climate station data.In the δ 18 O TRC series, the removal of the low-frequency variations which can be visualized through wavelets (Figure S3) results in slightly less strong, but still consistent correlation patterns.For the δ 13 C TRC series, the high-pass filtering removes the significance in the juniper-climate relationships almost completely.Contrasting, the correlations of the spruce δ 13 C TRC series with climate are enhanced in the positive (temperature) as well as negative (precipitation, rH, and vapor pressure) directions.Although correlations vary over time, this suggests that species-specific correlations with climate data since the 1960s are highly dependent on the chosen trend correction method.This is additionally corroborated by the difference displayed in moving correlations between the raw δ 13 C TRC series and the CO 2 -corrected δ 13 C TRC series (Figure S4).Such an influence of the chosen trend correction is in line with the results of [45].

Long-Term Trends in Carbon Isotope Ratios and their Implications for Recording Climate Signals
As displayed in Figure 4, even the corrected δ 13 C TRC time series do not indicate a uniform signal between species.In contrast to δ 18 O TRC , the δ 13 C TRC series reveal for both species a highly individual year-to-year variability as well as negative long-term trends.In addition, since 1900, a significantly more negative trend in spruce δ 13 C TRC values as compared to the juniper values is evident.By taking into account approximately comparable site conditions for both studied species, only certain internal physiological drivers can probably lead to the stronger depletion of the heavier 13 C isotopes in the spruce trees.Thus, we assume that species-specific ecophysiological responses to increasing atmospheric CO 2 or to increasing temperatures cause differences in the uptake of heavier 13 C in the tree-rings.The impact of this changing external forcing and the resulting different response of the two tree species become even more impressive when considering the δ 13 C TRC in a multi-centennial perspective.Global and regional studies show that trees tend to respond to an increase in atmospheric CO 2 through various ecophysiological adaptations, such as modulations in the response of photosynthesis, stomatal conductance, or a response in the water-use efficiency (iWUE) [49][50][51][52][53][54].It is likely that the substantial disparities and trends in our δ 13 C TRC series display species-specific ecophysiological reactions to a rising atmospheric CO 2 concentration and/or to a substantial altitude-dependent increase in temperature on the TP [55,56].Recently, it is reported that even within a sampling site such reactions can be species-specific [57].Consequently, such plant-internal adaptations will in turn lead to a blurred picture of climatic signals stored in the respective δ 13 C TRC values.The apparent negative δ 13 C TRC after correction for atmospheric CO 2 change can be therefore one reason for the weak correlations of δ 13 C TRC with climate parameters.In addition, the choice of the suitable correction factor can have an influence that should not be underestimated, as reported in [43].Probably our carbon isotope time-series indicate an individual response of the investigated tree species towards the changing CO 2 concentrations rather than recording an environmental signal.
Figure 8 shows the iWUE of the two investigated tree species for the period 1900-2007.With respect to the reported environmental change in terms of a general warming trend since the 1950s, we additionally split the investigation period into two individual segments with a threshold in 1955.It is apparent that the displayed individual iWUEs differ in their absolute amounts as well as in their long-term trends for the full period between the species.This points to different species-specific physiological strategies and adaptations either with regard to changes in climate and/or changes in the atmospheric CO 2 concentration.Whereas the iWUE of juniper trees is characterized by a more or less stable phase throughout the early period (−2% for the period 1900-1954), a steady and strong increase in the iWUE is apparent since the 1960s (+20%), reaching a maximum in 2005.With regard to the full period from 1900 to 2007, the relative increase in iWUE for juniper is 17%.In comparison, spruce tree iWUE is characterized by a slight increase (+1%) during 1900-1954.Similar to the juniper trees, spruce trees show a remarkably positive trend towards a higher iWUE for the recent six decades and therefore an adaptation for reducing water loss through transpiration.However, the absolute increase in iWUE of spruce since 1955 is markedly lower (+14%) and also evolves in its long-term rate during 1900-2007 generally slower than in juniper (+12%).Juniperus as a more drought tolerant species can probably better respond and react to the ongoing environmental changes with rising temperatures and lower humidity with a faster acceleration within their respective iWUE.On the other hand, more favorable microclimatic conditions with higher humidity at the spruce stands might help this species to better cope with the increase in temperature.The different trends and responses in iWUE indicate that plant physiological adaptations to changing external environmental parameters can be detected for both tree species.These plant physiological adaptations aggravate a more straightforward interpretation of our climate-proxy relationships.A necessary deciphering need in turn further analyses on the responsible plant-internal drivers.The future use of our δ 13 C TRC time series for climate reconstruction appears at this stage therefore challenging.Although admittedly reflecting an initial stage, our results on the iWUE are in good accordance to existing studies from the Tibetan Plateau dealing with Pinus and Abies species, where similar positive trends in iWUE were reported for the last decades [28,58].In addition, on a global scale, such responses can be observed in various tree species, e.g., from the tropics [59,60] or more temperate zones [53,[61][62][63][64][65][66].
1 We are aware that substantial ecophysiological investigations on and interpretations of iWUE need more knowledge about the corresponding changes in controlling plant physiological and environmental factors (e.g., c i , c a , and the carbon isotope discrimination ∆).Since this study mainly targeted a first estimation of the usability and the included potentials and limitations of a dual isotope and dual species approach at a single site, we focused on a first qualitative description of ecophysiological findings.Therefore, the identification and quantification of the proportion or the combination of, e.g., temperature rise and CO 2 concentration, on the two tree species clearly exceeds the scope of the presented manuscript.More complex and further comprehensive ecophysiological analyses of the obtained results will be subjects of an upcoming study.

Conclusions
This study represents the first substantial comparison of annually resolved, multi-centennial tree-ring δ 13 C and δ 18 O isotope time series from a high mountain site at the southeastern Tibetan Plateau.Although the external forcings and general site characteristics can be regarded as highly comparable, a thorough evaluation on responsible environmental factors causing the respective δ 13 C TRC and δ 18 O TRC variations in the tree species revealed in part substantial differences.
The displayed δ 18 O TRC variations of juniper and spruce trees show a congruent response to the analysed climate parameters, indicating a uniform response across species to a common forcing.As best results were obtained for various hydroclimatic parameters during the vegetation period, notably a summer monsoon signal was herein recorded.By combining and averaging the δ 18 O TRC time series of our investigated species, existing highly significant proxy-climate relationships on selected hydroclimatic variables could be in most cases retained.Microclimatic influences and/or species-specific differences in the individual tree species are presumably reduced in favor to a more representative and common site signal.Thus, when analyzing δ 18 O in tree-rings, the expressiveness of δ 18 O TRC -time series can be enhanced by studying two tree species at one location.
In contrast, only weak relationships occur between the δ 13 C TRC series and tested environmental parameters.This is all the more astonishing as the locations of both tree stands are close to each other.Consequently, mainly plant-internal strategies and species-specific effects or adaptations are biasing the recording of climate signals.In combination with a-still strongly preserved-long-term trend after atmospheric CO 2 trend correction in spruce, several questions remain open for an unambiguous interpretation of the δ 13 C TRC variations.Both studied species apparently reflect an individual response to environmental change.Juniper trees situated at higher elevations and southerly exposed react more sensitively with regards to the increasing irradiation and decrease in cloud cover and rH.This is for example evident from the changes in the iWUE.For spruce trees, the microclimatic advantage of a comparably higher humidity may be favorable, resulting in a delayed, lowered and lagged increase in iWUE.A verification of this assumption, however, requires a separate investigation, which should also include studies on seasonal isotope variations and coupled plant-internal effects.
Nevertheless, such species-specific investigations on similar settings might help to predict the influence and impact of future climate change on the adaptation ability of dominant tree species in the area of the treeline ecotone on the Tibetan plateau.For a quantification of such plant-internal responses, future studies should evaluate the complex interaction on the impact of a probable CO 2 -fertilization effect on iWUE, a general rise in temperature on the photosynthetic rate, and the interplay of growth and iWUE.(50-year windows) between the raw δ 13 C series of spruce and juniper (in violet) and the CO 2 -corrected δ 13 C series of the investigated two tree species (in orange).
local climate as recorded by the nearby climate station in the municipality of Qamdo (31 • 09 N/97 • 10 E; available period 1954-2012) is characterized by a mean annual temperature of 8.8 • C and an annual precipitation sum of 482 mm, indicating a sub-humid continental climate (Figure 1).The main summer season lasts from June to August and represents the period with highest temperatures (T JJA = 16.7 • C).Summer months with the maximum monthly temperatures are July (T July = 17.4 • C) and August (T Aug = 16.6 • C), respectively.The summer months July and August have the highest amounts of precipitation (Precip JJA = 296 mm; Precip July = 107 mm, Precip Aug = 101 mm).In contrast, winter months are cold and show little to no rainfall/snow (T Dec = 0 • C, Precip Dec = 2mm).Overall, local climate conditions display a typical summer monsoon climate with warm and humid summers contrasted by cold and dry winter conditions [4].

Figure 1 .
Figure 1.Overview of the investigation site: (A) location of the site Qamdo within China; (B) climate diagram of the climate station in Qamdo; and (C) digital elevation map of the Qamdo region with the investigation site (red box), location of the city and climate station of Qamdo (red dot).Yellow points indicate locations of the used 0.5° CRU TS 4.02 grid points.

Figure 1 .
Figure 1.Overview of the investigation site: (A) location of the site Qamdo within China; (B) climate diagram of the climate station in Qamdo; and (C) digital elevation map of the Qamdo region with the investigation site (red box), location of the city and climate station of Qamdo (red dot).Yellow points indicate locations of the used 0.5 • CRU TS 4.02 grid points.

Figure 2 .
Figure 2. Overview on the sampling site near Qamdo: (I) area of predominant Juniperus tibetica treestands on southerly expositions at 4400-4600 m asl.; and (II) vegetation belt with the Picea balfouriana forest on southeasterly expositions.

Figure 2 .
Figure 2. Overview on the sampling site near Qamdo: (I) area of predominant Juniperus tibetica tree-stands on southerly expositions at 4400-4600 m asl.; and (II) vegetation belt with the Picea balfouriana forest on southeasterly expositions.

3. 1 .
Figure3shows the resulting δ18 O TRC time series of the two studied tree species at Qamdo during 1685−2007.For a better comparability of the low-frequency variations between the δ 18 O TRC time series, the scale of the individual series was adjusted by an offset of 4‰.A first quantitative examination of the two time-series already reveals a high visual similarity (Table2).The mean values of the oxygen isotope series differ slightly within a range of about 0.6‰ (δ18 O TRC Juniper = 23.43‰,δ 18 O TRC Spruce = 22.98‰).Since both tree-species are situated within a short elevational range, this is assumed to be caused by slightly different plant-internal fractionations during cellulose synthesis and/or-on a microclimatic scale-a slightly differing moisture availability during the vegetation period.Probably caused by a higher proportion of earlywood[11,12], the juniper chronology is characterized by a higher autocorrelation (AC1 = 0.61) than the spruce time-series (AC1 = 0.26).It is apparent that the inter-annual variability of the individual time-series is considerably high (Figure3).In particular, during the recent decades, the number of positive and negative excursions in both time series exceeding a deviation of more than 1.5 SD increased compared to their long-term mean.

Figure 3 .
Figure 3. Annual δ 18 O TRC variations of Juniperus tibetica (blue) and Picea balfouriana (red) at Qamdo for the period 1685−2010.Bold lines represent the respective 11-year Fast Fourier Transform (FFT) filter for each time-series, dotted horizontal lines the individual species mean.Note: Individual y-axes are shifted against each other with an offset of 4‰ to better show the high synchronicity between both time-series.

Figure 4 .
Figure 4. Annual δ 13 C TRC variations of Juniperus tibetica (blue) and Picea balfouriana (red) at Qamdo for the period 1685−2011.Bold lines represent δ 13 C TRC CO 2 -corrected series and associated 11-year Fast Fourier Transform (FFT) filters, bright dotted time series indicate the species specific uncorrected δ 13 C TRC raw series, and dotted lines represent species-specific linear trends of the corrected δ 13 C TRC - series over the pre-industrial (1685−1849) and industrial period (1850−2007), respectively.

Figure 5 .
Figure 5. Climate-proxy relationships between δ 18 O TRC of: juniper trees (A); spruce trees (B); and the site-internal δ 18 O TRC mean (C) plotted against monthly and seasonal averages of different climate parameters (pre, precipitation; rH, relative humidity; sh, sunshine hours; tmn, minimum temperature; tmp, mean temperature; tmx, maximum temperature; vap, vapor pressure deficit; ws, wind speed) from Qamdo climate station for the period 1954-2007 (n = 54).The prefix "p" indicates relationships to the previous years' monthly means.Red/blue colors indicate positive/negative correlation coefficients (r).The respective level of significance is displayed by a plus sign (p < 0.05) or an asterisk (p < 0.01).

Figure 6 .
Figure 6.Climate-proxy relationships between δ 13 C TRC of juniper trees (A) and spruce trees (B) against different monthly and seasonal means of climate parameters (pre, precipitation; rH, relative humidity; sh, sunshine hours; tmn, minimum temperature; tmp, mean temperature; tmx, maximum temperature; vap, vapor pressure deficit; ws, wind speed) from the Qamdo climate station for the period 1954-2007 (n = 54)."p" shows relationships to the previous years' monthly means.Red/blue color indicates a positive/negative correlation coefficient (r).The respective level of significance is indicated by a plus sign (p < 0.05) or an asterisk (p < 0.01).

Figure 7 .
Figure 7. Spatial correlations (p < 0.01) of δ 18 O variations in juniper, spruce and the species' mean (vertically arranged) with climate variables using the Climate Explorer.(A-C) Selected climate variables for the individual species (juniper and spruce) and the calculated δ 18 O mean of our study site for the summer monsoon season (May to August).

Figure 7 .
Figure 7. Spatial correlations (p < 0.01) of δ 18 O variations in juniper, spruce and the species' mean (vertically arranged) with climate variables using the Climate Explorer.(A-C) Selected climate variables for the individual species (juniper and spruce) and the calculated δ 18 O mean of our study site for the summer monsoon season (May to August).

Figure 8 .
Figure 8. Intrinsic water-use efficiencies for Juniperus tibetica and Picea balfouriana at Qamdo.Dashed line indicates dividing line between the investigation periods, while old lines represent the individual linear trends of both species for the periods 1900-1954 and 1955-2007, respectively.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2076-3263/9/4/151/s1,FigureS1.Climate-proxy correlation of: juniper δ18 O TRC (A); spruce δ 18 O TRC (B); and the site-internal δ 18 O TRC mean (C) against monthly and seasonal means of different gridded climate parameters from the CRU TS4.01 dataset for the period 1954-2007 (cld, cloud cover; dtr, diurnal temperature range; frs, number of frost days; pet, potential evapotranspiration; tmn, minimum temperature; tmp, mean temperature; tmx, maximum temperature; wet, wet days frequency)."p" shows relationships to the previous years' monthly means.Red/blue colous indicate positive/negative correlation coefficients (r).The respective level of significance is indicated by a plus sign (p < 0.05) or an asterisk (p < 0.01); Figure S2.Climate-proxy correlation between juniper δ 13 C TRC (A) and spruce δ 13 C TRC (B) against different monthly and seasonal means of gridded climate parameters from the CRU TS4.01 dataset for the period 1954-2007 (cld, cloud cover; dtr, diurnal temperature range; frs, number of frost days; pet, potential evapotranspiration; tmn, minimum temperature; tmp, mean temperature; tmx, maximum temperature; wet, wet days frequency)."p" shows relationships to the previous years' monthly means.Red/blue colors indicate positive/negative correlation coefficients (r).The respective level of significance is indicated by a plus sign (p < 0.05) or an asterisk (p < 0.01); Figure S3.Wavelet analyses (Morlet) for both isotopes and species during the common period (1685-2007); (left) original values; and (right) high-pass filtered version.Periods of significant periodicities are framed by white lines; Figure S4.Moving correlations

Table 1 .
General site information and chronology statistics for the common period 1685-2007.

Table 1 .
General site information and chronology statistics for the common period 1685-2007.

Table 2 .
Correlation coefficients (r) and Gleichläufigkeit (GLK) between stable carbon and oxygen isotope time-series of Juniperus tibetica and Picea balfouriana on different frequency scales for the common period 1685-2007.All given values are highly significant at p < 0.01.