利用生物圈模型 (SiB2) 模拟青藏高原那曲草原近地面层能量收支

高志球, 卞林根, 程彦杰, 张雅斌, 王金星

高志球, 卞林根, 程彦杰, 等. 利用生物圈模型 (SiB2) 模拟青藏高原那曲草原近地面层能量收支. 应用气象学报, 2002, 13(2): 129-141.
引用本文: 高志球, 卞林根, 程彦杰, 等. 利用生物圈模型 (SiB2) 模拟青藏高原那曲草原近地面层能量收支. 应用气象学报, 2002, 13(2): 129-141.
Gao Zhiqiu, Bian Lingen, Cheng Yanjie, et al. Modeling of energy budget using simple biosphere model version 2(SiB2) over Tibetan Naqu prairie. J Appl Meteor Sci, 2002, 13(2): 129-141. .
Citation: Gao Zhiqiu, Bian Lingen, Cheng Yanjie, et al. Modeling of energy budget using simple biosphere model version 2(SiB2) over Tibetan Naqu prairie. J Appl Meteor Sci, 2002, 13(2): 129-141. .

利用生物圈模型 (SiB2) 模拟青藏高原那曲草原近地面层能量收支

资助项目: 

中国气象科学研究院博士后启动基金 

国家自然科学基金委员会预研究项目、国家教育部留学回国人员科研启动基金 

MODELINGOF ENERGY BUDGET USING SIMPLE BIOSPHERE MODEL VERSION 2(SiB2) OVER TIBETAN NAQU PRAIRIE

  • 摘要: 模式所需要的参数被合理地设置之后, 根据GAME/ Tibet (GEWEX亚洲季风试验/青藏高原试验) 那曲近地层观测站的资料, 将大气强迫变量代入SiB2(Simple Biosphere model version2), 文章模拟了该观测站地表能量收支。结果表明:SiB2能够较好地模拟青藏高原的能量收支情况, 净辐射、潜热通量和土壤热通量的模拟值和观测值吻合, 它们的相对误差分别为8% (低估)、6% (低估) 和3 %(低估)。同时, SiB2高估感热通量达40%。文章还给出了能量各分量的详细比较分析。
    Abstract: SiB2 (Simple Biosphere model version 2) is adopted to simulate surface energy budget in offline mode with certain atmospheric forcing values measured in the near surface observations taken over in the Tibetan short grass prairie as a part of the GAME/Tibet (DEWEX Asian Monsoon Experiment/Tibet) experiment. On the condition that all parameters are reasonably fixed, the simulated net radiation, latent heat flux and soil heat flux are shown to be basically in agreement with observations with their relative errors being 8% (underestimated), 6% (underestimated), and 3% (underestimated), respectively; in meantime, SiB2 overestimates sensible heat flux by 40%. The detail analysis about energy components is conducted.
  • 众所周知, 海洋和陆地是大气变化的两个重要外部强迫源。亚洲季风是由亚洲大陆和周边海域的热力、动力差异引起的。因此, 为了更好地研究季风和ENSO现象, 深入研究亚洲大陆尤其是青藏高原的能量、动量、水汽、二氧化碳和气溶胶的传输过程已经成为目前气候研究的重点之一[1~2]。过去的几十年, 随着青藏高原在单点和多点陆面过程研究的不断深入, 许多学者已取得了重要的研究成果。但是, 他们引用了传统公式:

    来计算感热和潜热, Cp为空气定压比热, ρ为空气密度, tsta分别为地表温度和某一高度上的气温, qs(ts) 和qa分别为地表饱和比湿 (地表温度为ts, 当然地表不总是处于饱和状态) 和该高度处的空气比湿, 则为对应高度处水平风速的绝对值, ChCq分别为湍流热量和水汽的交换系数。由于缺乏足够的测量数据, 同时, 关于大气和植被、植被和陆面、大气和土壤表层之间的能量、物质交换过程的物理细节缺乏更深入的认识, 难以取得一致的ChCq。这使得在关于青藏高原对亚洲季风和ENSO热力作用的认识上长期存在激烈争论[3]。因而, 有必要进一步应用模式来研究青藏高原的陆面过程。

    Sellers等的SiB2[4]以他们的SiB为基础, 应用能量、动量和质量守恒定律, 能够方便有效地模拟出土壤、大气和生物圈诸多参数以及陆面 (植被和裸地) 和大气之间的相互作用[5]。该模式已被嵌入其它区域和全球模式中[5]。在亚马逊河流域成功地模拟了陆面过程, 如感热、潜热通量, 地表温度, 空气动力学阻抗, 净同化率, 植被腔体传导率等等。

    为了估计亚洲季风期间青藏高原那曲短草草原表面的能量收支, 我们首先深入研究了SiB2所需要的全部参数, 进而利用1998年7月15~21日期间的观测资料, 将SiB2所需的6个大气强迫值 (向下短波辐射、向下长波辐射、水汽压、气温、风速和降水) 代入模式, 在此基础上, 模拟了净辐射、感热、潜热、地表热通量、地表有效辐射温度、土壤湿度等。

    Kim等对1998年GAME/Tibet (GEWEX亚洲季风试验/高原试验) 那曲近地层观测试验进行了详细的描述[6], 观测仪器的可靠性和方法在最近15年的试验中已得到了证实, 因而观测结果是可信的。该观测站位于那曲草原, 其地理位置为31°22.16′N、91°54.17′E, 那曲地处青藏高原中部, 属于高原亚寒带半湿润气候, 整个那曲地区植物高度一般小于5 cm。该观测站海拔高度为4496 m, 下垫面地势开阔平坦, 地表主要为沙粘土壤, 稀疏分布着细石块且不均匀地生长着高度为4~5 cm的草甸。测站具有良好的区域代表性。观测站东2 km处是连绵的山峰, 其平均海拔高度约为4650 m, 几乎没有植被覆盖; 南3 km处是青藏公路, 西南方向约6~7 km处是常年被积雪覆盖的群山; 西面则无大的起伏, 一望无际。包括三维超声风速仪在内的气象设备由韩国延世大学提供, 可以同时测量近地层大气4个高度 (3.5 m、2.2 m、1.3 m、0.7 m) 的风速、风向, 两层 (3.5 m、1.3m) 温度、湿度和水汽压; 地面以下两个深度 (0.015 m、0.04 m) 上土壤温度和深度为5 cm的土壤热通量, 近地层大气向上、向下的短波和长波辐射以及三维湍流风速资料。超声采样频率为20Hz。

    净辐射的4个分量, 即:向上、下短波辐射、向上、下长波辐射, 都是独立测量得到的, 一个CSAT3超声风速表结合五线热电偶和一个Krypton湿度表被安装在塔的2.85 m高处, 用涡旋相关法测量感热通量和潜热通量, 公式如下:

    (1)

    (2)

    我们用两个独立的土壤湿度表测量了地表到0.1 m深度的土壤水分; 测量了土壤5cm深处的热通量, 地面的土壤热通量计算是通过结合有限元方法和“强迫-恢复”方法来实现, 公式如下:

    (3)

    (4)

    发现两种方法计算结果一致。公式 (3) 和 (4) 中各符号的物理意义可见Garratt[7]和Stull[8]的文章。

    Sellers和他的合作者们用十多年时间设计了SiB2, 综合应用了地球生物圈的植被物理学、物候学和辐射传输的最新研究成果[4~5]

    用二流近似辐射传输模式[9], 计算了向上和向下的漫射辐射通量, 处理了通过最底层大气向下辐射通量的5个分量:可见光波段和近红外波段的直射和漫射辐射及热红外辐射。植被冠层、地表的反照率和净吸收辐射被参数化为太阳高度角和冠层与地表雪覆盖量的函数。

    SiB2中的净辐射通量由吸收的辐射减去发射的热辐射得到[5], SiB2把地球上的植被划分为9种类型, 并分别对其进行了详细述描, 已广泛使用在GCM研究中[10]。SiB2的基本方法是用类似电学上的欧姆定律来表示大气和地表之间的热量传输, 它分别给出了植被冠层之间, 冠层和地面之间空气感热、潜热通量的控制方程。

    (5)

    ΔXi的值为温度或水汽压, 分别对应于感热通量和潜热通量, 它们同时取决于植被冠层之间, 冠层和地面之间的总阻抗Σri, 如公式 (5) 所示。SiB2中计算感热通量的公式为:

    (6)

    其中Hg指裸土地面的感热通量, ; Hc是冠层部分的感热通量, , 因为观测站没有雪覆盖, 所以潜热通量λE在SiB2中可以表示为:

    (7)

    其中λEgs是裸土地表的潜热通量, , 且λEct是冠层的潜热通量, , SiB2中的净辐射计算公式表示为:

    (8)

    其中Rngs是裸土地表的净辐射, , 而Rnc是植被冠层的净辐射, 。在SiB2中, 用控制方程 (9~11) 分别计算植被、裸土地表和深土的温度 (Tc, Tg, Td),

    (9)

    (10)

    (11)

    在上述公式基础上, 植被热 (存储) 通量和裸土的热存储就可以计算出来, Sellers等[[4~5]和Randall等[10]详细描述了公式 (9~11) 中所有符号的物理意义。

    由公式 (9~11) 可知, 在SiB2中是用“强迫-恢复”方法计算地表、深土温度以及植被冠层的温度, 用冠层水储存计算降水截留。渗下的水储存到3层土壤中, 由上而下每层的厚度分别称为表面层、根区和补充层, 对应的厚度分别为0.02、0.13和0.85 m。水的传输在3层之间被描述为扩散过程, 植物的根穿透上两层, 水能够在补充层中发生重力排泄。

    SiB2用陆面-大气动力的相互作用, 将一层植被模型与植被微气象学、土壤消融过程、生长和生物量的新陈代谢联系在一起。

    SiB2中的大多数参数体现了植被的基本物理和生物特征, 如辐射传输模式中, 叶子的反射率和传输率, 土壤子模式中的水文参数等等。亚洲季风期间青藏高原那曲观测站的参数由实测和经验估计得到。考虑到实际测量的该区域土壤状况, 我们认为它属于类型6[5], 沙土占50 %, 卵石占10 %, 根占40 %, 土壤结构定为类型2[5]。一些近地层参数由测量得到, 因三维超声风速表 (CSAT 3) 装在高度为2.85 m处, 所以参考面也定在这一高度上。

    不随时间变化的植被参数包括:(i) 植被形态参数, 包括植物高度、冠层覆盖比、叶面角分布函数、根深度等等。本文所取的这些参数值都来自实际测量, 植被顶部和底部的高度分别为0.04 m和0.01 m。冠层覆盖比 (V) 和叶面指数 (LT) 在表层的能量和动量收支中也有重要的作用。植被覆盖比和叶面指数由相关的测量得到, 球形分布叶面角因子为0, 动量传输参数G1G4分别设置为1.449和11.785[5]。由于观测场地平坦, 平均地形坡设为0.04。(ii) 光特征参数是指叶子的传输率和反射率以及土壤的反射率。考虑到与波段和植物生命相关的4种情况:活着的植物在可见光波段, 活着的植物在近红外波段, 枯死的植物在可见光波段, 枯死的植物在近红外波段, 它们对应的叶层反射率和传输率分别设为0.105、0.58、0.36、0.58、0.07、0.25、0.22和0.38, 土壤的可见和近红外反射率分别设为0.11和0.225。所有的物理性质由参考文献[5]的表 5中6种类型得到。植物冠层项的最大核酮糖酶催化能力[5]为0.00003 mol·m-2·s-1, 光合作用内部量子效应参数为0.05 mol·m-2·s-1, 配对指数和叶子蒸腾分别为0.08, 0.95, 0.025, 根据植物腔体传导模型, 植物类型C4的腔体倾斜因子和最小腔体传导率分别设为4.0和0.4 mol·m-2·s-1, 温度阻碍因子S1S3S5分别表示光合作用中高温剪切因子、低温剪切因子和呼吸中高温剪切因子, 其值分别设为0.3、0.2和1.3 K-1。半阻碍温度参数S2S4S6也分别表示光合作用中高温剪切因子、低温剪切因子和呼吸中高温剪切因子, 它们的值分别设为313、288、328 K, 半阻碍水势参数为-200 m。

    随时间变化的参数通常由卫星资料得到[5], 我们由相应的测量计算并得到它们的估计值, 由绿色冠层从入射光中吸收的辐射比大约为0.95, 植被冠部的总叶面指数为0.4, 植被绿色部分 (N) 参数设为全绿色值 (1.0)。Z0为空气动力糙粗度, 取为0.013 m[11]。考虑到植被冠层高度, 植物的零平面位移 (D) 设为0.03 m。

    总体边界层阻抗系数C1, 土壤与植物层之间的阻抗系数C2的取值范围分别在1~16和10~1000之间。在计算感热和潜热通量时, 尤其是前者, C2的大小很重要。C1C2的取值, 对于大麦分别为9.3、278.0, 小麦分别为7.05、232.5, 玉米分别为15.1、131.0, 云杉分别为1.6、2450.0[5]。地表辐射温度的模拟表明, C1, C2分别为4.0和70.0。敏感性数值试验表明净辐射、感热通量、潜热通量和土壤热通量对C1, C2不敏感, 但是感热对C1C2的敏感程度比潜热大, C2C1重要, 为了单独研究C2对感热的作用, 在改变C2时, 其它量保持恒定。

    这个时期, 我们认为土壤组织结构为沙粘土。从参考文献[5]的表 4中查到土壤水势、饱和水传导率、土壤湿度指数和土壤空隙度, 分别为-0.07 m、35 m·s-1、4.9和0.44。

    总之, 我们结合卫星 (NOAA-14) 的资料和地面实测资料反演了地面的反射率, 利用实地的测量数据和美国全球植被覆盖数据, 确定了SiB2需要的所有参数。

    SiB2需要一些正确的初始值。它们是植被温度 (280.0 K), 地表温度 (282.0 K), 深土温度 (283.0 K), CAS (冠层空气) 温度 (282.0 K), 参照高度 (2.85 m) 上气温 (282.0 K)(测量得到), 土壤表面层湿度 (0.3), 根域湿度 (0.4), 交换域湿度 (0.4), 以上括号内值为同期内的观测值。

    SiB2需要6个强迫变量:向下短波辐射、向下长波辐射、水汽压、空气温度、测量高度上的水平风速和降水量。图 1给出了7月15~21日的强迫变量值, 向下短波辐射和向下长波辐射是重要的输入量。这一时期白天平均向下短波辐射为324.8 W·m-2, 夜间为-38 W·m-2, 最大峰值为1160 W ·m-2; 向下长波辐射几乎为一常数, 平均为328 W·m-2; 这一时期的平均水汽压为9.2 hPa, 并随天气发生变化; 平均空气温度为283.0 K, 也有较大的变化幅度。阴雨天气、晴好天气、整个时段的平均水平风速值分别为2.5 m·s-1、3.6 m·s-1、2.9 m·s-1, 在这7天, 降水亦很重要, 以上强迫值由测量得到, 如图 1所示。

    图  1  SiB2中的大气强迫变量时间变化曲线 (向下短波辐射 (实线, W·m -2), 地表向下长波辐射 (虚线, W·m -2), 水汽压 (hPa), 气温 (K) 及观测高度的水平风速 (m·s -1) 绝对值; 30 min时间步长内的降水量 (mm))

    我们选定7月15~21日一周时间为模拟周期, 其中16和17日是晴好天气, 18和19日是阴雨天气。将以上所述的所有参数、初始条件和强迫变量确定好之后, 我们用SiB2模拟了能量各分量, 并且分析了能量收支, 包括净辐射、感热通量、潜热通量、土壤热通量, 这些值也从观测资料中计算得到。同时模拟了地表温度、土壤湿度以及植被阻抗。

    模拟的净辐射有两个分量:植被冠层净辐射和地表净辐射, 它们在白天都有相同相位的明显变化, 图 3a给出了1998年7月15~21日那曲观测站植被冠层的净辐射通量和裸土地表的净辐射模拟结果。这一时期冠层平均净辐射在白天 (夜间) 为132.8 (-7.5) W·m-2, 白天 (夜间) 平均地表净辐射值为158.8(-23.0) W·m-2。平均模拟的白天 (夜间) 净辐射值为291.6(-29.0) W·m-2, 测量值为324.8(-37.0) W·m-2。由于夜间净辐射的主要部分为地表向上长波辐射, 而在夜间地表比植被冠层释放更多的热量, 图 2a表明, 夜间地表净辐射是净辐射分量中的重要部分是合理的。两个晴好的白天, 净辐射的两个峰值都出现在下午2 :30(北京时间, 下同), 白天 (从8 :00~20 :00[8]) 平均模拟峰值为777.0 W·m-2, 测量值为798.5 W·m-2。SiB2相对观测值低估了21.5 W·m-2, 占测量值的2.7 %, 这个值在仪器误差范围之内, 两个晴天的夜间 (从0 :00到8 :00, 从20 :00到24 :00), 净辐射的重要部分为地面上的长波辐射, 平均模拟净辐射为-43 W·m-2, 测量值为-54 W·m-2, 负号表示从地面向上, 二者绝对值差为11 W·m-2。对7天而言, 净辐射的平均绝对误差 (测量值减去模拟值) 为9.07 W·m-2

    图  2  1998年7月15~21日青藏高原那曲草原净辐射 (a)、感热通量 (b)、潜热通量 (c)、土壤热通量 (d) 及实测30 min时间步长内的降水量 (e) 模拟 (实线) 和测量 (点线) 结果
    图  3  模拟的1998年7月15~21日那曲观测站植被冠层 (实线) 和裸土地表 (虚线) 部分的净辐射 (a)、感热通量 (b)、潜热 (c)、土壤热通量 (d) 和实测30 min时间步长内的降水量 (e) 结果

    图 2b给出了模拟和计算的1998年7月15~21日那曲观测站的感热通量, 几乎所有模拟和测量的最大峰值都出现在下午2 :30, 7月16、17日两个晴天的平均最大模拟峰值为170.4 W·m-2, 观测峰值为132.4 W·m-2, 模拟值大于观测值, 其差额占观测值的28.7 %。在这两天的夜间, 平均模拟感热通量为-8.3 W · m-2, 测量值为-3.28 W·m-2, 负号表示低层温度低于参考高度的温度, 绝对值差为12.72 W·m-2。在7月17、18、19日3天, 两条曲线的形状和相位都非常一致, 几乎重合。SiB2高估7天中白天的感热通量而低估了夜间。图 2c中, 因为7月16、17日两个晴天模拟和测量的潜热通量较好的吻合, 所以我们认为主要冠层边界层的抗阻系数C1和地面到植被冠层腔体之间的空气动力学阻抗系数C2设置基本是正确的, 但图 2b表明SiB2高估了白天的感热通量。模拟感热通量有2个分量:植被感热通量和裸土地表感热通量见图 3b。他们白天有明显相同的相位变化, 在这7天中, 白天 (夜间) 平均植被冠层的感热通量为37.0(-16.6) W·m-2, 而白天 (夜间) 的平均裸土地表通量为40.3(8.5) W·m-2。白天 (夜间) 平均模拟感热通量值为77.4(-8.2) W·m-2, 测量的值为56.6(-3.1) W·m-2, 见图 2b图 3b中, 夜间植被冠层感热通量为负而地表感热通量为正值, 说明夜间植被冠层温度比空气温度低, 地面温度比空气温度高, 同时, 我们从图 3b中发现在总感热通量中, 冠层感热通量比裸土地表感热通量更重要。对7天而言, 感热通量的平均绝对误差 (测量值减去模拟值) 为-9.15 W·m-2

    同样地, 图 2c给出了1998年7月15~21日那曲观测站模拟和计算的潜热通量, 两条曲线的形状和相位都很一致, 二者峰值都在下午2 :30, 16、17日两个晴天的平均模拟峰值为307 W·m-2, 观测峰值为301.08 W·m-2。模拟值大于实测量值, 差值为6 W·m-2, 占测量值的2 %。两个晴天的晚间, 平均模拟潜热通量为-4.8 W·m-2, 测量值为-7.8W·m-2, 负值说明下层湿度小于参照高度上的湿度, 绝对差值为3.0 W·m-2。7月15、16、17和21日模拟和测量的潜热吻合较好, 7月19和20日, 白天模拟的值比测量的值大。在整个7天内, SiB2低估了夜间的潜热通量, 尽管测量值不是十分可信。观测的潜热分布 (图 2c) 体现了7月18日降水对潜热的影响, 但是模拟效果却没有体现该次降水对潜热的影响, 而7月19日降水对潜热的影响已经在模拟和测量的潜热分布中都得到体现。我们尚不明其原因。

    模拟潜热通量有2个分量:植被冠层潜热通量和裸土地表潜热通量 (见图 3c)。二者在白天都有明显同相位的变化, 白天 (夜间) 平均的冠层潜热通量为95.3(9.9) W·m-2, 平均裸土地表潜热通量为45.6(5.8) W·m-2, 平均模拟潜热通量为140.8 (15.7) W·m-2, 而观测值为145.1(20.1) W·m-2。由图 3c可以看到, 白天和夜间植被冠层潜热通量比裸土地表潜热通量更为重要, 虽然观测点植被冠层低矮稀疏, 但白天冠层潜热通量接近土壤潜热通量的2倍。对7天而言, 潜热通量的平均绝对误差 (测量值减去模拟值) 为4.73 W·m-2

    地表及浅层土壤热通量由许多因素决定, 包括太阳辐射、土壤物理性质、土壤水含量和植被覆盖率等。图 2d给出了1998年青藏高原模拟和测量的土壤热通量, 两条曲线的形状和相位十分相近, 16、17日两个晴天的峰值都出现在下午2 :30, 平均模拟峰值为246W·m-2, 测量峰值为257 W·m-2, 差值为11 W·m-2, 占观测值的4.3 %。两个晴好夜间 (从21 :00到07 :00), 平均模拟的土壤热通量为-66.9 W ·m-2, 测量值为-75.1 W·m-2, 负号表示热流量是从土壤深层到地表, 绝对差值为8.2 W·m-2, 占观测值的10.9 %。

    模拟土壤热通量有两个分量:裸土热储存和植被冠层的热储存 (见图 3d)。白天两者有同相位的显著变化, 白天 (夜间) 这一时期的平均植被冠层潜热通量为0.5(-0.8) W·m-2, 平均裸土潜热通量为72.8(-37.5) W·m-2图 2d中, 这7天白天 (夜间) 平均模拟的土壤热通量为73.36(-38.2) W·m-2, 测量值为61.2(-53.6) W·m-2。统计表明, 白天平均土壤热储存通量占土壤热通量的99 %以上。最大 (最小) 土壤热通量为353(-220) W·m-2, 最大 (最小) 土壤热储存通量为349(-212) W·m-2, 最大 (最小) 冠层热储存通量为11(-9) W·m-2。比较土壤热通量的两个部分表明, 裸土热储存通量占土壤热通量的99 %以上, 说明裸土在土壤热传输过程中有极其重要的作用。对7天而言, 土壤热通量的平均绝对误差 (测量值减去模拟值) 为-13.28 W·m-2

    能量收支描述为:

    (12)

    Rn为净辐射, H为感热通量, LE为潜热通量, G为土壤热通量, 观测资料很难使方程 (12) 闭合。图 4给出1998年7月中7天所观测的地表能量收支, HLE分别由涡旋相关法得到, G是地表处的土壤热 (存储) 通量, 净辐射中潜热通量、感热通量、土壤热通量和余差Rn-(H +LE +G) 的比例见图 4图 5给出了同期的地表能量收支的模拟结果。

    图  4  测量的1998年7月15~21日地表能量收支和余差Rn-(H +LE +G)
    图  5  1998年7月15~21日地表能量收支模拟结果

    总之, 在上述参数被合理确定之后, 将这7天由测量得到的大气强迫变量代入SiB2后, 该模式能较好地模拟出该测站的能量收支。我们进一步模拟了1998年7月2日~9月13日的能量收支, 7、8和9月中选定的3个7天时段平均能量收支, 如图 6所示。由图 6可以看出:3个7天时段的日平均最大测量净辐射值出现在8月; 9月净辐射观测和模拟值误差较大, 但模拟结果中, 7月份净辐射最大, 9月份最小, 该月模拟值误差达41 %。

    图  6  1998年7、8、9月3个7天周期的地表能量收支测量 (a) 和模拟 (b) 结果

    (1) 图 7给出了净辐射模拟值 (Rn, Model) 和同期向下短波辐射测量值 (Rswd) 的点聚图:其线性方程为Rn, Modeled=0.71 Rswd-31.0, 其相关系数为0.99, 表明可用该公式 (无需模式) 估算该地区净辐射, 模拟净辐射结果和这一时期测量净辐射结果的关系为Rn, Modeled=0.92Rn, Measured, 图略, 它表明SiB2低估了8 %的净辐射。从图 2a可以看出模式结果和测量值之间的误差在白天小于夜间, 夜间净辐射的主要部分是地表向上长波辐射, 这一误差在图 5中也可看到, 低估净辐射的原因可能来自SiB2处理通过最底层大气向下净辐射通量的5个主要部分时, 未充分考虑在青藏高原这种特殊的环境 (大气质量轻、大气透明度高) 中, 地形高度对辐射传输的影响。辐射方程中一些常数取决于大气中水汽和温度的分布, 不同季节测量值不相同[12]。事实上, 由于青藏高原的海拔高, 我们能够测量到非常大的辐射值, 有时总辐射超过太阳常数 (1367 W·m-2)。

    图  7  模拟的净辐射与测量得到向下短波辐射值间的拟合关系

    (2) 测量的鲍恩比是0.35, 而模拟值为0.46。

    (3) 图 4, 5分别给出了1998年7月15~21日能量收支的定量分析。平均测量净辐射为每天12.20 MJ·m-2, 模拟值为11.36 MJ·m-2, 绝对误差为0.84 MJ·m-2; 平均测量潜热通量为每天6.78 MJ·m-2, 模拟值为6.61 MJ·m-2, 绝对误差为0.17 MJ·m-2; 平均测量感热通量为每天2.31 MJ·m-2, 模拟值为3.01 MJ·m-2, 绝对误差为0.70 MJ·m-2; 土壤热通量为每天1.80 MJ·m-2, 模拟值为1.75 MJ·m-2, 绝对误差为0.05 MJ·m-2, 观测的余差Rn-(H +LE +G) 为2.63 MJ·m-2。造成误差的内在原因有两个方面:一是1km半径内平流和周围 (1 km半径) 山脉及山谷引起的局地环流作用可能引起10 %的误差, 二是仪器的精度允许10 %的误差。

    本文模拟的净辐射、潜热通量、土壤热通量和测量值吻合较好。他们的相对误差分别为8 %(低估), 6 %(低估) 和3 % (低估)。SiB2对感热高估达40 %。综上所述得出如下结论:

    (1) 在所有参数正确设置的情况下, SiB2能够很好地模拟出青藏高原那曲草原的净辐射、潜热通量和土壤热通量, 而对感热的模拟效果不好。

    (2) 在模拟能量分量及其相关物理量的过程中, 植被冠层的空气阻抗系数C2和土壤性质的确定至关重要。研究表明, 为了现实地估计能量收支, 对陆面参数进行深入研究是必要的。

    致谢: 与本文有关的研究工作得到周秀骥院士的指导, 特此致谢!
  • 图  1   SiB2中的大气强迫变量时间变化曲线 (向下短波辐射 (实线, W·m -2), 地表向下长波辐射 (虚线, W·m -2), 水汽压 (hPa), 气温 (K) 及观测高度的水平风速 (m·s -1) 绝对值; 30 min时间步长内的降水量 (mm))

    图  2   1998年7月15~21日青藏高原那曲草原净辐射 (a)、感热通量 (b)、潜热通量 (c)、土壤热通量 (d) 及实测30 min时间步长内的降水量 (e) 模拟 (实线) 和测量 (点线) 结果

    图  3   模拟的1998年7月15~21日那曲观测站植被冠层 (实线) 和裸土地表 (虚线) 部分的净辐射 (a)、感热通量 (b)、潜热 (c)、土壤热通量 (d) 和实测30 min时间步长内的降水量 (e) 结果

    图  4   测量的1998年7月15~21日地表能量收支和余差Rn-(H +LE +G)

    图  5   1998年7月15~21日地表能量收支模拟结果

    图  6   1998年7、8、9月3个7天周期的地表能量收支测量 (a) 和模拟 (b) 结果

    图  7   模拟的净辐射与测量得到向下短波辐射值间的拟合关系

  • Krishnamurti T N, Ramanathan Y.Sensitivity of the monsoon onset to different heating.Mon.Wea.Rev., 1982, 39:1090-1306.

    Murakami T, Chen L X, Xie A.Relationship among seasonal cycles, low frequency oscillations and transient disturbances as revealed from OLR data.Mon.Wea.Rev., 1986, 114:1456-1465. DOI: 10.1175/1520-0493(1986)114<1456:RASCLF>2.0.CO;2

    Chen L X, Reiter E, Feng Z Q.The atmospheric heat sources over the Tibetan plateau:May-Agugst 1979.Mon.Wea.Rev., 1985, 113:1771-1790. DOI: 10.1175/1520-0493(1985)113<1771:TAHSOT>2.0.CO;2

    Sellers P J, Randall D A, Collatz G J, et al.A revised land surface parameterization (SiB2) for atmospheric GCMs.Part I:Model formulation.Journal of Climate, 1996a, 9:676-705. DOI: 10.1175/1520-0442(1996)009<0676:ARLSPF>2.0.CO;2

    Sellers P J, Los S O, Tucker C J, et al.A revised land surface parameterization (SiB2) for atmospheric GCMs.Part II:The generation of global fields of terrestrial biophysical parameters from satellite data.Journal of Climate, 1996b, 9:706-737. DOI: 10.1175/1520-0442(1996)009<0706:ARLSPF>2.0.CO;2

    Kim J, Choi T, Lee H, et al.Energy partitioning and its imbalance over a prairie site in central Tibetan plateau during GAME-IOP 1998.Journal of Meteorological Society of Japan, 2001, (accepted).

    Garratt J R.Extreme maximum land surface temperatures.Journal of Applied Meteorology, 1992, 31:1097.

    Stull R B.Introduction to Boundary Layer Meteorology.London:Kluwer Academic Publishes, 1988.287.

    Dickinson R E.Land surface processes and climate-surface albedos and energy balance.Advances in Geophysics, 1983, 48.

    Randall D A, Dazlich D A, Zhang C.A revised land surface parameterization (SiB2) for atmospheric GCMs.Part III:The greening of the Colorado State University general circulation model.Journal of Climate, 1996, 9:738-763. DOI: 10.1175/1520-0442(1996)009<0738:ARLSPF>2.0.CO;2

    Gao Z, Wang J, Ma Y, et al.Study of roughness lengths and drag coefficients over Nansha Sea Region, Gobi, Desert, Oasis and Tibetan Plateau.Phys.Chem.Earth (B), 2000, 25 (2):141-145. DOI: 10.1016/S1464-1909(99)00126-4

    Monteith J L.Principles of Environmental Physics.New York:American Elsvier Publishing Company, Inc, 1973.39-52.

图(7)
计量
  • 文章访问数:  4534
  • HTML全文浏览量:  584
  • PDF下载量:  1936
  • 被引次数: 0
出版历程
  • 收稿日期:  2000-03-02
  • 修回日期:  2000-08-19
  • 纸刊出版:  2002-04-29

目录

/

返回文章
返回