黑龙江省(121.18°E~135.08°E，43.42°N~53.55°N)位于我国东北部，北隔黑龙江、东隔乌苏里江与俄罗斯相望，西临内蒙古，南邻吉林，地形以平原、丘陵为主(图 1)；属温带大陆性季风气候，夏季凉爽短促，冬季寒冷而漫长；年降水量400~650 mm，年均温-5 ℃~5 ℃，无霜期100~150 d，研究区的降水形式主要为降雨和降雪，夏季降水量约占全年降水量的60%，集中降雪期为每年11月至次年1月(图 2).黑龙江地区位于我国东部季风区的最北缘，纬度较高，不同季节水汽来源复杂，不仅受到来自太平洋的暖湿水汽，还受到西风带的水汽来源以及内陆水汽再蒸发.因此探究该地区大气降水δ¹⁸O的变化特征及水汽来源，对于中高纬度季风区大气环流和水循环机制的理解具有重要意义.

根据GNIP黑龙江地区的2个站点数据：哈尔滨站(监测时段1986-1998年，样品数量30个，采样单位：石家庄水文与工程地质研究所)、齐齐哈尔站(监测时段1988-1992年，样品数量50个，采样单位：未知)，由于部分样品数据缺失，可用的完整数据样品有80个.所有监测数据记录均为月平均值，根据月平均气温差异，研究区所有样品分为降雨样品和降雪样品，降雪样品采集过程中在室温下将雪融化再收集起来.研究区样品的降水形式包括降雨、降雪、雨夹雪等.其中，齐齐哈尔站有记录降水样品类型，降雨样品32个，雨夹雪样品2个，降雪样品16个；哈尔滨站降水样品类型的记录缺失，本文根据月平均气温分别为-4.6 ℃和-13.8 ℃的两个月为降雪样品，判断哈尔滨站有降雨样品28个，降雪样品2个.

本文采用的气象资料为美国国家环境预报中心(NCEP/NCAR)提供的逐月再分析格点资料，空间分辨率为水平2.5×2.5网格点，下载地址为http://www.esrl.noaa.gov.用到的数据资料有：地面到300 hPa高度场经向风速vwnd、纬向风速uwnd、比湿shum、通量气压pres.sfc. GrADS气象软件建立模型模拟水汽通量，追踪降水气团输送影响降水中δ¹⁸O变化. GrADS使用过程大致包括以下环节：下载数据资料文件及前期处理；编写及检查数据描述文件(*.ctl)；设定编写脚本文件(*.gs)；进行图形绘制及展示；保存描述文件及图像.

降水中稳定同位素比值(¹⁸O/¹⁶O或²H/¹H)用相对于维也纳标准平均海洋水(V-SMOW)的千分差来表示：

式中R_水样和R_V-SMOW分别代表水样中氧或氢的稳定同位素比值和V-SMOW中的稳定同位素比值.降水中稳定同位素比值的加权平均值计算公式为

式中，δ为加权平均值，P_i为降水量，δ_i为其相应的稳定同位素值.

大气降水中δ¹⁸O与δD之间的关系被定义为全球大气降水线(GMWL)：δD=8δ¹⁸O+10，其斜率和截距受到水汽来源、输送方式、水汽凝结温度以及降水期间的温度、湿度等的影响^[33]，使得局地大气降水线(LMWL)的斜率和截距有所偏离.根据GNIP在研究区两站点监测的稳定同位素数据，计算得出两站点大气降水线方程，齐齐哈尔站点：δD=7.58δ¹⁸O-0.14，R=0.98，n=50，哈尔滨站点：δD=5.52δ¹⁸O-19.4，R=0.83，n=30如图 3(a).综合所有样品数据得到的研究区域大气降水方程为δD=7.49δ¹⁸O-0.86，R=0.97，n=80如图 3(b).在大气降水线方程中，斜率a表现出较低值，而斜率值也较低，主要是由于哈尔滨站点较齐齐哈尔纬度偏低，年平均气温高，且位置更靠近海洋，在夏半年季风影响下有更多湿润水汽输送.

黑龙江地区大气降水线方程δD=7.49δ¹⁸O-0.86与柳鉴容^[21]通过GNIP和CHNIP的928个样品得出的中国东部季风区大气降水线方程δD=7.48δ¹⁸O+ 1.01相比，斜率和截距十分接近.并与CHNIP的三江站δD=7.29δ¹⁸O-6.71，海伦站δD=7.71δ¹⁸O+ 2.58对比^[29]，斜率非常相近.首先，降水线方程的相关系数能够说明大气降水中δD、δ¹⁸O具有很好的相关性；其次，该方程的斜率和截距与全国大气降水线接近，这说明东部季风区大气降水过程中的水汽源地基本都来自于海洋蒸发.黑龙江地区位于东部季风区内纬度较高地区，仅考虑海表面水的蒸发和凝结物降落时2个非平衡蒸发过程对δD和δ¹⁸O的影响时，由于凝结物降落过程中的蒸发速率与温度成正比，δD和δ¹⁸O与温度增大的方向相同，即大气水线将沿着稳定同位素比率增大的方向逐渐偏离^[22]；因此温度越高，湿度越小，偏离程度就越大，大气水线的截距值也将随着偏离的程度而偏向小值^[23].因此，黑龙江地区大气水线的斜率和截距值可以较好的反映中高纬度东部季风区的自然地理状况和气象条件.

为了深入探究研究区不同降水类型对降水中稳定同位素变化的影响，本文根据降水类型的差异将两站点的样品分为5组，试图分别进行δD、δ¹⁸O值相关性分析，结果如表 1.

由于哈尔滨站降雪样品和齐齐哈尔站雨夹雪样品数量不足，没有进行相关性分析，根据3组计算结果，16个齐齐哈尔降雪样品与同站点降水样品相比，斜率值较高(7.54)，截距值较高(-3.15).与哈尔滨站点降雨样品值对比，差距较大并且都偏高值^[6].这和水汽来源和相变过程中的气象条件有很大关系，降雪主要集中在冬半年，降水水汽干燥，降水过程中经过了多次再蒸发，使得降水中δ值偏低，斜率值较高^[24].与哈尔滨站点的差异变化则主要是受到纬度较高导致的月平均气温较低和水汽输送距离的不同^[24].

GNIP监测时间段内，黑龙江地区降水中δ¹⁸O、气温、降水量的年际变化特征如图 2所示，其中哈尔滨站年平均气温变化范围是4.17~10.97 ℃，年降水量变化范围是347~655 mm，δ¹⁸O的变化范围是-8.26‰~-10.97‰.齐齐哈尔站年平均气温变化范围是3.63~5.94 ℃，年降水量变化范围是350~705 mm，δ¹⁸O的变化范围是-12.27‰~-17.58‰.可以看出，两监测站降水量变化不大，但年平均气温有0.5~5 ℃差异，δ¹⁸O、δD均在全球降水的变化范围之间，但δ¹⁸O、δD值总体偏低，且纬度越高，δ值越低，齐齐哈尔站较哈尔滨站δ值偏低，体现了明显的纬度效应^[5]，即重同位素随着纬度的升高而逐渐贫化的现象，这主要是瑞利分馏在温度变化和水汽输送过程的体现.

根据监测时段2个站点的δ¹⁸O值与气象资料，对δ¹⁸O与该月气温及降水量分别进行回归检验，得出δ¹⁸O=0.35T-16.26，R=0.53，n=80；δ¹⁸O=0.04P-16.04，R=0.19，n=80.可见，在年际尺度上，δ¹⁸O表现出明显的温度效应，即温度是在中高纬度地区影响大气降水中稳定同位素组成的主要因素，δ¹⁸O表现为与温度正相关，随着纬度升高，相关性越明显^{[4, 26]}；但δ¹⁸O与降水量相关性不显著，研究表明，降水量效应往往出现在低纬度的沿海地区，这可能与水汽来源有关，但在中高纬度地区，温度效应大大掩盖了降水量效应^[18]，使得黑龙江地区δ¹⁸O与降水量没有表现出显著的相关性.

为了探究年际尺度上研究区水汽相变及与气象要素的关系，本文分别对两监测站点的不同降水类型样品与气温、降水量的相关性做回归分析.齐齐哈尔站32个降雨样品与降水量和气温的相关关系为：δ¹⁸O=0.01P-12.18，R=0.01，n=32，δ¹⁸O=0.39T-16.85，R=0.32，n=32.齐齐哈尔站16个降雪样品与降水量和气温的相关关系如下：δ¹⁸O=0.20P-23.39，R=0.04，n=16，δ¹⁸O=0.17T-19.86，R=-0.03，n=16.哈尔滨站28个降雨样品与降水量、气温的回归分析如下：δ¹⁸O=-0.003P-9.33，R=-0.03，n=28；δ¹⁸O=0.05T-10.58，R=-0.02，n=28.其中，降雨和降雪样品和降水量均没有体现出相关性，监测样品和气温的相关性也不显著，只有齐齐哈尔站点的降雨样品和气温表现出相关性，这可能是由于高纬度地区降雨量效应不明显，温度效应在降雪样品中体现也不显著.

黑龙江地区全年内δ¹⁸O的变化范围是-26.40‰~-8.02‰，平均值为-14.58‰；δD的变化范围是-207.40‰~-62.95‰，平均值为-110.35‰，表现出夏高冬低的季节变化特征.其中哈尔滨站月平均气温最高月为7月(23.7 ℃)，月平均气温最低月为1月(-18.2 ℃)，降水量变化范围为3.4~151.4 mm，分别出现在1月和7月. δ¹⁸O的变化范围为-5.88‰~-17.81‰，和月均温呈现出显著的正相关，δ¹⁸O最低值-17.81‰出现在1月如图 4(a).齐齐哈尔站最低月气温为1月(-18 ℃)，月均温最高月为7月(21.9 ℃)和8月(21.8 ℃)，月降水量变化范围为1.4~175.6 mm，分别为1月和7月. δ¹⁸O的变化范围为-7.70‰~-26.42‰，最低δ¹⁸O值在2月(-26.42‰)和12月(-24.11‰)，如图 4(b).与哈尔滨站相比，齐齐哈尔站的降水中δ¹⁸O值显著偏低，月平均气温较低，且降水量季节变化大.

根据黑龙江地区两个站点年内δ¹⁸O值与相关气象监测数据，分别计算δ¹⁸O与月均温(T)、月降水量(P)的回归方程：哈尔滨站：δ¹⁸O=0.17T-13.09，R=0.27，n=10；δ¹⁸O=0.03P-13.28，R=0.15，n=10.齐齐哈尔站：δ¹⁸O=0.42T-17.40，R=0.87，n=12；δ¹⁸O=0.08P-19.34，R=0.42，n=12.可以看出，在年内尺度上，黑龙江地区存在显著的温度效应，δ¹⁸O与气温表现出明显的正相关；降水量效应却不显著，但是δ¹⁸O表现出与月降水量正相关，即夏半年，降水量偏多，δ¹⁸O偏高.

我们利用NCEP/NCAR逐月再分析资料计算了1990-1991年的夏季(a-c)、冬季(d-f)逐月整层水汽输送场，结果见图 5，可以看出，6-8月，影响我国季风区的水汽主要来源于印度洋的西南季风和太平洋的东南季风，但影响黑龙江地区的大气降水则主要以东南季风带来的太平洋水汽为主^[27]，水汽来源于较近的海洋性气团，降水中δ¹⁸O表现为高值；12月至次年2月，受蒙古-西伯利亚高压控制，降水较少，且主要是降雪为主，降水中δ¹⁸O表现为低值.

研究基于黑龙江地区GNIP两站点降水中稳定同位素以及美国NECP/NCAR逐月再分析数据，分析了黑龙江地区年际尺度和季节尺度下降水中稳定同位素的变化特征，计算两站点不同降水类型下δ¹⁸O-δD相关性，以及雨样、雪样中的稳定同位素效应，并利用水汽通量模型追踪季节上影响研究区的降水气团差异，分析了水汽输送过程中两站点间的关系，得出以下结论：

1) 黑龙江地区大气降水线为δD=7.49δ¹⁸O-0.86，R=0.97，n=80，其中哈尔滨站点：δD=5.52δ¹⁸O-19.42，R=0.83，n=30，齐齐哈尔站点：δD=7.58δ¹⁸O-0.14，R=0.98，n=50.哈尔滨站点斜率a表现出较低值，而斜率值也较低，主要是由于哈尔滨站点较齐齐哈尔纬度偏低，年平均气温高，且位置更靠近海洋，在夏半年季风影响下有更多湿润水汽输送，反映了黑龙江地区的自然地理状况和气象条件.

2) 研究区降水中δ¹⁸O、δD值总体偏低，且纬度越高δ值越低，齐齐哈尔站较哈尔滨站δ值偏低，年内尺度上δ¹⁸O表现出夏高冬低的季节变化特征.通过分析两站点不同降水类型中δ¹⁸O-δD相关性，降雪样品与降水样品相比，斜率值较高(7.54)，截距值较高(-3.15)，这和水汽来源和相变过程中的气象条件有很大关系，降雪主要集中在冬半年，降水水汽干燥，降水过程中经过了多次再蒸发，使得降水中δ值偏低，斜率值较高.

3) 年际尺度上，δ¹⁸O与温度线性关系显著(δ¹⁸O=0.35T-16.26，R=0.53，n=80)，表现出明显的温度效应，但降水量效应不显著(δ¹⁸O=0.04P-16.04，R=0.19，n=80)，通过两监测站点的不同降水类型样品做回归分析，发现温度效应和降水量效应均不显著.在年内尺度上，哈尔滨站：δ¹⁸O=0.17T-13.09，R=0.27，n=10；δ¹⁸O=0.03P-13.28，R=0.15，n=10，齐齐哈尔站：δ¹⁸O=0.42T-17.40，R=0.87，n=12；δ¹⁸O=0.08P-19.34，R=0.42，n=12，存在显著的温度效应，δ¹⁸O与气温表现出明显的正相关；降水量效应却不显著，但是δ¹⁸O表现出与月降水量正相关，即夏半年，降水量偏多，δ¹⁸O偏高.

4) 运用逐月再分析资料模拟水汽输送场，结果表明，夏季黑龙江地区降水中δ¹⁸O呈现出高值，大气降水主要来源于低纬度太平洋的暖湿水汽和局地水汽；冬季风期间降水形式以降雪为主，δ¹⁸O值偏低，主要是受到蒙古-西伯利亚高压和高纬度水汽及云团内水汽再蒸发的影响.