Characteristics of δD and δ¹⁸O in Precipitation and Water Vapor Sources in Zigui Section of the Three Gorges Reservoir

三峡库区秭归段（106°50′～110°50′ E，29°16′～31°25′ N）位于湖北省西部，长江两岸，三峡工程坝上库首。属亚热带季风气候区，年均气温14.0～18.8 ℃，年均降水量1 000～1 200 mm，干季（11月—次年4月）湿季（5—10月）分明，湿季降水量占全年降水量的70%左右。境内山脉起伏，垂直气候变化明显。库区属于我国暴雨中心之一，强降水过程导致水土流失严重，并伴有滑坡、泥石流等自然灾害发生。研究区植物资源丰富，森林覆盖率高，植被以人工林为主，其中马尾松（Pinus massoniana Lamb.）林分布在库区酸性或者中性土壤上，柏木（Cupressus funebris Endl.）林是钙质土壤上的主要森林植被类型，这两种森林植被类型在库区森林中面积最大，伴生有木姜子（Litsea cubeba (Lour.) Pers.）、漆树（Toxicodendron vernicifluum (Franch.) F. A. Barkl.）等植物，林下植被主要有胡枝子（Lespedeza bicolor Turcz.）、狗脊（Woodwardia japonica (L. f.) Sm.）等。

2016年9月至2019年8月，在湖北秭归长江三峡库区森林生态系统国家定位观测研究站（110°18′ E，30°37′ N）和秭归县九岭头林场（110°48′ E，30°58′ N）的空旷地各放置1个雨量筒，在每个雨量筒上放置1个漏斗，漏斗内放置乒乓球以防止水分蒸发，分别采集大气降水样品，并将两个地点每天采集的降水样品混合（即日降水样品）。每天上午7:00—8:00用采样瓶收集大气降水，并用Parafilm封口膜密封。采集的样品置于−5～0 ℃冰箱低温保存，以备测试。降水量、温度等气象数据由秭归长江三峡库区森林生态站提供和野外实测获得，即在各采样点的空旷地放置自制雨量计和温度计，每天早7:00和晚19:00，记录降水量和温度等气象数据。

大气降水氢氧同位素（δD、δ¹⁸O）由清华大学稳定同位素实验室的MAT 253同位素比率质谱仪和Flash 2000 HT元素分析仪测定（δD测定精度为 ± 1‰，δ¹⁸O测定精度为 ± 0.2‰）。测定结果以相对于维也纳标准平均海水（VSOMW）的千分差（‰），即：

式中：δ为氢（氧）同位素的比值，R_s为样品中²H / ¹H（或¹⁸O / ¹⁶O），R_d为VSOMW中²H / ¹H（或¹⁸O / ¹⁶O）。

大气气团传输途径运用拉格朗日积分轨迹模型（HYSPLIT）的后向轨迹法进行模拟，选取垂直方向500、1 000、1 500 m作为模拟的初始高度，以研究区作为终点，模拟每次降水前144 h内水汽运动轨迹。

采用SPSS 24.0对数据进行处理，用相关分析确定大气降水δD（δ¹⁸O）与降水量、温度的相关关系，线性回归分析拟合大气降水δD与δ¹⁸O关系式。用Origin 9.1软件制图。

由2016年9月至2019年8月三峡库区气象资料可知（图1），该库区秭归段总降水量为2 484.9 mm，其中湿季（5—10月）降水量为1 734.3 mm，约占全年降水总量的70%，干季（11月—翌年4月）降水量为757.5 mm。气温变化同降水量变化趋势一致，表现为夏季高温多雨、冬季寒冷干燥的气候特征。

Figure 1.  The monthly variation of precipitation and temperature in Zigui section of the Three Gorges Reservoir Area

三峡库区秭归段大气降水δD的变化范围为−137.34‰～−3.66‰，均值为−51.81‰；δ¹⁸O的变化范围为−18.38‰～−2.90‰，均值为−8.44‰（图2，表1）。可见，该库区秭归段大气降水δD（δ¹⁸O）变化范围均介于全球^[10]与中国^[21]大气降水氢氧同位素变化范围内，均值低于全球及中国大气降水氢氧同位素均值。

Figure 2.  The monthly variation of precipitation δD and δ¹⁸O in Zigui section of the Three Gorges Reservoir Area

湿季大气降水的δD和δ¹⁸O变化范围分别为−137.34‰～−3.66‰和−18.38‰～−3.24‰，均值分别为−64.90‰和−9.74‰；干季大气降水的δD和δ¹⁸O变化范围分别为−131.95‰～−4.10‰和−18.17‰～−2.90‰，均值分别为−30.61‰和−6.33‰（表1）。由此可知，三峡库区秭归段干季大气降水δD和δ¹⁸O富集，而湿季大气降水δD和δ¹⁸O贫化。

根据三峡库区秭归段3个水文年大气降水δD（δ¹⁸O）实测值，采用最小二乘法拟合该地区大气降水δD-δ¹⁸O的关系式为：δD =（8.52 ± 0.12）δ¹⁸O +（20.11 ± 1.14）（R² = 0.96，p < 0.01，n = 186）（图3）。

Figure 3.  Relationship between δD and δ¹⁸O of precipitation in Zigui section of the Three Gorges Reservoir Area

研究区干季大气降水δD-δ¹⁸O关系式为：δD =（8.43 ± 0.34）δ¹⁸O +（22.76 ± 2.29）（R² = 0.90，p < 0.01，n = 71）；湿季大气降水δD-δ¹⁸O关系式为：δD =（8.09 ± 0.11）δ¹⁸O +（13.95 ± 0.19）（R² = 0.98，p < 0.01，n = 115）。湿季大气降水线的斜率和截距均小于干季（图3）。

三峡库区秭归段大气降水过量氘（d-excess）变化范围为1.52‰～37.76‰，均值为15.72‰，明显高于全球平均过量氘值（10‰）。研究区过量氘（d）季节变化显著，干季过量氘（d）均值为20.06‰，湿季过量氘均值为13.05‰，湿季过量氘值明显低于干季（图4，表1）。降水过量氘（d）的这种季节变化是季风区气候的特点之一。

Figure 4.  Monthly variation of d -excess of precipitation in Zigui section of the Three Gorges Reservoir Area

将三峡库区秭归段大气降水中δD（δ¹⁸O）分别与降水量（P）建立线性回归方程，得到P-δD的关系为：δD =（−0.88 ± 0.24）P +（−38.41 ± 4.47）（R² = 0.05，p < 0.01，n = 186）；P-δ¹⁸O的关系为：δ¹⁸O =（−0.09 ± 0.02）P +（−7.04 ± 0.52）（R² = 0.04，p < 0.05，n = 186）。由表2可知，在年尺度上，三峡库区秭归段大气降水δD（δ¹⁸O）与降水量（P）呈显著负相关关系，表现出降水量效应；在季节尺度上，湿季降水量效应显著，而干季降水量效应不显著。

将秭归段大气降水中δD（δ¹⁸O）分别与温度（T）建立线性回归方程，得到T-δD的关系为：δD =（−1.80 ± 0.28）T +（−24.23 ± 4.87）（R² = 0.18，p < 0.01，n = 186）；T-δ¹⁸O的关系为：δ¹⁸O =（−0.17 ± 0.03）T +（−5.84 ± 0.58）（R² = 0.12，p < 0.01，n = 186）。由表2可知，在年尺度上，大气降水δD（δ¹⁸O）与温度（T）呈显著负相关关系，表现出反温度效应；湿季和干季δD、δ¹⁸O与温度相关性不显著。

为进一步探究三峡库区秭归段大气降水水汽来源，选取该地区3个水文年内干季和湿季共6场典型降水事件（干季和湿季各选择3场不同量级降水事件），利用HYSPLIT水气团轨迹模型对降水气团的运移路径进行模拟。HYSPLIT模拟结果显示（图5），湿季大气降水气团主要来自于西太平洋、印度洋及我国南海，干季大气降水气团主要来自于亚欧大陆内部及局地的水汽蒸发。

Figure 5.  Trajectory map of air mass in typical rain events in Zigui section of the Three Gorges Reservoir Area

三峡库区秭归段大气降水δD与δ¹⁸O的线性关系及长江流域季风区多个城市（长沙^[8]、武汉^[22]、重庆^[23]、南京^[24]）大气降水线的斜率和截距大于中国大气降水线（δD = 7.9 δ¹⁸O + 8.2）^[21]，反映出这些地区大气降水水汽输送类型相似，且在一定程度上受到海洋水汽的影响。同时，本研究得到的三峡库区秭归段大气降水δD和δ¹⁸O的关系式的斜率和截距略大于长江流域季风区多个城市^[22-24]，产生这种现象的原因可能有以下两方面：（1）该研究采样点位于湖北省秭归县三峡库区，受该库区长江水分的蒸发和长江两岸较高森林覆盖率的影响，空气相对湿度较大，雨滴降落过程中受云下二次蒸发效应影响较小，斜率和截距偏大^[25]；（2）三峡库区秭归段地理位置及地形差异导致斜率和截距与其他城市相比存在一定程度的偏移，研究区海拔变化较大，区内山区、丘陵相间，极易形成局地性降水^[26]。

大气降水中的氢氧同位素与其蒸发过程条件、分馏凝结过程、水汽传输形式与局地降水量、温度等气象要素密切相关^[27]。我国降水氢氧同位素组成主要受降水量效应和温度效应控制^[28]。在年尺度上，三峡库区秭归段降水量效应显著；在季节尺度上，只有湿季有降水量效应，这可能与秭归段全年大部分降水均发生在湿季有关。研究区降水量效应与同处在我国长江流域季风区的武汉^[22]、重庆^[23]等地研究结果一致。这可能由以下两方面原因造成：首先，库区秭归段湿季降水充沛，空气饱和差很容易补偿，雨水下落过程中不易蒸发，出现较低的δ值，从而引起显著的降水量效应。其次，中纬度季风区降水中氢氧同位素具有明显的降水量效应^[29]，而三峡库区秭归段、武汉及重庆等均属于中纬度季风区。

温度也是影响降水氢氧同位素组成的主要因素之一^[30-31]。在年尺度上，研究区大气降水与温度呈负相关关系，即反温度效应。这种现象在长江流域^[18]等地均有报道。已有研究表明，反温度效应主要与研究区所处纬度有关^[11]。研究区位于中低纬度季风区，湿热的气候使雨滴降落过程中蒸发分馏不显著，长江蒸发水汽与地面蒸发水汽增加了降水量，在某种程度上，降水量效应抑制和掩盖了温度效应，使其呈现出反温度效应。

大气降水中过量氘（d）与水汽源地的气象条件密切相关，因此，d值可以作为指示地区大气降水水汽来源有效指标^[2]。以往的研究表明，水汽源地相对湿度是过量氘（d）季节变化的主要驱动因素^[32]。研究区过量氘存在明显的季节变化，干季d均值（19.47‰）显著高于湿季d均值（12.76‰），这与李维杰^[33]在我国西南地区的贵阳、重庆、成都等地研究结果一致，符合我国季风影响区域的降水特点。这可能由于库区湿季期间降水主要来源于海洋水汽凝结，空气湿度大，因此d值较小；而干季降水主要来自于内陆水汽蒸发，空气湿度小，水汽源区在不平衡条件下经历了快速蒸发过程，因此d值较大。此外，d值的变化规律与HYSPLIT模拟的水汽来源较为一致，可见，地区大气降水氢氧同位素特征能为识别其水汽来源提供有效信息。

HYSPLIT模拟结果显示，三峡库区秭归段干季和湿季水汽来源存在明显差异。在干季，该库区降水主要受大陆性季风控制，降水气团主要来自高纬度西风带输送的亚欧大陆内部及局地水汽的蒸发补给，降水量小，空气湿度较低，因此，降水δD和δ¹⁸O富集。而在湿季，随着东亚季风和南亚季风于5月中旬在我国沿海地区登陆，大量湿润、蒸发较弱的海洋气流向内陆推进，降水量明显增加，大气降水δD和δ¹⁸O表现出明显的下降趋势；远距离海洋水汽输送到三峡库区过程中，δD和δ¹⁸O不断冷凝而优先从水汽中分离出来，致使δD和δ¹⁸O逐渐贫化。由此可见，本研究区大气降水δD（δ¹⁸O）的周期性季节变化与水汽来源密切相关。此结果与同处长江流域季风区的武汉^[22]、重庆^[23]、长沙^[8]等地结果一致。

（1）三峡库区秭归段大气降水氢氧同位素关系式为：δD = 8.52 δ¹⁸O + 20.11（R² = 0.96，p < 0.01，n = 186），反映出三峡库区秭归段大气降水在一定程度上受到海洋水汽的影响。

（2）三峡库区秭归段大气降水过量氘（d）具有明显的季节变化，主要与水汽来源季节性差异密切相关。湿季水汽主要受东南季风和西南季风所携带海洋水汽的影响，过量氘（d）值较低；干季大气降水气团主要来自于亚欧大陆内部及局地的水汽蒸发，过量氘（d）值较高。

（3）不同时间尺度下，影响三峡库区秭归段降水氢氧同位素变化的环境因子存在差异。年尺度上，主要受“降水量效应”和“反温度效应”影响。在季节尺度上，降水量是湿季降水氢氧同位素值变化的主要影响因素。