气象学报  2021, Vol. 79 Issue (4): 659-673   PDF    
http://dx.doi.org/10.11676/qxxb2021.037
中国气象学会主办。
0

文章信息

张璐, 黄倩, 张宏昇, 张强, 田红瑛. 2021.
ZHANG Lu, HUANG Qian, ZHANG Hongsheng, ZHANG Qiang, TIAN Hongying. 2021.
干湿地表的湍流特征及其对深对流影响的大涡模拟
Large eddy simulation of turbulence effects on deep-convection triggering over dry and wet surfaces
气象学报, 79(4): 659-673.
Acta Meteorologica Sinica, 79(4): 659-673.
http://dx.doi.org/10.11676/qxxb2021.037

文章历史

2020-06-29 收稿
2021-04-09 改回
干湿地表的湍流特征及其对深对流影响的大涡模拟
张璐1,2 , 黄倩1 , 张宏昇2 , 张强3,4 , 田红瑛1     
1. 兰州大学大气科学学院,半干旱气候变化教育部重点实验室,兰州,730000;
2. 北京大学物理学院大气与海洋科学系,气候与海-气实验室,北京,100871;
3. 中国气象局兰州干旱气象研究所/甘肃省干旱气候变化与减灾重点实验室/中国气象局干旱气候变化与减灾重点开放实验室,兰州,730020;
4. 甘肃省气象局,兰州,730020
摘要: 利用大涡模式模拟了对流边界层结构演变以及深对流触发过程。通过改变鲍恩比的敏感性试验研究不同大气初始状况下湿润和干旱下垫面湍流特征及其对深对流触发过程的影响。结果表明:干旱下垫面的混合层干而暖,厚度较大;湿润下垫面相反。由于地表感热通量对热力湍流形成的作用更大,干旱下垫面上湍流混合和夹卷作用更强,使得水汽和相当位温在边界层内分布更均一,而在边界层顶有较大的负扰动;干旱下垫面上对流强度较湿润下垫面大,但均表现为泡状对流,水平方向上呈网状结构。不同下垫面上深对流的发生与大气初始状况有关,当初始时刻1—3 km的逆温强度较弱时(0.15 K/(100 m)),边界层内湍流迅速发展,深对流首先在干旱下垫面发生,但因对流有效位能较小,云层厚度小于湿润下垫面。当1—3 km的逆温强度增加到0.55 K/(100 m)时,云层形成时间较晚,云层厚度明显减小,仅当边界层顶的比湿较大时,有深对流发生,但仍首先发生在干旱下垫面,考虑贯穿对流在边界层顶引起的较强冷却作用,云层厚度大于湿润下垫面。
关键词: 干、湿地表    初始温、湿度廓线    湍流特征    深对流    大涡模拟    
Large eddy simulation of turbulence effects on deep-convection triggering over dry and wet surfaces
ZHANG Lu1,2 , HUANG Qian1 , ZHANG Hongsheng2 , ZHANG Qiang3,4 , TIAN Hongying1     
1. Key Laboratory for Semi-Arid Climate Change of the Ministry of Education,College of Atmospheric Science,Lanzhou University,Lanzhou 730000,China;
2. Laboratory for Climate and Ocean-Atmosphere Studies,Department of Atmospheric and Oceanic Sciences,School of Physics,Peking University,Beijing 100871,China;
3. Institute of Arid Meteorology,China Meteorological Administration (CMA);Key Laboratory of Arid Climatic Change and Reducing Disaster of Gansu Province;Key Open Laboratory of Arid Climatic Change and Disaster Reduction of CMA,Lanzhou 730020,China;
4. Gansu Meteorological Bureau,Lanzhou 730020,China
Abstract: Using the Met Office large eddy model (LEM), a series of boundary-layer convection and deep convection experiments are performed to examine the characteristics of turbulence and deep convection triggering mechanism under different initial potential temperature and specific humidity conditions over wet and dry surfaces. Model results indicate that the mixed layer is warmer and drier over dry surface with a larger depth, and the opposite is true over wet surface. Because the surface sensible heat flux is more efficient in the formation of thermal turbulence, turbulent mixing and entrainment are stronger over dry surface, making the water vapor mixing ratio and equivalent potential temperature more uniform within the convective boundary layer (CBL) but a greater negative disturbance at the top of the CBL. The convection structure is the same as the bubble-like convection with a horizontal net-like structure over both wet and dry surfaces, but more vigorous over dry surface. The triggering of deep convection over different surfaces varies dramatically with the initial atmospheric condition. With weak atmospheric inversion (0.15 K/(100 m)) in the 1-3 km layer, deep convection first occurs over dry surface, but the accompanied clouds are thinner than those over wet surface because of the smaller convective available potential energy. Increasing the inversion intensity in the 1-3 km to 0.55 K/(100 m) not only delays the formation of clouds, but also decreases the cloud thickness. Only when combined with high initial specific humidity at the top of the CBL, can deep convections occur, also first over dry surface, and the associated clouds are thicker as a result of a significant cooling of the CBL top.
Key words: Wet and dry surfaces    Initial thermal and humidity profiles    Characteristics of turbulence    Deep convection    Large eddy simulation    
1 引 言

土壤湿度与降水的相互作用是气候系统中的一个重要过程(张人禾等,2016张述文等,2016),它不仅与突发、多发的中尺度灾害天气过程相联系,而且还影响到不同空间尺度的水循环,尤其是在干旱半干旱地区的表现更加显著(Taylor,et al,2012张强等,2019)。影响土壤湿度与降水相互作用的关键在于土壤湿度如何通过地表能量分配、边界层过程以及与云的耦合影响深对流触发(Dirmeyer,et al,2014Santanello,et al,2018)。

目前,中外许多学者认为深对流触发主要存在两种可能的作用机制。第一种机制主要考虑中尺度环流的作用。模式研究(Adler,et al,2011Maronga,et al,2013Rieck,et al,2014Kang,2016Lee,et al,2019)表明:在适当的大尺度流场条件下,土壤湿度的非均匀分布可引发二级环流,低层大气的水平辐合使上升运动加强,空气块很容易被抬升到自由对流高度并释放对流有效位能,有利于深对流的触发。这种情况的深对流多形成于相对干旱的地区,并被许多观测研究证实(Taylor,et al,20102012Taylor,2015)。第二种机制侧重于边界层热力学演变。土壤湿度的不同导致地表感热和潜热通量分配的不同,使地表向边界层输送的热量和水汽存在差异,影响边界层温、湿度结构和演变规律,进而影响深对流触发。利用一维边界层模式,发现深对流首先发生于湿润或是干旱下垫面很大程度上取决于模式设置的大气初始状况(Ek,et al,1994Findell,et al,2003Gentine,et al,2013Tawfik,et al,2015Yin,et al,2015Bhowmick,et al,2018)。然而,由于“土壤湿度—边界层过程—深对流发生”物理过程十分复杂(Ford,et al,2015Guillod,et al,2015Tuttle,et al,2016Santanello,et al,2018栾澜等,2018),受大气初始状况影响的深对流触发机制尚未得到系统理论解释。以往的研究多集中于土壤湿度的非均匀分布对深对流的影响(Maronga,et al,2013Kang,2016Cao,et al,2018),侧重于二级环流的动力作用。较少探讨湍流作用下大气边界层如何响应土壤湿度变化、如何传递并导致深对流触发,以及初始大气温湿状况对干、湿地表上深对流形成和发展的影响等。

利用大涡模式模拟研究水平均匀的湿润和干旱下垫面上湍流的结构特征以及深对流的形成和演变过程,通过改变不同大气初始状况下地表鲍恩比的敏感性数值试验,研究不同初始大气状况下湿润和干旱下垫面对对流边界层结构、湍流特征、对流形式和强度以及各气象要素水平分布的影响,并进一步探讨湿润和干旱下垫面上深对流的触发机制和演变过程。鉴于深对流触发一直是数值天气预报的瓶颈问题,该研究结果有助于对灾害对流天气预报能力的提高,以及区域尺度到全球尺度的边界层参数化方案和积云对流参数化方案的改进。

2 大涡模式及数值试验

采用的大涡模式是英国气象局的Large Eddy Model(LEM,Version 2.4)(Gray,et al,2001)。LEM是一种高分辨率、非静力平衡大气数值模式,它利用谱滤波方式对Navier-Stokes方程组求网格体积平均,从而准确求解大尺度湍流涡旋的运动,而对于小尺度湍流通量和湍流能量采用Smagorinsky-Lilly次网格模型进行模拟(黄倩等,2014)。LEM的三相云微物理参数化方案中包含了34种水在不同相态之间转化的微物理过程。LEM可以模拟湍流和边界层干、湿对流过程(包括浅层对流和层积云)(Huang,et al,2010王蓉等,2015任燕等,2018Huang,et al,2018张蕴帅等,2021),也可用于研究中尺度对流系统的降水对流云问题(Sang,et al,2018)。

本研究设计了6组采用不同初始温、湿度廓线的湿润和干旱下垫面试验进行对比,利用大涡模式模拟边界层对流和深对流(表1),试验名称的第1个字母表示对流类型:B为边界层对流,D为深对流;第2个字母表示下垫面类型:W为湿润下垫面,D为干旱下垫面;第3个字母表示初始位温;第4个字母表示初始比湿。所有试验均采用Arakawa-C交错网格,水平方向上为50×50个等距网格,垂直方向为120层上疏下密的张弛网格,结合边界层对流和深对流的时空尺度和模式特点(Agee,et al,1999Findell,et al,2003张强等,20082019Kang,2016Cioni,et al,2017Wang,et al,2019王蓉等,2019张蕴帅等,2021),两类试验的模拟区域、模拟高度、水平分辨率和模拟时长的设置如表2所示;根据实际大气在清晨的温湿分布特征,图1给出了理想初始位温和比湿廓线,在0.9 km以上设置3个不同强度的逆温层(表3)。本研究只考虑热力湍流的作用,初始风场设置为0。根据地表能量平衡理论,有

表 1  不同初始位温廓线(θ1θ2)和比湿廓线(q1q2q3)条件下湿润和干旱下垫面的边界层对流和深对流试验 Table 1  Experiments for boundary-layer convection and deep convection cases initialized by different potential temperature(θ1θ2)and specific humidity(q1q2q3)profiles over wet and dry surfaces
对流类型 下垫面 θ1q1 θ1q2 θ1q3 θ2q1 θ2q2 θ2q3
边界层对流(B)
湿润(W) BWθ1q1 BWθ1q2 BWθ1q3 BWθ2q1 BWθ2q2 BWθ2q3
干旱(D) BDθ1q1 BDθ1q2 BDθ1q3 BDθ2q1 BDθ2q2 BDθ2q3
深对流
(D)
湿润(W) DWθ1q1 DWθ1q2 DWθ1q3 DWθ2q1 DWθ2q2 DWθ2q3
干旱(D) DDθ1q1 DDθ1q2 DDθ1q3 DDθ2q1 DDθ2q2 DDθ2q3
表 2  边界层对流和深对流试验的模拟范围、水平分辨率和模拟时间长度 Table 2  Model domain,horizontal resolution and simulation period for boundary-layer convection and deep convection cases
对流类型 模拟区域 模拟高度 水平分辨率 模拟时长
边界层对流 10 km×10 km 6 km 0.2 km 8 h
深对流 100 km×100 km 10 km 2 km 12 h
图 1  初始位温 (θ1θ2) 廓线和比湿 (q1q2q3) 廓线 Fig. 1  Initial profiles of potential temperature (θ1θ2) and specific humidity (q1q2q3
表 3  不同高度处初始位温(θ1θ2)的逆温强度 Table 3  Initial inversion intensity at different heights
初始位温
高度(km)
0-0.9 0.9-1 1-3 3-5
θ1 0 4 K/(100 m) 0.15 K/(100 m) 0.65 K/(100 m)
θ2 0 2 K/(100 m) 0.55 K/(100 m) 0.35 K/(100 m)
$ {R}_{\rm n}=H+LE+G+\Delta Q $ (1)

式中,Rn为地表净辐射,H为地表感热通量,LE为地表潜热通量,G为土壤热通量,ΔQ为地表能量残余。不考虑地表能量残余并忽略土壤热通量,由于土壤湿度的差异会引起地表感热和潜热通量分配的差异,综合夏季午后干旱和湿润下垫面的观测和模拟结果(吕世华等,1995张强等,2003Zhang,et al,2011周甘霖等,2019),所有试验的地表感热和潜热通量之和均取500 W/m2,干旱下垫面的地表感热和潜热通量分别设定为400 和100 W/m2,即地表鲍恩比为4,而湿润下垫面与此相反,地表鲍恩比为0.25。

模式的侧边界采用周期侧边界条件,上、下边界采用刚性边界条件,并在距模式顶约2/3高度以上加入Newtonian 阻尼层,以吸收模式上边界反射的重力内波。模式的下边界条件由Businger-Dyer方案给出。为了激发湍流运动,在初始时刻对1 km高度以下的位温和比湿场加入±0.2 K和±0.05 g/kg的随机扰动。模式对时间和空间均采用中央差分方案进行数值计算,积分时间步长为0.01 s,模拟1 h 后输出第一次数据,以后每隔 15 min 输出一次数据。

3 模拟结果及分析 3.1 湿润和干旱下垫面边界层结构

从不同大气初始条件下湿润和干旱下垫面上平均虚位温和水汽混合比廓线2 和6 h的模拟结果(图2)可以看出,由于湍流作用将地表热量和水汽向上输送,从2—6 h平均虚位温升高,混合层向上发展。干旱下垫面上混合层厚度较湿润下垫面大,但后者的混合层较前者湿润。初始位温θ1在0.9—1 km是强度较大的逆温层(图1),逆温强度为初始位温θ2的2倍(表3)。对比图2ab发现,2 h后,初始位温为θ2时,湿润下垫面上混合层已经发展到1 km以上,而初始位温为θ1时,湿润下垫面上混合层的垂直发展受0.9—1 km的逆温层中较大负浮力的影响而被限制在1 km以下。6 h后,该逆温层对湿润下垫面上混合层的影响仍然存在,而初始位温为θ1的干旱下垫面上较强的热力湍流使混合层在2 h后已经穿透该逆温层并迅速向上发展,6 h后混合层高度甚至接近3 km(图2c)。然而,同一时间初始位温为θ2的干旱下垫面上混合层高度约为2.4 km(图2d),这是因为初始位温θ1在1 km以上的逆温强度较小,为边界层的发展提供了有利条件,虽然在3 km以上稳定度再次增大,但是至模式运行结束(8 h),边界层并没有到达该高度,所以3 km以上的逆温层对边界层发展的影响很小。另外,所有试验中湍流均能贯通0.9—1 km高度上的逆温层,因此,影响边界层发展高度的关键因素为1—3 km的逆温强度。混合层内水汽混合比随高度升高而略减小,尤其是混合层发展到1 km附近时,由于其上水汽混合比较小,干空气夹卷进入混合层使水汽混合比降低。干旱下垫面上混合层因垂直增长而变干,而湿润下垫面上混合层则因较大的地表潜热通量而进一步变湿。

图 2  湿润和干旱下垫面平均虚位温 (实线) 和水汽混合比 (虚线) 垂直廓线2 h (a、b) 和6 h (c、d) 的模拟结果 (a、c. θ 1,b、d. θ 2 Fig. 2  Vertical profiles of averaged virtual potential temperature (solid lines) and water vapor mixing ratio (dashed lines) over wet and dry surfaces after 2 h (a,b) and 6 h (c,d) of integration (a,c. θ1,b,d. θ2

图3给出了不同大气初始状况下湿润和干旱下垫面上垂直速度、水汽混合比以及相当位温在边界层不同高度上概率密度函数(PDF)的6 h模拟结果。从垂直速度来看(图3adg),其PDF呈正倾斜(峰值在左侧)分布,说明上升气流强于下沉气流但下沉气流多于上升气流。到 $ {z}_{i} $ 高度处,正倾斜分布更加明显,这是因为在边界层低层,地表加热使气块受到浮力作用而上升,但是随高度的升高,浮力作用减弱,只有少量发展很强的上升气流才能到达(甚至进入)自由大气。水汽混合比的PDF在边界层内呈负倾斜分布(图3eh),Couvreux等(2005)认为这主要与由夹卷层向下延伸的干舌有关。干舌的侵入使混合层中空气变干,PDF表现为长尾拖向左侧,在 $\dfrac{\,2\,}{\,3\,} {z}_{i} $ 高度上最显著。而在边界层顶,较强的上升运动将边界层低层湿润的空气带到较干的夹卷层,上升运动和干空气夹卷的共同作用使水汽混合比的PDF趋于平缓(图3b),随着高度的升高,上升运动对水汽场产生的正扰动增大,使得水汽混合比的PDF在1.1 $ {z}_{i} $ 高度上变为正倾斜分布(图略)。在边界层内,相当位温受水汽影响较大,PDF呈负倾斜分布(图3fi);而在边界层顶,受自由大气中干暖空气的夹卷作用影响,相当位温的PDF分布更加平缓(图3c)。

图 3  不同初始位温条件下湿润和干旱下垫面垂直速度 (a、d、g)、水汽混合比 (b、e、h) 和相当位温 (c、f、i) 在 $\dfrac{\,1\,}{\,3\,} {z}_{i} $ (g—i)、 $\dfrac{\,2\,}{\,3\,} {z}_{i} $ (d—f) 和 $ {z}_{i} $ (a—c) 高度处概率密度函数的6 h模拟结果 Fig. 3  Probability density functions of vertical velocity (a,d,g),water vapor mixing ratio (b,e,h) and equivalent potential temperature (c,f,i) at the heights of $\dfrac{\,1\,}{\,3\,} {z}_{i} $ (g—i), $ \dfrac{\,2\,}{\,3\,}{z}_{i} $ (d—f) and $ {z}_{i} $ (a—c) over wet and dry surfaces under different initial potential temperatures after 6 h of integration

图3还可以看出,湿润下垫面上垂直速度的分布范围比干旱下垫面窄,说明上升和下沉气流的强度均弱于干旱下垫面。此外,在边界层内,干旱下垫面上由于强烈的湍流混合作用,水汽混合比和相当位温的分布更加集中;而在边界层顶,PDF的偏度值大于湿润下垫面,表明干旱下垫面上较强的夹卷作用对水汽混合比和相当位温场产生较大的负扰动。进一步考虑初始位温的影响,在相同的下垫面上,边界层顶的逆温强度越弱,越有利于边界层增长,从而为垂直运动提供有利的发展条件,使得边界层内较强上升和下沉速度的概率密度越大(图3dg);而在边界层顶,由于直接受逆温层的阻挡,随着逆温强度的增大,较强上升和下沉运动的概率密度均减小(图3a)。同时,初始位温通过影响边界层的高度进而影响边界层的温、湿状况,即当地表感热(潜热)通量一定时,边界层越低,热量(水汽)在越薄的气层中混合,边界层越暖(湿),反之亦然。初始比湿对气象要素水平分布的影响较小(图略)。

3.2 湿润和干旱下垫面边界层对流特征

边界层内存在多种形式的对流,如泡状对流、卷状对流等,不同形式的边界层对流发展对深对流的触发过程也不同(Weckwerth,2000Huang,et al,2009Salesky,et al,2017Lee,et al,2019)。图4给出了湿润(图4a、c)和干旱(图4b、d)下垫面上边界层对流2和6 h的模拟结果。仅在热力湍流的作用下,湿润和干旱下垫面上对流形式基本一致,均为泡状对流(王敏仲等,2018)。2 h后,湿润下垫面对流边界层高度为0.9 km,上升和下沉气流排列规则。在水平方向上,上升气流呈网状结构,中间为下沉气流(图略)。同一时间,干旱下垫面上升气流和下沉气流均强于湿润下垫面,但分布不规则,下沉气流范围较大。随着对流发展,湿润下垫面上热泡发生倾斜和合并,而干旱下垫面上由于地表持续的加热作用,边界层内形成两个强上升区,垂直尺度达2.8 km,水平尺度约1.3 km,热泡内最大上升速度为9.4 m/s。部分热泡上升进入夹卷层后,受到负浮力的作用,上升运动减弱,并返回混合层,在边界层顶附近形成独立的下沉运动中心,如y=−1.3 km(图4a)、y=3 km(图4d)等。

图 4  模拟的湿润 (a、c,试验BWθ1q1) 和干旱 (b、d,试验BDθ1q1) 下垫面垂直速度沿x = 0 km的垂直剖面( a、b. 2 h, c、d. 6 h;虚线为边界层顶高度) Fig. 4  The y-z cross sections of vertical velocity w along the line of x = 0 km over wet (a,c for case BWθ1q1) and dry (b,d for case BDθ1q1) surfaces after 2 h (a,b) and 6 h (c,d) of integration (The height of the boundary layer is plotted by dashed line)

改变两种下垫面的初始位温和初始比湿,进一步验证了以上结论,并发现初始比湿对湿润和干旱下垫面上边界层对流差异的影响较小,初始位温的影响主要体现在当1—3 km的逆温强度增大时,热泡的垂直发展受到抑制,水平尺度增大,且水平方向上的网状结构随时间消散缓慢,6 h后仍可分辨(图略)。

3.3 湿润和干旱下垫面湍流特征

图5给出了平均水平和垂直脉动速度方差(以下简称速度方差) $ \dfrac{1}{2}({{\sigma }_{u}^{2}}+{{\sigma }_{v}^{2}}) $ $ {{\sigma }_{w}^{2}} $ 和脉动虚位温方差(以下简称虚位温方差) $ {{\sigma }_{{\theta }_{v}}^{2}} $ 廓线的6 h模拟结果,其中速度方差和虚位温方差分别采用对流速度尺度( $ {w}_{*} $ )和温度尺度( $ {\theta }_{{v}_{*}} $ )进行归一化(Stull,1988),即

图 5  模拟的6 h不同初始状况下湿润和干旱下垫面无量纲水平速度方差 (a、d)、无量纲垂直速度方差 (b、e) 和无量纲虚位温方差 (c、f) 廓线 (a、b、c. θ1,d、e、f. θ2 Fig. 5  Vertical profiles of normalized horizontal velocity variance (a,d),vertical velocity variance (b,e) and virtual potential temperature variance (c,f) over wet and dry surfaces after 6 h of integration under different initial potential temperature and specific humidity conditions (a,b,c. θ1,d,e,f. θ2
$ {w}_{*}={\left[\frac{{g}{z}_{i}}{{\theta }_{v}}\overline{{\left({w'}{\theta }_{v}'\right)}_{\rm s}}\right]}^{1/3} $ (2)
$ {\theta }_{{v}_{*}}= -\frac{\overline{{\left({w}'{{\theta }_{v}}'\right)}_{\rm s}}}{{w}_{*}} $ (3)

式中,g为重力加速度,文中取9.8 m/s2 $ {z}_{i} $ 为对流边界层顶的高度, $ {\theta }_{v} $ 为混合层内平均虚位温, $ \overline{{\left({w}'{\theta }_{v}'\right)}_{\rm s}} $ 为地表动力学感热通量。图5ad显示,水平速度方差在地表受摩擦的影响有最大值,随高度的升高而迅速减小,在(0.3—0.7) $ {z}_{i} $ 随高度的变化较小,而在混合层顶附近,上升运动受逆温层的阻挡转变为水平运动,使水平速度方差又略微增大,1—3 km的逆温强度越大(θ2),增幅越明显(图5d),然后水平速度方差再次减小。垂直速度方差从地表向上迅速增大,在 $\dfrac{\,1\,}{\,3\,} {z}_{i} $ 处达到最大值,之后随高度的升高而减小。与水平速度方差相比,干旱下垫面上垂直速度方差较湿润下垫面增大更明显,说明干旱下垫面更有利于热力湍流的发生和发展。初始位温为θ2时,两种下垫面上方差差异减小(图5df),说明逆温强度的增大使湍流强度对下垫面的热力响应降低。

由于对流泡在近地面形成,虚位温方差在近地面出现极大值,在边界层中部因较强的湍流混合,数值较小且几乎不随高度变化,在夹卷层中受冷空气上冲和暖空气夹卷的共同影响出现极大值。对比图5cf发现,初始位温为θ2时,夹卷层中虚位温方差较大,这与1—3 km高度较强的逆温有关,从图2d也可以看出,初始位温为θ2时,6 h后边界层顶附近的虚位温梯度较大,这将在夹卷层中产生较大的虚位温扰动。此外,初始位温为θ1时,受初始时刻0.9—1 km高度强逆温层的影响(图2c),湿润下垫面上虚位温方差在夹卷层的极大值略大于干旱下垫面(图5c);而初始位温为θ2时,6 h后湿润和干旱下垫面上虚位温方差基本一致(图5f)。初始比湿对湍流统计量的影响较小,可以忽略。

通过对比分析湿润和干旱下垫面上边界层对流试验的结果,发现湿润和干旱下垫面上对流边界层的结构和湍流特征有较显著的差异,而湍流动能(TKE)是湍流强度的度量。图6给出了湿润(试验BWθ2q1)和干旱(试验BDθ2q1)下垫面上TKE的源(浮力项)和汇(耗散项)随高度分布的6 h模拟结果,并考虑了大(网格)尺度湍涡和小(次网格)尺度湍涡对湍流动能产生的贡献。浮力项在近地面最大,随高度线性递减,到夹卷层减为负值并消耗湍流动能。尽管图6中两个试验采用了相同的初始温湿廓线且地表感热和潜热通量之和一致,但干旱下垫面上湍流动能的产生却强于湿润下垫面,这种差异是由地表能量分配的差异造成的:干旱下垫面上较大的浮力项来源于较大的地表感热通量,而湿润下垫面上地表潜热通量的增加虽然补偿了地表感热通量的减小,但是地表感热通量对于热力湍流形成的作用更强。在近地面受地表摩擦的影响,湍涡尺度较小;在较大高度上,湍涡的发展不再受地表的限制,能够形成较大尺度的湍涡,因此耗散项的最大值出现在近地层,并随高度减小。大尺度湍涡在近地层之上的混合层对浮力项有主要贡献,而小尺度湍涡在近地层对浮力项的贡献不能忽略,约占浮力项的40%。

图 6  湿润 (蓝线,试验BWθ2q1) 和干旱 (红线,试验BDθ2q1) 下垫面湍流动能收支方程中浮力项 (Buo) 和耗散项 (Diss) 随高度分布的6 h模拟结果 (BuoCN和BuoSG分别代表网格尺度和次网格尺度湍涡对浮力项的贡献) Fig. 6  Vertical profiles of the buoyancy term (Buo) and dissipation term (Diss) in the turbulence kinetic energy budget equation over wet (blue lines,case BWθ2q1) and dry (red lines,case BDθ2q1) surfaces after 6 h of integration (BuoCN and BuoSG represent the contributions of grid-scale and subgrid-scale turbulence eddies to buoyancy term respectively)
3.4 湿润和干旱下垫面深对流触发过程

未饱和的气块在地表受热后沿干绝热线上升,其位温和水汽混合比保持不变,但温度却以干绝热递减率降低,使气块中的水汽逐渐趋于饱和。达到抬升凝结高度后,水汽发生凝结并释放潜热,位温增大,温度递减率降低,气块温度按湿绝热递减率降低。在整个过程中,相当位温始终保持不变。当气块达到自由对流高度后,由于温度高于周围环境大气,受浮力作用将产生向上的加速度而持续上升。当其他条件不变时,增大低层大气的水汽含量,会使抬升凝结高度降低,相当位温增大,进而自由对流高度降低;而升高低层大气的温度,则会使抬升凝结高度升高,但相当位温升高,自由对流高度降低。

图7来看,初始大气在910 hPa高度附近均存在较强的逆温。初始位温θ1在900—705 hPa高度的温度递减率接近于干绝热递减率,而在705—545 hPa高度的温度递减率更接近于湿绝热递减率,在这种层结中,较小的相当位温增量都将引起对流有效位能的显著增加;而初始位温θ2则在900—705 hPa高度的温度递减率接近于湿绝热递减率,在705—545 hPa高度的温度递减率接近于干绝热递减率,由于900—705 hPa高度的层结曲线斜率较小,状态曲线与层结曲线相交于该气层时,较小的相当位温增量都将引起自由对流高度明显下降。从图7还可以发现,初始大气温度随着高度升高而降低,温度露点差也随之减小,空气更容易达到饱和,在745 hPa高度之上,初始比湿q1q2q3的温度露点差依次减小,且初始比湿为q3时,547 hPa(4.7 km)高度以上的空气处于饱和状态(图7ef),更易形成较厚的云层。

图 7  6种初始大气状况 (a. θ1 q1,b. θ2 q1,c. θ1 q2,d. θ2 q2,e.θ1 q3,f. θ2 q3) 的斜温图 (黑实线为层结曲线,蓝实线为露压曲线,红虚线为状态曲线) Fig. 7  Thermodynamic profiles (i.e.,skewT-lgp diagrams) of six atmospheric conditions combined freely by the potential temperature θ1θ2 and specific humidity q1q2q3 at the initial time (a. θ1 q1,b. θ2 q1,c. θ1 q2,d. θ2 q2,e. θ1 q3,f. θ2 q3;Temperature is represented by black solid line and dew point temperature by blue line,the red dashed line represents the state curve profiles)

为了进一步研究不同初始位温和比湿条件下湿润和干旱下垫面上湍流特征对深对流的发生和发展的影响,图8给出了深对流试验模拟的云水凝结物(液态水和冰水)混合比、边界层高度、抬升凝结高度、自由对流高度随时间的变化。初始位温为θ1时(图8af),随着地表热量的输送,相当位温增大,引起对流有效位能持续增加(图略)。湿润下垫面上对流有效位能(最大为360.3 J/kg)明显大于干旱下垫面(最大值为62.3 J/kg),Yin等(2015)在采用美国南部大平原实测资料初始化的模拟研究中也得出了相似的结论,Taylor等(2010)在对Sahel地区中尺度对流系统的个例分析中指出由于较大的地表潜热通量,湿润下垫面上相当位温较大,使对流有效位能大于干旱下垫面。由于1—3 km的逆温强度较弱,干旱下垫面上边界层对流穿透0.9—1 km的初始逆温层后迅速向上发展,并于3 h后与抬升凝结高度相交(图8df)。虽然湿润下垫面上边界层增长缓慢,但是湍流向上输送了较多水汽,使边界层增湿,因此抬升凝结高度较低且基本不随时间变化,边界层顶与抬升凝结高度在7 h后相交(图8ac)。在抬升凝结高度之上,上升气流中水汽发生凝结,云层开始形成。自由对流高度则随着地表热量的输入而不断降低,与抬升凝结高度几乎同时和边界层顶相交。当边界层顶发展到自由对流高度后,对流有效位能转化为动能,垂直速度迅速增大(湿润下垫面上由2 m/s增大到12 m/s,干旱下垫面上由4 m/s增大到7 m/s),云层急剧增厚,上升运动不再受到抑制,因此认为边界层顶与自由对流高度相交时,深对流被触发(Findell,et al,2003Bhowmick,et al,2018)。从图8af可以看出,当初始大气在1—3 km的逆温强度较弱时,深对流先在干旱下垫面上发生,这与Findell等(2003)得出的结论一致。由于对流有效位能较小,干旱下垫面上的云层厚度小于湿润下垫面。随着深对流的发展,边界层顶以上初始比湿较大时湿润下垫面(试验DWθ1q3)在10 h后对流有效位能出现次峰值(196.4 J/kg),Cioni等(2017)研究中也发现当深对流优先在干旱下垫面上发生时,湿润下垫面上对流有效位能形成两个峰值。这是因为在湿润下垫面上较厚云层的形成过程中,水汽的凝结释放大量潜热,使对流有效位能增加并形成极大值,而在初始位温为θ1的其他两组对比试验中,由于745 hPa高度以上水汽含量较低(图7ac),对流有效位能仅形成一个峰值。

图 8  所有深对流试验中液态水和冰水混合比之和 (色阶)、边界层高度 (BLH,黑实线)、抬升凝结高度 (LCL,黑点划线)、自由对流高度 (LFC,黑虚线) 随时间的变化 (a—f. 采用初始位温θ1,g—l. 采用初始位温θ2;a、d、g、j. 采用初始比湿q1,b、e、h、k. 采用初始比湿q2,c、f、i、l. 采用初始比湿q3;边界层高度由水汽梯度极小值法计算得出 (张宏昇等,2020);云层由0.01 g/kg等值线表示) Fig. 8  Profiles of the sum of liquid and ice water mixing ratios (color shaded),boundary layer height (BLH,black solid lines),lifting condensation level (LCL,black dash-dotted lines) and the level of free convection (LFC,black dashed lines) as functions of time (a-f. using the initial temperature θ1 and g-l using the initial temperature θ2;a,d,g,j. using the initial specific humidity q1;b,e,h,k. using the initial specific humidity q2 and c,f,i,l. using the initial specific humidity q3;The height of the boundary layer is computed using the minimum value of water vapor gradient approach;the cloud layer is indicated by the 0.01 g/kg mixing ratio contour)

增大1—3 km的逆温强度(图8gl),云层的形成滞后,云层厚度也明显减小,仅当边界层顶以上初始比湿较大时(试验DWθ2q3和DDθ2q3)有深对流发生,深对流仍然先发生于干旱下垫面(图8il)。而Findell等(2003)则认为当1—3 km的温度递减率接近湿绝热递减率时,干旱下垫面(鲍恩比为1.4—1.6)上边界层的发展受到抑制,而湿润下垫面(鲍恩比为0.3—0.4)上边界层的增湿将有利于自由对流高度降低,并与边界层顶相交,从而触发深对流。在本研究中,虽然湿润下垫面(试验DWθ2q3)自由对流高度在10 h后明显降低(图8i),但是干旱下垫面(试验DDθ2q3)边界层的增长并没有受到抑制,并且由于3 km以上较低的逆温强度(表3),边界层高度甚至略高于试验DDθ2q3。造成这种差异的原因可能是本研究中干旱下垫面上较大的鲍恩比(值为4)更有利于边界层的发展,削弱了覆盖逆温增强对边界层增长的抑制作用,导致深对流先在干旱下垫面上发生。虽然湿润下垫面上对流有效位能和边界层以上的水汽含量均高于干旱下垫面,但是后者云水凝结物的含量和云层厚度均大于前者。云水凝结物主要有两个来源:(1)由湍流将地表蒸发的水汽输送到抬升凝结高度发生凝结(用水汽垂直输送最大量表征其强度),(2)自由大气因增湿降温达到饱和而发生凝结(用温度露点差表征其强度)。由于深对流发生前湿润下垫面上最大垂直速度较小(小于3 m/s),湿润下垫面抬升凝结高度以上的水汽垂直输送最大量明显小于干旱下垫面(图9a),由湍流输送引起的水汽凝结受到抑制;另外,虽然湿润下垫面上边界层顶附近的平均水汽垂直输送量较大(图9b),但是干旱下垫面上较强的上升气流(贯穿对流)一方面增加了自由大气中的暖空气向下夹卷混合(图9c),另一方面较冷的混合层空气穿透边界层顶进入自由大气,使边界层顶以上的空气显著冷却(Ek,et al,1994Westra,et al,2012Gentine,et al,2013),温度露点差减小且明显低于湿润下垫面,在1.2 $ {z}_{i} $ 以上甚至接近0(图9d),水汽向液态水转化的效率增大,因此干旱下垫面上云层厚度大于湿润下垫面。但当深对流发生后,随着最大垂直速度的增大,湿润下垫面上云层开始增厚(图8i)。

图 9  湿润 (蓝线,试验DWθ2q3) 和干旱 (红线,试验DDθ2q3) 下垫面 (a) 水汽垂直输送最大值、(b) 平均水汽垂直输送量、(c) 平均热量垂直输送量和 (d) 温度露点差随高度分布的10 h模拟结果 (虚线为抬升凝结高度,点划线为自由对流高度) Fig. 9  Profiles of the maximum (a) and the mean (b) vertical transport of water vapor,mean vertical transport of heat (c) and mean dew-point depression (d) after 10 h of integration over wet (blue lines,case DWθ2q3) and dry (red lines,case DDθ2q3) surfaces (The lifting condensation level is plotted by dashed line and the level of free convection by is denoted by dash-dotted line)
4 结论和讨论

利用大涡模式模拟了对流边界层发展以及深对流的触发过程。通过改变不同大气初始条件下地表鲍恩比的敏感性试验,研究湿润和干旱下垫面边界层结构、湍流特征及其对深对流触发的影响。结果表明:

(1)干旱下垫面的混合层干而暖,厚度较大,湿润下垫面则相反。边界层内,干旱下垫面的上升和下沉运动均强于湿润下垫面,且干旱下垫面上强烈的湍流混合作用使水汽混合比和相当位温更加均一;而在边界层顶,干旱下垫面上较强的夹卷作用对水汽混合比和相当位温造成了较大的负扰动。1—3 km高度的逆温强度是影响边界层高度的关键因素,且初始位温主要通过影响边界层的高度间接影响边界层内水汽混合比、相当位温等变量的分布。

(2)仅考虑热力湍流作用,湿润和干旱下垫面的对流运动均表现为泡状对流,水平方向上呈网状结构。由于地表感热通量比地表潜热通量对热力湍流形成的作用更强,干旱下垫面对流强度和对流尺度均比湿润下垫面大。1—3 km的逆温强度增大,下垫面对湍流强度的影响减弱,湿润和干旱下垫面热泡的垂直发展受到抑制,水平尺度增大,网状结构维持时间增加;初始比湿对湿润和干旱下垫面边界层对流差异的影响较小。

(3)不同下垫面深对流的触发和发展与大气初始温湿状况有关。当初始大气在1—3 km的逆温强度较弱时(0.15 K/(100 m)),边界层内湍流迅速发展,干旱下垫面首先发生深对流,但因对流有效位能较小,云层厚度小于湿润下垫面。1—3 km的逆温强度增加到0.55 K/(100 m),对流云形成时间滞后,云层厚度明显减小,仅在边界层顶水汽饱和度较高的初始比湿条件下有深对流发生,且深对流仍然首先发生在干旱下垫面,考虑热泡的贯穿对流过程在边界层顶引起的较强冷却,云层厚度大于湿润下垫面。

以往关于深对流触发机制的研究或侧重于土壤湿度非均匀分布诱发的次级环流上升支的抬升作用,或仅采用一维边界层模式获取边界层的热力演变,而没有给出湍流过程对它的影响。文中针对热力湍流在不同大气温湿状况下对水汽、热量的垂直输送在“土壤湿度—边界层热力演变—深对流”物理过程中的作用等方面开展研究。研究结果有助于深入理解深对流的触发过程,改进积云对流参数化方案,进而提高对灾害对流天气的预报能力。由于本研究使用的初始温、湿廓度线均是理想的廓线,结果具有一定的局限性。下一步的工作将利用实测的探空资料对模式初始化,进行个例模拟,对该研究结果进一步充实。

参考文献
黄倩, 王蓉, 田文寿等. 2014. 风切变对边界层对流影响的大涡模拟研究. 气象学报, 72(1): 100-115. Huang Q, Wang R, Tian W S, et al. 2014. Study of the impacts of wind shear on boundary layer convection based on the large eddy simulation. Acta Meteor Sinica, 72(1): 100-115. (in Chinese)
栾澜, 孟宪红, 吕世华等. 2018. 青藏高原土壤湿度触发午后对流降水模拟试验研究. 高原气象, 37(4): 873-885. Luan L, Meng X H, Lü S H, et al. 2018. Simulation on afternoon convective precipitation triggered by soil moisture over the Qinghai-Tibetan Plateau. Plateau Meteor, 37(4): 873-885. (in Chinese)
吕世华, 陈玉春. 1995. 绿洲和沙漠下垫面状态对大气边界层特征影响的数值模拟. 中国沙漠, 15(2): 116-123. Lü S H, Chen Y C. 1995. The numerical simulation of the features of the planetary boundary layer of the oasis and the Gobi desert in the arid region. J Desert Res, 15(2): 116-123. (in Chinese)
任燕, 黄倩, 张君霞等. 2018. 大涡模式分辨率对海洋信风区大气边界层结构和演变模拟的影响. 热带气象学报, 34(1): 23-33. Ren Y, Huang Q, Zhang J X, et al. 2018. Effect of large eddy resolution on simulation of structure and evolution of the marine trade-wind boundary layer. J Trop Meteor, 34(1): 23-33. (in Chinese)
王敏仲, 徐洪雄, 王寅钧等. 2018. 塔克拉玛干沙漠夏季晴空对流边界层大涡模拟. 中国沙漠, 38(6): 1275-1286. Wang M Z, Xu H X, Wang Y J, et al. 2018. Large-Eddy simulation of summer clear sky convective boundary layer in the Taklimakan desert. J Desert Res, 38(6): 1275-1286. (in Chinese)
王蓉, 黄倩, 田文寿等. 2015. 边界层对流对示踪物抬升和传输影响的大涡模拟研究. 大气科学, 39(4): 731-746. Wang R, Huang Q, Tian W S, et al. 2015. Study of a large eddy simulation of the effects of boundary layer convection tracer uplift and transmission. Chinese J Atmos Sci, 39(4): 731-746. DOI:10.3878/j.issn.1006-9895.1502.14155 (in Chinese)
王蓉, 黄倩, 岳平. 2019. 大涡模式水平分辨率对边界层夹卷过程及示踪物垂直传输的影响. 干旱气象, 37(1): 48-56. Wang R, Huang Q, Yue P. 2019. Effect of large eddy horizontal resolution on simulation of entrainment and tracer vertical transport. J Arid Meteor, 37(1): 48-56. (in Chinese)
张宏昇, 张小曳, 李倩惠等. 2020. 大气边界层高度确定及应用研究进展. 气象学报, 78(3): 522-536. Zhang H S, Zhang X Y, Li Q H, et al. 2020. Research progress on estimation of atmospheric boundary layer height. Acta Meteor Sinica, 78(3): 522-536. (in Chinese)
张强, 曹晓彦. 2003. 敦煌地区荒漠戈壁地表热量和辐射平衡特征的研究. 大气科学, 27(2): 245-254. Zhang Q, Cao X Y. 2003. The influence of synoptic conditions on the averaged surface heat and radiation budget energy over desert or Gobi. Chinese J Atmos Sci, 27(2): 245-254. DOI:10.3878/j.issn.1006-9895.2003.02.11 (in Chinese)
张强, 王胜. 2008. 西北干旱区夏季大气边界层结构及其陆面过程特征. 气象学报, 66(4): 599-608. Zhang Q, Wang S. 2008. A study on atmospheric boundary layer structure on a clear day in the arid region in northwest China. Acta Meteor Sinica, 66(4): 599-608. DOI:10.3321/j.issn:0577-6619.2008.04.013 (in Chinese)
张强, 岳平, 张良等. 2019. 夏季风过渡区的陆-气相互作用: 述评与展望. 气象学报, 77(4): 758-773. Zhang Q, Yue P, Zhang L, et al. 2019. land-atmosphere interaction over the summer monsoon transition zone in China: A review and prospects. Acta Meteor Sinica, 77(4): 758-773. (in Chinese)
张人禾, 刘栗, 左志燕. 2016. 中国土壤湿度的变异及其对中国气候的影响. 自然杂志, 38(5): 313-319. Zhang R H, Liu L, Zuo Z Y. 2016. Variations of soil moisture over China and their influences on Chinese climate. Chinese J Nat, 38(5): 313-319. (in Chinese)
张述文, 刘源, 曹帮军等. 2016. GLDAS和CMIP5产品的中国土壤湿度—降水耦合分析及变化趋势. 气候与环境研究, 21(2): 188-196. Zhang S W, Liu Y, Cao B J, et al. 2016. Soil moisture–precipitation coupling and trends in China, based on GLDAS and CMIP5 products. Climatic Environ Res, 21(2): 188-196. (in Chinese)
张蕴帅, 黄倩, 马耀明等. 2021. 黄河源区鄂陵湖湖面和湖边草地对流边界层湍流结构特征的大涡模拟研究. 大气科学, 45(2): 435-455. Zhang Y S, Huang Q, Ma Y M, et al. 2021. Large eddy simulation study of the turbulent structure characteristics of the convective boundary layer over Ngoring Lake and surrounding grassland in the Source Region of the Yellow River. Chinese J Atmos Sci, 45(2): 435-455. (in Chinese)
周甘霖, 李耀辉, 孙旭映等. 2019. 我国北方不同下垫面地表能量通量的变化特征. 干旱气象, 37(4): 577-585. Zhou G L, Li Y H, Sun X Y, et al. 2019. Characteristics of surface energy fluxes over different types of underlying surfaces in North China. J Arid Meteor, 37(4): 577-585. (in Chinese)
Adler B, Kalthoff N, Gantner L. 2011. Initiation of deep convection caused by land-surface inhomogeneities in West Africa: A modelled case study. Meteor Atmos Phys, 112(1-2): 15-27. DOI:10.1007/s00703-011-0131-2
Agee E, Gluhovsky A. 1999. LES model sensitivities to domains, grids, and large-eddy timescales. J Atmos Sci, 56(4): 599-604. DOI:10.1175/1520-0469(1999)056<0599:LMSTDG>2.0.CO;2
Bhowmick M, Parker D J. 2018. Analytical solution to a thermodynamic model for the sensitivity of afternoon deep convective initiation to the surface Bowen ratio. Quart J Roy Meteor Soc, 144(716): 2216-2229. DOI:10.1002/qj.3340
Cao B J, Zhang S W, Li D Q, et al. 2018. Effect of mesoscale land use change on characteristics of convective boundary layer: Semi-idealized large eddy simulations over northwest China. J Meteor Res, 32(3): 421-432. DOI:10.1007/s13351-018-7185-8
Cioni G, Hohenegger C. 2017. Effect of soil moisture on diurnal convection and precipitation in Large-Eddy simulations. J Hydrometeorol, 18(7): 1885-1903. DOI:10.1175/JHM-D-16-0241.1
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. DOI:10.1256/qj.04.167
Dirmeyer P A, Wang Z Y, Mbuh M J, et al. 2014. Intensified land surface control on boundary layer growth in a changing climate. Geophys Res Lett, 41(4): 1290-1294. DOI:10.1002/2013GL058826
Ek M, Mahrt L. 1994. Daytime evolution of relative humidity at the boundary layer top. Mon Wea Rev, 122(12): 2709-2721. DOI:10.1175/1520-0493(1994)122<2709:DEORHA>2.0.CO;2
Findell K L, Eltahir E A B. 2003. Atmospheric controls on soil moisture–boundary layer interactions. Part Ⅰ: Framework development. J Hydrometeorol, 4(3): 552-569. DOI:10.1175/1525-7541(2003)004<0552:ACOSML>2.0.CO;2
Ford T W, Rapp A D, Quiring S M. 2015. Does afternoon precipitation occur preferentially over dry or wet soils in Oklahoma?. J Hydrometeorol, 16(2): 874-888. DOI:10.1175/JHM-D-14-0005.1
Gentine P, Holtslag A A M, D'Andrea F, et al. 2013. Surface and atmospheric controls on the onset of moist convection over land. J Hydrometeorol, 14(5): 1443-1462. DOI:10.1175/JHM-D-12-0137.1
Gray M E B, Petch J, Derbyshire S H, et al. 2001. Version 2.3 of the Met Office Large Eddy Model. Exeter, UK: The Met Office
Guillod B P, Orlowsky B, Miralles D G, et al. 2015. Reconciling spatial and temporal soil moisture effects on afternoon rainfall. Nat Commun, 6: 6443. DOI:10.1038/ncomms7443
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. DOI:10.1175/2008MWR2450.1
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 Saharan boundary layer. J Geophys Res, 115: D05201.
Huang Q, Marsham J H, Tian W S, et al. 2018. Large-eddy simulation of dust-uplift by a haboob density current. Atmos Environ, 179: 31-39. DOI:10.1016/j.atmosenv.2018.01.048
Kang S L. 2016. Regional Bowen ratio controls on afternoon moist convection: A large eddy simulation study. J Geophys Res Atmos, 121(23): 14056-14083. DOI:10.1002/2016JD025567
Lee J M, Zhang Y Y, Klein S A. 2019. The effect of land surface heterogeneity and background wind on shallow cumulus clouds and the transition to deeper convection. J Atmos Sci, 76(2): 401-419. DOI:10.1175/JAS-D-18-0196.1
Maronga B, Raasch S. 2013. Large-eddy simulations of surface heterogeneity effects on the convective boundary layer during the LITFASS-2003 experiment. Bound Layer Meteor, 146(1): 17-44. DOI:10.1007/s10546-012-9748-z
Rieck M, Hohenegger C, van Heerwaarden C C. 2014. The influence of land surface heterogeneities on cloud size development. Mon Wea Rev, 142(10): 3830-3846. DOI:10.1175/MWR-D-13-00354.1
Salesky S T, Chamecki M, Bou-Zeid E. 2017. On the nature of the transition between roll and cellular organization in the convective boundary layer. Bound Layer Meteor, 163(1): 41-68. DOI:10.1007/s10546-016-0220-3
Sang W J, Huang Q, Tian W S, et al. 2018. A large eddy model study on the effect of overshooting convection on lower stratospheric water vapor. J Geophys Res Atmos, 123(18): 10023-10036.
Santanello J A Jr, Dirmeyer P A, Ferguson C R, et al. 2018. Land-atmosphere interactions: The LoCo perspective. Bull Amer Meteor Soc, 99(6): 1253-1272. DOI:10.1175/BAMS-D-17-0001.1
Stull R B. 1988. An Introduction to Boundary Layer Meteorology. Dordrecht: Kluwer Academic Publishers, 118-120, 164
Tawfik A B, Dirmeyer P A, Santanello J A Jr. 2015. The heated condensation framework. Part I: Description and southern Great Plains case study. J Hydrometeorol, 16(5): 1929-1945. DOI:10.1175/JHM-D-14-0117.1
Taylor C M, Harris P P, Parker D J. 2010. Impact of soil moisture on the development of a Sahelian mesoscale convective system: A case-study from the AMMA Special Observing Period. Quart J Roy Meteor Soc, 136(S1): 456-470. DOI:10.1002/qj.465
Taylor C M, de Jeu R A M, Guichard F, et al. 2012. Afternoon rain more likely over drier soils. Nature, 489(7416): 423-426. DOI:10.1038/nature11377
Taylor C M. 2015. Detecting soil moisture impacts on convective initiation in Europe. Geophys Res Lett, 42(11): 4631-4638. DOI:10.1002/2015GL064030
Tuttle S, Salvucci G. 2016. Empirical evidence of contrasting soil moisture-precipitation feedbacks across the United States. Science, 352(6287): 825-828. DOI:10.1126/science.aaa7185
Wang M Z, Xu X D, Xu H X, et al. 2019. Features of the deep atmospheric boundary layer over the Taklimakan Desert in the summertime and its influence on regional circulation. J Geophys Res Atmos, 124(23): 12755-12772. DOI:10.1029/2019JD030714
Weckwerth T M. 2000. The effect of small-scale moisture variability on thunderstorm initiation. Mon Wea Rev, 128(12): 4017-4030. DOI:10.1175/1520-0493(2000)129<4017:TEOSSM>2.0.CO;2
Westra D, Steeneveld G J, Holtslag A A M. 2012. Some observational evidence for dry soils supporting enhanced relative humidity at the convective boundary layer top. J Hydrometeorol, 13(4): 1347-1358. DOI:10.1175/JHM-D-11-0136.1
Yin J, Albertson J D, Rigby J R, et al. 2015. Land and atmospheric controls on initiation and intensity of moist convection: CAPE dynamics and LCL crossings. Water Resour Res, 51(10): 8476-8493. DOI:10.1002/2015WR017286
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): 1586-1594. DOI:10.1007/s11430-011-4207-0