石家庄地区大气降水中稳定同位素的环境效应及其不同算法下的大气水线方程数学专业论文设计

1、题 目：石家庄地区大气降水中稳定同位素的环境效应及其不同算法下的大气水线方程目录摘要2Abstract21引言32研究区概况及研究方法42.1研究区概况42.2研究方法43降水同位素数据的大气水线53.1大气降水线53.2不同季节6种不同算法下大气水线的差异63.3不同季节不同算法下的大气水线比较与选择84石家庄地区大气降水的氢氧同位素特征分析94.1大气降水中氢氧同位素的月均变化特征94.2大气降水氘盈余的变化125降水同位素的环境效应分析135.1降水中18O和D的温度效应135.2降水中18O和D的降雨量效应155.3降水中18O和D与水汽压的相关关系166结论与讨论17参考文献18致谢

2、22石家庄地区大气降水中稳定同位素的环境效应及其不同算法下的大气水线方程摘要：根据GNIP中位于水文地质工程地质研究所（河北石家庄）1985年到2003年测得石家庄多年的大水中氢氧同位素的量和气象资料，用6种不同算法计算大气水线方程，对不同的算法进行比较,再结合地理环境特征和气象因素分析表明：春夏秋季PWRMA算法下的rmSSEav最小，因而选择算法PWRMA计算石家庄春夏秋季的大气水线。冬季PWMA的rmSSEav最小，故选择PWMA来计算石家庄冬季的大气水线方程。大气降水中d的值呈7月到次年3月高，4-6月低的特点。之所以全国大气水线方程的截距和斜率都比该地区大，主要是石家庄气候干燥的缘故

3、。并且得到了石家庄地区D和18O月均变化表现为后半年(411月)高、前半年(12月次年3月)低、且过量氘受水汽压、温度、降水量的综合环境效应。其中温度效应占主导，呈正相关关系且相关性明显。关键词：石家庄地区；不同算法；大气水线方程；环境效应；Environmental Effects of Stable Isotopes in Atmospheric Precipitation in Shijiazhuang Area and Atmospheric Water Line Equation under Different AlgorithmsAbstract: This paper calcu

4、lates the atmospheric waterline equation with six different algorithms based on the quantity and meteorological data of hydrogen and oxygen isotopes in the flood water of Shijiazhuang from 1985 to 2003 from the Institute of Hydrogeological Engineering Geology (Shijiazhuang, Hebei) in GNIP. The rmSSE

5、av of this algorithm is the smallest in spring, summer and autumn PWRMA. The rmSSEav of winter PWMA is the smallest, so choose PWMA to calculate the atmospheric waterline equation in Shijiazhuang winter. d values in precipitation are high from July to March and low from April to June. The reason why

6、 the intercept and slope of the national atmospheric waterline equation are larger than that of the region is mainly due to the dry climate in Shijiazhuang. The D and 18O monthly changes in shijiazhuang area were shown to be high in the second half (april-oct), low in the first half (november-march)

7、, and excessive deuterium was subjected to vapor pressure, temperature and precipitation. Among them, temperature effect is dominant, positive correlation and obvious correlation.Key words:Shijiazhuangarea; different algorithms; atmospheric waterline equation; Environmental effects;1引言首次国际大气降水同位素研究开

8、始于20世纪50年代1，在20世纪60年代全球氢氧同位素监测网(GNIP)由世界气象组织(WMO)和国际原子能机构(IAEA)共同创建，为监测全世界大气降水同位素奠定了基础。大气水线方程是进行区域研究的基准，是降水中D和18O采用最小二乘法进行回归分析得到的，它可以表示为D=a18O +b2，其中a反映18O和D分馏快慢，b表示D偏离于平衡状态的程度3。Craig(1963)4通过对北美大陆不同水体中氢氧同位素的研究得出D和18O间存在D=818O +10的关系，从此以后它们的这种关系被称作全球大气水线方程(GMWL)。降水同位素监测在我国研究的比较晚，在1983年前我国仅香港一个观测点5，后

9、来创建银川、昆明等好几个站点6。至今，我国在上述站点的数据基础上得到多有效的理论支撑，对降水在局部地区，甚至在大尺度，大空间7-18的稳定同位素分布特征、水汽输送过程、变化规律与机制、同位素效应等已获得了一系列重要的认识。郑淑蕙等19最先研究得出我国大气水线方程，将其表示为D=7.918O.+8.2;经过对不同大气水线分析，其不同不仅与采集降水样本年限、地点和数量有关，与算法的选取密不可分。大气水线计算多用 (OLSR)，其缺点是对样品给予同等重要的地位(IAEA，1992)。王圣杰20基于6种不同的算法研究过天山地区的大气水线。发现在OLSR、MA、RMA、PWLSR、PWMA、PWRMA中

10、，该地区的大气水线计算方法中最适合的是PWLSR算法。因此，本文结合GNIP数据，用OLSR、MA、RMA、PWLSR、PWMA、PWRMA这6种不同计算方法应用到石家庄地区的大气水线中，并对计算结果进行比较。结合降水中同位素的水汽压、降雨量、温度以及大气降水中氘盈余同石家庄的大气水线研究，为以后该地区的大气水线特征提供数据基础。2研究区概况及研究方法2.1研究区概况石家庄市地处中纬度河北省中南部，温带大陆性季风气候。太阳辐射的季节性变化显著，地面的高低气压活动频繁，春温夏热秋凉冬冷，四季温度有明显区别，寒暑悬殊。该地区雨量分布不均，夏冬季长，春秋季短。夏季受海洋流的影响，降雨占全年65%左右

11、，空气湿度在夏季中期接近100%。石家庄地区属于东北和南部季风过渡区间，这些多重条件对大气降水氢氧同位素变化产生深刻的影响21。2.2研究方法本文引用GNIP中石家庄站点气象数据资料，均记录了每月15号18O、D值以及水汽压、降水量和温度等监测数据，它表示本月上半月和上个月下半月的同位素数据平均数值。本文应用大气降水中与算术平均值和加权平均值计算，公式为:式中:;；。氢氧同位素的算术平均数反映该地区的气候特征，氢氧同位素降水量加权平均数则表示降水对同位素浓度大小的影响。除最常用的最小二乘回归(OLSR)来研究大气水线方程，还有主轴回归(MA)、简化主轴回归(RMA)、以及对应加权PWLSR、P

12、WMA和PWRMA;由rmSSEv rmSSEt.和rmSSEp的算术平均数得到平均标准化均方根误差(rmSSEav)，如果该值越靠近1，表示回归效果越好越稳定。三种加权的计算中rmSSEav也采用求算术平均数，上述6种算法的大气水线斜率、截距以及rmSSEav等主要参数采用澳大利亚核科学技术组织(ANSTO)研发的软件Local Meteoric Water Line Freeware. 计算得到。3降水同位素数据的大气水线3.1大气降水线大气降水线是指全球或特定区域中进行区域研究的基准，是降水中D和18O采用最小二乘法进行回归分析得到的，某种程度上说明研究区域的自然状况，同时也对水文地质和

13、水循环中稳定同位素的研究意义重大。Craig22对全球范围内的大气降水样本进行分析，得到全球大气水线(GMWL)方程: 。氢氧同位素间分馏的比率可以用水线方程的斜率表示，海洋表面蒸发时的动力同位素效应所导致氘偏离平衡状态的程度，用截距表示23。全球大气水线(GMWL)作为参照，因为各地气象条件、地理因素和降水过程的不同，各个站点大气水线都会与GMWL有幅度不同的偏移。部分专家调查研究了我国及各局地大气降水线方程，郑淑蕙等19得到的我国大气降水线方程。柳鉴容等24得出为东部季风区局地大气降水线方程，截距和斜率比我国大气降水线方程的小。本文跟据GNIP所提供的D和18O气象资料，用最常用的方法OL

14、SR得出的石家庄地区大气降水线方(与全国大气降水线方程相比，该方程的截距和斜率均表现得较小)。这是因为石家庄地区不同水汽来源导致不同的氢氧稳定同位素比率，并且该地区降水过程中也进行一定程度的二次蒸发，其蒸发模式与全球大气降水线有明显的不同。图1最小二乘法得出的大气降水线Fig.1 Precipitation line from least square method3.2不同季节6种不同算法下大气水线的差异表1春季不同算法的斜率截距和平均标准化均方根误差分布Table 1 Slope Intercept and Mean Standardized Root Mean Square Error

15、Distribution of Different Algorithms in SpringabrmSSEavOLSRRMAMAPWLSRPWRMAPWMA3.66657.297614.32204.31326.784810.5357-28.3061-7.905131.5604-29.7673-16.29534.15021.38251.09991.39231.25651.09351.1720表1是基于石家庄GNIP数据下6种不同算法春季的大气水线。从表可以看出，对比OLSR和PWLSR,后者的斜率大于前者,截距小于前者;RMA和加权简化主轴回归的对比，加权的斜率和截距比前者小；主轴回归和加权主轴

16、回归相比，后者的斜率和截距表现为后者小于前者。表2夏季不同算法的斜率截距和平均标准化均方根误差分布Table 2 Slope Intercept and Mean Standardized Root Mean Square Error Distribution of Different Algorithms in SummerabrmSSEavOLSRRMAMAPWLSRPWRMAPWMA5.42098.119612.06055.25457.628610.9761-14.11656.905737.6029-14.06225.006631.89421.19851.06191.19941.2117

17、1.06301.1442表2是基于石家庄GNIP数据下6种不同算法夏季的大气水线。从表可以看出，对比OLSR和PWLSR,后者的斜率小于前者,截距大于前者;RMA和加权简化主轴回归的对比，加权的截距和斜率比前者小；对主轴回归和加权主轴回归相比，后者的斜率和截距表现为后者小于前者。表3秋季不同算法的斜率截距和平均标准化均方根误差分布Table 3 Slope Intercept and Mean Standardized Root Mean Square Error Distribution of Different Algorithms in AutumnabrmSSEavOLSRRMAMAP

18、WLSRPWRMAPWMA2.76837.269118.77965.25888.338313.1071-28.25325.415791.5209-10.921012.711349.30821.60651.13251.63971.15601.10791.2805表3是基于石家庄GNIP数据下6种不同算法秋季的大气水线。从表可以看出，对比OLSR和PWLSR, 后者的截距和斜率大于前者; RMA和加权简化主轴回归的对比，加权的截距和斜率比前者小；对主轴回归和加权主轴回归相比，后者的斜率和截距表现为后者小于前者。表4冬季不同算法的斜率截距和平均标准化均方根误差分布Table 4 Slope Inte

19、rcept and Mean Standardized Root Mean Square Error Distribution of Different Algorithms in WinterabrmSSEavOLSRRMAMAPWLSRPWRMAPWMA4.57766.46369.01915.63916.13516.6477-26.0222-6.463420.0382-9.47234.39020.86311.16201.05261.16241.04291.02541.0205表4是基于石家庄GNIP数据下6种不同算法冬季的大气水线。从表可以看出，对比OLSR和PWLSR, 后者的截距和斜率

20、大于前者;RMA和加权简化主轴回归的对比，加权的截距和斜率比前者小；对主轴回归和加权主轴回归相比，后者的斜率和截距表现为后者小于前者。表1-4是在石家庄GNIP数据6种算法下四季的大气水线。可以看出，前3种未加权的算法中OLSR算法的截距、斜率最小，MA算法的最大。加权后斜率和截距与上述趋势一致。3.3不同季节不同算法下的大气水线比较与选择在6种不同算法的大气水线方程中，各算法的mSSEav都较小，相比之下在春夏秋季PWRMA算法的rmSSEav最小，最接近1，故最为稳定，因而选择算法PWRMA计算石家庄春夏秋季的大气水线。冬季PWMA算法的rmSSEav最小，最接近1，故定，因而选择算法PW

21、MA计算石家庄冬季的大气水线。4石家庄地区大气降水的氢氧同位素特征分析4.1大气降水中氢氧同位素的月均变化特征表5石家庄大气降水不同月份18O的差异及均值()Table 5 vary and average value of 18Oin different months of precipitation in Shijiazhuang (). 月份1 2 3 4 5 6 7 8 9 10 11 12 1985 -7.3 -7.6 -11.2 -10.0 -7.5 -15.4 1986 -5.4 -12.4 -5.3 -1.5 -2.1 -10.2 -7.5 -7.7 -7.7 -6.3 -4.

22、7 -10.0 1987 -9.9 -14.9 -12.2 -4.1 -4.1 -8.7 -4.7 -8.7 -5.6 -8.2 -7.2 1988 -16.0 -9.4 -5.0 -3.0 -8.1 -6.1 -5.8 1989 -8.3 -12.5 -4.8 -3.0 -9.5 -7.8 -9.8 -7.2 -4.8 -6.7 -6.4 1990 -12.8 -12.9 -12.0 -7.5 -3.9 -6.5 -11.9 -5.4 -8.5 -6.4 -10.2 -5.8 1991 -8.5 -9.1 -6.8 -9.8 -5.3 -9.5 -5.2 -5.2 -6.1 -3.0 -11

23、.6 1992 -4.1 -10.3 -8.3 -6.3 -5.1 -7.9 1995 -9.4 -12.2 -12.6 -6.2 1996 -14.2 -7.9 -1.1 -2.1 -5.1 -10.0 -8.8 -2.6 1997 -4.2 -9.1 -3.8 -9.4 -5.2 -7.1 -6.3 -4.8 -12.2 1998 -14.1 -8.4 -4.5 -2.4 -5.0 -8.2 -8.5 -8.3 -6.2 1999 -0.3 -3.5 -3.6 -9.0 -8.0 -5.3 -15.3 2000 -14.0 -8.6 -7.2 -10.1 -5.1 -7.8 -9.0 20

24、01 -13.3 -3.2 -9.9 -7.4 -7.6 -10.6 -8.2 2002 -2.6 -0.6 -5.5 -9.7 -7.0 -4.2 -7.9 -20.6 -11.5 2003 -11.5 -2.5 -5.8 -2.0 -7.4 -8.9 -7.5 -10.8 -8.0 -4.0 均值-11.9 -4.2 -7.6 -2.3 -4.4 -6.7 -8.8 -7.3 -7.5 -8.8 -10.5 -11.8 表6石家庄地区大气降水不同月份D的变化范围及平均值()Table 6 variation range and average value of D in different

25、 months of precipitation in Shijiazhuang area ()月份1234567891011121985-48-47.6-83.5-61.6-41.8-97.21986-45.1-82.2-30.2-15-9.5-71.1-58.6-54.9-53.6-39.6-32.4-69.91987-57.1-111.3-103.7-14.6-21.5-58.9-29.5-60.2-41.7-55-44.91988-100.5-60.6-21.5-23.4-53.6-37.8-33.21989-53.4-90.1-23.1-34.2-67.7-55.7-69.7-55.

26、5-30.3-55.9-471990-91.9-80.2-82.1-50.3-45.7-56.3-95.1-37.6-57.3-47.7-60.3-36.41991-54.2-63.1-64.7-71.9-46-82.4-42.6-29.4-43-26.9-65.11992-30.1-58.8-30.2-66.1-46.1-31.4-50.4H31995-48.1-47.7-83.6-61.7-41.9-97.3-45.21996-82.3-30.3-15.1-9.6-71.2-58.7-54.1-53.7-39.7-32.5-69.1-57.21997-111.4-103.8-14.7-21

27、.6-58.1-29.6-60.3-41.8-55.1-44.1-100.61998-60.7-21.6-23.5-53.7-37.9-33.3-53.51999-90.2-23.2-34.3-67.8-55.8-69.8-55.6-30.4-55.1-47.1-91.12000-80.3-82.2-50.4-45.8-56.4-95.2-37.7-57.4-47.8-60.4-36.5-54.32001-63.2-64.8-71.1-46.1-82.5-42.7-29.5-43.1-26.1-65.2-30.22002-58.9-30.3-66.2-46.2-31.5-50.5H42003-

28、48.2-47.8-83.7-61.8-41.1-97.4-45.3均值-91.1 -66.5 -33.6 -33.5 -53.0 -60.5 -55.1 -50.3 -44.7 -51.1 -60.1 -61.7 根据所得到的资料数据，石家庄地区大气降水的18O范围为20.620.3，均值为7.74;D变化范围为111.39.5，平均值为54.37，虽然变化幅度大，但都在郑淑慧等19研究的我国降水的变化范围中。按月份分析了该地区18O和D的特征。发现石家庄地区大气降水中18O和D月变化特征均表现为(411月)高、(12月次年3月)低。图2石家庄大气降水18O平均值随月份的变化Fig. 2 v

29、ariation of 18Omean value of precipitation in Shijiazhuang with month图3石家庄大气降水D平均值随月份的变化Fig. 3 variation of D mean value of precipitation in Shijiazhuang with month图4石家庄大气降水18O和D的加权平均值随月份的变化Fig.4 variation of weighted average of 18Oand D in Shijiazhuang with month由图4研究表明，石家庄地区大气降水中氧氢同位素月均值变化整体趋势一致。1

30、8O的变化在水循环中受水汽蒸发作用十分明显，所以以18O为例来研究氢氧稳定同位素的特征变化。总体来说石家庄地区降水中的18O表现为12月到次年3月低、4月到11月高的特征，表现出温带季风气候的特征25。5月份加权平均值最高是同位素最为富集的时期，;1月降水中氢氧稳定同位素月平均值出现最低值;降水中氢氧稳定同位素在7-10月份较为稳定。之所以导致降水中稳定同位素不同可能有:(1)在运输，凝结，降落过程中水汽受该地气象及地理因素的作用，同位素发生的富集和耗尽程度不同，从而18O值出现了月均差异;(2)研究区受降水中水汽源差异和不同季风气候的影响。4.2大气降水氘盈余的变化由全球大气降水线GMWL，

31、Dansgaard1提出了氘盈余(d)的定义，并把它定义为d=D818O。各地区的氘盈余不同，是由于全球大气降水线和各地大气降水线方程都存在一定程度的偏离。大气降水主要来自海水，当氘盈余d值为0时，海水的水汽凝结与蒸发达到动态平衡速率相等。因为全球大气降水的氘盈余是10。所以d值的月差异和地域不同实质是大气降水在蒸发凝结中18O和D分馏速度的不同，这除了与形成该地区降水的气候条件和水汽平衡条件有关，还与降水源的风速、相对湿度、温度等因素密不可分26,27。有学者研究22，当水汽源地相对湿度较低时，d值较高，反之d值较低;在干旱情况下，强烈的蒸发会引起动力分馏系数提高，从而使d值增加。图5石家庄

32、地区d的月变化Fig. 5 monthly variation of d in Shijiazhuang area由图5石家庄地区氘盈余（d）值月均变化，从图中可以看出石家庄地区氘盈余(d)月均变化在-4.1815.59范围中，47月总体是显著小于其它月份，夏秋两季d值大部分小于10，春冬两季d值大部分在10以上。5降水同位素的环境效应分析降水中稳定同位素的环境效应是同位素和各气象因子之间的相关性，因为各个地域的地理和气候差异，所以影响同位素变化的主导因子也会不同。已有专家发现温度效应和降水量效应是引起我国大气降水稳定同位素时空变化主要因素，非季风区受温度效应影响显著，季风区受降水量效应作用更

33、加突出8,28。温度、降雨量、风速风向和水汽压等是影响我国东部季风区降水18O的主导因素24，其中，降雨量和温度的影响尤其显著。5.1降水中18O和D的温度效应温度效应是18O和D与地面平均气温往往呈线性关系。温度对降水同位素组成的影响是各因素中最为重要的因素29，温度效应主要是由于蒸发过程中分馏作用随温度的升高而减弱造成的。进而使18O的值增高。但仅有同位素含量趋于稳定时，气温和同位素正相关十分显著30。温度效应常在中高纬度内陆地区更为突出，越到内陆地区，正相关性越显著31，这种现象在我国主要表现在季节温度变化比较大的地区。如图6，图7得到18O与温度t的线性方程，D与温度t的线性方程。图6

34、石家庄地区18O与月均温度的关系Fig.6 relationship between 18O and monthly mean temperature in Shijiazhuang area图7石家庄地区D与月均温度的关系Fig7. The relationship between D and the monthly temperature in Shijiazhuang area如图6，图7所示，18O、D与温度之间的线性关系反映了温度效应，从图中得出，在我国石家庄地区温度效应十分显著，但在很多东部季风区和南方季风区的温度效应表现很弱32，可以看出虽然石家庄地区与北方内陆地区降水同位素温度

35、效应特征接近，但不是十分明显，没有西北地区典型。5.2降水中18O和D的降雨量效应一般地，雨量越大，降水18O和D值越低，这种效应称为降雨量效应。主要发生在中低纬沿海地区，我国内陆表现通常不显著，并且它的产生与大气强烈的对流现象相关。尤其在暴雨时期该效应表现的越显著33。Dansgaard1研究分析全球各站台所观测的大气降水的同位素数据资料，得到降水量效应。已有学者分析研究，在潮湿的空气等温情况下，水汽凝结为液态水将服从瑞利的分馏模式。因此产生降雨量效应可能与雨滴下落过程中的水蒸气交换和水蒸气蒸发有关34。蔡明刚等27提出研究区中降水的水汽来源较为单一是“降雨量效应”的重要条件之一，如果在季风

36、区，季风也能引起该效应。图8石家庄地区18O值与月均降雨量的关系Fig. 8relationship between 18Ovalue and monthly average rainfall in Shijiazhuang area图9石家庄地区D值与月均降雨量的关系Fig.9 relationship between D value and monthly average rainfall in Shijiazhuang area如图8，图9，对石家庄地区1985-2003年（1992-1994年数据缺失）的18O和D与降雨量P的关系，得出石家庄18O与降雨量的线性方程:，D与降雨量的线性

37、方程: 。(见图8与9)。5.3降水中18O和D与水汽压的相关关系图10石家庄地区18O值与月均水汽压的关系Fig.10 Relationship between 18O value and monthly water vapor pressure in Shijiazhuang area图11石家庄地区D值与月均水汽压的关系Fig.11 nexus between D value and monthly water vapor pressure n Shijiazhuang area结果表明石家庄地区大气降水中氢氧同位素不但与温度、降水量线性相关，而且与水汽压相关。6结论与讨论石家庄位于我国

38、的中高纬东部季风地区，受到严重的季风气候作用。对石家庄地区大气降水中氢氧同位素的特征变化和组成的分析，研究它们与气象要素以及水汽来源的关系，得出如下结论:（1）不同算法存在着不同程度的差异。对最小二乘回归(OLSR)、主轴回归(MA)、简化主轴回归(RMA)3种计算方法比较，OLSR算法的斜率和截距最小，RMA算法的居中，MA算法的最大。对加权简化主轴回归(PWRMA)、加权主轴回归(PWMA)和加权最小二乘回归(PWLSR)算法的斜率和截距也是PWLSR算法的最小，PWRMA算法的居中，PWMA算法的最大。大气水线方程的截距和斜率在春、夏、秋、冬4季存在着差异性。对比OLSR和PWLSR,的

39、对比发现，秋冬两季后者的斜率和截距比前者大，春季后者的斜率比前者大但截距后者比前者小，夏季斜率和截距与春季相反RMA和PWRMA的比较发现，春夏两季后者的斜率和截距比前者小，在秋季后者的斜率和截距比前者大，冬季后者的斜率比前者小但截距后者比前者大对MA和PWMA，在春夏秋冬四季，加权的斜率和截距表现为后者小于前者。（2）在6种算法的大气水线中，算法的平均标准化均方根误差mSSEav中春夏秋3季PWRMA算法的rmSSEav最小，故选择算法PWRMA计算石家庄春夏秋季的大气水线。冬季PWMA算法的rmSSEav最小，故最为稳定，因而选择算法PWMA计算石家庄冬季的大气水线。(3)石家庄地区气象要

40、素及水汽来源的差异对降水中氢氧同位素含量都有影响。并对石家庄地区大气降水的氢氧同位素特征分析得到该地区的18O表现为12月到次年3月低、4到11月高的特征。(4)石家庄地区气候偏干旱，降水蒸发效应十分明显，其水线方程为，斜率和截距都小于我国乃至全球的降水线方程。d值变化特征呈47月总体是显著小于其它月份，夏秋两季d值大多小于10，春冬两季d值大多在10以上。反映了不同季风区降水水汽来源地及蒸发条件有差异。(5)石家庄大气降水的氢氧稳定同位素的温度效应显著，降雨量效应微弱，全年尺度下18O与温度t的线性方程为18O=0.1101t9.311(2=0.2541)，D与温度t的线性方程为D=0.56

41、63t62.44(2=0.2315）18O与降雨量的线性关系:18O=0.0062P8.263(2=0.0.0197)，D与降雨量的线性关系:D=0.0403P56.233(2=0.0 291)。石家庄地区大气降水的氢氧稳定同位素在不同月份均具有温度效应，结果得出石家庄地区同位素比率受水汽压，降水量和温度共同作用，但温度的影响更显著。在6种算法的大气水线中，春夏秋季3个季节选择PWRMA算法来计算大气水线。冬季PWMA算法的rmSSEav最小更接近1，因而选择算法PWMA来计算石家庄冬季的大气水线。参考文献：1Dansgaard W. The abundance of 18 O in atmo

42、spheric water and water vapour. J. Tellus, 1953, 5(4) : 461-469.2Friedman I. Deuterium content of natural waters and other substancesJ. Geochimica et Cosmochimica Acta, 1953, 4(1-2) : 89-103.3黄一民,章新平,孙葭.长沙大气水线及与局地气象要素的关系J.长江流域资源与环境,2014,23(10):1412-1417.4Craig.H,GordonLI,HoribeY.Isotopic.exchangeeff

43、ects.in.the.evaporation.of.water:Low-temperature.experimental.resultsJ.Journal.of.Geophysical.Research,1963,68(17): 5079-5087.5李正红,张发旺.全球降水氢氧同位素研究进展J.勘察科学技术,2004(1):1-6.6赵琦经,顾慰祖,顾文燕,等.中国降水同位素站网J.水文,1995(5):25-27.7陈中笑,程军,郭品文,等.中国降水稳定同位素的分布特点及其影响因素J.大气科学学报,2010,33(6):667- 679.8Liu Z F, Tian L D, Chai

44、X R, et al. A model-based determination of spatial variation of precipitation 18O over China J. Chemical Geology, 2008, 249(1-2) : 203-212.9Zhang X P, Nakawo M, Yao T D, et al. Variations of stable isotopic compositions in precipitation on the Tibetan Plateau and its adjacent regionsJ. Science in Ch

45、ina ( Series D), 2002, 45 (6) : 481-493.10Tian L D, Yao T D, Schuster P F, et al. Oxygen- 18 concentrations in recent precipitation and ice cores on the Tibetan PlateauJ. Journal of Geophysical Research, 2003, 108( D9) : 4293-4302.11陈粉丽,张明军,马潜,李小飞,王圣杰,李菲.兰州及其周边区域大气降水(18)O特征及其水汽来源J.环境科学,2013,34(10):3

46、755-3763.12Liu J R，Song X F，Sun X M，et al. Isotopic composition of precipitation over arid Northwestern China and its implications for the water vapor originJ. Journal of Geographical Sciences, 2009, 19(2) : 164-174.13章新平,姚檀栋.我国降水中18O的分布特点J.地理学报,1998,53(7):356-364.14Zhang X P, Liu J M, Sun W Z, et a

47、l. Relations between oxygen stable isotopic ratios in precipitation and relevant meteorological factors in southwest ChinaJ. Science in China ( Series D), 2007, 50( 4) : 571-581.15Liu J R, Song X F, Yuan G F, et al. Characteristics of 18O in precipitation over eastern monsoon China and the water vap

48、or sourcesJ. Chinese Science Bulletin, 2010, 55( 2) : 200-211.16于津生,虞福基,刘德平.中国东部大气降水氢、氧同位素组成J.地球化学,1987,16(1):22-26.17李小飞,张明军,马潜,李亚举,王圣杰,汪宝龙.我国东北地区大气降水稳定同位素特征及其水汽来源J.环境科学,2012,33(09):2924-2931.18刘鑫,宋献方,夏军,等.黄土高原岔巴沟流域降水氢氧同位素特征及水汽来源初探J.资源科学,2007,29(3):59-66.19郑淑蕙,侯发高,倪葆龄.我国大气降水的氢氧稳定同位素研究J.科学通报,1983,28(13):801806.20王圣杰.天山地区降水稳定氢氧同位素特