环境科学  2019, Vol. 40 Issue (12): 5272-5285   PDF    
祁连山北坡中段降水稳定同位素特征及水汽来源分析
张百娟1,2, 李宗省2, 王昱1, 李永格2, 吕越敏2, 袁瑞丰2, 桂娟2     
1. 兰州理工大学能源与动力工程学院, 兰州 730000;
2. 中国科学院西北生态环境资源研究院(筹), 内陆河流域生态水文重点实验室, 甘肃省水文水资源工程技术中心, 兰州 730000
摘要: 依据祁连山北坡中段8个站点的降水样品,结合同期气象数据,从降水同位素特征、影响因素以及水汽来源等方面进行分析,结果表明:①研究区降水稳定同位素季节变化明显,表现为夏半年富集,冬半年贫化的特征;在空间尺度上δ18O值随海拔的升高而减小,年降水δ18O的海拔效应为-0.19‰/100 m;②各站点的局地大气降水线的斜率和截距表现为随海拔的增加而增加的趋势,表明2000 m以上的高海拔山区受到更强烈的局地再循环水汽的影响;③研究区降水中稳定同位素温度效应显著,δ18O的温度效应为0.64‰,且仅在夏季存在微弱的降水量效应;④研究区云下蒸发作用显著.在5、6、7和8月,研究区降水δ18O的平均雨滴蒸发率分别为23%、11%、12%和16%,云下蒸发富集率46%、27%、38%和32%;⑤在夏季连续降雨条件下,研究区降水的水汽来源主要为西风水汽,同时受到局地蒸发水汽的影响.本研究结果有助于进一步了解内陆河水文循环过程,为进一步开展干旱区同位素水文研究奠定基础.
关键词: 祁连山北坡      大气降水      稳定同位素      水汽来源      影响因素     
Characteristics of Stable Isotopes and Analysis of Water Vapor Sources of Precipitation at the Northern Slope of the Qilian Mountains
ZHANG Bai-juan1,2 , LI Zong-xing2 , WANG Yu1 , LI Yong-ge2 , Lü Yue-min2 , YUAN Rui-feng2 , GUI Juan2     
1. College of Energy and Power Engineering, Lanzhou University of Technology, Lanzhou 730000, China;
2. Key Laboratory of Ecohydrology of In-land River Basin/Gansu Hydrology and Water Resources Engineering Research Center, Northwest Institute of Eco-environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
Abstract: This study is based on precipitation samples from eight sites at the northern slope of the Qilian Mountains, combined with meteorological factors over the same period. Precipitation isotope characteristics, influence factors and the vapor sources of precipitation were analyzed, and the results show that:① The stable isotopes of precipitation in the study area show obvious seasonal changes, which are characterized by enrichment in the summer half-year and depletion in the winter half-year. The spatial precipitation δ18O value shows a significant downward trend with increasing altitude, and the altitude effect of the annual precipitation δ18O is -0.19‰/100 m, respectively; ② At all stations, the slope and intercept of local meteoric water lines show an increasing trend from low altitude to high altitude. The high-altitude mountains above 2000 m are affected by local water vapor recirculation; ③ The temperature effect is more significant and the temperature effect of δ18O is 0.64‰, and there is only a weak precipitation effect in summer; ④ The results indicate that sub-cloud evaporation has a great influence on the δ18O of precipitation; the average raindrop evaporation rate of δ18O is 23%, 11%, 12%, and 16%, and the δ18O composition has been enriched by 46%, 27%, 38%, and 32% in May, June, July, and August from cloud base to ground, respectively.⑤ Under the condition of continuous rainfall in summer, the vapor sources of precipitation mainly come from the west and are affected by local evaporation of water vapor. The study enhances knowledge of isotopic evolution of precipitation and provides a basis for further study of isotopic hydrology in arid regions.
Key words: Northern slope of Qilian Mountains      precipitation      stable isotopes      vapour sources of precipitation      influence factors     

降水是大气圈和水圈物质与能量交换过程中最积极的因子, 降水稳定同位素由氢氧同位素组成, 区别于其他溶于水的同位素, 它们本身就是水分子的组成部分, 是描述水循环过程的天然示踪剂, 因而在水循环和水文过程研究中具有重要的意义[1~4].

目前, 国内外学者对全球各地区降水同位素做了大量研究. Craig[3]指出降水中氢氧稳定同位素中存在显著的线性关系(δD=δ18O+10);Dansgaard[2]分析了降水稳定同位素的时空分布特征, 并提出同位素效应和过量氘等基本概念. Giustini等[5]的研究表明, 海拔和纬度是影响意大利降水同位素特征的主要因素;Li等[6]初步建立全球山地高海拔区降水中δ18O的垂直变化梯度为-0.40‰/100 m(5000 m以上).国内学者对降水同位素的研究主要集中在降水稳定同位素特征、地区水汽来源分析等方面[7~11]. Li等[12]综述了我国降水中稳定同位素的组成和分布, 阐明影响同位素比值的因素;刘忠方等[13]利用BW模型建立我国降水中δ18O与纬度和海拔定量关系的模型;Yao等[14, 15]的研究结果表明西宁和德令哈降水δ18O表现出温度效应;Tian等[16, 17]在青藏高原的聂拉木和拉萨地区研究降水δ18O仅在夏季季风活动时期偏负;同时已有研究表明在玉树、那曲、沱沱河、改则和狮泉河地区降水δ18O变化特征极其相似[17~20];汪少勇等[21]研究了长江源区降水氢氧稳定同位素特征及水汽来源;Yang等[22]分析了祁连山马粪沟流域降水和各水体稳定同位素组成的月尺度变化;Feng等[23]的研究表明温度是控制祁连山排露沟流域降水事件同位素变化的主要因素;袁瑞丰等[24]讨论祁连山石羊河流域降水同位素季节变化特征以及同位素效应;有学者系统地分析了祁连山中段黑河中上游地区降水氢氧同位素的变化特征以及与温度的相关关系、确定黑河流域大气降水水汽来源[25~29].

尽管前人对祁连山降水稳定同位素已经做了大量研究, 但研究区主要集中在中低海拔地区.由于高海拔地区观测数据及采样难度大, 使得其研究相对缺乏, 因此本文根据祁连山中段8个站点降水样品, 结合同期气象因子, 从降水同位素特征、局地大气降水线和云下蒸发等方面进行分析, 对祁连山高海拔地区降水同位素进一步研究, 深入了解干旱区降水同位素演化过程, 以期为进一步开展干旱区同位素水文研究奠定基础.

1 研究区概况

研究区位于祁连山北坡中段(图 1), 祁连山(93.5°~103.4°E, 35.5°~40°N, 海拔4000~5000 m)位于亚欧大陆中部, 总面积181689.1 km2, 是青藏高原东北部最大的边缘山系, 其地形地貌较为复杂.南北以柴达木、茶卡盆地及河西走廊为界, 东西以秦岭、六盘山及阿尔金山为界[30, 31].祁连山区远离海洋, 长期受西风气流控制, 具有典型的大陆性气候和高原气候的特征, 年均降水量300~700 mm左右, 是一座天然“高山水塔”, 同时也是我国河西内陆水系、青海内陆水系以及黄河上游水系的发源地[32, 33].祁连山北坡中段属于黑河流域发源地, 年均气温0.34℃, 年降水量约446 mm, 夏季降水量约占全年的86%, 同时降水量从东向西呈减少趋势, 并随海拔升高而增加[34~36].祁连山中段具有明显的垂直地带性, 土地利用类型主要有冰雪、草甸、草原和少量灌耕地等[37].

图 1 研究区位置及降水采样点分布示意 Fig. 1 Geographical position of study area and sampling in precipitation

2 材料与方法 2.1 样品采集

选取祁连山北坡中段8个站点为降水样品采集点(图 1), 采样点位置与降水样品明细见表 1.降水样品都是以次降水为单位进行采集, 为减小建筑物其他可能污染源对样品的影响, 须将样品收集器远离地面.每次采样结束, 立即将样品装在提前净化好的聚乙烯瓶中并用封口膜封住, 同时贴好标签存放至冰箱冷藏, 以防止蒸发, 并且由气象观测人员记录同期气象数据.

表 1 研究区采样点状况 Table 1 Sampling stations in study region

2.2 分析方法

降水样品在中国科学院寒区旱区环境与工程研究所内陆河流域生态水文重点实验室完成测试, 采用稳定同位素红外光谱技术的液态水稳定同位素分析仪(Model DLT-100, Los Gatos Research, Inc, MountainView, CA)测定降水中的δ18O和δD值, 18O/16O和D/H的测试精度分别为0.2‰和0.5‰.同位素测试的标准样选择为5个, 分别为:LGR1A(δ18O=-19.5‰, δD=-154.5‰)、LGR2A(δ18O=-15.5‰, δD=-116.5‰)、LGR3A(δ18O=-13.1‰, δD=-96.4‰)、LGR4A(δ18O=-7.7‰, δD=-51.1‰)和LGR5A(δ18O=-2.8‰, δD=-9.5‰), 测定数据的有效位数符合张自超等[40]对稳定同位素分析结果的有效位数规则, 即所有数据均保留一位有效数字.测定结果通过校准, 结果以Vienna“标准平均海水”(V-SMOW)的千分差来表示:

(1)

式中, RsampleRV-SMOW分别代表降水样与国际通用标准样中的氧或氢稳定同位素比率;R为同位素的比值(18O/16O和D/H).为了量化不同地区的大气降水线与全球大气降水线在斜率和截距上的偏移及比较这种差异, Dansgaard[2]在前人基础上首次定义了过量氘:d-excess=δD-8δ18O, 可用于判断降水来源地气候特征.

3 结果与分析 3.1 降水稳定同位素的时间变化

图 2所示, 研究区8个站点降水稳定同位素的季节变化明显, 均表现为夏半年富集, 冬半年贫化的特征, 与我国西南地区以及太湖流域氢氧同位素组成所表现的季节效应相反[41, 42], 与干旱区大多数站点季节变化过程相似[17].

图 2 研究区各站点降水同位素的时间变化 Fig. 2 Temporal variation of precipitation isotopes in the study area

采样期间, 研究区降水δ18O(δD)值的变化范围为从-29.7‰~8.3‰(-219.2‰~60.9‰), 年降雨量加权值为-6.1‰(-31.7‰), d-excess的年降雨量加权值为17.1‰, 变化范围为-23.5‰~40.4‰.降水稳定同位素值年际变化较大, δ18O值冬、夏相差高达30‰之多, d-excess值相差超过50‰.原因是冬、春季降水水汽主要由西风输送, 而且长的传输路径使水汽重同位素不断贫化, 而夏季、秋季稳定同位素主要受强烈蒸发影响, 对降水量较小的降水事件影响更为显著, 已有研究证实在黑河流域存在降水量小, 但是降水中重同位素含量高的现象[43, 44].如峨堡站降水中δ18O值为6.8‰(降水量为0.8 mm, 2013-09-16), 祁连站降水中δ18O值为8.3‰(降水量为1.3 mm, 2013-07-04).进一步研究发现, 研究区祁连站和野牛沟站降水δ18O和d-excess值在时间尺度上变化趋势一致, 其余6个站点降水δ18O和d-excess值则表现出反向变化趋势, 其原因主要是d-excess受水汽源地水体蒸发时周围空气相对湿度的影响[45], 陆地蒸发性水汽补充到原始水汽中使得降水中d-excess增加, 此类似现象也发生在托来河流域及古浪河流域[40, 46].

3.2 降水稳定同位素的空间变化

图 3(a)所示, 在区域尺度上降水中δ18O(δD)值随海拔的升高而呈明显的下降趋势, 相对较高值主要分布在走廊平原, 特别是在张掖地区, 而低值主要分布在祁连山区, 从海拔最低的张掖站到海拔最高的冰沟站, δ18O(δD)年降雨量加权值变化范围为-2.3‰(-4.4‰)~-7.9‰(-42.7‰).而d-excess值在区域尺度上与δ18O值相反, 较高值分布在祁连山区, 各站点的年均值随海拔的升高呈增加趋势[图 3(b)]. d-excess年降雨量加权值的变化范围为14.1‰~20.5‰, 造成这种趋势的原因可由以下两点解释:①张掖地区所在的走廊平原地区受到强烈的云下二次蒸发作用;②在祁连山区局地水汽再循环作用更为强烈, 反映出在空间上降水稳定同位素随海拔和局地气候状况的变化.进一步计算结果表明, 研究区年δ18O、δD和d-excess的海拔效应分别为-0.19‰/100 m、-0.72‰/100 m和0.77‰/100 m, Poage等[47]的研究表明除同位素递减率较高的极端纬度地区外, 全球降水中δ18O垂直变化梯度的平均值为-0.28‰/100 m(海拔5000 m以下地区);吴锦奎等[43]研究黑河流域δ18O的海拔梯度为-0.47‰/100 m, 而在乌鲁木齐河流域存在降水同位素正高程梯度0.12‰/100 m, 此现象是由水汽再循环和云下蒸发过程影响所致[48].

图 3 δ18O和d-excess的空间变化 Fig. 3 Spatial change of δ18O and d-excess

3.3 局地大气降水线(LMWL)

图 4(a)所示为研究区局地大气降水线(δD=7.82δ18O+13.77), 降水线斜率小于全球降水线的斜率, 但截距高于10‰, 表明研究区在云下蒸发作用的同时受到强烈局地蒸发影响, 导致降水线斜率偏小, Friedman等[49]首先提出雨滴受到蒸发作用使得大气降水线斜率小于8.进一步分析得到研究区固态降水及液态降水局地大气降水线存在明显的差异, 固态降水大气水线为δD=8.29δ18O+25.70, 液态降水大气水线为δD=7.84δ18O+12.78, 两者斜率和截距都相差较大, 造成如此大差异的原因是固态降水大多发生在冬、春季, 温度较低, 而液态降水发生在夏、秋季, 气温高且蒸发强烈, 固态降水除了受到干燥的环境和云下蒸发作用之外, 更多地受到局地再循环水汽的影响.研究区各站点的局地大气降水线:冰沟(δD=8.13δ18O+22.49)、祁连(δD=7.83δ18O+10.33)、民乐(δD=8.18δ18O+13.29)、张掖(δD=6.43δ18O+8.25)、峨堡(δD=8.16δ18O+15.48)、野牛沟(δD=8.12δ18O+19.21)、托勒站(δD=7.69δ18O+15.53)和嘉峪关站(δD=6.86δ18O-17.7).同一区域8个站点局地降水线斜率以及截距差异如此之大的原因是二次蒸发的强弱程度和局地水汽再循环作用不同.在祁连山高海拔区(冰沟、峨堡、祁连、野牛沟和托勒站)主要受到强烈的水汽再循环作用, 使得斜率和截距升高, 而在走廊平原区(包括民乐、张掖和嘉峪关站)较低的斜率主要是较为干燥的大气环境和较强的蒸发过程所导致.同时, 将研究区分为2000 m以上高海拔山区(冰沟、祁连、民乐、峨堡、野牛沟、托勒)以及2000 m以下走廊平原区(嘉峪关站和张掖站), 分析发现2000 m以上的高海拔山区各站点降水线方程截距和斜率都相对较高, 表明高海拔山区除了受到云下蒸发作用之外, 还受到局地再循环水汽以及水汽源地的影响;而2000 m以下走廊平原区(嘉峪关站和张掖站)降水线斜率和截距都相对较低, 可能受到云下二次蒸发的影响.

图 4 研究区局地大气降水线和各站点大气降水线斜率、截距与海拔的关系 Fig. 4 Local meteoric water line for study area and the relationship between the slope, intercept, and altitude of the local precipitation line

图 4(b)所示, 研究区各站点大气降水线的斜率和截距随着海拔的增加而增加, 反映出高海拔地区局地再循环水汽对降水稳定同位素浓度的影响较大, 而低海拔地区主要受云下蒸发过程的影响, 此现象在托来河流域[40]、古浪河流域[46]和乌鲁木齐河流域[50]也得到证实.进一步分析发现研究区降水线截距整体上呈现增加的趋势, 而在2000 m以下1500 m左右截距相差较大, 主要体现在张掖站以及嘉峪关站, 相较于张掖站, 嘉峪关站受到强烈云下蒸发过程的影响, 导致降水线截距较小, 而2250 m左右及以上, 各站点大气降水线斜率变化不大且都接近于8, 其原因主要是受到平衡分馏作用的影响.

3.4 降水稳定同位素与局地气象因子的相关关系

降水同位素值变化是一个蒸发和凝结的物理过程, 影响其分布的因素主要有区域环境背景和区域地理特征(温度、降水量、相对湿度等)[12].如图 5所示, 降水同位素与相对湿度、降水量没有明显的相关关系, 而表现出显著的温度效应, 因此本文着重分析了降水的温度效应以及降水量效应.

图 5 研究区各站点降水同位素与气象因子的相关关系 Fig. 5 Relationship between precipitation isotope and meteorological factors in the study area

3.4.1 降水稳定同位素与温度的关系

表 2为研究区降水稳定同位素数据与温度的关系, 以气温0℃和10℃为界, 据此可将所有数据划分为3类(其中张掖站在观测期所采数据都在10℃以上, 冰沟站都在10℃以下):气温0℃以下为第一组, 该组数据中所有站点降水中δ18O与日均温度呈明显的正相关, d-excess和日均温度展现出显著的负相关, 其原因是当气温在0℃以下时, 降水过程中主要受平衡分馏影响[51, 52], 而且δ18O和气温的相关系数随台站海拔的升高而增大, 反映了水汽随海拔抬升而发生的绝热冷却过程.当气温在0~10℃之间时, 除峨堡站、祁连、野牛沟站外, δ18O值随气温变化无明显变化趋势, 展现出稳定状态, 而d-excess值除峨堡站外都展现出稳定状态, 反映出气温上升引起的δ18O富集效应以及降水稳定同位素受到的云下蒸发作用被局地水汽再循环过程的影响所抵消.第三组数据δ18O随气温的升高而升高, 呈现出不明显的温度效应, 而d-excess则随着气温的上升显著降低, 这些降水事件主要发生在夏季, 受到了强烈的云下蒸发影响.此前, Pang等[52]在新疆的乌鲁木齐河流域也发现了类似的现象, 已有研究也认为这是强烈云下蒸发影响的明显表征[53].

表 2 研究区各站点不同温度范围内的温度效应 Table 2 Temperature effects in different temperature range of each site in the research area

研究区降水中稳定同位素温度效应显著, δ18O的温度效应为0.64‰, 即温度每升高1℃, δ18O增加0.64‰.峨堡站、野牛沟站、祁连站、民乐站、张掖站、托勒站、嘉峪关站和冰沟站所有次降水中δ18O表现出的温度效应依次为0.67‰、0.28‰、0.78‰、0.65‰、0.18‰、0.60‰、0.37‰和0.62‰, 除张掖站外降水中δ18O与日平均气温都具有明显的正相关关系, 同时随着海拔的降低, 降水稳定同位素表现的温度效应也减小.而在观测期内d-excess除了峨堡站表现出显著的温度效应(0.30‰)外, 其他各站点与日均温度相关关系不明显.这些研究结果表明, 降水稳定同位素演化主要受云下蒸发过程的影响, 而且局地水汽再循环过程引起的稳定同位素贫化效应被强烈的云下蒸发影响所主导, 同时反映出内陆干旱区不同海拔降水稳定同位素温度效应的差异性及其降水同位素演化过程的复杂性.

3.4.2 降水稳定同位素与降水量的关系

研究区降水量主要集中在夏季, 约占全年降水量的71%, 整个采样期间降水中δ18O不存在降水量效应, 因此本文重点分析研究区夏季降水稳定同位素的特征.如表 3所示, 各站点降水稳定同位素在夏季都表现出一定的降水量效应, 其中在峨堡站和托勒站8月、祁连站和民乐站7月的降水量效应较为显著, 而张掖站8月以及托勒站、嘉峪关站和冰沟站7月与降水量呈正相关关系, 不存在温度效应.造成此现象可能有以下两种原因:①降水水汽来源的不同;②在夏季较高的气温下, 张掖站受到更为强烈的蒸发作用导致同位素分馏, 使得降水量大但稳定同位素值仍较高.

表 3 研究区夏季降水稳定同位素的特征 Table 3 Characteristics of stable isotopes of summer precipitation in the study area

4 讨论 4.1 云下蒸发对降水同位素的影响

通常, 降水受到云下二次蒸发的影响将导致δ18O、δD富集以及d-excess贫化, 并且大气降水线的斜率和截距也会发生变化[54].在不同气象条件下降水稳定同位素值的变化以及大气降水线斜率和截距的变化可用来反映云下二次蒸发的存在[55~57].如表 4所示, 就降水时气温在0℃以上的降水事件而言, 随着气温升高, 降水稳定同位素值逐步偏正, 而δD-δ18O关系线的斜率和截距展现出明显的降低趋势. δ18O值从-8.6‰(当0 < T≤5℃)升高到-0.7‰(当T>20℃), 相应地δD-δ18O关系线的斜率从8.59降低到6.28, 其截距从15.6‰降低到5.7‰.这些事实表明, 当降水时的气温在0℃以上时, 云下蒸发对降水稳定同位素产生明显的富集作用.更为显著的是, 随着气温升高蒸发过程增强, 降低的斜率和截距印证了雨滴降落过程中云下蒸发不断加强的事实.需要强调的是, 云下蒸发对降雪事件或者0℃以下降水事件无明显影响.随着降水量的增加, δ18O值展现出一个降低趋势, 而δD-δ18O关系线的斜率和截距呈现出增加趋势(表 4).就较大降水量事件或连续降水事件而言, 大气饱和水汽压相对较高, 因此相对湿度也不断增大, 从而引起云下蒸发强度逐步减弱, 进而稳定同位素浓度富集程度较低.综上, 云下蒸发过程对降水稳定同位素浓度产生了明显影响, 特别是对小降水量降水事件的影响更大, 而且这一影响随着降水量的增加而逐步减弱.

表 4 研究时段内不同气温、降水和相对湿度下δD-δ18O关系的斜率、截距和R2变化 Table 4 Values of δ18O and δD-δ18O correlation for parameters for different air temperatures, precipitation, and relative humidity intervals

水汽从源地运移到研究区的过程中, 相对湿度对降水稳定同位素的影响最为明显, 当相对湿度是100%时, 稳定同位素平衡分馏发生, 其值遵循瑞利分馏规律[4].如表 4所示, 随着相对湿度降低, δ18O值从-5.5‰(70% < H)增加到-5.2‰(H≤50%), 而且δD-δ18O关系线中的斜率和截距明显降低:斜率从8.07降低到7.83, 截距从14.1‰降低到11.9‰, 再次反映较高的相对湿度下较弱的云下蒸发影响.

4.1.1 降落雨滴的云下蒸发率(f)

以上分析证实, 云下蒸发对降水稳定同位素浓度的影响是显著的.因此, 计算了降落雨滴的云下蒸发率(f), 即雨滴从云底到地面被蒸发的水量, 蒸发掉的部分用百分数来表示.依据之前的研究[58, 59], 祁连山和河西走廊地区雨滴的平均半径为0.4 mm.计算得研究区5、6、7和8月的平均雨滴蒸发率分别为23%、11%、12%和16%.就空间分布而言(图 6), 除6月份外, 平均雨滴蒸发率从西北到东南展现出明显的降低趋势, 较低值主要分布在高海拔的祁连山区, 而高值主要分布在走廊平原, 特别是嘉峪关地区.不同地形区的平均雨滴蒸发率值也印证了这一特征, 在祁连山高海拔区(冰沟、峨堡、祁连、野牛沟和托勒站)5、6、7和8月的平均雨滴蒸发率分别为16%、11%、9%和9%, 在走廊平原区(包括民乐、张掖和嘉峪关站)相应的值分别为35%、11%、17%和27%.

图 6 研究区降落雨滴的云下蒸发率的空间变化 Fig. 6 Spatial pattern on the evaporation rate of falling drops from cloud base to ground

4.1.2 δ18O云下蒸发富集率(E)

Froehlich等[53]探讨了云下蒸发过程对降水同位素的影响, 发现云下蒸发可以改变降水同位素达20%, 因此进一步评估了云下蒸发对降水稳定同位素浓度的影响, 并将其定义为δ18O云下蒸发富集率(E).计算表明, 云底降水的稳定同位素受云下蒸发影响, 在5、6、7和8月分别富集了46%、27%、38%和32%.从空间分布来看(图 7), 除6月份外, δ18O云下蒸发富集率(E)从祁连山到河西走廊方向展现出明显的降低趋势, 区域尺度上从西北到东南方向逐步降低.这一特征也可由不同地形区相应的均值反映出来, 在祁连山高海拔区(冰沟、峨堡、祁连、野牛沟和托勒站)5、6、7和8月的δ18O云下蒸发富集率(E)分别为30%、34%、33%和26%, 在走廊平原区(包括民乐、张掖和嘉峪关站)相应的值分别为72%、15%、46%和41%.就季节变化而言, 研究区内各站点δ18O云下蒸发富集率(E)在5月均为最高值, 6月较低, 而后6~8月其值波动增高, 而且走廊平原的站点其年内波动幅度显著大于祁连山区(图 7).

图 7 研究区δ18O云下蒸发富集率空间变化 Fig. 7 Spatial pattern of the relative influence of sub-cloud evaporation on the δ18O composition of precipitation

4.2 水汽来源对降水同位素的影响

本研究表明, 大气降水同位素值不仅与局地气象因子存在相关关系, 同时也与大气水气输送过程紧密联系.降水的水汽气团来源及其性质是决定降水同位素组成的重要因素[60], 有研究表明德令哈地区水汽来源主要为局地蒸发[61];陕甘宁地区夏半年主要有来自孟加拉湾的水汽、东南季风水汽和西风带水汽, 冬半年以西风带水汽为主[62];滇南蒙自地区干季降水的水汽主要来自于西风带输送及局地再蒸发水汽, 而湿季降水的水汽主要来源于远源海洋水汽的输送[63].通过前人大量的研究, 对我国西北干旱区外来水汽对降水的影响已有了清晰的认识[17, 64, 65], 降水的水汽来源主要为北极的极地气团、中亚的大陆气团、西太平洋赤道地区的海洋气团以及印度洋的海洋气团.在所有水汽源中, 祁连山地区水汽来源以西风水汽为主, 在夏季季风最强时, 在天气尺度下也可到达研究区, 并有部分局地再循环水汽补给作用[62].

冰沟站在整个观测期内d-excess值的年际波动非常大, 其值的变化范围为从-34.16‰~50.42‰, d-excess的年降雨量加权值为20.5‰, 说明降水水汽来源以大陆性局地水汽气团为主, 并且水源区气候干燥, 夏季d-excess的年降雨量加权值是20.4‰, 冬季为23.7‰.冬半年多为固态降水, 温度较低, 云下二次蒸发的影响十分微弱, 因此呈现出较高的d-excess值.因此本文以冰沟为例, 在夏季连续降雨条件下运用HYSPLIT模型模拟了1500、2500和5000 m高度水汽来源轨迹(图 8).在2014年6、7和8月连续降水条件下, 水汽来源主要是西风水汽.其中在6月26日, 1500、2500和5000 m高度下水汽来源分别为西南季风携带的海洋性水汽、西风水汽以及北极的极地气团, 并随着降水时间的持续, 水汽来源逐渐向北偏移, 在6月28日, 3个高度的水汽来源在中亚地区重合;在7月16、17和18日, 1500和2500 m高度下水汽来源都为局地蒸发水汽, 而5000 m高度下7月16和17日水汽来源为西风水汽, 7月18日变为西南季风水汽;8月6日, 1500 m高度水汽来源为东南季风水汽, 随着降水时间的持续, 在8月7日为北方大陆气团, 8月8日变为西风水汽;而2500和5000 m高度水汽来源在8月的6、7和8日都为西风水汽.综上, 在夏季连续降雨条件下, 冰沟站降水的水汽来源主要为西风水汽, 同时受到局地蒸发水汽的影响.

图 8 冰沟站夏季连续降水水汽来源后向轨迹图 Fig. 8 Back trajectories for continuous summer precipitation events at Binggou station

5 结论

(1) 研究区降水稳定同位素季节变化明显, 表现为夏半年富集, 冬半年贫化的特征;在时间尺度上除祁连站和野牛沟站, 其余站点降水δ18O和d-excess值均表现出反向变化趋势.在空间尺度上δ18O(δD)值随海拔的升高而减小, 年降水δ18O、δD和d-excess的海拔效应分别为-0.19‰/100 m、-0.72‰/100 m和0.77‰/100 m.

(2) 建立研究区局地大气降水线(δD=7.82δ18O+13.77)以及各个站点局地大气降水线, 其斜率和截距都随海拔的增加而增加, 同时2000 m以上的高海拔山区受到更强烈的局地再循环水汽的影响.

(3) 降水同位素温度效应显著且仅在夏季存在微弱的降水量效应, 全年降水δ18O的温度效应为0.64‰.峨堡站、野牛沟站、祁连站、民乐站、张掖站、托勒站、嘉峪关站和冰沟站降水中δ18O表现出的温度效应依次为0.67‰、0.28‰、0.78‰、0.65‰、0.18‰、0.60‰、0.37‰和0.62‰.

(4) 研究区云下蒸发对降水稳定同位素浓度的影响显著.在5、6、7和8月, 研究区降水δ18O的平均雨滴蒸发率分别为23%、11%、12%和16%, 云下蒸发富集率46%、27%、38%和32%.在空间尺度上, 除6月份外平均雨滴蒸发率以及云下蒸发富集率都展现为从西北到东南明显地降低趋势.

(5) 运用HYSPLIT模型模拟发现, 在夏季连续降雨条件下, 研究区降水的水汽来源主要为西风水汽, 同时受到局地蒸发水汽的影响.

参考文献
[1] 顾慰祖. 同位素水文学[M]. 北京: 科学出版社, 2011.
Gu W Z. Isotope hydrology[M]. Beijing: Science Press, 2011.
[2] Dansgaard W. Stable isotopes in precipitation[J]. Tellus, 1964, 16(4): 436-468. DOI:10.3402/tellusa.v16i4.8993
[3] Craig H. Isotopic variations in meteoric waters[J]. Science, 1961, 133(3465): 1702-1703. DOI:10.1126/science.133.3465.1702
[4] Craig H, Gordon L I. Deuterium and oxygen-18 variations in the ocean and the marine atmosphere[M]. Pisa: Laboratorio di geologianucleare, 1965.
[5] Giustini F, Brilli M, Patera A. Mapping oxygen stable isotopes of precipitation in Italy[J]. Journal of Hydrology:Regional Studies, 2016, 8: 162-181. DOI:10.1016/j.ejrh.2016.04.001
[6] Li Z, Yao T D, Tian L D, et al. Variations of δ18O in precipitation from the Muztagata Glacier, East Pamirs[J]. Science in China Series D:Earth Sciences, 2006, 49(1): 36-42. DOI:10.1007/s11430-004-5090-8
[7] Wu J K, Ding Y, Ye B, et al. Spatio-temporal variation of stable isotopes in precipitation in the Heihe River Basin, Northwestern China[J]. Environmental Earth Sciences, 2010, 61(6): 1123-1134. DOI:10.1007/s12665-009-0432-7
[8] 赵良菊, 尹力, 肖洪浪, 等. 黑河源区水汽来源及地表径流组成的稳定同位素证据[J]. 科学通报, 2011, 56(1): 58-67.
Zhao L J, Yin L, Xiao H L, et al. Isotopic evidence for the moisture origin and composition of surface runoff in the headwaters of the Heihe River basin[J]. Chinese Science Bulletin, 2011, 56(4-5): 406-416.
[9] 吴军年, 王红. 张掖大气降水的δ18O特征及水汽来源[J]. 安徽农业科学, 2011, 39(3): 1601-1604.
Wu J N, Wang H. Characteristics of the δ18O in precipitation and water vapor source in Zhangye[J]. Journal of Anhui Agricultural Sciences, 2011, 39(3): 1601-1604. DOI:10.3969/j.issn.0517-6611.2011.03.127
[10] 王宁练, 张世彪, 蒲健辰, 等. 黑河上游河水中δ18O季节变化特征及其影响因素研究[J]. 冰川冻土, 2008, 30(6): 914-920.
Wang N L, Zhang S B, Pu J C, et al. Seasonal variation of δ18O in river water in the upper reaches of Heihe River basin and its influence factors[J]. Journal of Glaciology and Geocryology, 2008, 30(6): 914-920.
[11] 王宁练, 张世彪, 贺建桥, 等. 祁连山中段黑河上游山区地表径流水资源主要形成区域的同位素示踪研究[J]. 科学通报, 2009, 54(15): 2148-2152.
Wang N L, Zhang S B, He J Q, et al. Tracing the major source area of the mountainous runoff generation of the Heihe River in northwest China using stable isotope technique[J]. Chinese Science Bulletin, 2009, 54(16): 2751-2757.
[12] Li Y J, Zhang M J, Wang S J, et al. Stable isotope in precipitation in China:a review[J]. Sciences in Cold and Arid Regions, 2012, 4(1): 83-90. DOI:10.3724/SP.J.1226.2012.00083
[13] 刘忠方, 田立德, 姚檀栋, 等. 中国大气降水中δ18O的空间分布[J]. 科学通报, 2009, 54(6): 804-811.
Liu Z F, Tian L D, Yao T D, et al. Spatial distribution of δ18O in precipitation over China[J]. Chinese Science Bulletin, 2009, 54(6): 804-811.
[14] Yao T D, Thompson L G, Jiao K Q, et al. Recent warming as recorded in the Qinghai-Tibetan cryosphere[J]. Annals of Glaciology, 1995, 21: 196-200. DOI:10.3189/S0260305500015810
[15] Yao T D, Thompson LG, Mosley-Thompson E, et al. Climatological significance of δ18O in North Tibetan ice cores[J]. Journal of Geophysical Research:Atmospheres, 1996, 101(D23): 29531-29537. DOI:10.1029/96JD02683
[16] Tian L, Yao T, Schuster P F, et al. Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau[J]. Journal of Geophysical Research:Atmospheres, 2003, 108(D9): 4293.
[17] Tian L D, Yao T D, MacClune K, et al. Stable isotopic variations in west China:a consideration of moisture sources[J]. Journal of Geophysical Research:Atmospheres, 2007, 112(D10): D10112.
[18] Yu W S, Yao T D, Tian L D, et al. Oxygen-18 isotopes in precipitation on the eastern Tibetan Plateau[J]. Annals of Glaciology, 2006, 43: 263-268. DOI:10.3189/172756406781812447
[19] 田立德, 姚檀栋, 沈永平, 等. 青藏高原那曲河流域降水及河流水体中氧稳定同位素研究[J]. 水科学进展, 2002, 13(2): 206-210.
Tian L D, Yao T D, Shen Y P, et al. Study on stable isotope in river water and precipitation in Naqu River basin, Tibetan Plateau[J]. Advances in Water Science, 2002, 13(2): 206-210. DOI:10.3321/j.issn:1001-6791.2002.02.013
[20] 田立德, 马凌龙, 余武生, 等. 青藏高原东部玉树降水中稳定同位素季节变化与水汽输送[J]. 中国科学D辑:地球科学, 2008, 51(8): 986-992.
Tian L D, Ma L L, Yu W S, et al. Seasonal variations of stable isotope in precipitation and moisture transport at Yushu, eastern Tibetan Plateau[J]. Science in China Series D:Earth Sciences, 2008, 51(8): 1121-1128.
[21] 汪少勇, 王巧丽, 吴锦奎, 等. 长江源区降水氢氧稳定同位素特征及水汽来源[J]. 环境科学, 2019, 40(6): 2615-2623.
Wang S Y, Wang Q L, Wu J K, et al. Characteristics of stable isotopes in precipitation and moisture sources in the headwaters of the Yangtze River[J]. Environmental Science, 2019, 40(6): 2615-2623.
[22] Yang Y G, Xiao H L, Qin Z D, et al. Hydrogen and oxygen isotopic records in monthly scales variations of hydrological characteristics in the different landscape zones of alpine cold regions[J]. Journal of Hydrology, 2013, 499: 124-131. DOI:10.1016/j.jhydrol.2013.06.025
[23] Feng F, Feng Q, Liu X D, et al. Stable isotopes in precipitation and atmospheric moisture of pailugou catchment in Northwestern China's Qilian Mountains[J]. Chinese Geographical Science, 2017, 27(1): 97-109. DOI:10.1007/s11769-017-0849-y
[24] 袁瑞丰, 李宗省, 蔡玉琴, 等. 干旱内陆河流域降水稳定同位素的时空特征及环境意义[J]. 环境科学, 2019, 40(5): 2122-2131.
Yuan R F, Li Z X, Cai Y Q, et al. Space-time characteristics and environmental significance of stable isotopes in precipitation at an arid inland river Basin[J]. Environmental Science, 2019, 40(5): 2122-2131.
[25] 张应华, 仵彦卿. 黑河流域中上游地区降水δ18O变化特征[J]. 冰川冻土, 2007, 29(3): 440-445.
Zhang Y H, Wu Y Q. Characteristics of the δ18O in precipitation in the upper and middle reaches of Heihe River[J]. Journal of Glaciology and Geocryology, 2007, 29(3): 440-445. DOI:10.3969/j.issn.1000-0240.2007.03.015
[26] 张应华, 仵彦卿. 黑河流域中上游地区降水中氢氧同位素与温度关系研究[J]. 干旱区地理, 2007, 30(1): 16-21.
Zhang Y H, Wu Y Q. Relation between oxygen and hydrogen isotopes in precipitation and temperature in Heihe River Basin, China[J]. Arid Land Geography, 2007, 30(1): 16-21. DOI:10.3321/j.issn:1000-6060.2007.01.003
[27] 张应华, 仵彦卿. 黑河流域大气降水水汽来源分析[J]. 干旱区地理, 2008, 31(3): 403-408.
Zhang Y H, Wu Y Q. Precipitation from different water vapor sources in the Heihe River Basin, China[J]. Arid Land Geography, 2008, 31(3): 403-408.
[28] Zhang Y H, Song X F, Wu Y Q. Use of oxygen-18 isotope to quantify flows in the upriver and middle reaches of the Heihe River, Northwestern China[J]. Environmental Geology, 2009, 58(3): 645-653. DOI:10.1007/s00254-008-1539-y
[29] 张应华, 仵彦卿. 黑河流域中游盆地地下水补给机理分析[J]. 中国沙漠, 2009, 29(2): 370-375.
Zhang Y H, Wu Y Q. Analysis of groundwater replenishment in the middle reaches of Heihe River[J]. Journal of Desert Research, 2009, 29(2): 370-375.
[30] 刘雪梅, 张明军, 王圣杰, 等. 2008-2014年祁连山区夏季降水的日变化特征及其影响因素[J]. 地理学报, 2016, 71(5): 754-767.
Liu X M, Zhang M J, Wang S J, et al. Diurnal variation of summer precipitation and its influencing factors of the Qilian Mountains during 2008-2014[J]. Acta Geographica Sinica, 2016, 71(5): 754-767.
[31] 王希强, 陈仁升, 刘俊峰. 气候变化背景下祁连山区负积温时空变化特征分析[J]. 高原气象, 2017, 36(5): 1267-1275.
Wang X Q, Chen R S, Liu J F. Spatial and temporal variation characteristics of accumulated negative temperature in Qilian Mountains under climate change[J]. Plateau Meteorology, 2017, 36(5): 1267-1275.
[32] 张耀宗, 张勃, 刘艳艳, 等. 近半个世纪以来祁连山区气温与降水变化的时空特征分析[J]. 干旱区资源与环境, 2009, 23(4): 125-130.
Zhang Y Z, Zhang B, Liu Y Y, et al. Variation characteristics of air temperature and precipitation in Mt. Qilian region in recent half century[J]. Journal of Arid Land Resources and Environment, 2009, 23(4): 125-130.
[33] 郑勤, 陈仁升, 韩春坛, 等. 祁连山TRwS204与中国标准雨量筒降水观测对比研究[J]. 高原气象, 2018, 37(3): 747-756.
Zheng Q, Chen R S, Han C T, et al. Intercomparison study on precipitation observations of TRwS204 and Chinese standard precipitation gauge in the Qilian Mountains[J]. Plateau Meteorology, 2018, 37(3): 747-756.
[34] Chen R S, Song Y X, Kang E S, et al. A cryosphere-hydrology observation system in a small alpine watershed in the Qilian Mountains of China and its meteorological gradient[J]. Arctic, Antarctic, and Alpine Research, 2014, 46(2): 505-523. DOI:10.1657/1938-4246-46.2.505
[35] Chen R, Liu J, Kang E, et al. Precipitation measurement intercomparison in the Qilian Mountains, Northeastern Tibetan Plateau[J]. The Cryosphere Discussions, 2015, 9(2): 2201-2230. DOI:10.5194/tcd-9-2201-2015
[36] 丁永健, 叶佰生, 周文娟. 黑河流域过去40a来降水时空分布特征[J]. 冰川冻土, 1999, 21(1): 42-48.
Ding Y J, Ye B S, Zhou W J. Temporal and spatial precipitation distribution in the Heihe Catchment, northwest China, during the past 40 a[J]. Journal of Glaciology and Geocryology, 1999, 21(1): 42-48.
[37] 王金叶, 常学向, 葛双兰, 等. 祁连山(北坡)水热状况与植被垂直分布[J]. 西北林学院学报, 2001, 16(S1): 1-3.
Wang J Y, Chang X X, Ge S L, et al. Vertical distribution of the vegetation and water and heat conditions of Qilian Mountain (northern slope)[J]. Journal of Northwest Forestry University, 2001, 16(S1): 1-3.
[38] Li Z X, Feng Q, Wang Q J, et al. The influence from the shrinking cryosphere and strengthening evopotranspiration on hydrologic process in a cold basin, Qilian Mountains[J]. Global and Planetary Change, 2016, 144: 119-128. DOI:10.1016/j.gloplacha.2016.06.017
[39] 李永格, 李宗省, 冯起, 等. 托来河流域不同海拔降水稳定同位素的环境意义[J]. 环境科学, 2018, 39(6): 2661-2672.
Li Y G, Li Z S, Feng Q, et al. Environmental significance of the stable isotopes in precipitation at different altitudes in the Tuolai River Basin[J]. Environmental Science, 2018, 39(6): 2661-2672.
[40] 张自超, 丁悌平. 关于同位素地质测试数据的数据处理及结果表示[J]. 岩矿测试, 2000, 19(1): 77-79.
Zhang Z C, Ding T P. Data processing and expression for analytical results in isotope geology[J]. Rock and Mineral Analysis, 2000, 19(1): 77-79. DOI:10.3969/j.issn.0254-5357.2000.01.015
[41] 张贵玲, 角媛梅, 何礼平, 等. 中国西南地区降水氢氧同位素研究进展与展望[J]. 冰川冻土, 2015, 37(4): 1094-1103.
Zhang G L, Jiao Y M, He L P, et al. Hydrogen and oxygen isotopes in precipitation in Southwest China:progress and prospects[J]. Journal of Glaciology and Geocryology, 2015, 37(4): 1094-1103.
[42] 胡勇博, 肖薇, 钱雨妃, 等. 水汽源地和局地蒸发对大气降水氢氧稳定同位素组分的影响[J]. 环境科学, 2019, 40(2): 573-581.
Hu Y B, Xiao W, Qian Y F, et al. Effects of water vapor source and local evaporation on the stable hydrogen and Oxygen isotopic compositions of precipitation[J]. Environmental Science, 2019, 40(2): 573-581.
[43] 吴锦奎, 杨淇越, 丁永建, 等. 黑河流域大气降水稳定同位素变化及模拟[J]. 环境科学, 2011, 32(7): 1857-1866.
Wu J K, Yang Q Y, Ding Y J, et al. Variations and simulation of stable isotopes in precipitation in the Heihe River basin[J]. Environmental Science, 2011, 32(7): 1857-1866.
[44] Zhou S Q, Nakawo M, Sakai A, et al. Water isotope variations in the snow pack and summer precipitation at July 1 Glacier, Qilian Mountains in northwest China[J]. Chinese Science Bulletin, 2007, 52(21): 2963-2972. DOI:10.1007/s11434-007-0401-z
[45] Merlivat L, Jouzel J. Global climatic interpretation of the deuterium-oxygen 18 relationship for precipitation[J]. Journal of Geophysical Research:Oceans, 1979, 84(C8): 5029-5033. DOI:10.1029/JC084iC08p05029
[46] 桂娟, 李宗省, 冯起, 等. 古浪河流域大气降水稳定同位素的时空特征及其环境意义[J]. 环境科学, 2019, 40(1): 149-156.
Gui J, Li Z X, Feng Q, et al. Space-time characteristics and environmental significance of the stable isotopes in precipitation in the Gulang River Basin[J]. Environmental Science, 2019, 40(1): 149-156.
[47] Poage M A, Chamberlain C P. Empirical relationships between elevation and the stable isotope composition of precipitation and surface waters:Considerations for studies of paleoelevation change[J]. American Journal of Science, 2001, 301(1): 1-15. DOI:10.2475/ajs.301.1.1
[48] Kong Y L, Pang Z H. A positive altitude gradient of isotopes in the precipitation over the Tianshan Mountains:effects of moisture recycling and sub-cloud evaporation[J]. Journal of Hydrology, 2016, 542: 222-230. DOI:10.1016/j.jhydrol.2016.09.007
[49] Friedman I, Machta L, Soller R. Water-vapor exchange between a water droplet and its environment[J]. Journal of Geophysical Research, 1962, 67(7): 2761-2766. DOI:10.1029/JZ067i007p02761
[50] 冯芳, 李忠勤, 金爽, 等. 天山乌鲁木齐河流域山区降水δ18O和δD特征及水汽来源分析[J]. 水科学进展, 2013, 24(5): 634-641.
Feng F, Li Z Q, Jin S, et al. Characteristics of δ18O and δD in precipitation and its water vapor sources in the upper Urumqi River basin, Eastern Tianshan[J]. Advances in Water Science, 2013, 24(5): 634-641.
[51] Clark I D, Fritz P. Environmental isotopes in hydrogeology[M]. Boca Raton: CRC Press, 1997.
[52] Pang Z H, Kong Y L, Froehlich K, et al. Processes affecting isotopes in precipitation of an arid region[J]. Tellus B:Chemical and Physical Meteorology, 2011, 63(3): 352-359. DOI:10.1111/j.1600-0889.2011.00532.x
[53] Froehlich K, Kralik M, Papesch W, et al. Deuterium excess in precipitation of Alpine regionsmoisture recycling[J]. Isotopes in Environmental and Health Studies, 2008, 44(1): 61-70. DOI:10.1080/10256010801887208
[54] Zhang M J, Wang S J. A review of precipitation isotope studies in China:basic pattern and hydrological process[J]. Journal of Geographical Sciences, 2016, 26(7): 921-938. DOI:10.1007/s11442-016-1307-y
[55] 孟玉川, 刘国东. 长江流域降水稳定同位素的云下二次蒸发效应[J]. 水科学进展, 2010, 21(3): 327-334.
Meng Y C, Liu G D. Effect of below-cloud secondary evaporation on the stable isotopes in precipitation over the Yangtze River basin[J]. Advances in Water Science, 2010, 21(3): 327-334.
[56] Wu H W, Zhang X P, Li X Y, et al. Seasonal variations of deuterium and oxygen-18 isotopes and their response to moisture source for precipitation events in the subtropical monsoon region[J]. Hydrological Processes, 2015, 29(1): 90-102. DOI:10.1002/hyp.10132
[57] 赵诗坤, 庞朔光, 文蓉, 等. 海河流域降水稳定同位素的云底二次蒸发效应[J]. 地理科学进展, 2015, 34(8): 1031-1038.
Zhao S K, Pang S G, Wen R, et al. Influence of below-cloud secondary evaporation on stable isotope composition in precipitation in the Haihe River Basin, China[J]. Progress in Geography, 2015, 34(8): 1031-1038.
[58] Kong Y L, Pang Z H, Froehlich K. Quantifying recycled moisture fraction in precipitation of an arid region using deuterium excess[J]. Tellus B:Chemical and Physical Meteorology, 2013, 65(1): 19251. DOI:10.3402/tellusb.v65i0.19251
[59] Li Z X, Feng Q, Yong S, et al. Stable isotope composition of precipitation in the south and north slopes of Wushaoling Mountain, northwestern China[J]. Atmospheric Research, 2016, 182: 87-101. DOI:10.1016/j.atmosres.2016.07.023
[60] 侯典炯, 秦翔, 吴锦奎, 等. 乌鲁木齐大气降水稳定同位素与水汽来源关系研究[J]. 干旱区资源与环境, 2011, 25(10): 136-142.
Hou D J, Qin X, Wu J K, et al. Characteristics of stable isotopes in precipitation and the water vapor sources in Urumqi[J]. Journal of Arid Land Resources and Environment, 2011, 25(10): 136-142.
[61] 朱建佳, 陈辉, 巩国丽. 柴达木盆地东部降水氢氧同位素特征与水汽来源[J]. 环境科学, 2015, 36(8): 2784-2790.
Zhu J J, Chen H, Gong G L, et al. Hydrogen and oxygen isotopic compositions of precipitation and its water vapor sources in Eastern Qaidam Basin[J]. Environmental Science, 2015, 36(8): 2784-2790.
[62] 刘洁遥, 张福平, 冯起, 等. 陕甘宁地区降水稳定同位素特征及水汽来源[J]. 应用生态学报, 2019, 30(7): 2191-2200.
Liu J Y, Zhang F P, Feng Q, et al. Characteristics of stable isotopes in precipitation over Shaanxi-Gansu-Ningxia and its water vapor sources[J]. Chinese Journal of Applied Ecology, 2019, 30(7): 2191-2200.
[63] 李广, 章新平, 许有鹏, 等. 滇南蒙自地区降水稳定同位素特征及其水汽来源[J]. 环境科学, 2016, 37(4): 1313-1320.
Li G, Zhang X P, Xu Y P, et al. Characteristics of stable isotopes in precipitation and their moisture sources in Mengzi Region, Southern Yunnan[J]. Environmental Science, 2016, 37(4): 1313-1320.
[64] Aizen V, Aizen E, Melack J, et al. Isotopic measurements of precipitation on central Asian glaciers (Southeastern Tibet, northern Himalayas, central Tien Shan)[J]. Journal of Geophysical Research:Atmospheres, 1996, 101(D4): 9185-9196. DOI:10.1029/96JD00061
[65] Li Z X, Gao Y, Wang Y M, et al. Can monsoon moisture arrive in the Qilian Mountains in summer?[J]. Quaternary International, 2015, 358: 113-125. DOI:10.1016/j.quaint.2014.08.046