青海省同德县乔木状甘蒙柽柳的年龄及其生长对环境的响应
方欧娅1,*, 贾恒锋2, 邱红岩1, 任海保1
1中国科学院植物研究所植被与环境变化国家重点实验室, 北京 100093
2三峡大学生物与制药学院, 湖北宜昌 443002

* 通信作者Author for correspondence(E-mail:oyfang@ibcas.ac.cn)

摘要

青海省同德县高原河漫滩地区存在乔木状的甘蒙柽柳( Tamarix austromongolica), 且在即将修建的水电站库区成林分布, 引发了学者和公众在其是否属于古树而需要保护以及如何保护问题上的广泛关注。该研究利用树木年轮学方法, 调查了水库淹没区和非淹没区中胸径较大的甘蒙柽柳的树龄, 分析了甘蒙柽柳的生长历史及其对环境因素的响应。研究发现胸径较大的甘蒙柽柳可由多株合生形成, 测量胸径与植株个体年龄之间无直接联系。该地区甘蒙柽柳径向生长量在20世纪70-80年代迅速提升, 80年代后生长趋于平稳。甘蒙柽柳的径向生长响应于7、8月份的水分条件, 径流的增加在一定程度上促进了甘蒙柽柳的生长。未成熟植株的径向生长与生长季的风速有着显著的负相关关系。研究结果对甘蒙柽柳特殊形态的形成认知具有理论意义, 对该树种保护方案的评估具有实践意义。

关键词: 甘蒙柽柳; 乔木状; 树轮; 径流; 风速
doi: 10.17521/cjpe.2017.0088
Age of arboreous Tamarix austromongolica and its growth response to environment in Tongde County of Qinghai, China
FANG Ou-Ya1,*, JIA Heng-Feng2, QIU Hong-Yan1, REN Hai-Bao1
1State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, the Chinese Academy of Sciences, Beijing 100093, China
and2College of Biology and Pharmacy, China Three Gorges University, Yichang, Hubei 443002, China
Abstract

Aims The objectives were to identify the age of the arboreous Tamarix austromongolica in the flood plain area of the Qinghai Plateau and clarify the response patterns of T. austromongolica’s growth to the environmental factors. We focused on social issues about whether the T. austromongolica should be protected and how to protect in the reservoir area of a hydropower station. Methods In this study, arboreous T. austromongolicain both reservoir submerged and non-submerged areas were sampled and measured based on the dendrochronology method. The ages were estimated based on the geometrical characteristics of the pith and the identified age of the inner ring. The correlation and response analysis showed the relationship between T. austromongolica’s growth and environmental factors. Important findings We accurately determined the age and historical growth dynamics of the T. austromongolica with large diameter at breast height ( DBH). The results showed a special accretion phenomenon in arboreous T. austromongolica, which accelerated the DBH increasing, i.e. no direct relationship existed between the plants’ DBH and ages of the individuals. Radial growth of T. austromongolica, increased rapidly in the 1970s and 1980s and began to stabilize in the late 1980s, and mainly responded to the runoff in July and August of the Yellow River. Increasing runoff would promote the radial growth of T. austromongolica. The growth of the immaturate plant showed significant negative correlation with the wind speed in the growing season. The results will be of theoretical significance to the formation of the special morphology of the T. austromongolica, and will provide scientific practical guidance in designing the protection schemes.

Key words: Tamarix austromongolica; arboreous; tree ring; runoff; wind speed

在青海省同德县巴沟乡然果村黄河岸边海拔约2 660 m的河漫滩沙地上, 存在约10 hm2乔木状且具较大树干胸径的甘蒙柽柳(Tamarix austromongolica)林, 位于即将建设的羊曲水电站水利工程淹没区。这片甘蒙柽柳林是否存在古树而应当被保护, 以及如何实施保护, 受到了学者和公众的广泛关注, 且引起了极大的争议(吴玉虎, 2012; 新华社, 2016; 中国周刊, 2016), 从而直接影响了水利工程的进程和生态保护工程的开展。我们以水利工程即将淹没地带的然果村甘蒙柽柳林为主要研究区, 结合以工程淹没区外的班多村甘蒙柽柳林为辅助研究区, 针对社会关注问题开展了树木年轮学研究, 旨在回答以下两个重点问题: 1)甘蒙柽柳林中是否存在年龄较大的古树; 2)该地区甘蒙柽柳特殊的乔木状形态是如何形成的, 甘蒙柽柳生长历史及其与环境条件的关系是怎样的?

树木年轮学在历史定年(邵雪梅等, 2007; Li et al., 2015)以及树木生长对环境响应的评估(Babst et al., 2013; Gou et al., 2015; Zhang et al., 2015)等方面有着突出的贡献。树木的生长受到环境因素的限制作用, 其长期的轮宽变化与气候、环境变化之间存在一定的相关关系(Li et al., 2016; Yin et al., 2016), 从而使得树木的径向生长对气候、环境等因素具有较为明确的响应(Fang et al., 2016; Qin et al., 2016; Wang et al., 2017)。明确这种响应关系有助于评估气候变化、环境变迁对树木生长的影响(Rohner et al., 2016; Gao et al., 2017), 为未来气候变化或者人为改变(如移植等)的植物生境条件下的植物生长评估提供依据。

目前针对柽柳属(Tamarix)植被的研究对象通常为幼龄树, 成熟甘蒙柽柳分枝结构复杂, 干形较差而较难开展相关研究(Brotherson et al., 1984; 肖生春等, 2005; Xia et al., 2016)。以柽柳属植物为研究对象的树木年轮学工作主要集中在树轮水文学方面, 附近河流湖泊或地下水的水位变化与柽柳生长通常具有显著关系(肖生春等, 2004; Xiao et al., 2014), 从而使得在气候变化影响下, 柽柳生长状况、森林群落结构发生改变(Horton & Clark, 2001)。另一方面, 柽柳对水分条件的响应有助于分析历史时期水文条件变化(Xiao et al., 2016)。由于柽柳独特的生长环境, 其埋藏亚化石样本可以应用于古环境的探究(Lipp et al., 1996; Yang et al., 1996; Frumkin, 2009)。以柽柳为对象的树木年轮学研究在多个学科中具有一定的应用价值, 如研究者利用柽柳年轮对堆积物的响应, 判断河漫滩沉积年, 沙丘形成演化, 从而进行地貌学相关研究(Friedman et al., 2005; Birken & Cooper, 2006; Wang et al., 2010; 张锦春等, 2014)。基于前期科研工作者的成果, 我们提出假设, 认为位于高原河漫滩的甘蒙柽柳的生长同样也响应于河流水文条件等环境因子的变化。然而, 针对近年来发现的研究区特殊形态的甘蒙柽柳, 其历史生长过程及环境响应关系尚未见报道。

本研究以青海省同德县乔木状甘蒙柽柳为对象, 调查胸径较大的植株年龄, 分析甘蒙柽柳径向生长与流域径流等环境要素的关系。对该甘蒙柽柳林大树年龄的准确判断是评估其生态价值所亟待解决的重要问题, 同时, 甘蒙柽柳的生长对环境的响应也是评估其在移植环境下生长状况的重要依据。

1 数据和方法
1.1 研究区概况

研究区属于秦岭地槽褶皱系的青海南山冒地槽带, 断裂构造发育。黄河对高原面不断侵蚀, 形成阶地和较为宽大的河谷, 并在河流弯曲凸岸形成河漫滩。河漫滩以石砾质冲积土和沙质冲积土为主, 土层较厚, 土壤发育程度差。在此之上分布有少量耐旱耐水抗风且耐碱土的植被, 典型的有甘蒙柽柳和小叶杨(Populus simonii)。该区域为温带大陆性干旱气候, 常年低温干旱, 据兴海气象站1981-2010年数据统计, 全年平均气温1.7 ℃, 最暖月在7月, 平均气温12.7 ℃, 最冷月在1月, 平均气温-11.0 ℃。降水稀少, 且主要集中在5-9月, 最大降水量出现在7月, 平均7月降水总量约86.6 mm, 年降水总量约为377.8 mm。2-4月为该地大风发生频率最高的季节, 受到山地地形的影响, 东西向山谷风偏大。

本研究主要针对羊曲水电站淹没区然果村甘蒙柽柳林(35.54° N, 100.16° E, 海拔2 678 m)开展研究, 同时也涉及了库区附近非淹没区班多村甘蒙柽柳林地(35.33° N, 100.26° E, 海拔2 730 m)(图1)。然果村林地甘蒙柽柳种群呈斑块状分布, 径级呈现较强的正偏态和尖顶分布特征, 胸径集中分布在20-30 cm。林地具有大量胸径超过30 cm的植株, 其中测量胸径值最大为132.1 cm, 该植株被媒体称为“ 柽柳王” 。胸径20 cm以上的植株平均树高约为8 m。班多村林地甘蒙柽柳种群径级呈对称、常峰态分布, 较然果村具有更多的胸径超过30 cm的植株, 然而超过50 cm的植株在该区域十分罕见。两地甘蒙柽柳平均胸径和树高无显著差异, 且林下均缺乏幼苗更新。

图1 青海同德县甘蒙柽柳采样点。Fig. 1 Sampling sites of Tamarix austromongolica in Tongde, Qinghai.

1.2 树轮样本采集及测量

我们于2016年11月下旬赴青海省同德县开展了实地调查, 利用直径5.15 mm的生长锥着重针对胸径较大的甘蒙柽柳进行采样并同时对其胸径进行测量。在然果村以及班多村分别采集了37株和15株活立木的树芯样本, 同时分别获取了4株和1株倒木圆盘样本(表1)。针对个别胸径较大且疑似多株合生树木获取多根样芯, 如“ 柽柳王” (编号RG-X05), 我们进行了5个方向的树芯采集。由于然果村甘蒙柽柳林普遍受到泥石流等灾害的影响, 树木个体存在损伤, 采样时尽量避免树木损伤部位, 对较为明确的竖直干材进行样本采集。依照同样的方法, 在班多村进行了树轮样本的采集。取样中, 尽量钻取到树木髓心位置或者能够在树芯上显示出树轮圆弧的近树木髓心位置, 为后续的树龄准确判断提供依据。

表1 同德县甘蒙柽柳采样树木年龄 Table 1 Ages of sampling Tamarix austromongolicain Tongde County

在实验室, 将树芯样本固定后连同圆盘样本用砂纸进行精细打磨, 直至在显微镜下能够清晰显示木质部细胞。继而利用Lintab 6树轮宽度测量仪(Rinntech, Heidelberg, Germany)对所有样本的树轮宽度进行测量, 测量精度为0.001 mm。

1.3 定年及建立年表

采样时间2016年11月甘蒙柽柳当年生长已完成, 树皮内部木质部最外轮判断为2016年生长, 依此对树轮进行初步定年。对比研究区所有样本的树轮宽度变化, 结合骨架图和COFECHA程序(Holmes, 1983)进行准确的交叉定年, 从而判断树轮样本年轮的确切形成年份。在对甘蒙柽柳历史生长状况进行初步了解后, 通过负指数函数或直线拟合并去除其生长趋势, 建立标准化树轮年表。

1.4 树龄判断

具有髓心的树木均以最内轮年份为树木年龄的判断依据, 部分样芯采集时获取到近髓心位置的圆弧形特征, 我们依据其几何结构判断最内轮距离髓心长度从而估算树木的实际年龄(Norton et al., 1987)。圆盘样本均具有髓心, 经过交叉定年判断倒伏年份, 获取实际树龄值。需要说明的是, 这里所述的实际树龄均为树干处年龄, 植株生长到采样高度年份为生长起始年。在树芯样本中存在22株22根样芯上不存在髓心圆弧特征, 无法准确获知树木年龄, 仅记录样芯最内轮年份, 判断其最小年龄。

1.5 相关及响应分析

为了探究甘蒙柽柳径向生长与气候、环境要素之间的相关及响应关系, 本研究将树轮标准年表与气候环境因子月值数据集的数据关系进行分析。所涉及的气候数据包括唐乃亥水文站1956-2007年黄河年径流系列统计以及兴海气象站1960-2015年月平均温度、降水、风速数据集。相关及响应分析均利用Dendro 2002软件(Biondi & Waikul, 2004)完成, 其中响应分析采用自举法随机执行500次。由于前一年环境因子对树木当年生长有着一定的作用, 故而研究选用前一年生长季末至当年生长季末的月值数据进行分析。

2 结果
2.1 甘蒙柽柳的年龄及多株合生现象

在52株活立木的57根树芯样本中, 有8株9根树芯样本具有树木髓心部分, 该类样本经过交叉定年即可获取树木的实际年龄值。样芯中, 24株树的26根树芯最内侧树轮显示出有接近髓心的圆弧形特征, 间接获取其树木的年龄值。5个圆盘样本均获取准确年龄。结果表明, 最大年龄115年, 源于一株倒木圆盘样本RG-P41, 另外, RG-X15A样本未采集到近髓心部位, 然而实测样本最内轮形成年份已达公元1931年, 立木实际年龄可能接近或超过100年。

研究着重对然果村“ 柽柳王” (实测胸径132.1 cm)进行了树龄分析, 5个位置取样分析的结果显示, 每个树芯样本在距树皮12-22 cm处显示了各自的髓心。交叉定年工作发现各树芯的年龄具有较大差异, 最小年龄为31年, 最大年龄为86年, 证明了该“ 柽柳王” 的确为多株合生。此外, 该多株合生的植株其中一根样芯中包含2段树芯样本, 每段样本经交叉定年发现它们具有同时生长的时段, 故而进一步判断该甘蒙柽柳为至少包含2层的多株合生植株。

2.2 甘蒙柽柳生长动态

1960年以前, 样本量较小, 年表受到个体特征影响较大。20世纪70年代, 甘蒙柽柳生长逐渐加快, 且综合同期同地区甘蒙柽柳区域每木检尺工作, 甘蒙柽柳林植株年龄集中分布在20-40年之间, 即在该时期存在大量新生个体。20世纪80年代初以来, 甘蒙柽柳径向生长趋于稳定。

研究区甘蒙柽柳年径向生长量个体间差异较大, 在取样树木样本中, 树木平均年生长速率在1.8 mm· a-1至8.3 mm· a-1之间。然果村与班多村的树木历史生长过程没有显著差异, 两地的树轮宽度年表(图2)显示, 高频及低频趋势均基本吻合。其中, 班多村极端低值波动幅度略小, 表现出更为稳定的生长动态。

图2 甘蒙柽柳树轮标准年表。BD, 班多村; RG, 然果村。Fig. 2 Tree ring standard chronologies of Tamarix austromongolicas. BD, Banduo Village; RG, Ranguo Village.

两地年表均显示在20世纪70年代甘蒙柽柳径向生长多次出现极低值, 并在1988年、1992年、1995年、2002年、2008年以及2015年同时具有极低值, 该时期或年份甘蒙柽柳的径向生长量较低。

2.3 甘蒙柽柳径向生长对环境的响应

同德县甘蒙柽柳的径向生长主要与黄河径流流量以及风速存在较为显著的相关性。甘蒙柽柳径向生长与7、8月黄河径流流量之间存在正相关关系, 显著与7月径流正相关及响应关系(p< 0.05)体现在班多村研究区(图3), 在然果村则与8月的径流量间的相关性更强(图4)。

图3 班多村甘蒙柽柳生长与径流量及风速的相关及响应。左侧为甘蒙柽柳径向生长与逐月径流流量的相关及响应, 右侧为径向生长与月平均风速的相关及响应。横坐标p表示前一年的月份, 误差条表示95%置信区间。Fig. 3 Correlation and response relationship ofTamarix austromongolica growth to monthly runoff and wind in Banduo Village. Left for the relationship between growth and runoff; Right for the relationship between growth and wind. p in the abscissa indicates the month of the previous year, and the error bars show the 95% confidence intervals.

图4 然果村甘蒙柽柳生长与径流量及风速的相关及响应。左侧为甘蒙柽柳径向生长与逐月径流流量的相关及响应, 右侧为径向生长与月平均风速的相关及响应。横坐标p表示前一年的月份, 误差条表示95%置信区间。Fig. 4 Correlation and response relationship ofTamarix austromongolicagrowth to monthly runoff and wind in Ranguo Village. Left for the relationship between growth and runoff; Right for the relationship between growth and wind. p in the abscissa indicates the month of the previous year, and the error bars show the 95% confidence intervals.

甘蒙柽柳径向生长与风速之间的显著负相关主要发生在植被生长季, 在班多村, 甘蒙柽柳生长与当年4-8月的风速之间存在显著的负相关(p < 0.05), 明确的显著负响应发生在当年的5-6月(图3)。而在然果村, 甘蒙柽柳径向生长与风速之间的显著负相关关系发生在生长季早期或是生长季前, 同时响应分析结果表明, 径向生长与各月平均风速主成分向量的回归函数并没有达到95%显著性水平, 然果村甘蒙柽柳径向生长对月平均风速不存在显著的响应过程(图4)。

主要研究区年表与环境因子的滑动相关的结果表明, 甘蒙柽柳生长与环境因子的相关关系在时间尺度上存在较为剧烈的变化。近50年, 该流域径流7月径流量与然果村甘蒙柽柳生长的相关性逐渐减弱, 与此同步, 8月份径流流量与甘蒙柽柳生长的相关逐渐增强(图5)。

图5 然果村甘蒙柽柳生长与径流流量的滑动相关。不同颜色所对应的数值表示不同的相关系数。Fig. 5 Moving correlation betweenTamarix austromongolica growth and monthly runoff. The values represented by the different colors show the different correlation coefficients.

然果村甘蒙柽柳生长与多个月份的月平均风速的负相关在近50年来显著减弱, 且在近年来出现反向特征。强风对生长季甘蒙柽柳生长的不利影响逐渐减弱(图6)。

图6 然果村甘蒙柽柳与风速的滑动相关。不同颜色所对应的数值表示不同的相关系数。Fig. 6 Moving correlation betweenTamarix austromongolica growth and monthly wind speed. The values represented by the different colors show the different correlation coefficients.

唐乃亥水文站记录在1971年、1988年, 2001年及2002年, 7-8月黄河径流总量出现极低值(低于4.0 × 109m3), 远小于1956-2007年7-8月黄河径流的平均值(6.27 × 109m3)。同期两地树轮指数平均值分别为0.71、0.91、0.93和0.89, 均为低生长年。

3 讨论
3.1 研究区乔木状甘蒙柽柳的最大年龄估计及多株合生现象

早期, 有研究者(Brotherson et al., 1984; Haigh, 1998)曾对柽柳属植被的年龄-胸径关系进行了探究, 初步建立了年龄-胸径关系方程, 认为在一定范围内随着年龄的增长, 柽柳胸径呈现线性增加趋势。然而, 我们在对研究区乔木状甘蒙柽柳的调查中发现, 树龄确切超过百年的甘蒙柽柳一株(倒木), 疑似超过百年的甘蒙柽柳一株, 早期人们对该地区甘蒙柽柳年龄的直观判断存在高估。两株甘蒙柽柳均位于然果村, 相对于班多村林区, 然果村甘蒙柽柳平均年龄(仅包含采样样本)较大, 存在较多年龄> 50年的植株。

柽柳属植物干形复杂, 多数研究着重于其幼体的生长发育, 而对成年个体年龄研究较少。本研究所涉及的甘蒙柽柳个体胸径较大, 年龄相对偏大, 与传统柽柳研究中的研究对象存在区别, 具有其独特之处。

最为典型的独特之处在于其多株合生现象的存在。多株合生现象在羊曲水电站库区及周边地区甘蒙柽柳林中十分普遍, 从而使得作为灌木或者小乔木的甘蒙柽柳能够具有高大乔木的外形特征。我们发现, 组成多株合生树的个体具有不同的年龄, 因此, 多株合生树的形成是由于不同时期生长出的植株由于分布的聚集, 在径向生长过程中逐渐发生挤压合生在一起。在研究区, 泥石流灾害频发, 日平均风速大, 成年植株通常具有较为明显的损伤。我们认为, 生长在成年植株周围或甚至生长在成年植株损伤而产生的空洞或裂隙中的幼苗可能更易在恶劣的环境中存活, 从而促使了多株合生现象的发生。

3.2 甘蒙柽柳的生长与生态环境因子的关系

由于立地条件、树木健康状况等多种因素的影响, 同一地区同一树种树木个体间的生长差异在树木学研究中普遍存在(Bowman et al., 2013)。本研究甘蒙柽柳林所在研究区位于高原河漫滩, 常年大风, 生存环境恶劣, 从而导致了局地微环境变动剧烈, 甘蒙柽柳个体之间的生长速率存在较大的差异。

水分对柽柳属植物分布及生长的影响是显著的。在我国, 蒋进和高海峰(1992)对中国主要柽柳属植物的耐旱性进行了研究, 证明了水分条件对柽柳属植物生长的显著影响。杨维康等(2002)对新疆柽柳属植物分布进行了较为全面的分析, 验证了水分条件是决定该地区柽柳属植物现代地理分布的关键环境因子。本研究的结果表明, 高原河漫滩的乔木状甘蒙柽柳与黄河季节性径流量间存在关联, 显著的响应关系发生在班多村甘蒙柽柳径向生长与生长当年7月径流之间。7-8月黄河径流流量的增加对甘蒙柽柳的生长具有积极的作用, 体现了夏季充足的水分条件在甘蒙柽柳生长过程中的重要性。类似的, 利用树木年轮学方法, 肖生春和肖洪浪(2006)分析了极端干旱区湖岸柽柳径向生长对水环境演变的响应。侯佳文等(2015)分析了博斯腾湖滨多枝柽柳(Tamarix ramosissima)生长对湖水水位的敏感性, 表明了柽柳径向生长与当年以及前一年的湖泊水位之间存在显著的相关性。本研究中, 甘蒙柽柳径向生长与径流之间的相关性仍然存在时间上的变化, 在然果村, 甘蒙柽柳径向生长与7月径流量正相关性减弱的同时, 与8月径流量正相关性增强, 一定程度上反映了甘蒙柽柳生长与径流显著相关的时间段在近50年发生推迟。

据当地水文站点及文献记载(Wang et al., 2006; Gou et al., 2007), 20世纪70年代黄河自然径流量显著较低, 与本研究所建立的甘蒙柽柳树轮宽度年表低值年频发时期相吻合; 进一步单年分析的结果同时表明了在黄河7-8月径流量出现低值的年份, 树轮径向生长量较低, 一定程度上证实了水文径流与河漫滩甘蒙柽柳生长的关系。由于水分条件的关键作用(郭京衡等, 2014), 使得我国柽柳属植被主要分布在干旱、半干旱地区河岸及湖滨沙地, 河流及湖泊水文条件通常被认为是影响该地区柽柳属植物生长的重要环境因子(侯佳文等, 2015; 朱绪超等, 2015); 同时, 沙漠腹地地区的柽柳也与水文条件有着密切的联系, 其通常分布在古河道或旧河床地区, 地下水位是影响柽柳生长的关键因素(袁国富等, 2012)。

柽柳是中国西北地区重要的防风阻沙植被, 在然果村, 甘蒙柽柳林的防风效应体现明显, 本研究发现, 研究区植被生长当年5-6月平均风速是影响甘蒙柽柳生长的最主要的环境限制因子, 突出表现在两地研究区甘蒙柽柳径向生长与风速的显著负相关关系。刘国军等(2008)对塔克拉玛干沙漠中的沙生柽柳(Tamarix taklamakanensis)的研究表明, 多风环境中沙生柽柳的根系形态与沙漠主导风向有着密切的关系。在本研究中, 班多村甘蒙柽柳生长显著负响应于5-6月风速, 然果村甘蒙柽柳生长与风速的关系较为复杂, 体现在近50年来负相关关系的显著减弱甚至反向。1971年甘蒙柽柳生长出现极低值, 结合区域甘蒙柽柳调查数据后判断林龄结构, 研究认为20世纪70年代以来存在大量新生甘蒙柽柳个体。从而认为, 逐渐成熟的甘蒙柽柳林对强风的抵抗性增强, 负相关性减弱。

4 结论

青海省同德县羊曲水电站库区及附近甘蒙柽柳林中存在极少量树龄超过百年的植株, 其较大的胸径以及乔木状形态的形成, 主要是由于多株合生现象的存在。该地区现存甘蒙柽柳多生于20世纪70-80年代, 20世纪80年代后, 甘蒙柽柳生长较为平稳。甘蒙柽柳径向生长主要响应于7-8月黄河径流以及生长季风速。7-8月径流量的增加将有利于甘蒙柽柳生长, 但水电站库区然果村甘蒙柽柳林对径流响应月份表现出逐渐推迟的过程。生长季的强风对未成熟的甘蒙柽柳的生长是不利的, 随着甘蒙柽柳林林龄的增大, 径向生长与风速的负相关性减弱。青海省高原河漫滩乔木状甘蒙柽柳多株合生现象的发现在植被特殊形态的研究方面具有理论意义, 其年龄的准确判断在对其生态价值的评估方面具有实践意义。甘蒙柽柳历史生长及其对环境、气候因子的响应关系的明确, 为其在气候变化下或是人为改变的环境下的生长评估提供了科学依据, 对高原干旱区、半干旱区生态可持续发展具有建设性的科学价值。

致谢 感谢国家林业局组织该项研究工作, 感谢青海省林业厅的野外协助以及数据支持。非常感谢中国科学院植物研究所张齐兵研究员对该研究工作的大力指导, 以及植物研究所树轮实验室工作人员张启对树轮样本的协助处理及测量。

The authors have declared that no competing interests exist.

作者声明没有竞争性利益冲突.

参考文献
[1] Babst F, Poulter B, Trouet V, Tan K, Neuwirth B, Wilson R, Carrer M, Grabner M, Tegel W, Levanic T (2013). Site- and species-specific responses of forest growth to climate across the European continent. Global Ecology and Biogeography, 22, 706-717. [本文引用:1]
[2] Biondi F, Waikul K (2004). DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree- ring chronologies. Computers & Geosciences, 30, 303-311. [本文引用:1]
[3] Birken AS, Cooper DJ (2006). Processes of Tamarix invasion and floodplain development along the lower Green River, Utah. Ecological Applications, 16, 1103-1120. [本文引用:1]
[4] Bowman DMJS, Brienen RJW, Gloor E, Phillips OL, Prior LD (2013). Detecting trends in tree growth: Not so simple. Trends in Plant Science, 18, 11-17. [本文引用:1]
[5] Brotherson JD, Carman JG, Szyska LA (1984). Stem-diameter age relationships of Tamarix ramosissima in central Utah. Journal of Range Management, 37, 362-364. [本文引用:2]
[6] China Weekly (2016). Life and death robbery of Qinghai ancient Tamarix forest. http: //www. fx361. com/page/2016/1217/400445. shtml. Cited: 2017-03-12. (in Chinese)
[中国周刊 (2016). 青海古柽柳林生死劫. http://www.fx361.com/page/2016/1217/400445.shtml. Cited: 2017-03-12. ] [本文引用:1]
[7] Fang OY, Wang Y, Shao XM (2016). The effect of climate on the net primary productivity (NPP) of Pinus koraiensis in the Changbai Mountains over the past 50 years. Trees, 30, 281-294. [本文引用:1]
[8] Friedman JM, Vincent KR, Shafroth PB (2005). Dating floodplain sediments using tree-ring response to burial. Earth Surface Processes and Land forms, 30, 1077-1091. [本文引用:1]
[9] Frumkin A (2009). Stable isotopes of a subfossil tamarix tree from the dead sea region, israel, and their implications for the intermediate bronze age environmental crisis. Quaternary Research, 71, 319-328. [本文引用:1]
[10] Gao LL, Gou XH, Deng Y, Yang MX, Zhang F (2017). Assessing the influences of tree species, elevation and climate on tree-ring growth in the Qilian Mountains of northwest China. Trees, 31, 393-404. [本文引用:1]
[11] Gou XH, Chen FH, Cook E, Jacoby G, Yang MX, Li JB (2007). Streamflow variations of the Yellow River over the past 593 years in western China reconstructed from tree rings. Water Resources Research, 43, W06434. doi: DOI:10.1029/2006WR005705. [本文引用:1]
[12] Gou XH, Deng Y, Gao LL, Chen FH, Cook E, Yang MX, Zhang F (2015). Millennium tree-ring reconstruction of drought variability in the eastern Qilian Mountains, northwest China. Climate Dynamics, 45, 1761-1770. [本文引用:1]
[13] Guo JH, Zeng FJ, Li CJ, Zhang B (2014). Root architecture and ecological adaptation strategies in three shelterbelt plant species in the southern Taklimakan Desert. Chinese Journal of Plant Ecology, 38, 36-44. (in Chinese with English abstract)
[郭京衡, 曾凡江, 李尝君, 张波 (2014). 塔克拉玛干沙漠南缘三种防护林植物根系构型及其生态适应策略. 植物生态学报, 38, 36-44. ] [本文引用:1]
[14] Haigh SL (1998). Stem diameter-age relationships of Tamarix ramosissima on lake shore and stream sites in southern Nevada. The Southwestern Naturalist, 43, 425-429. [本文引用:1]
[15] Holmes RL (1983). Computer-assisted quality control in tree- ring dating and measurement. Tree-Ring Bulletin, 43, 69-78. [本文引用:1]
[16] Horton JL, Clark JL (2001). Water table decline alters growth and survival of Salix gooddingii and Tamarix chinensis seedlings. Forest Ecology and Management, 140, 239-247. [本文引用:1]
[17] Hou JW, Haimit Y, Ye M (2015). The sensitivity of ring width of Tamarix ramosissina to lake level of the Bosten Lake. Journal of Desert Research, 35, 667-673. (in Chinese with English abstract)
[侯佳文, 海米提·依米提, 叶茂 (2015). 博斯腾湖滨柽柳(Tamarix ramosissina)生长对湖水位的敏感性. 中国沙漠, 35, 667-673. ] [本文引用:1]
[18] Jiang J, Gao HF (1992). A study on drought resistant ordination of Tamarix L. Arid Zone Research, 9(4), 41-45. (in Chinese with English abstract)
[蒋进, 高海峰 (1992). 柽柳属植物抗旱性排序研究. 干旱区研究, 9(4), 41-45. ] [本文引用:1]
[19] Li MQ, Shao XM, Yin ZY, Xu XG (2015). Tree-ring dating of the reshui-1 tomb in Dulan County, Qinghai Province, north-west China. PLOS ONE, 10, e0133438. doi: DOI:10.1371/journal.pone.0133438. [本文引用:1]
[20] Li YJ, Fang KY, Cao CF, Li DW, Zhou FF, Dong ZP, Zhang Y, Gan ZF (2016). A tree-ring chronology spanning 210 years in the coastal area of southeastern China, and its relationship with climate change. Climate Research, 67, 209-220. [本文引用:1]
[21] Lipp J, Trimborn P, Edwards T, Waisel Y, Yakir D (1996). Climatic effects on the δ18O and δ13C of cellulose in the desert tree Tamarix jordanis. Geochimica et Cosmochimica Acta, 60, 3305-3309. [本文引用:1]
[22] Liu GJ, Zhang XM, Li XR, Wei J, Shan LS (2008). The adaptation of the Tamarix root to the windy environment. Chinese Science Bulletin, 53, 147-150. (in Chinese)
[刘国军, 张希明, 李晓荣, 魏疆, 单立山 (2008). 塔克拉玛干柽柳的根对多风环境的适应. 科学通报, 53, 147-150. ] [本文引用:1]
[23] Norton D, Palmer J, Ogden J (1987). Dendroecological studies in new zealand 1. An evaluation of tree age estimates based on increment cores. New Zealand Journal of Botany, 25, 373-383. [本文引用:1]
[24] Qin L, Yuan YJ, Zhang RB, Wei WS, Yu SL, Fan Z, Chen F, Zhang TW, Shang HM (2016). Tree-ring response to snow cover and reconstruction of century annual maximum snow depth for northern Tianshan Mountains, China. Geochronometria, 43, 9-17. [本文引用:1]
[25] Rohner B, Weber P, Thürig E (2016). Bridging tree rings and forest inventories: How climate effects on spruce and beech growth aggregate over time. Forest Ecology and Management, 360, 159-169. [本文引用:1]
[26] Shao XM, Wang SZ, Xu Y, Zhu HF, Xu XG, Xiao YM (2007). A 3500-year master tree-ring dating chronology from the northeastern part of the Qaidam Basin. Quaternary Sciences, 27, 477-485. (in Chinese with English abstract)
[邵雪梅, 王树芝, 徐岩, 朱海峰, 许新国, 肖永民 (2007). 柴达木盆地东北部3500年树轮定年年表的初步建立. 第四纪研究, 27, 477-485. ] [本文引用:1]
[27] Wang H, Yang Z, Saito Y, Liu JP, Sun X (2006). Interannual and seasonal variation of the Huanghe (Yellow River) water discharge over the past 50 years: Connections to impacts from ENSO events and dams. Global and Planetary Change, 50, 212-225. [本文引用:1]
[28] Wang XC, Zhang MH, Ji Y, Li Zongshan, Li M, Zhang YD (2017). Temperature signals in tree-ring width and divergent growth of Korean pine response to recent climate warming in northeast Asia. Trees, 31, 415-427. [本文引用:1]
[29] Wang XM, Zhang CX, Zhang JW, Hua T, Lang LL, Zhang XY, Wang L (2010). Nebkha formation: Implications for reconstructing environmental changes over the past several centuries in the Ala Shan Plateau, China. Palaeogeography Palaeoclimatology Palaeoecology, 297, 697-706. [本文引用:1]
[30] Wu YH (2012). I am sad for the wild Tamarix in Qinghai. http: //blog. sciencenet. cn/home. php?mod=space&uid=38998&do=blog&id=495736. Cited: 2017-03-12. (in Chinese)
[吴玉虎 (2012). 我为青海野生古柽柳悲哀. http://blog.sciencenet.cn/home.php?mod=space&uid=38998&do=blog&id=495736. Cited: 2017. 3. 12. ] [本文引用:1]
[31] Xia JB, Zhao XM, Chen YP, Fang Y, Zhao ZG (2016). Responses of water and salt parameters to groundwater levels for soil columns planted with Tamarix chinensis. PLOS ONE, 11, e0145828. doi: DOI:10.1371/journal.pone.0145828. [本文引用:1]
[32] Xiao SC, Peng XM, Tian QY (2016). Climatic and human drivers of recent lake-level change in east Juyan Lake, China. Regional Environmental Change, 16, 1063-1073. [本文引用:1]
[33] Xiao SC, Xiao HL (2006). Response of radial growth of lakeshore tamarisk to changes in the water environment in extreme and arid regions. Journal of Beijing Forestry University, 28(2), 39-45. (in Chinese with English abstract)
[肖生春, 肖洪浪 (2006). 极端干旱区湖岸柽柳径向生长对水环境演变的响应. 北京林业大学学报, 28(2), 39-45. ] [本文引用:1]
[34] Xiao SC, Xiao HL, Peng XM, Tian QY (2014). Intra-annual stem diameter growth of Tamarix ramosissima and association with hydroclimatic factors in the lower reaches of China’s Heihe River. Journal of Arid Land , 6, 498-510. [本文引用:1]
[35] Xiao SC, Xiao HL, Si JH, Liu FM (2005). Growth characteristics of Tamarix ramosissima in arid regions of China. Acta Botanica Boreali-Occidentalia Sinica, 25, 1012-1016. (in Chinese with English abstract)
[肖生春, 肖洪浪, 司建华, 刘发民 (2005). 干旱区多枝柽柳的生长特性. 西北植物学报, 25, 1012-1016. ] [本文引用:1]
[36] Xiao SC, Xiao HL, Zhou MS, Si JH, Zhang XY (2004). Water level change of the West Juyan Lake in the past 10years recorded in the tree ring of the shrubs in the lake banks. Journal of Glaciology and Geocryology, 26, 557-562. (in Chinese with English abstract)
[肖生春, 肖洪浪, 周茂先, 司建华, 张小由 (2004). 近百年来西居延海湖泊水位变化的湖岸林树轮记录. 冰川冻土, 26, 557-562. ] [本文引用:1]
[37] Xinhua News Agency (2016). Wild Tamarix forest facing submerged in Tongde Qinghai, site protection or transplant? http: //www. chinanews. com/sh/2016/11-09/8058287. shtml. Cited: 2017-03-01. (in Chinese)
[新华社 (2016). 青海同德野生柽柳林面临被淹没, 原址保护还是移植? http://www.chinanews.com/sh/2016/11-09/8058287.shtml. Cited: 2017-03-01. ] [本文引用:1]
[38] Yang W, Spencer RJ, Krouse HR (1996). Stable sulfur isotope hydrogeochemical studies using desert shrubs and tree rings, Death Valley, California, USA. Geochimica et Cosmochimica Acta, 60, 3015-3022. [本文引用:1]
[39] Yang WK, Zhang DY, Yin LK, Zhang LY (2002). Distribution and cluster analysis on the similarity of the Tamarix communities in Xinjiang. Arid Zone Research, 19(3), 6-11. (in Chinese with English abstract)
[杨维康, 张道远, 尹林克, 张立运 (2002). 新疆柽柳属植物(Tamarix L. )的分布与群落相似性聚类分析. 干旱区研究, 19(3), 6-11. ] [本文引用:1]
[40] Yin Z-Y, Li MQ, Zhang Y, Shao XM (2016). Growth-climate relationships along an elevation gradient on a southeast- facing mountain slope in the semi-arid eastern Qaidam Basin, northeastern Tibetan Plateau. Trees, 30, 1095-1109. [本文引用:1]
[41] Yuan GF, Zhang P, Xue SS, Zhuang W (2012). Change characteristics in soil water content in root zone and evidence of root hydraulic lift in Tamarix ramossissima thickets on sand dunes. Chinese Journal of Plant Ecology, 36, 1033-1042. (in Chinese with English abstract)
[袁国富, 张佩, 薛沙沙, 庄伟 (2012). 沙丘多枝柽柳灌丛根层土壤含水量变化特征与根系水力提升证据. 植物生态学报, 36, 1033-1042. ] [本文引用:1]
[42] Zhang JC, Yao T, Liu CZ, Liu SZ, Sun T, Yuan HB, Tang JN, Ding F, Li XM, Liu R (2014). Climate environmental change and stable carbon isotopes in age layers of Tamarix sand -hillocks in Kumtag Desert. Acta Ecologica Sinica, 34, 943-952. (in Chinese with English abstract)
[张锦春, 姚拓, 刘长仲, 刘世增, 孙涛, 袁宏波, 唐进年, 丁峰, 李学敏, 刘瑞 (2014). 库姆塔格柽柳沙包年层稳定碳同位素与气候环境变化. 生态学报, 34, 943-952. ] [本文引用:1]
[43] Zhang QB, Evans MN, Lyu LX (2015). Moisture dipole over the Tibetan Plateau during the past five and a half centuries. Nature Communications, 6, 8062. [本文引用:1]
[44] Zhu XC, Yuan GF, Shao MA, Yi XB, Du T (2015). Spatial pattern of riparian vegetation in desert of the lower Tarim River Basin. Chinese Journal of Plant Ecology, 39, 1053-1061. (in Chinese with English abstract)
[朱绪超, 袁国富, 邵明安, 易小波, 杜涛 (2015). 塔里木河下游河岸带植被的空间结构特征. 植物生态学报, 39, 1053-1061. ] [本文引用:1]