中国气象学会主办。
文章信息
- 黄倩, 王蓉, 田文寿, 左洪超, 张强. 2014.
- HUANG Qian, WANG Rong, TIAN Wenshou, ZUO Hongchao, ZHANG Qiang. 2014.
- 风切变对边界层对流影响的大涡模拟研究
- Study of the impacts of wind shear on boundary layer convection based on the large eddy simulation
- 气象学报, 72(1): 100-115
- Acta Meteorologica Sinica, 72(1): 100-115.
- http://dx.doi.org/10.11676/qxxb2014.007
-
文章历史
- 收稿日期:2013-4-10
- 改回日期:2013-9-29
2. 中国气象局干旱气象研究所, 甘肃省干旱气候变化与减灾重点实验室, 兰州, 730020
2. Key Laboratory of Arid Climate Change and Reducing Disaster of Gansu Province, Institute of Arid Meteorology China Meteorological Administration, Lanzhou 730020, China
1 引言
对流边界层也称为混合层,它是对流层中直接受地表影响的最低一层大气。对流边界层中的空气由于热泡对流和风切变作用而充分混合。浮力和风切变是影响对流边界层湍流发展和物理量时空分布的重要机制,但是浮力和风切变对边界层发展的影响是不同的。Deardorff(1972)指出在极不稳定时,摩擦速度不再对湍流强度产生影响,并且给出了分别由浮力和风切变主导边界层气流运动的不同尺度参数。虽然在极不稳定时,磨擦速度不再影响湍流强度,但它对边界层对流形式的影响不能忽略。早期的研究(LeMone,1973)曾指出,地表热通量决定边界层对流的强度,而风切变有助于将边界层对流组织成不同的形式。Sykes等(1989)用摩擦速度(u*)和对流速度尺度(w*)之比的大小说明风切变对有组织的边界层对流形成的重要作用。风切变对边界层对流形式的影响也得到了其他研究者的证实(Grossman,1982;Weckwerth,et al,1997;Huang,et al,2009)。
对流边界层的发展对一些天气现象,如云的形成、沙尘暴的产生以及污染物的传输等有直接影响。除了地表热量向上传输加热边界层大气从而使边界层向上发展以外,对流边界层之上自由大气中浮力较强的暖空气向下混合进入边界层,以及来自边界层的上冲热泡形成的夹卷过程对对流边界层增长的贡献也不容忽略。尽管自由对流的大涡模拟试验结果表明,由于热泡与对流边界层之上自由大气的相互作用,在热泡中形成较强的涡旋,从而将自由大气界面向下折卷,导致逆温层中的暖空气卷入对流边界层是形成夹卷作用的重要机理(Sullivan,et al,1998),但是风切变对夹卷过程的影响也受到广泛地关注。Kim等(2003)认为在强风切变情况下,有组织的边界层对流增强了夹卷层中开尔文-亥姆霍兹波的形成,从而加强了夹卷层与自由大气的湍流交换,而且增大风切变使夹卷层显著增厚。Pino等(2003)和Conzemius等(2006)也都认为风切变能增强夹卷通量,从而有利于对流边界层厚度的抬升。近年来随着对风切变对夹卷热通量(Kim,et al,2006;Pino,et al,2006;Sun,et al,2009)、夹卷层厚度(袁仁明,2005;林恒等,2010)的参数化研究,对风切变对夹卷过程的影响有了更深入的了解。
通过以往的研究,虽然对风切变对有组织的边界层对流、边界层增厚及对流夹卷过程等的影响有一定的认识,但是对风切变对对流边界层以及夹卷层中不同高度的物理量,如位温和垂直速度等水平分布的影响,以及风切变对对流边界层中示踪物传输高度的影响尚不明确。本研究以敦煌干旱区为研究区域,进一步深入分析风切变对对流边界层物理量(位温、垂直速度和示踪物)水平分布的影响,研究风切变对示踪物传输高度的影响。这些研究结果将为今后进一步研究干旱和半干旱区边界层对流对沙尘抬升和传输高度的影响补充知识。 2 大涡模式及数值试验
利用“西北干旱区陆-气相互作用野外观测实验”2000年5—6月在甘肃敦煌进行的加强观测期间6月3日敦煌站的小球探空资料,以及6月3日敦煌双墩子戈壁站的地表热通量观测资料,采用英国气象局提供的大涡模式LEM version2.4(Gray,et al,2001)模拟研究干旱区边界层结构及演变特征,并利用敏感性数值试验诊断分析风切变对对流边界层发展及示踪物垂直传输的影响。
LEM version2.4是一个用于模拟湍流尺度和云尺度问题的三维高分辨率数值模式,其基本方程组为
其中,u是速度矢量,θ是位温,qn代表张量变量,p是气压,ρ是空气密度,B′是浮力,τ是次网格切应力,hθ是θ的次网格通量,hqn是qn的次网格通量,δi3是克罗内克函数,Ω是地球的角速度,εijk为交变单位张量,是与微物理过程有关的θ源项,代表与辐射有关的θ源项,是与微物理过程有关的qn源项。如果假设则式(3)可以写成 这里平均参考状态是由ρs、θs和ps来定义的。该模式利用滤波方法对纳维-斯托克斯方程组求网格体积平均,从而准确计算大湍流涡旋的运动,而对小尺度湍流通量和能量采用一阶湍流闭合的次网格模式模拟。模式中次网格切应力τij和标量通量hn,i可以由以下两式给出
υ和υh是湍流粘性系数,Sij是形变应力,其具体形式如下其中,Rip是理查逊数,S2=SijSij,λ是混合长,定义为这里,λ0是一个与粘滞耗散率有关的长度尺度,z0是粗糙度长度,κ是冯·卡曼常数。模式水平方向的侧边界条件是周期性的。上下边界采用钢性边界层条件,为了减少由于模式上边界引起的重力波反射,在距模式顶约2/3高度以上加入牛顿阻尼吸收层。模式下边界采用莫宁-奥布霍夫相似理论计算地表切应力及地表感热和潜热通量。另外,为了使边界层对流在模式中发展,在模式1 km高度以下对位温和水汽场分别进行了±0.2 K和±0.05 g/kg的随机扰动(Huang,et al,2010)。模式水平方向的模拟区域是5 km×5 km,分辨率是100 m,垂直方向的模拟高度是5 km,垂直方向采用张弛网格,最小垂直格距在近地面,取为2.7 m,最大格距在模式4 km高度以上,为150 m。模式中使用的地表地转风由美国国家环境预报中心/国家大气研究中心(NCEP/NCAR)2.5°×2.5°再分析资料计算得到,地转风切变由小球探空资料1 km高度的风速和地表地转风资料求得。模式在计算时采用有限差分方法,模拟时间是10 h,每隔0.5 h输出一次数据。
为了模拟风切变对对流边界层演变以及对位温、垂直风速等物理量分布的影响,设计了4个不同风切变的数值试验,试验的初始位温和水汽混合比廓线都取2000年6月3日07时(北京时,下同)的小球探空观测廓线,所用地表热通量都是6月3日敦煌双墩子戈壁站的观测资料。试验B是仅有浮力驱动边界层发展的试验,此试验中初始风场和地转风速都为0。试验BG、BT、BGT是浮力和风切变共同驱动边界层发展的试验,这些试验中的地表热通量与试验B中使用的相同,所不同的是BG试验中的地转风切变存在于地面向上100 m厚度的近地层内;而试验BT的地转风切变存在于1 km高度向下的100 m厚度的气层内;试验BGT中地转风切变存在于整个对流边界层内(图 1b中分别给出了试验BG、BT和BGT的地转风廓线)。为了研究风切变对对流边界层物质传输的影响,在所设计的4个试验中,从地面向上到100 m厚度的气层内加入绝对浓度为100的被动示踪物。
3 模拟结果及分析 3.1 风切变对边界层物理量水平分布的影响 3.1.1 风切变对对流边界层发展的影响图 1是2000年6月3日不同时次实测和模拟的位温廓线。为了考虑背景场的影响,大涡模式可以在模拟过程中将水平平均位温、水汽混合比和风速廓线以2 h的时间间隔收敛到实测探空廓线(Marsham,et al,2008;Huang,et al,2009)。图 1中红线代表以2 h的时间间隔把位温和水汽混合比收敛到实测值的模拟结果,而蓝线代表以2 h的时间间隔把位温、水汽混合比和风速同时收敛到实测值的模拟结果。另外,通过分析模拟的湍流动能随时间的演变特征,发现模拟初期模式的调整时间约为2.5 h,所以分析结果从10时开始。图 1的模拟结果和实测结果都显示,在12时约为500 m厚的对流边界层之上是厚度约1 km的逆温层,逆温层以上的残留层清晰可辨;14时由于混合层中的对流向上发展已逐渐侵蚀对流层之上逆温层,到16时对流边界层扩展到4 km的高度,对流边界层与残留层已经混合成一层深厚的大气边界层(张强等,2004)。Zhang等(2006)通过对敦煌干旱区实测探空资料的分析,指出辐射、地表特征和风切变等热力和动力因子是引起敦煌地区超厚边界层的主要原因,并深入研究了地表过程的热力特征对敦煌干旱区超厚边界层形成的贡献(Zhang,et al,2011)。本研究将从动力因子角度出发,进一步解释风切变对敦煌干旱区超厚边界层的影响。另外,模式同时考虑背景动力和热力强迫影响(将位温、水汽混合比和风速收敛到实测值)的模拟结果与实测值更接近,这说明大涡模式能够模拟出与实际观测一致的对流边界层的平均结构及演变特征。
为了便于比较有风切变和无风切变的敏感性试验结果,试验B、BG、BT和BGT都没有把模拟的位温、水汽混合比和风速平均廓线值收敛到实测探空廓线。对比试验B、BG、BT和BGT模拟的不同时次的平均位温和示踪物浓度随高度的变化(图 2)可以看出,13时浮力驱动的对流边界层厚度约为600 m,浮力和风切变共同驱动的对流边界层厚度约为700 m,随着时间推移对流边界层不断变暖,对流边界层顶的高度也逐渐升高,到16时试验B中的对流边界层的厚度为1.6 km,而试验BGT中的对流边界层的厚度达到2.5 km。通过分析得到在相同条件下,浮力驱动的对流边界层厚度比浮力和风切变共同驱动的对流边界层厚度小,而浮力和风切变共同驱动的对流边界层比浮力驱动的对流边界层升温更快,其中试验BGT中的对流边界层升温最快。通常边界层之上覆盖有逆温层,白天一方面由湍流运动将地表的热量向上输送加热边界层大气,另一方面混合层顶以上的夹卷过程将上部逆温层中的暖空气向下夹卷进入混合层也加热边界层大气。由于试验BG、BT和BGT与试验B的唯一区别是试验BG、BT、BGT中考虑了风切变对边界层发展的影响,而强风切变可以增加对流边界层的夹卷热通量(Moeng,et al,1994),因此有风切变的对流边界层(试验BG、BT、BGT)比仅有浮力驱动的对流边界层升温快,厚度也较大,而且对流边界层中的垂直混合运动更容易打通逆温层而与残留层混合成一层大气(张强等,2004)。
从图 2给出的平均示踪物浓度随高度变化的分布中看出,示踪物被传输的高度与对流边界层的发展密切相关,随着边界层内对流向上发展,对流边界层顶的高度不断升高,近地层的示踪物也被传输到较高的高度。其中,有切变影响的对流边界层中示踪物传输的高度比只有浮力驱动的对流边界层中示踪物传输得高,而且从图 2可以看出,随着风切变增大(试验BGT),夹卷过程增强,示踪物传输的高度更高。3.1.2 风切变对物理量水平分布的影响
图 3给出了模式模拟的边界层不同高度示踪物浓度、位温及垂直速度的概率密度函数(PDFs)分布。垂直速度的概率密度函数(图 3c、f、i)分布显示,无论是浮力驱动的边界层对流还是浮力和风切变共同驱动的边界层对流中,在边界层低层(0.3Zi,其中Zi代表对流边界层顶的高度),上升气流和下沉气流的强度基本一致,只是下沉气流略多于上升气流;而在边界层中上部(0.7Zi),下沉气流比上升气流多,但是上升气流强度比下沉气流的强;在边界层顶部(1.0Zi),这种下沉气流多而弱、上升气流少而强的分布特征更加明显。垂直速度的这些特征与Moeng等(1990)和LeMone(1990)的研究结果一致,即由地表热力驱动的边界层大气中,边界层底部具有不规则的上升运动,而只有少数较强的上升气流能达到边界层顶部。另外,从图 3c、f、i还能看出,在有风切变的个例中,边界层中上层的下沉气流增多而且其速度增大,而上升气流的强度减弱,这种现象在边界层顶部更明显。由位温的概率密度函数分布看出(图 3b、e、h),从0.3Zi到0.7Zi为正倾斜(概率密度函数的峰值在左侧)分布,即上升气流较下沉气流暖,到1.0Zi位温的概率密度函数分布趋于平缓,尤其是在考虑风切变的影响后,概率密度函数分布变得更平缓。这主要是因为受地表加热影响,边界层中下层上升气流为暖空气,而下沉气流是相对较冷的空气,到边界层顶部,由于夹卷作用,上部逆温层的暖空气被夹卷进入边界层从而也加热了边界层大气,因此靠近边界层顶的大气位温的概率密度函数分布相对更平缓(Couvreux,et al,2005;Huang,et al,2009)。示踪物绝对浓度的概率密度函数分布显示(图 3a、d、g),从0.3Zi到0.7Zi为负倾斜(概率密度函数的峰值在右侧),即上升气流对应的示踪物浓度较小,在1.0Zi处,由于上部不含示踪物空气的夹卷,其概率密度函数的分布也变得更平缓,这一特征在有风切变的试验(分别为试验BT和BGT)中更为突出。
从边界层顶以上1.1Zi、1.2Zi和1.3Zi处示踪物浓度、位温及垂直速度的概率密度函数分布(图 4)中可以看出(图 4c、f、i),浮力和风切变共同驱动的边界层对流中,上升气流和下沉气流的强度变化范围比浮力驱动的边界层对流的变化范围大,而且风切变越大上升气流和下沉气流的强度变化范围也越大。在1.1Zi,试验B模拟位温的概率密度函数呈正倾斜分布,而从1.2Zi到1.3Zi,位温的概率密度函数分布变为负倾斜(图 4b、e、h),也就是说在1.1Zi高度,上升气流比下沉气流暖,而从1.2Zi到1.3Zi高度下沉气流要比上升气流暖,这说明越靠近上部夹卷层,向下卷入的气流越暖。从图中容易看出,如果考虑风切变的影响,试验BG、BT和BGT中位温的概率密度函数分布更平缓些,其中,整层边界层都存在风切变时(试验BGT),位温的概率密度函数分布从1.1Zi到1.3Zi都较平坦。图 4a、d、g也显示考虑风切变的影响,边界层顶以上的示踪物浓度的概率密度函数分布相对平缓。
3.2 风切变对夹卷层的影响 3.2.1 热通量及湍流特征量随高度的变化
从不同模拟试验得到的动力学热通量廓线(图 5)可以看出,热通量在地面有最大值,约为0.33 K·m/s,随高度增加热通量线性递减,到混合层顶附近减小为0,再往上热通量出现负值区,而且在边界层顶热通量达到负的最大值。热通量的负值区通常称为夹卷层,夹卷层之上热通量为0。图 5还显示边界层风切变对夹卷层厚度和夹卷强度(用负的最大热通量值表示)有较大影响。在仅有浮力驱动的试验B中,夹卷层的厚度只有0.6Zi,其强度也只有-0.05 K·m/s左右;当近地层或边界层顶有风切变时(试验BG和BT),夹卷层的厚度和夹卷强度分别增大到1.0Zi和-0.1 K·m/s左右;当边界层整层存在风切变时(试验BGT),夹卷层厚度增大到1.6Zi,而夹卷强度增大到-0.18 K·m/s。另外,还可以注意到,当边界层整层存在风切变时,其混合层厚度比无风切变、近地层有风切变或边界层顶有风切变模拟试验中的都小。这是由于增大边界层风切变使边界层对流的强度减弱(Huang,et al,2009),即上冲热泡的速度减小,从而抑制混合层向上发展,导致其厚度减小。图 5显示随着风切变增大,夹卷层厚度也增大,这是因为风切变使机械作用加强,进而导致夹卷作用增强和夹卷层增厚,图 5中显示的特征与Kim等(2003)的研究结果一致。
图 6a、b、c分别是无量纲速度方差(、,以下简称为水平速度方差,以下简称为垂直速度方差)的垂直分布,图 6d是无量纲位温方差,以下简称为位温方差)的垂直分布。其中w*和θ*分别是对流速度尺度和温度尺度,即
式中,g为重力加速度,Zi是对流边界层顶的高度,θv是虚位温,是近地层虚位温通量。图 6a、b显示,只有浮力驱动的模拟试验中(实线),无量纲水平速度方差在近地面有最大值(约0.6),而且随高度增加而递减;到边界层中部,水平速度方差随高度变化较小;到边界层顶部又略有增大,然后很快减小。当边界层中有风切变时(点线、虚线、点划线),模拟的水平速度方差明显大于无风切变时的水平速度方差。从图中还能看出,水平速度方差除了在近地面有极大值外,在边界层顶部又出现了一个峰值,然后随高度增加逐渐减小为0。由于受地表摩擦作用的影响,近地面的风切变往往较大,从而使水平速度方差在近地面达到最大。而在边界层顶附近,由于上部逆温层的覆盖,热泡的上升运动受到限制而转化为水平速度分量,因此水平速度方差在边界层 顶附近增大,这一特征在有风切变的试验中尤其明显。另外还注意到,在边界层顶附近,试验BGT中的水平速度方差比试验BG和BT中的大,而且顺风方向(x方向)的水平速度方差增加得更大。图 6c不同试验中的垂直速度方差随高度的变化基本一致,垂直速度方差都是从地表向上逐渐增大,到1/3Zi处达到最大,然后随高度增加而逐渐减小。其中有风切变的试验中垂直速度方差值比无风切变试验中的大。另外,与水平速度方差相比较,加入风切变后,水平速度方差增大的幅度更大,这也说明了风切变更有利于水平方向湍流的产生(Stull,1988)。从无量纲位温方差的垂直分布(图 6d)可以看出,位温方差随高度的增加而减小,在边界层中部达到最小值,之后又随着高度的增加而增大,在夹卷层出现极大值后又随高度减小。夹卷层中位温方差的增大是由于自由大气的暖夹卷空气与冷的上冲热泡之间的较大差别造成的。图 6d还显示试验B夹卷层的位温方差比试验BG、BT或BGT的小,这是因为试验B中夹卷热通量明显小于其他3个试验的夹卷热通量。另外,从试验BG、BT和BGT夹卷层的位温方差看出,虽然试验BGT的夹卷热通量最大,但是BGT夹卷层的位温方差却小于试验BG或BT的值,这也许与夹卷层中平均位温梯度的变化有关,如图 2b、c显示,试验BG与BT夹卷层的位温梯度比试验BGT的大。3.2.2 湍流动能的收支
通过对比分析浮力驱动边界层对流和浮力与风切变共同驱动边界层对流的试验结果,可以发现风切变对对流边界层物理属性有较显著的影响,特别是对夹卷层厚度及夹卷强度有重要影响,因此有必要进一步分析风切变对决定边界层发展的湍流动能收支及其垂直分布有何影响。一般在水平均匀、略去下沉的条件下,湍流动能收支方程可以表示为
其中,u′、v′、w′是脉动速度分量;U、V是平均风速的水平分量;p是气压; 0是标准密度; v是背景虚位温;e是湍流动能,根据Stull(1988)的定义可以表示为湍流动能收支方程式(15)中等号左边的项表示湍流动能的局地变化项,等号右边的项依次表示湍流动能的切变产生项(S)、浮力项(B)、湍流输送项(T)、气压传输项(P)和耗散项(D),其中,切变产生项和浮力项是湍流动能的源,耗散项是汇,而湍流传输项和气压传输项既不产生湍流动能也不消耗湍流动能,只在垂直方向上对湍流动能起再分配的作用。从图 7a看出,试验B中浮力项是湍流动能的源项,由式(15)知在准平衡状态下,当切变产生项为0时,浮力项、湍流传输项和气压传输项之和与耗散项是相平衡的。图 7a还显示,近地层以上湍流动能的耗散项随高度的变化较小,而且湍流传输项有随高度线性递增的变化趋势,这与浮力产生项随高度线性递减有很好的对应。这一结果与以往的观测和模拟研究结果都有较好的一致性(Lenschow,1974;Moeng,et al,1994;Pino,et al,2003)。图 7b是试验BG中湍流动能收支随高度的变化,可以看出,近地层的风切变不仅使湍流动能的切变产生项增加,而且还影响其他项的大小,如风切变对浮力项在边界层顶部的影响比底部和中部的影响更明显。从图 7b还能看出,切变产生项在近地面有最大值,然后随高度的增加而减少,在边界层中部减少到最小后基本不随高度变化,到边界层顶部又逐渐随高度增加,在夹卷层增大到最大,最后又随高度的增加而减小到0。由于模式的研究区域是敦煌干旱沙漠地区,较大的地表热通量形成尺度较大的热泡,热泡在垂直方向充分混合,从而减小了边界层中部的风切变,因此切变产生项在边界层中部较小。另外,与无风切变试验中的湍流动能的耗散项比较,当近地层加入风切变时耗散项在近地层和夹卷层都有所增大。因为湍流动能的耗散与小尺度湍涡有关(Stull,1988),这也进一步说明由于下垫面的阻挡,风切变在近地层产生的湍涡尺度较小,而且近地层风切变还能增强夹卷层内的小尺度湍流运动。图 7c显示,在1 km高度加入风切变后,湍流动能收支中各项的垂直分布与近地层加入风切变所对应的湍流动能收支基本一致,这说明近地层风切变对夹卷层湍流动能的贡献不容忽视。图 7d显示当边界层整层存在风切变时湍流动能的切变产生项在夹卷层明显增大,而且为了平衡湍流动能切变产生项的增大,耗散项也增大了。另外,湍流输送项和气压传输项在垂直方向有相互补偿的特点。通过比较有风切变和无风切变敏感性试验中湍流动能收支在垂直方向的分布,发现增大边界层风切变,夹卷层中湍流动能的产生也增强,而夹卷层中湍流活动的增强又有利于上部逆温层中暖空气的向下夹卷,因此在研究夹卷通量的参数化问题时正确引入风切变的影响有重要意义。
3.3 风切变对示踪物垂直传输的影响 3.3.1 垂直速度与示踪物浓度从不同模拟试验在不同时次的垂直速度及示踪物浓度的垂直分布结果(图 8)明显看出,仅有浮力驱动的对流边界层,对流泡在垂直方向发展的高度较高,有的对流泡可以从地面向上发展到整个对流边界层,如图 8a中水平方向-2 km处、图 8b中水平方向-0.5 km处、图 8c中水平方向-1.5 km和1.1 km处的对流泡。而在浮力和风切变共同驱动的对流边界层中,对流泡在垂直方向上发展时,很难从地表向上发展到整个对流边界层,尤其当边界层风切变增大时(如图 8j)这一现象更明显。另外,还注意到,有风切变时向上发展的对流泡发生倾斜,如图 8l中水平方向-1 km处的对流泡和水平方向0.5 km处的对流泡,这里对流泡的不同倾斜方向也许与对流泡中的不同涡旋方向有关。对流泡的倾斜与风切变有关,由图 6a、b的分析可知,风切变加强了边界层顶对流泡的垂直动量向水平动量的转化,因此倾斜的对流泡在边界层顶更容易将夹卷层空气向下混合,从而使对流边界层进一步增暖、变厚。这与Kim等(2003)研究风切变对夹卷层涡旋结构影响的结果一致,即风切变越大,边界层顶的涡旋运动也越强,从而导致较强的夹卷作用。
从图 8示踪物浓度的垂直分布看出,随着边界层对流的向上发展,示踪物在垂直方向上传输的高度也在不断增加。比较试验B和试验BG中示踪物传输的高度发现,尽管风切变仅存在于近地层中,但示踪物传输的高度明显升高。而试验BT中示踪物传输的高度与试验BG的差别不大。当增大边界层风切变时(试验BGT),在12时30分,混合层顶之上就有一层明显增厚的存在示踪物的气层,如图 8j中垂直方向0.7—2 km,而在相同时刻试验B中这一层的厚度不足150 m。这也进一步证明了风切变有利于夹卷运动,从而促进了混合层与残留层的贯通,使示踪物在垂直方向被传输到更高的高度。3.3.2 示踪物在垂直方向的概率分布
从示踪物浓度在垂直方向的概率密度函数(PDFs)分布(图 9)可以看出,在只有浮力驱动的边界层对流试验中(图 9a),示踪物浓度在垂直方向的概率密度函数分布较对称,即随着高度增加示踪物浓度的概率密度函数没有发生倾斜,浓度的变化范围为7.5—12.1,到边界层顶附近浓度迅速递减,而且示踪物传输的最大高度约为1.3 km。如果在边界层底部加入风切变(图 9b),示踪物浓度的概率密度函数从地面到500 m高度为垂直递减的分布,从500 m到边界层顶,示踪物浓度的概率密度函数分布不仅随高度递减,还随高度发生倾斜,即浓度随高度逐渐减小,而且示踪物传输的最大高度大约是1.7 km。图 9c显示的示踪物传输的最大高度与图 9b的相近,只是概率密度函数分布从近地层向上就开始倾斜,而且在混合层内同一高度浓度的变化也略小于图 9b的结果。当边界层整层存在风切变时(图 9d),示踪物传输的最大高度可以达到2 km,其浓度的概率密度函数分布在垂直方向倾斜最明显,即示踪物浓度随高度递减更迅速,而且在同一高度浓度变化范围也较大。通过比较图 9中仅有浮力驱动的边界层对流与浮力和风切变共同驱动的边界层对流中示踪物浓度的概率密度函数垂直分布可以发现,浮力作用使示踪物浓度的概率密度函数分布在垂直方向不发生倾斜,这与浮力产生的较大尺度对流泡的上冲密切相关。
为了进一步解释增大风切变示踪物浓度随高度递减的原因,又计算了边界层顶示踪物浓度的垂直通量(图 10),可以看出,试验BG、BT和BGT边界层顶示踪物浓度的垂直通量都比试验B的大,即风切变加强了示踪物在边界层顶的向上输送,从而导致边界层中上部示踪物浓度随高度递减。试验BGT显示示踪物浓度的概率密度函数分布在垂直方向倾斜最明显,随高度递减最快,这与图 10d中边界层顶示踪物浓度向上通量最大有很好的一致性。
4 结 论本研究利用大涡模式模拟研究了干旱区下垫面上空对流边界层的演变,通过对比分析了浮力驱动的与浮力和风切变共同驱动的边界层结构、物理属性的差异,给出了风切变对边界层对流及示踪物垂直传输的影响特征及机理。模拟和诊断分析结果表明:
(1)浮力和风切变共同驱动的对流边界层较仅有浮力驱动的对流边界层增暖快,厚度也较大,而且有风切变的边界层对流更容易打通逆温层之上的残留层。在边界层下部,上升气流和下沉气流分布较对称;而在边界层中上层,上升气流的区域虽然少于下沉气流,但是上升气流的强度较大。风切变不仅使下沉气流的区域变大,还使下沉气流的强度增强。另外,由于边界层顶的夹卷作用,位温和示踪物浓度的概率密度函数分布趋于平缓,这一特征在加入风切变的试验中更突出。
(2)在仅有浮力驱动边界层发展的条件下,湍流动能收支中湍流传输项随高度线性递增,浮力产生项随高度线性递减。在相同地表热通量的条件下,增大风切变,夹卷层厚度和夹卷强度都增大。在边界层顶,增大风切变有利于热泡的垂直动量转化为水平动量。另外,增大风切变使夹卷层中湍流动能增强,同时夹卷层中湍流动能的耗散也增大,说明夹卷作用主要是由小尺度湍涡运动所致。
(3)在仅有浮力驱动的对流边界层中,边界层对流的强度较相同条件下有风切变的大,而且对流泡发展的高度也较高。有风切变的情况下,边界层对流泡发生倾斜,这种倾斜的对流泡与垂直的对流泡相比,在边界层顶更容易将上部逆温层中的暖空气向下夹卷混合,从而使边界层进一步增暖,增厚。因此示踪物在浮力和风切变共同驱动的边界层中比浮力驱动的边界层中传输的高度更高,但是示踪物浓度随高度递减,而浮力驱动的边界层中对流泡的强烈垂直混合及夹卷活动相对较弱,示踪物浓度的概率密度函数分布在垂直方向不发生倾斜。
林恒, 孙鉴泞, 卢伟. 2010. 有切变对流边界层夹卷厚度参数化的大涡模拟研究. 南京大学学报(自然科学版), 46(6): 616-624 |
袁仁明. 2005. 白天混合层顶部夹卷层厚度的特征研究. 地球物理学报, 48(1): 19-24 |
张强, 卫国安, 侯平. 2004. 初夏敦煌荒漠戈壁大气边界结构特征的一次观测研究. 高原气象, 23(5): 587-596 |
Conzemius R, Fedorovich E. 2006. Dynamics of sheared convective boundary layer entrainment. Part I: Methodological background and large-eddy simulations. J Atmos Sci, 63(4): 1179-1199 |
Couvreux F, Guichard F, Redelsperger J L, et al. 2005. Water-vapour variability within a convective boundary-layer assessed by large-eddy simulations and IHOP_2002 observations. Quart J Roy Meteor Soc, 131(611): 2665-2693 |
Deardorff J W. 1972. Numerical investigation of neutral and unstable planetary boundary layers. J Atmos Sci, 29(1): 91-115 |
Gray M E B, Petch J, Derbyshire S H, et al. 2001. Version 2.3 of the Met Office large eddy model. The Met Office, Exeter, UK |
Grossman R L. 1982. An analysis of vertical velocity spectra obtained in the BOMEX fair-weather, trade-wind boundary layer. Bound Layer Meteor, 23(3): 323-357 |
Huang Q, Marsham J H, Parker D J, et al. 2009. A comparison of roll and nonroll convection and the subsequent deepening moist convection: An LEM case study based on SCMS data. Mon Wea Rev, 137(1): 350-365 |
Huang Q, Marsham J H, Parker D J, et al. 2010. Simulations of the effects of surface heat flux anomalies on stratification, convective growth, and vertical transport within the Sharan boundary layer. J Geophys Res, 115(D5): doi:10.1029JD012689 |
Kim S W, Park S U, Moeng C H. 2003. Entrainment processes in the convective boundary layer with varying wind shear. Bound Layer Meteor, 108(2): 221-245 |
Kim S W, Park S U, Pino D, et al. 2006. Parameterization of entrainment in a sheared convective boundary layer using a first-order jump model. Bound Layer Meteor, 120(3): 455-475 |
LeMone M A. 1973. The structure and dynamics of horizontal roll vortices in the planetary boundary layer. J Atmos Sci, 30(6): 1077-1091 |
LeMone M A. 1990. Some observations of vertical velocity skewness in the convective planetary boundary layer. J Atmos Sci, 47(9): 1163-1169 |
Lenschow D H. 1974. Model of the height variation of the turbulence kinetic energy budget in the unstable planetary boundary layer. J Atmos Sci, 31(2): 465-474 |
Marsham J H, Parker D J, Grams C M, et al. 2008. Observations of mesoscale and boundary-layer scale circulations affecting dust transport and uplift over the Sahara. Atmos Chem Phys, 8(23): 6979-6993 |
Moeng C H, Rotunno R. 1990. Vertical-velocity skewness in the buoyancy-driven boundary layer. J Atmos Sci, 47(9): 1147-1162 |
Moeng C H, Sullivan P P. 1994. A comparison of shear-and buoyancy-driven planetary boundary layer flows. J Atmos Sci, 51(7): 999-1022 |
Pino D, Vilà-Guerau de Arellano J, Duynkerke P G. 2003. The contribution of shear to the evolution of a convective boundary layer. J Atmos Sci, 60(16): 1913-1926 |
Pino D, Vilà-Guerau de Arellano J, Kim S W. 2006. Representing sheared convective boundary layer by zeroth-and first-order-jump mixed-layer models: Large eddy simulation verification. J Appl Meteor Climatol, 45(9): 1224-1243 |
Sullivan P, Moeng C H, Stevens B, et al. 1998. Structure of the entrainment zone capping the convective atmospheric boundary layer. J Atmos Sci, 55(19): 3042-3064 |
Sun J N, Xu Q J. 2009. Parameterization of sheared convective entrainment in the first-order jump model: Evaluation through large-eddy simulation. Bound Layer Meteor, 132(2): 279-288 |
Stull R B. 1988. An Introduction to Boundary Layer Meteorology. Netherlands: Kluwer Academic Publisher, 666pp |
Sykes R I, Henn D S. 1989. Large-eddy simulation of turbulent sheared convection. J Atmos Sci, 46(8): 1106-1118 |
Weckwerth T M, Wilson J W, Wakimoto R M, et al. 1997. Horizontal convective rolls: Determining the environmental conditions supporting their existence and characteristics. Mon Wea Rev, 125(4): 505-526 |
Zhang Q, Wang S, Li Y Y. 2006. Study on physical mechanism of influence on atmospheric boundary layer depth in arid regions of northwest China. Acta Meteor Sinica, 20(5): 1-12 |
Zhang Q, Zhang J, Qiao J, et al. 2011. Relationship of atmospheric boundary layer depth with thermodynamic processes at the land surface in arid regions of China. Sci China: Earth Sci, 54(10): 1585-1594 |