半干旱区不同下垫面大气湍流通量比较分析
文章快速检索     高级检索
  高原气象  2017, Vol. 36 Issue (5): 1325-1335  DOI: 10.7522/j.issn.1000-0534.2016.00101
0

引用本文 [复制中英文]

陈丽晶, 张镭, 梁捷宁, 等. 2017. 半干旱区不同下垫面大气湍流通量比较分析[J]. 高原气象, 36(5): 1325-1335. DOI: 10.7522/j.issn.1000-0534.2016.00101
[复制中文]
Chen Lijing, Zhang Lei, Liang Jiening, et al. 2017. Analysis of Turbulent Flux on Different Underlying Surfaces over Semi-arid Area[J]. Plateau Meteorology, 36(5): 1325-1335. DOI: 10.7522/j.issn.1000-0534.2016.00101.
[复制英文]

资助项目

国家自然科学基金项目(41475008);国家重大科学研究计划项目(2012CB955302);兰州大学中央高校基本科研业务费专项资金(lzujbky-2015-9)

通讯作者

张镭.E-mail:zhanglei@lzu.edu.cn

作者简介

陈丽晶(1991-), 女, 福建福州人, 硕士研究生, 主要从事大气边界层和大气湍流研究, E-mail:chenlj10@lzu.edu.cn

文章历史

收稿日期: 2016-03-29
定稿日期: 2016-09-27
半干旱区不同下垫面大气湍流通量比较分析
陈丽晶, 张镭, 梁捷宁, 周旭     
兰州大学大气科学学院, 半干旱气候变化教育部重点实验室, 兰州 730000
摘要: 为了更深入地认识不同条件下的陆-气相互作用,选取2007年1月至2008年12月兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL)与中科院吉林通榆观测站的湍流资料,并结合常规气象观测资料,利用Swinbank于1951年提出的涡动相关法分析三种下垫面上的湍流通量,从而定量地描述和认识陆-气相互作用,并比较分析半干旱区不同下垫面热通量和动量通量的日变化与季节变化特征。结果表明:(1)SACOL、草原站和农田站上感、潜热通量的日变化均表现为单峰型,其中感热通量Hs日峰值可分别达到140.3,157.3和144.8 W·m-2;潜热通量LvE日峰值为Hs的40%~75%。(2)夜间热通量主要集中于-22.6~24.2 W·m-2,夜间Hs一年四季均表现为负值,而LvE几乎保持正值。(3)各测站上动量通量τ具有明显日变化,但由于风向、风速差异,通榆地区τ的日峰值可达0.25 kg·m-1·s-2,而SACOL仅有0.08 kg·m-1·s-2。(4)通榆站上τ的日变化具有显著季节差异,春、秋季的日较差远大于其他季。(5)排除热力因子影响,黄土高原下垫面对湍流活动的削弱作用更强,在该下垫面上近地层大气更容易也更快趋于稳定。
关键词: 陆-气相互作用    湍流通量    涡动相关    SACOL    吉林通榆    
1 引言

陆地作为地球大气的下边界, 与其在不同的时空尺度上相互影响、相互作用(鲍婧, 2012)。研究陆-气相互作用, 对更合理地解释陆-气间的各种物质、能量交换过程, 更深入地认识地球大气, 更全面地把握其运动规律, 进而更好地进行天气预报, 更准确地分析全球气候的变化趋势具有重要意义(Yang et al, 2012; Zuo et al, 2016)。为充分认识全球陆-气相互作用过程, 研究不同区域内陆-气相互作用是非常必要的。干旱半干旱区面积广大, 约占全球陆地总面积的41%(李宏宇等, 2015), 且该区域下垫面的不均一性和特殊的水热交换特征(Wang and Huang, 2010), 使其对气候变化的响应较为敏感, 从而成为全球气候变化研究中的关键地区(陈继伟, 2013; Narisma et al, 2007), 该区域的陆-气相互作用过程在全球陆地能量、水汽循环和气候变化中具有举足轻重的作用。陆地与大气之间通过边界层所进行的物质、能量交换与输送作为陆-气相互作用过程中非常重要的一环, 已受到广泛关注(Jiang et al, 2013), 在各种尺度的数值模拟研究中, 湍流通量的模拟也是数值模式物理框架中极为关键的组成部分(Garratt, 1993)。通过分析大气湍流通量可以加深对不同类型下垫面与大气间相互作用的定量理解, 进一步改进模式中的参数化方案, 完善陆面模式、发展气候模式, 提高天气预报和气候变化预测的准确率。因此, 研究分析大气湍流通量具有重要的理论意义和应用价值。

早在20世纪80年代中后期, 国际地圈-生物圈计划(International Geosphere-Biosphere Programme, IGBP)和世界气候研究计划(World Climate Research Programme, WCRP)等组织已在全球范围内开始进行陆面过程的研究, 目前相关实验已超过50项(王介民, 1999)。最近几十年, 我国也相继开展了大量陆-气相互作用的观测试验, 如“黑河地区地气相互作用野外观测试验研究”(HEIFE)(胡隐樵等, 1994; 张强和胡隐樵, 1995; Wang and Mitsuta, 1991, 1992)、“亚洲季风试验-青藏高原试验”(GAME-Tibet)(刘辉志和洪钟祥, 2000; 张强等, 2009)、“淮河流域水分和能量循环实验”(HUBEX)(徐敏等, 2002)、“内蒙古半干旱草原土壤-植被-大气相互作用”(IMGRASS)(吕达仁等, 2002)以及“中国西北干旱区陆-气相互作用”(NWC-ALIEX)(鲍艳和吕世华, 2006; Zhang et al, 2005)等。这一系列的实验获得了相当有价值的数据资料和研究成果, 为我国陆面过程的研究打下了坚实的理论和实践基础。王介民等(1990)胡隐樵等(1990)发现干旱区戈壁或沙漠下垫面能量输送主要以感热为主, 潜热可忽略不计; 牛国跃和王介民(1992)在SiB框架的基础上, 改进系统主要过程的参数化, 模拟了黑河实验区各种能量交换过程; 黄荣辉等(2013)总结了NWC-ALIEX经过连续12年的观测和多次加强观测所取得的干旱区陆面过程参数、边界层和地-气相互作用特征等方面的研究, 并对干旱区陆面过程参数化方案的发展、陆面过程模式的改进、中国西北干旱区感热输送特征以及西北干旱区地-气相互作用对中国东部气候的影响及其机理进行了总结, 并揭示了中国西北干旱区春、夏季具有高感热输送特征及此高感热对中国东部夏季气候变异的重要影响。

目前, 对于半干旱区陆-气相互作用的认识仍有不足, 现有的陆面过程模式和各参数化方案在半干旱地区的模拟效果误差较大。且已开展的多项陆面过程的观测试验, 多集中于湿润、半湿润、干旱区和青藏高原地区, 半干旱区的观测站很少, 能提供连续观测的站点更少。兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University, SACOL)和中科院吉林通榆观测站是位于我国半干旱区的两个重要测站(涂钢, 2007), 为揭示半干旱区陆-气交换特征、改进陆面过程参数化方案等方面研究提供了珍贵的观测资料。刘辉志等(2006)Jiang et al(2007)利用通榆站资料, 研究了退化草地和农田上的地-气相互作用, 发现在非生长季两种不同下垫面上的有效能量分配均以感热通量为主, 在生长季潜热通量与感热通量相当。鲍婧(2012)使用SACOL的湍流数据, 分析了黄土高原上近地层的陆-气相互作用, 发现动量通量和感热通量均有明显的日变化, 且动量通量始终由上向下输送, 感热通量在白天是由下向上输送, 到了夜间则由上向下输送。张强和赵鸣(1998)对西北干旱区非均匀下垫面的陆-气相互作用做了大量研究; 左金清等(2010)Zuo et al(2009)梁捷宁等(2013a, 2013b, 2014)和Liang et al(2014)深入研究了黄土高原复杂地形上大气边界层的湍流特征, 分析了低空急流对湍流活动的作用, 并从能量不闭合角度讨论了其对CO2通量的影响。

然而, 已有的工作主要对半干旱区各测站进行单站资料分析, 缺乏对不同类型下垫面间地-气作用特征的比较及其物理机制的解析, 而这种比较分析是深入认识半干旱区陆-气相互作用所必须的。本文利用2007年1月至2008年12月SACOL、通榆草原站和农田站的湍流资料, 研究半干旱区不同下垫面上动量通量与热通量(感热、潜热)的日变化与季节变化, 深入分析不同下垫面上大气湍流对物质、能量输送的物理机制, 讨论不同植被类型对地面湍流通量交换的影响, 进一步深入研究我国半干旱区的陆-气相互作用, 为改进模式中的参数化方案, 完善陆面模式提供参考。

2 观测站与仪器 2.1 SACOL

SACOL始建于2005年, 设在兰州大学榆中校区萃英山顶(35.946°N, 104.137°E), 海拔1965.8 m。地处中国黄土高原半干旱区, 属温带大陆性气候, 年平均气温6.7 ℃, 年平均降雨量381.8 mm, 相对湿度为63%, 全年盛行西北风和东南风, 平均风速1.6 m·s-1, 全年日照时数2607.2 h。下垫面为典型的黄土高原残塬地貌, 沟壑纵横, 土壤为第四纪黄土风蚀形成的灰钙土。植被基本为原生类, 以长芒草、冷蒿和赖草等短小植被为主, 冬季地表植被高约0.1 m, 夏、秋季覆盖率在80%左右。SACOL拥有国际先进的观测仪器, 是由我国自主建设的半干旱区长期观测站, 已被批准加入CEOP, 并作为此计划的全球协同加强观测站之一(Huang et al, 2008)。SACOL地面通量相关的观测项目主要有近地层气象要素、地面辐射通量、地面湍流通量以及土壤热通量等。本文所用资料主要是涡动相关通量。

2.2 吉林通榆观测站

中科院吉林通榆观测站建于2002年10月, 位于吉林省白城市通榆县新华乡境内(44°34′N, 122°53′E), 平均海拔184 m, 地处冲积平原, 地形平坦, 下垫面主要为草原植被, 现已基本被开垦、沙化, 土地退化和盐碱化严重。土壤主要有风沙土、淡黑钙土、盐碱土和草甸土, 土壤质地较粗, 有机质含量低, 风蚀严重。该地属温带大陆性季风气候, 年平均气温为5.2 ℃, 年平均降水量为404.3 mm, 相对湿度为57%, 全年主导西南风, 平均风速4.2 m·s-1(刘辉志等, 2000, 2006; 涂刚, 2007)。

观测站由分别在农田和草原的两个观测点组成(下称农田站、草原站), 相距约5 km。农田站下垫面在非生长季为裸土, 生长季主要作物为玉米, 作物最高达1.8 m。草原站下垫面为退化草地, 夏季草最高不超过0.1 m, 覆盖率约60%;冬季草高约0.5 m, 覆盖率约40%。与SACOL相同, 通榆观测站的湍流观测系统由超声风速温度仪(CSAT3)、湿度和CO2脉动仪(LI-COR7500) 组成, 可测动量、热量、水汽及CO2通量, 样本的取样时间均为10 Hz, 且所有仪器均采用128 W太阳能电池供电。不同的是仪器的架设高度, SACOL上仪器架设高度为2.88 m, 农田站为2 m, 草原站为2.5 m(刘辉志等, 2000, 2006; 涂刚, 2007)。

3 数据分析方法 3.1 涡动相关法

涡动相关法(Eddy Covariance Technique, EC), 由澳大利亚科学家Swinbank(1951)于1951提出, 通过直接测量大气各种属性的湍流脉动值(如脉动风速、脉动温度等)来计算湍流通量。该法建立在其所依据的物理原理上, 而不是建立在经验之上, 或通过其它物理量间接推论而来。

根据涡动相关法, 某一物理变量S的垂直通量F可以用下式来表述:

$ F = \frac{1}{T}\int_0^T {ws{\rm{d}}t}, $ (1)

式中: w为垂直风速; s为物理量S在平均时间T内的脉动。采用Reynolds分解, 将$w = \bar w + w', s = \bar s + s'$, 其中$\bar w$$\bar s$分别表示垂直风速和物理量S的平均态, w′、s′分别表示垂直风速和物理量S的瞬时扰动值。则:

$ F = \overline {ws} + \overline {w's'}, $ (2)

假设下垫面均匀平坦, 则$\bar w$, 垂直平流项$\overline {ws} = 0$, 湍流通量简化为:

$ F = \overline {w's'}, $ (3)

根据上式, 感热通量Hs、潜热通量LvE和动量通量τ可表示为:

$ Hs = \rho {c_p}\overline {w'{{T'}_s}}, $ (4)
$ {L_v}E = {L_v}\overline {w'{{\rho '}_v}}, $ (5)
$ \tau = \rho u_ * ^2, $ (6)

式中: ${u_ * } = {\left({{{\overline {u'w'} }^2} + {{\overline {v'w'} }^2}} \right)^{\frac{1}{4}}}$为摩擦速度; ρ为空气密度(单位: kg·m-3); ρv为水汽浓度(单位: kg·m-3); u′、v′、w′为三维脉动风速; Ts(K)为实际空气的脉动超声虚温; cp=1004.64×(1+0.84·q)为湿空气的定压比热容(单位: J·kg-1·K-1); Lv=(2.5008-0.0023668T)·106为凝结潜热系数(单位: J·kg-1)。

3.2 数据处理方法

尽管涡动相关通量观测仪器(特别是传感器结构)及其野外应用已有了较大的进步, 通量的计算仍需进行一些修正。

3.2.1 对原始记录进行预处理

u′, v′, w′, Ts, q′等10~20 Hz的原始观测数据分别记为x, 以5 min为移动窗口, 计算逐一窗口的平均值x及标准偏差SDx。如果某时刻x满足$\left| {x - \bar x} \right| \ge 3.25 \cdot S{D_x}$, 且连续满足此条件点的个数小于4, 则将其视为“野点”, 并剔除, 再用其相邻两侧的测量值做线性内插进行插补(Vickers and Mahrt, 1997)。

3.2.2 超声虚温修正

超声风速温度仪测得的温度为超声虚温Ts, 受湿度影响, 计算温度和感热通量时应考虑湿度订正。Schotanus et al(1983)提出超声虚温的修正方法, Foken(2011)对其进行了简化, 结果如下:

$ {T_s} = T \cdot \left({1 + 0.51q} \right), $ (7)
$ \overline {w'T'} = \overline {w'{{T'}_s}} - 0.51 \cdot \bar T \cdot \overline {w'q'} . $ (8)
3.2.3 坐标旋转

三维超声风速温度仪观测的坐标系为超声坐标系(Sonic coordinates), 当下垫面为倾斜地形或者仪器安装发生倾斜时, 水平风速或者平行于倾斜地面的风速便会在垂直方向上产生分量, 严重影响垂直风速的观测(王国印, 2013)。为了消除这种影响, 需要将超声坐标系转换为平均流线型坐标系(Tilt correction)(Finnigan, 1983; Kaimal and Finnigan, 1994; Tanner and Thurtell, 1969; Wilczak, 2001), 即坐标旋转修正。方法主要有二次坐标旋转(Double Roation, DR)、三次坐标旋转(Triple Roation, TR)和平面拟合(Planner Fit, PF)。Kainmal and Finnigan(1994)指出TR在计算应力时误差较大, 故不推荐使用。王国印(2013)利用SACOL涡动相关通量观测数据, 发现平面拟合坐标旋转对感热通量和潜热通量的影响较小, 故本文使用的是平面拟合法。

4 结果分析 4.1 热通量

在SACOL、草原站和农田站所处的三种不同类型下垫面(黄土高原、退化草地和农田)上, 热通量的日变化受太阳辐射影响, 一年四季均表现为单峰型, 且感热通量Hs和潜热通量LvE呈同位相变化(图 1~5)。日出后(06:00-08:00, 地方时, 下同), 太阳辐射增强, 近地层产生较大的垂直温度梯度, 促进湍流发展, 热通量不断加强, 至午时达到最大。SACOL、草原站和农田站上Hs的年平均日峰值可分别达到140.3, 157.3和144.8 W·m-2, 为LvE的1~2.5倍, 与李宏宇等(2010)利用2008年7月SACOL的观测资料发现“从感热和潜热积分值之比得到的波恩比平均是1.95”的结论相近, 说明白天能量分配以Hs为主。之后, 由于太阳辐射的减少以及湍流运动产生的热量扩散交换, 垂直温度梯度减小, 热量传递效率降低。16:00-18:00, 温度随高度的递增抑制湍流活动, 通量不断减小。夜间, HsLvE均稳定在零值附近。从季节上看, 热通量的日变化幅度和日峰值随季节变化明显, SACOL和农田站上Hs的日变化幅度表现为春季>夏季>秋季>冬季(图 2), 此种季节变化差异与总辐射日变化和净辐射日变化随季节的变化情况一致(图 5)。相较于2007年, 2008年春季SACOL上净辐射的增加和相对湿度减少, 导致近地层垂直温度梯度增大, 使得Hs的日峰值有明显增加; 而2008年夏季和秋季通榆草原站和农田站Hs日峰值的减少, 则是由净辐射的减少和相对湿度的增加引起的, 说明太阳辐射和相对湿度对Hs的影响密切。由于LvE受降水影响显著, 因此在多雨的夏季, SACOL、草原站和农田站上LvE的日变化幅度最大(图 3, 图 6), 与涂钢(2007)对通榆地区研究所得出“潜热通量日变化夏季最强”的结论一致, SACOL。草原站和农田站上夏季LvE的日峰值分别可达195.4, 211.4和167.5 W·m-2。通榆草原站和农田站的降水主要集中于68月, 而SACOL主要集中于810月(图 6), 因此SACOL秋季多雨的特征导致该季节SACOL上的LvE日峰值远大于后两者(图 3, 图 6)。与2007年同期相比, 2008年春季SACOL上冻土消融、积雪融化导致LvE出现显著增加; 而2008年夏季SACOL上净辐射的减少导致水分的蒸发量减少、植被的蒸腾作用减弱, 从而使得LvE比2007年夏季有所减少。2008年夏、秋季, 草原和农田上的降水量均大于2007年同期, 但受下垫面影响, 农田上茂盛的农作物能增加土壤含水量, 因而导致2008年夏、秋季农田上的LvE出现显著增加, 而退化草地上有所减少。

图 1 2007年(a)和2008年(b)SACOL、通榆草原站和农田站热通量的年平均日变化 Figure 1 Annual averagely diurnal variation of heat fluxes at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 2 2007年(a)和2008年(b)SACOL、通榆草原站和农田站感热通量的季平均日变化 Figure 2 Seasonal averagely diurnal variation of sensible heat fluxes at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 3 2007年(a)和2008年(b)SACOL、通榆草原站和农田站潜热通量的季平均日变化 Figure 3 Seasonal averagely diurnal variation of latent heat fluxes at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 4 2007年(a)和2008(b)年SACOL、通榆草原站和农田站各辐射量的年平均日变化 Figure 4 Annual averagely diurnal variation of radiation at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 5 2007年(a, c)和2008年(b, d)SACOL、通榆草原站和农田站总辐射(a, b)和净辐射(c, d)的季平均日变化 Figure 5 Seasonal averagely diurnal variation of global radiation (a, b) and net radiation (c, d) at SACOL, grassland and cropland station of Tongyu in 2007 (a, c) and 2008 (b, d)
图 6 2007-2008年SACOL、通榆草原站和农田站月平均气温(a)、月平均相对湿度(b)和月累计降水量(c) Figure 6 Monthly average of temperature (a), relative humid-ity (b) and monthly accumulation of precipitation (c) at SACOL, grassland and cropland station of Tongyu from 2007 to 2008

从SACOL、草原站和农田站日间热通量的季均值上看, 降水分布的季节性差异导致LvE受季节变化影响比Hs更显著(图 6, 图 7)。在多雨的夏季, 降水时释放的凝结潜热使各测站上的LvE均出现明显增加。相较于2007年, 2008年SACOL上降水的季节分布更为均匀, 因此LvE的季节变化也更平缓。而SACOL所在的榆中地区秋季多雨, 且通榆地区全年降雨的89%集中在69月(图 6), 导致SACOL上秋季LvE平均值约为通榆地区的2倍(图 3, 图 7)。此外, 草原站和农田站地表植被的分布影响局地的降水分布, 通榆地区的降雨具有很强的局地性(刘辉志等, 2006), 导致两测站的LvE在夏、秋季出现较大差异(图 7)。

图 7 2007年(a)和2008年(b)SACOL、通榆草原站和农田站日间热通量的季均值 Figure 7 Seasonal average of heat fluxes at SACOL, grassland and cropland station of Tongyu during the daytime in 2007 (a) and 2008 (b)

夜间, 大气层结稳定, 湍流活动不明显, 热通量主要集中于-22.6~24.2 W·m-2(图 1, 图 8)。其中, Hs在各季均表现为负值(图 8), 说明夜间贴地层逆温普遍存在, 热量始终由大气向地表传输, 且地球大气对陆面的保温作用在退化草地上更明显, 并在秋季达到最强。这是因为在无植被覆盖的裸土地上, 夜间地表温度下降得更快, 从而在近地面形成更大的垂直温度梯度, 导致热量更多地由高空向下传递。而LvE很大程度上受制于地表水分条件(李宏宇等, 2010), 由于生长季雨水丰富, 植物生长旺盛, 且茂密的植被增强了下垫面的储水能力, 夜间水分通过植物蒸腾和表层土壤水分蒸发以水汽的形态由地表向上输送, 导致夜间LvE在生长季均为正值且数值较大。而在非生长季, 降水减少, 植物蒸腾作用减弱甚至消失, 夜间蒸散作用减弱, 夜间LvE数值有所减少(图 8)。冬、春季夜间通榆地区LvE出现负值, 可能与地面的霜冻现象有关(涂刚, 2007)。

图 8 2007年(a)和2008年(b)SACOL、通榆草原站和农田站夜间热通量季均值 Figure 8 Seasonal average of heat fluxes at SACOL, grassland and cropland station of Tongyu at nighttime in 2007 (a) and 2008 (b)
4.2 动量通量

与SACOL相比, 草原站和农田站上动量通量τ的日变化更明显, 其日峰值和日较差远大于SACOL(图 9, 图 10)。观测期内, 通榆地区τ的年平均日峰值可达0.25 kg·m-1·s-2, 而SACOL仅有0.08 kg·m-1·s-2; τ的年平均日较差在通榆达0.17 kg·m-1·s-2, SACOL为0.05 kg·m-1·s-2。从前述热通量的结果看, 两地的Hs相近, 说明在湍流发展过程中热力作用相近, 动力作用是造成两地τ差异的主要原因。而在与涂钢(2007)对通榆地区风速的研究结果对比中发现, 通榆地区τ与风速具有同位相的日变化与季变化, 说明近地面风速对τ有显著影响。通过对2007-2008年SACOL、通榆草原站和农田站各季平均风向风速的比较分析发现, 三个测站的年平均风速多集中在2~4 m·s-1, 但草原站和农田站上大于4 m·s-1的风速分别占37.90%和33.88%, 而SACOL上4 m·s-1以上风速仅占21.17%; 6 m·s-1以上风速在草原站和农田站分别占16.45%和12.60%, 在SACOL仅为6.15%(图 11)。可见, 通榆地区的大风频率远高于SACOL。另一方面, 虽然通榆地区的主导风向为西南风, 但西北风、西风、南风、东南风和东风的出现频率较为接近, 风速风向的改变有利于湍流的形成及湍流通量的输送; 而SACOL以西北风为主, 其他方向风较少, 这种情形下更有利于物质和能量的平流输送(图 11)。因此, 风向风速上的巨大差异是导致两地τ产生如此大差异的主要原因之一。此外, 下垫面类型、植被类型和空气动力学粗糙度亦是影响τ的重要因素, 其对τ的具体影响程度需要做进一步分析。通榆地区τ的日变化季节性差异也更显著, 春、秋、冬季多变的风向导致这三个季节里τ的日变化幅度远大于夏季; 且其日峰值在秋季最大, 高达0.31 kg·m-1·s-2(图 10)。

图 9 2007年(a)和2008年(b)SACOL、通榆草原站和农田站动量通量的年平均日变化 Figure 9 Annual averagely diurnal variation of momentum at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 10 2007年(a)和2008年(b)SACOL、通榆草原站和农田站动量通量的季平均日变化 Figure 10 Seasonal averagely diurnal variation of momentum at SACOL, grassland and cropland station of Tongyu in 2007 (a) and 2008 (b)
图 11 2007-2008年草原站(左)、农田站(中)和SACOL站(右)各季风向风速图 (a)春季, (b)夏季, (c)秋季, (d)冬季 Figure 11 Seasonal average of wind direction and wind speed in grassland station (left), cropland station (medium), and SACOL (right) from 2007 to 2008. (a) spring, (b) summer, (c) autumn, (d) winter

夜间, 各测站上的τ均下降至日间的40%~50%, 但仍表现为退化草地>农田>SACOL, 且草原站和农田站上τ在夜间保持平稳, 而在SACOL上则不断减小(图 9, 图 10, 图 12)。由于白天, τ主要受热力因子控制, 夜间则受动力因子影响, 说明排除热力因子作用, 黄土高原下垫面对湍流活动的削弱作用较强, 在该种类型下垫面上近地层大气更容易也更快趋于稳定。τ在SACOL、草原站和农田站上具有显著的季节变化, 由于秋、冬季SACOL地面风速较小, τ在SACOL上表现为春季>夏季>秋季>冬季; 而秋、冬季通榆地区地面风速较大、风向变化剧烈, 草原站和农田站上τ较夏季出现大幅增加(图 11, 图 12)。因此, SACOL与通榆站上的τ在秋、冬季差别最大。而退化草地上的τ高于农田, 则是因为农田区周围零星的防护林导致退化草原和农田下垫面近地面层的风特性不同, 年平均风速可差2 m·s-1以上(刘辉志等, 2006)。但受风速限制, SACOL上τ的季节变化远小于通榆(图 10, 图 12)。

图 12 2007年(a)和2008年(b)SACOL、通榆草原站和农田站日间和夜间动量通量季平均 Figure 12 Seasonal average of momentum at SACOL, grassland and cropland station of Tongyu during the daytime and nighttime in 2007 (a) and 2008 (b)
5 结论

利用2007年1月至2008年12月SACOL和吉林通榆草地站、农田站的湍流资料, 并结合常规气象观测资料, 比较分析了半干旱区不同下垫面上感热通量Hs、潜热通量LvE和动量通量τ的日变化及季节变化特征。得到以下主要结论:

(1) SACOL、草原站和农田站上HsLvE的日变化在一年四季中均表现为典型的单峰结构。其中, Hs日峰值可分别达到140.3, 157.3和144.8 W·m-2; LvE的日峰值为Hs的40%~75%。夜间, 热通量较日间均大幅减少, 并平稳地维持在-22.6~24.2 W·m-2, 且Hs为负值, LvE为正值。

(2) 降水分布的季节性差异导致LvE受季节变化影响比Hs更显著。夏季, 随着降水的增加, 各测站上LvE也明显增加。而榆中地区秋季多雨, 导致了SACOL上秋季LvE平均值约为通榆地区的2倍。

(3) 受地面风向和风速的影响, 三种下垫面上τ的日变化幅度表现为退化草地>农田>SACOL, 通榆地区τ的日峰值可达0.25 kg·m-1·s-2。而SACOL以西北风为主, 年平均风速小、且年变化幅度不大, 一年四季τ均保持在0.03~0.06 kg·m-1·s-2间。通榆τ的日变化季节性差异更大, 春、秋季的日较差远大于其他季节。且相较于通榆站夜间τ的平稳变化, SACOL上τ的不断减小说明黄土高原下垫面对近地层大气的湍流运动的削弱作用更明显。

致谢 衷心感谢兰州大学半干旱气候与环境观测站(SACOL)与吉林通榆“干旱化和有序人类活动”长期观测野外试验站提供的宝贵数据资料, 以及各位审稿专家提出的宝贵意见和建议, 特此也对各位编辑的认真校对表示感谢。
参考文献
Finnigan J J. 1983. A streamline coordinate system for distorted turbulent shear flows[J]. Journal of Fluid Mechanics, 130: 241–258. DOI:10.1017/S002211208300107X
Foken T. 2011. Eddy covariance[M]. New York: Springer.
Garratt J R. 1993. Sensitivity of climate simulations to land surface and atmospheric boundary-layer treatments-A Review[J]. J Climate, 6(3): 419–449. DOI:10.1175/1520-0442(1993)006<0419:SOCSTL>2.0.CO;2
Huang J P, Zhang W, Zuo J Q, et al. 2008. An overview of the semi-arid climate and environment research observatory over the Loess Plateau[J]. Adv Atmos Sci, 25(6): 906–921. DOI:10.1007/s00376-008-0906-7
Jiang H M, Liu S H, Liu H P. 2013. Spectra and cospectra of turbulence in an internal boundary layer over a heterogeneously irrigated cotton field[J]. Acta Meteor Sinica, 27(2): 233–248. DOI:10.1007/s13351-013-0208-6
Jiang J F, Yan X D, Huang Y, et al. 2007. Simulation of CO2 and sensible/latent heat fluxes exchange between land surface and atmosphere over cropland and grassland in semi-arid region, China[J]. Journal of Forestry Research, 18(2): 114–118. DOI:10.1007/s11676-007-0022-z
Kaimal J C, Finnigan J J. 1994. Atmospheric boundary layer flows:Their structure and measurement[M]. Oxford: Oxford University Press.
Liang J N, Zhang L, Wang Y, et al. 2014. Turbulence regimes and the validity of similarity theory in the stable boundary layer over complex terrain of the Loess Plateau, China[J]. J Geophys Res, 119: 6009–6021. DOI:10.1002/2014JD021510
Narisma G T, Foley J A, Licker R, et al. 2007. Abrupt changes in rainfall during the twentieth century[M]. American Geophysical Union: , L06710.
Schotanus P, Nieuwstadt F T M, De Bruin H A R. 1983. Temperature measurement with a sonic anemometer and its application to heat and moisture fluxes[J]. Bound Layer Meteorology, 26(1): 81–93. DOI:10.1007/BF00164332
Swinbank W C. 1951. The measurement of vertical transport of heat and water vapor by eddies in the lower atmosphere[J]. J Meteor, 8: 135–145. DOI:10.1175/1520-0469(1951)008<0135:TMOVTO>2.0.CO;2
Tanner C B, Thurtell G W.1969.Anemoclinometer measurements of reynolds stress and heat transport in the atmospheric surface layer[M].University of Wisconsin Tech.Rep.ECOM-66-G22-F, 82.
Vickers D, Mahrt L. 1997. Quality control and flux sampling problems for tower and aircraft data[J]. J Atmos Ocean Technol, 14: 512–526. DOI:10.1175/1520-0426(1997)014<0512:QCAFSP>2.0.CO;2
Wang G Y, Huang J P. 2010. Observation analysis of land-atmosphere interactions over the Loess Plateau of northwest China[J]. J Geophys Res, 115(17): 1–15. DOI:10.1029/2009JD013372
Wang J M, Mitsuta Y. 1991. Turbulence structure and transfer characteristics in the surface layer of the HEIFE Gobi area[J]. J Meteor Soc Japan(69): 587–593.
Wang J M, Mitsuta Y. 1992. Evaporation from the desert-A preliminary result of HEIFE[J]. Bound Lay Meteor, 59(4): 413–418. DOI:10.1007/BF02215461
Wilczak J M, Oncley S P, Stage S A. 2001. Sonic anemometer tilt correction algorithms[J]. Bound Lay Meteor, 99: 127–150. DOI:10.1023/A:1018966204465
Yang Q D, Zuo H C, Xiao X, et al. 2012. Modelling the effects of plastic mulch on water, heat and CO2 fluxes over cropland in an arid region[J]. J Hydrol, 452-453: 102–118. DOI:10.1016/j.jhydrol.2012.05.041
Zhang Q, Huang R H, Wang S, et al. 2005. NWC-ALIEX and its research advances[J]. Adv Earth Sci, 20(4): 60–74.
Zuo H C, Chen B L, Wang S X, et al. 2016. Observational study on complementary relationship between pan evaporation and actual evapotranspiration and its variation with pan type[J]. Agricultural and Forest Meteorology, 222: 1–9. DOI:10.1016/j.agrformet.2016.03.002
Zuo J Q, Huang J P, Wang J M, et al. 2009. Surface turbulent flux measurements over the Loess Plateau for a semi-arid climate change study[J]. Adv Atmos Sci, 26(4): 679–691. DOI:10.1007/s00376-009-8188-2
鲍婧. 2012. 黄土高原半干旱区近地层陆-气相互作用研究[D]. 兰州: 兰州大学, 1. Bao Jing.2012.Study on land-atmosphere interaction in surface layer over semi-arid Loess Plateau of China[D].Lanzhou:Lanzhou University, 1.
鲍艳, 吕世华. 2006. 干旱、半干旱区陆-气相互作用的研究进展[J]. 中国沙漠, 26(3): 134–140. Bao Yan, Lü Shihua. 2006. Review of land-atmosphere interaction research in Arid and Semi-arid regions[J]. Journal of Desert Research, 26(3): 134–140.
陈继伟. 2013. 人类活动对干旱区农田地气相互作用的影响[D]. 兰州: 兰州大学, 1-2. Chen Jiwei.2013.The impact of human activities on land and atmosphere interaction over cropland in the arid area[D].Lanzhou:Lanzhou University, 1-2. http://www.cnki.com.cn/Article/CJFDTOTAL-HJKB201208016.htm
胡隐樵, 高由禧, 王介民, 等. 1994. 黑河实验(HEIFE)的一些研究成果[J]. 高原气象, 13(3): 225–236. Hu Yinqiao, Gao Youxi, Wang Jiemin, et al. 1994. Some achievents in scientific research during HEIFE[J]. Plateau Meteor, 13(3): 225–236.
胡隐樵, 奇跃进, 杨选利. 1990. 河西戈壁小气候和热量平衡特征的初步分析[J]. 高原气象, 9(2): 113–119. Hu Yinqiao, Qi Yuejin, Yang Xuanli. 1990. Preliminary analysis about characteristics of microclimate and heat energy budget in Hexi Gobi(Huayin)[J]. Plateau Meteor, 9(2): 113–119.
黄荣辉, 周德刚, 陈文, 等. 2013. 关于中国西北干旱区陆-气相互作用及其对气候影响研究的最近进展[J]. 大气科学, 37(2): 189–210. DOI:10.3878/j.issn.1006-9895.2012.12303 Huang Ronghui, Zhou Degang, Chen Wen, et al. 2013. Recent prograss in studies of air-land interaction over the arid area of northwest China and its impact on climate[J]. Chinese J Atmos Sci, 37(2): 189–210. DOI:10.3878/j.issn.1006-9895.2012.12303
李宏宇, 符淙斌, 郭维栋, 等. 2015. 干旱区不同下垫面能量分配机理及对微气候反馈的研究[J]. 物理学报, 64(5): 059201(1-5). DOI:10.7498/aps.64.059201 Li Hongyu, Fu Zongbin, Guo Weidong, et al. 2015. Study of energy partitioning and its feedback on the micro climate over different surfaces in an arid zone[J]. Acta Phys Sinica, 64(5): 059201(1-5). DOI:10.7498/aps.64.059201
李宏宇, 张强, 王胜. 2010. 陇中黄土高原夏季陆面辐射和热量特征研究[J]. 地球科学进展, 25(10): 1070–1081. Li Hongyu, Zhang Qiang, Wang Sheng. 2010. Research on characteristics of land-surface radiation and heat budget over the Loess Plateau of central Gansu in summer[J]. Adv Earth Sci, 25(10): 1070–1081.
梁捷宁, 张镭, 鲍婧, 等. 2013a. 黄土高原复杂地形受中尺度运动影响的稳定边界层湍流特[J]. 大气科学, 37(1): 113–123. DOI:10.3878/j.issn.1006-9895.2012.12004 Liang Jiening, Zhang Lei, Bao Jing, et al. 2013a. Turbulence influenced by mesoscale motions in the stable boundary layer over complex terrain of the Loess Plateau[J]. Chinese J Atmos Sci, 37(1): 113–123. DOI:10.3878/j.issn.1006-9895.2012.12004
梁捷宁, 张镭, 张武, 等. 2013b. 黄土高原半干旱区地表能量不闭合及其对二氧化碳通量的影响[J]. 物理学报, 62(9): 1–12. DOI:10.7498/aps.62.099203 Liang Jiening, Zhang Lei, Zhang Wu, et al. 2013b. Energy balance analysis over Loess Plateau and the consequences for carbon dioxide flux[J]. Acta Phys Sinica, 62(9): 1–12. DOI:10.7498/aps.62.099203
梁捷宁, 张镭, 田鹏飞, 等. 2014. 黄土高原复杂地形上边界层低空急流对近地层湍流的影响[J]. 地球物理学报, 57(5): 1387–1398. DOI:10.6038/cig20140504 Liang Jiening, Zhang Lei, Tian Pengfei, et al. 2014. Impact of low-level jets on turbulent in nocturnal boundary layer over complex terrain of the Loess Plateau[J]. Chinese J Geophys, 57(5): 1387–1398. DOI:10.6038/cig20140504
刘辉志, 洪钟祥. 2000. 青藏高原改则地区近地层湍流特征[J]. 大气科学, 24(3): 289–300. DOI:10.3878/j.issn.1006-9895.2000.03.01 Liu Huizhi, Hong Zhongxiang. 2000. Turbulent characteristics in the surface layer over Gerze area in the Tibetan Plateau[J]. Chinese J Atmos Sci, 24(3): 289–300. DOI:10.3878/j.issn.1006-9895.2000.03.01
刘辉志, 涂钢, 董文杰. 2006. 半干旱地区地气界面水汽和二氧化碳通量的日变化及季节变[J]. 大气科学, 30(1): 109–110. DOI:10.3878/j.issn.1006-9895.2006.01.09 Liu Huizhi, Tu Gang, Dong Wenjie. 2006. Seasonal and diurnal variations of the exchange of water vapor and CO2 between the land surface and atmosphere in the semi-arid area[J]. Chinese J Atmos Sci, 30(1): 109–110. DOI:10.3878/j.issn.1006-9895.2006.01.09
吕达仁, 陈佐忠, 陈家宜, 等. 2002. 内蒙古半干旱草原土壤-植被-大气相互作用(Imgrass)综合研究[J]. 地学前缘, 9(2): 295–306. Lü Daren, Chen Zuozhong, Chen Jiayi, et al. 2002. Composite study on inner Mongolia semi-arid grassland soil-vegetation-atmosphere interaction(IMAGRASS)[J]. Earth Science Frontiers, 9(2): 295–306.
牛国跃, 王介民. 1992. 简易一维陆面过程的数值模拟[J]. 高原气象, 11(4): 411–422. Niu Guoyue, Wang Jiemin. 1992. A numerical simulation of one-dimensional land-surface process model[J]. Plateau Meteor, 11(4): 411–422.
涂钢. 2007. 半干旱区不同下垫面地气相互作用研究[D]. 北京: 中国科学院大气物理研究所, 1-12. Tu Gang.2007.Study of land-atmosphere interaction on different underlying surface in semi-arid area[D].Beijing:The Institute of Atmospheric Physics, 1-12.
王国印. 2013. 半干旱区陆-气相互作用的观测与研究[D]. 兰州: 兰州大学, 42-48. Wang Guoyin.2013.Observation of land-atmosphere interaction over semi-arid regions[D].Lanzhou:Lanzhou University, 42-48.
王介民. 1999. 陆面过程试验和地气相互作用研究-从HEIFE到IMGRASS和GAME-Tibet/TIPEX[J]. 高原气象, 18(3): 280–294. Wang Jiemin. 1999. Land surface process experiments and interaction study in China-from HEIFE to IMGRASS and GAME-Tibet/TIPEX[J]. Plateau Meteor, 18(3): 280–294.
王介民, 刘晓虎, 祁勇强. 1990. 应用涡旋相关方法对戈壁地区湍流输送特征的初步研究[J]. 高原气象, 9(2): 120–129. Wang Jiemin, Liu Xiaohu, Qi Yongqiang. 1990. A preliminary study of turbulence transfer characteristics in Gobi area with an eddy correlation technique[J]. Plateau Meteor, 9(2): 120–129.
徐敏, 蒋维楣, 胡非, 等. 2002. 湍流湿度脉动特性研究-根据淮河流域实验资料[J]. 地球物理学报, 45(1): 17–25. Xu Min, Jiang Weimei, Hu Fei, et al. 2002. Study on the characteristics of turbulent moisture-based on the data obtained from Huaihe River basin experiment[J]. Chinese J Geophys, 45(1): 17–25.
张强, 胡隐樵. 1995. 热平流影响下湿润地表的通量-廓线关系[J]. 大气科学, 19(1): 8–20. DOI:10.3878/j.issn.1006-9895.1995.01.02 Zhang Qiang, Hu Yinqiao. 1995. The flux-profile relationship under the condition of heat advection over moist surface[J]. Chinese J Atmos Sci, 19(1): 8–20. DOI:10.3878/j.issn.1006-9895.1995.01.02
张强, 王胜, 张杰, 等. 2009. 干旱区陆面过程和大气边界层研究进展[J]. 地球科学进展, 24(11): 1185–1194. DOI:10.11867/j.issn.1001-8166.2009.11.1185 Zhang Qiang, Wang Sheng, Zhang Jie, et al. 2009. The prograss on land surface process and atmospheric boundary layer in arid regions[J]. Adv Earth Sci, 24(11): 1185–1194. DOI:10.11867/j.issn.1001-8166.2009.11.1185
张强, 赵鸣. 1998. 干旱区绿洲与荒漠相互作用下陆面特征的数值模拟[J]. 高原气象, 17(4): 1–8. Zhang Qiang, Zhao Ming. 1998. A numerical simulation of characteristics of land-surface process under interaction of oasis and desert in arid region[J]. Plateau Meteor, 17(4): 1–8.
左金清, 王介民, 黄建平, 等. 2010. 半干旱草地地表土壤热通量的计算及其对能量平衡的影响[J]. 高原气象, 29(4): 840–848. Zuo Jinqing, Wang Jiemin, Huang Jianping, et al. 2010. Estimaion of ground heat flux for a semi-arid grassland and its impact on the surface energy budget[J]. Plateau Meteor, 29(4): 840–848.
Analysis of Turbulent Flux on Different Underlying Surfaces over Semi-arid Area
CHEN Lijing , ZHANG Lei , LIANG Jiening , ZHOU Xu     
College of Atmospheric Science, Lanzhou University, Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, Lanzhou 730000, China
Abstract: In order to understand the land-atmosphere interaction between different conditions more deeply, we take eddy covariance technique proposed by Swinbank in 1951 primarily to analyse the turbulent flux which could describe the land-atmosphere interaction quantificationally on three types of underlying surfaces.We contrast the diurnal and seasonal variations of heat flux and momentum flux on loess plateau, grassland and cropland in semi-arid area by using turbulent data from January 2007 to December 2008, which comes from Semi-arid Climate and Environment Observatory of Lanzhou University (SACOL), the grassland and the cropland station in Tongyu.Combined with the conventional meteorological elements, the results show that:(1) At the SACOL, grassland station and cropland station of Tongyu, the diurnal variations of sensible and latent heat flux both show a typical "single-peak" structure.The peak value of sensible heat flux (Hs) in SACOL, grassland and cropland is 157.3 W·m-2, 144.8 W·m-2 and 140.3 W·m-2, which is the double of latent heat flux (LvE).(2) The heat flux mainly locate in -22.6~24.2 W·m-2 at night, and Hs is negative, LvE is positive.(3) The diurnal variation of the momentum flux(τ) shows the same law in the three station, the daily peak of τ in Tongyu is 0.25 kg·m-1·s-2, which in SACOL is 0.08 kg·m-1·s-2.(4) The diurnal variation of τ has an obvious difference in difference seasons, the daily range in spring and autumn is bigger than other seasons.(5) Except thermal factors, loess plateau has a stronger extinction effect on turbulent, atmosphere will be stable easier and quicker.
Key Words: Land-atmosphere interaction    Turbulent Flux    Eddy Covariance Technique    SACOL    Tongyu