2.2试验方案设计
黄河源区草地退化试验中的初始场、积分时间、模拟区域和水平分辨率等与柳媛普等[34]的试验设计相同,不同的是本文采用的模式为WRF(Version3.6.1),物理过程参数化方案为:
RRTM长波辐射、Dudhia短波辐、Noah陆面过程、Kain-Fritsch积云对流参数化方案、YSU行星边界层方案以及WSM6微物理过程[32]。
表2.1实验方案设计
表2.2植被参数
为了证明分段积分方法模拟源区草地退化前后的准确度更高,进而得到更可信的结果,本文设计两组试验。
试验一包括控制试验及敏感性试验,分别在源区草地退化前后的条件下采用连续积分方法进行模拟,试验中模拟范围、敏感性试验中修改植被类型的方案、区域以及修改土壤湿度的方案与柳媛普等[34]的设计方案相同;试验二也包括控制试验与敏感性试验,采用分段积分方法进行与试验一相同的模拟研究。
试验方案的具体设计见表2.1。
两组试验的控制试验采用模式原有的植被类型--草地,将NCEP提供的大尺度再分析资料作为初始值;两组试验的敏感性试验中将源头区及其东北方向兴海和唐乃亥一带植被改为稀疏植被,称为草原退化区,其它地区仍为草地,同时,还将草原退化区地下10cm和200cm深度处土壤湿度改为实际观测值。
控制试验模拟区域地形高度、土地覆盖情况及源区范围如图2.1(a)所示,敏感性试验的土地覆盖情况如图2.1(b)所示。
模拟区域中涉及到的主要植被参数见表2.2。
图2.1[34]控制试验及敏感性试验地形、植被分布(a)控制试验地形高度、植被分布及源区范围,(b)敏感性试验植被分布,图中粗线为长江、黄河
2.3分段积分方法验证
首先分别将两组试验中控制试验模拟的地面温度、2m高度处气温和NCEP提供的再分析资料(1°×1°)中对应量进行对比分析,检验分段积分方法在黄河源区的模拟能力。
两组试验中控制试验模拟得到的以及由再分析资料得到的源区2003年7月地面2m高度处气温的月平均日变化见图2.2(a),从中可以看出分段积分模拟得到的月平均日变化曲线很接近再分析资料得到的月平均日变化曲线,只是数值稍微偏低。
对于连续积分得到的结果,在14:
00(北京时下同)数值比再分析资料中的值小约2℃,然而在02:
00,连续积分得到的数值比再分析资料得到的数值高。
图2.2(b)为两组试验中控制试验模拟得到的以及由再分析资料得到的源区地面气温月平均日变化曲线。
从图中可以看出分段积分方法模拟得到的曲线基本上与再分析资料得到的曲线吻合,然而对于连续积分得到的结果,在14:
00数值比再分析资料得到的数值高,02:
00的数值又比再分析资料得到的相应数值低约2℃。
总体来看,对温度的模拟分段积分方法的准确度明显高于连续积分方法,由此可以相信,分段积分方法在源区草地退化对该区域气候影响的模拟的准确性比传统的积分方法高,模拟得到的结果更可信。
图2.2两组控制试验与再分析资料得到的2m高度处温度月平均日变化(a)和地面温度月平均日变化(b)(单位:
°C)
2.4结果与分析
植被类型改变时,下垫面的反照率、土壤有效湿度、粗糙度等相应的发生变化,对地—气系统的能量造成影响,进一步引起近地层内各气候环境要素发生变化,下面通过地表温度、2m高度处气温、湿度和能量变化等方面分析源区草地退化对该区域气候的影响并验证分段积分方法得到的结果与连续积分方法得到的结果存在差异。
2.4.1气温变化
2.4.1.1气温的月平均日变化
气温日变化受太阳辐射强度的影响,一天当中在12:
00左右达到最高值。
图2.3(c)与图2.4(c)为柳媛普等将MM5模式得到的源区草地面积退化前后离地面2m高度处气温及地面温度分别求区域平均得到区域平均2m高度气温及地面温度7月份变化,然后对7月份区域平均2m高度气温及地面温度求平均得到2m高度气温及地面温度月平均日变化情况[34]。
采用相同的处理方法,得到分段积分与连续积分2m高度处气温及地表温度7月平均日变化情况(图2.3(a、b)、图2.4(a、b))。
对于2m高度处气温,由图2.3(c)可知,MM5模式模拟得到的源区草地退化前月平均2m高度处气温最高值出现在一天当中的12:
00为14.0℃左右,最低气温出现在凌晨02:
00低于4℃。
源区草原退化后2m高度气温白天升高比较明显,尤其在11:
00-17:
00时段内,相反,20:
00到次日清晨在降低,进一步加剧了昼夜温差。
然而柳媛普等曾证明利用MM5模式模拟源区草地退化时温度的模拟结果偏低[35]即源区草地退化前后的实际值应高于图2.3(c)的结果。
从连续积分方法模拟得到的源区草地退化前后7月2m高度处气温月平均日变化(图2.3(a))可以看出,源区草地退化后,2m高度处气温在08:
00-17:
00时间段内小于退化前,在17:
00-02:
00时间段内反而高于退化前的温度,从总体来看,白天温度减小夜晚温度增加,这说明源区草地退化后该区域的昼夜温差将减小,这与MM5模式模拟得到的结果相反;从分段积分模拟实验所得结果(图2.3(b))可以看出,与源区草地退化前相比,源区草地退化后2m高度处温度在08:
00-20:
00时间段内增加,且在14时增加值达到最大为1.1℃,晚上20时到次日凌晨虽有升高但不明显,由此可知源区草地退化后该区域的2m高度处温度有不同程度的升高,昼夜温差增大,这与MM5模式模拟得到的结论相同,且分段积分得到的数值比MM5模式模拟的数值高。
因此分段积分方法与连续积分方法模拟得到的结果存在差异,分段积分模拟得到的结果可信度更高。
08:
0014:
0020:
0002:
0008:
0014:
0020:
0002:
00
(北京时)(北京时)
(c)
图2.3连续积分方法(a)、分段积分方法(b)和柳媛普等[34](c)模拟得到的源区草地退化前后7月2m高度处气温月平均日变化(单位:
°C)(注:
图(c)中E1代表草地退化前试验即控制试验,E2代表草地退化后试验即敏感性试验,下同)
对于地表温度,由图2.4(c)可知,柳媛普等利用MM5模式模拟的源区草地退化后地表温度在白天升高非常明显在08:
00开始升高,14:
00升高幅度最大,约为3℃,20:
00开始降低,进一步加剧了昼夜温差。
从连续积分方法模拟得到的源区草地退化前后7月地面温度月平均日变化(图2.4(a))可以看出,源区草地退化后模拟的地面温度与退化前相比,退化后得到的地面温度在白天升高幅度较小,夜晚温度增加,最终使地面温度增加,但不会导致昼夜温差增大。
分段积分模拟得到的结果见图2.4(b),将源区草地退化前后的结果相比,可以看出相比于源区草地退化前,退化后的地表温度从08:
00-20:
00增加,20:
00到次日凌晨降低,这将导致该区域昼夜温差增大,与MM5模式模拟得到的结论相符且模拟的数值高于MM5模式得到的数值。
由此看来传统的积分方法与分段积分方法对地面温度的模拟结果存在差异,分段积分方法模拟得到的结果更准确。
因此,当源区草地退化后,2m高度处温度与地面温度均有不同程度的增加,且地面温度增加的幅度更大。
同时由于温度在白天增加,夜晚降低,这将导致黄河源区昼夜温差增大,对生物生存更加不利。
(a)
(b)
08:
0014:
0020:
0002:
0008:
0014:
0020:
0002:
00
(北京时)(北京时)
(c)
图2.4连续积分方法(a)和分段积分方法(b)和柳媛普等利用MM5模式[34](c)模拟得到的源区草地退化前后7月地面温度月平均日变化(单位:
℃)
2.3.1.2气温的空间变化
由气温的月平均日变化可以看出源区草地退化导致该区域地表温度的增加幅度较2m高度处气温的增加幅度大。
为了明显突出分段积分与连续积分的差异,下面比较分段积分方法与连续积分方法所得源区退化前后地面温度差值场并分析源区草地退化后对该区域地面温度的影响。
其中地面温度差值场是用源区草地退化后的结果求得源区7月平均地面温度场减去源区草地退化前相应的平均地面温度场得到的。
柳媛普等研究表明,源区草原退化后对地面的影响范围主要在草原退化区[34]。
因此下面针对草地退化区来分析源区草地退化后对该区域的影响。
本文中草地退化区有两部分,区域一(32.2°~35.5°N,95.8°~99.5°E)和区域二(34.8°~36.1°N,99.8°~102.1°E)。
对于退化区域一,从连续积分方法模拟的地表温度差值场(图2.5(a))可以看出地表温度改变范围在0~4℃,且在(33°~34°N,97°~98°E)范围内有两个等值线密集区且东北侧的密集区的值为1.6℃,西南侧密集区的值为2.4℃。
同时,在(34°~34°30’N,98°30’~99°30’E)范围内有等值线密集区,值为1.6℃。
从图2.5(b)中可知分段积分模拟退化区域一的地表温度改变范围是0~3.2℃,同样在(33°~34°N,97°~98°E)范围内,有两个等值线密集区且东北侧的密集区的值为2.4℃,西南侧密集区的值为1.6℃,与连续积分密集区的值相反,在(34°~34°30’N,98°30’~99°30’E)范围内有等值线密集区,值为0.8℃,比连续积分模拟得到的结果小0.8℃。
图2.5连续积分(a)与分段积分(b)得到的退化区域一地表温度月平均变化(单位:
℃)
图2.6连续积分(a)与分段积分(b)得到的退化区域二地表温度月平均变化(单位:
℃)
对于退化区域二,连续积分模拟得到的地表温度改变范围由图2.6(a)可知为3.8~5.0℃,在(35°45’N,100°30’E)附近存在最大值。
分段积分模拟的地表温度(图3.6(b))改变范围在3.8~5.0℃,与连续积分模拟的结果相差不多,但最大值位于(35°15’N~101E)附近。
总之,分段积分模拟得到的源区草地退化前后在退化区地表温度月平均变化与连续积分得到相应结果存在差异。
源区草地退化会导致退化区的地面温度升高,且对退化区域二的影响大于退化区域一的影响。
2.3.2空气湿度变化
由下垫面植被变化引起的地—气间热通量变化,除了影响地—气系统温度发生变化,也会导致大气湿度状况的改变。
与地表温度变化区域相同,2m高度处湿度降低区主要位于草地退化区,并且草原退化区中心降低幅度最大,下面针对草地退化区来比较分段积分与传统积分方法在敏感性试验中对2m高度处湿度模拟结果的差异并分析源区草地退化后对退化区空气湿度的影响。
图2.7连续积分(a)与分段积分(b)模拟得到的退化区域一2m高度处湿度月平均变化
(单位:
g·kg-1)
图2.8连续积分(a)与分段积分(b)模拟得到的退化区域二2m高度处湿度月平均变化
(单位:
g·kg-1)
对于退化区域一,由连续积分模拟得到的源区草地退化前后2m高度处湿度的月平均变化(图2.7(a))可知中,湿度的减小幅度在0~1.2g·kg-1范围内并且在中心位置附近出现湿度降低最大值0.8g·kg-1;对于分段积分模拟得到的2m高度处湿度(图2.7(b))减小幅度为0.2~1.0g·kg-1。
对于退化区域二,连续积分得到的空气湿度减小幅度0.6~1.6g·kg-1(图2.8(a))且最大值1.6g·kg-1位于退化区域二的西南区域。
分段积分模拟得到的空气湿度减小幅度在1.0~1.4g·kg-1且中心存在最大值1.4g·kg-1。
总之,连续积分比分段积分模拟的空气湿度减小幅度相差0.2g·kg-1,虽然数值上相差不多,但出现最大值的位置相差很大,这证明分段积分方法与连续积分方法的模拟结果是有较大差异的。
另外,由2m高度处湿度差值场可知源区草地退化后会使源区的空气湿度减小,空气变的干燥,不利于植被的生长,也不利于牲畜生活。
2.3.3能量变化
图2.9连续积分(a)和分段积分(b)模拟得到的源区草地退化前后区域一的感热通量月平均空间变化(单位:
W·m-2)
图2.10连续积分(a)和分段积分(b)模拟得到的源区草地退化前后区域二的感热通量月平均空间变化(单位:
W·m-2)
对于退化区域一的感热通量,连续积分得到的源区草地退化前后感热通量的空间变化(见图2.9(a))改变幅度在-16~8W·m-2其中,大部分地区感热通量减小,东部某些区域感热通量有小幅度增加;由分段积分得到的源区草地退化后的感热通量空间变化(见图2.9(b))可以看出,感热通量改变幅度在0~32W·m-2,其中西部地区增加的幅度较小,东南区域增加幅度大,增加的最大值为32W·m-2。
对于退化区域二的感热通量,从连续积分得到的源区草地退化前后感热通量空间变化(见图2.10(a))可以看出,感热通量增加幅度在8~16W·m-2;分段积分得到的源区草地退化后该区域的感热通量空间变化(见图3.10(b))范围在26~38W·m-2,增加幅度明显高于连续积分得到的结果。
图2.11连续积分(a)和分段积分(b)模拟得到的源区草地退化前后区域一的潜热通量月平均空间变化(单位:
W·m-2)
图2.12连续积分(a)和分段积分(b)模拟得到的源区草地退化前后区域二的潜热通量月平均空间变化(单位:
W·m-2)
对于退化区域一的潜热通量,由连续积分模拟得到的该区域潜热通量月平均空间变化减少,减小量在0~60W·m-2(见图2.11(a))。
分段积分模拟得到的该区域潜热通量明显减小,减小量在0~80W·m-2(见图2.11(b)),其中在东南区域减小幅度达到80W·m-2,大于连续积分得到的结果。
对于退化区域二的潜热通量,连续积分模拟得到的该区域潜热通量月平均空间变化幅度在64~80W·m-2(见图2.12(a)),分段积分模拟得到的潜热通量月平均空间变化范围在92~104W·m-2(见图2.12(b))大于连续积分模拟得到的结果。
由此可以看出,在能量通量的模拟方面,分段积分方法与连续积分方法模拟得到的结果存在较大的差异,有时甚至出现相反方向的变化。
从分段积分的模拟结果可知,源区草原退化后潜热通量的改变幅度比感热通量的改变幅度大,说明草原退化对潜热通量的影响高于对感热通量的影响,即对地气之间由于水发生相变引起的热量交换的影响要高于对地—气之间因湍流运动引起的热量交换的影响。
第三章结论与讨论
本文利用中尺度天气预报模式WRF,在设定的源区草地退化前后的条件下,分别利用连续积分方法与分段积分方法进行1个月的模拟,得到以下几点结论。
1)对于源区草地退化前的模拟,与传统积分方法相比,分段积分方法模拟得到的结果更接近真实值,由此可以确信在敏感性试验中采用分段积分方法得到的结论更可信。
2)分别采用分段积分方法和连续积分方法模拟得到的黄河源区草地退化前后2m高度处温度、地面温度、2m高度处湿度、以及潜热通量和感热通量等气象要素的变化存在差异,有些甚至出现相反的变化。
3)在全球增暖的大背景下黄河源区草地退化引起该区域平均2m高度处温度和地表温度升高,其中地表温度升高的幅度更大。
4)源区草地退化引起该区域空气湿度明显减小,空气变的干燥,不利于植被生长和牲畜生活。
5)源区草地退化引起该区域感热通量增加,潜热通量减小;感热通量增加量比潜热通量减小量小,源区草地退化对该区域的潜热影响更大。
总之,分段积分方法相比于传统积分方法能够有效地克服初始场误差和模式误差的积累,提高模式对气候变化的模拟能力,使得到的结果准确度更高,不确定性更小,保证了分段积分方法作敏感性试验比传统方法更可信。
分析结果表明,在敏感性试验中,分段积分与连续积分得到的结果存在较大的差异,由于分段积分方法对气候的模拟比连续积分方法更准确,我们有理由相信在敏感性试验中分段积分所得的结果可信度更高。
同