大气边界层研究进展

车军辉， 赵平， 史茜， 杨秋彦

1 中国气象科学研究院灾害天气国家重点实验室, 北京 100081 2 南京信息工程大学气象灾害预报预警与评估协同创新中心， 南京 210044 3 山东省气象服务中心， 济南 250031 4 山东省气象台， 济南 250031

0 引言

大气边界层通常指大气底部直接受地球表面影响的一层，其高度约为1～1.5 km，响应时间尺度小于1小时(Stull，1988)，是地球各个圈层相互作用的关键区域.受地面热力与动力影响，该层大气运动具有明显的湍流性质，并且湍流过程对热量、动量和水汽的垂直输送导致气象要素呈现显著的日变化，因此多尺度大气边界层过程在中尺度气象模式、大气环流模式、天气预报模式、气候预测模式以及大气环境质量预报模式中都具有十分重要的作用(刘树华等，2013).大气边界层高度也是判断湍流混合、垂直扰动、对流传输、云带、大气污染物扩散以及分析大气环境容量的主要特征量(Therry and Lacarrère，1983；Holtslag and Boville，1993；Hong and Pan，1996；Beyrich，1997；Collier et al.，2005；刘树华等，2013)，对云、对流以及大气污染的发展、演变有重要影响，因此数值预报模式对边界层高度的精细参数化不够已经成为制约天气、气候数值预报模式和环境质量预报模式发展的瓶颈之一(Yanai and Li，1994；Garratt，1994；Baklanov et al.，2011；刘树华等，2013；Miao et al.，2015).同时，由于大气边界层是人类活动空间区域所在地，发生在大气边界层内的相关环境问题直接关系人类的健康与生存，因此大气边界层的研究长期受到国际科学界的高度重视.

较早的边界层研究主要集中在近地层(图1)，其中湍流理论和大气探测技术的发展推动了大气边界层研究(Monin and Obukhov，1954).从20世纪70年代开始，飞机、雷达、系留气球、卫星遥感等大气探测技术的发展以及数值模拟技术的出现，促使大气边界层研究逐渐从近地层扩展至整个边界层.大气边界层一般可分为对流、稳定和中性三类，它们的物理性质有着明显的差异.Clarke等(1971)和Kaimal等(1976)分别给出了稳定和对流边界层气象要素的垂直分布结构；Deardorff(1972)首次开展了对中性和对流边界层的大涡数值模拟.这些研究极大地加深了对大气边界层特征的认识，特别是关于大气边界层的区域性特征.青藏高原作为全球最大的复杂地形，其边界层过程的独特性吸引着国内外学者的广泛关注(如，叶笃正等，1958；Tao and Ding，1981；Xu et al.，2002；左洪超等，2004；Li and Gao, 2007；Sun et al.，2007；吕雅琼等，2008；陈学龙等，2010；李茂善等，2011；赵平等，2018；Zhao et al.，2019).

当前非均匀或复杂下垫面、特殊地区和气象条件下的大气边界层特征已经成为国际边界层气象学研究的重点，而外场观测和数值模拟已经成为研究大气边界层特征和形成发展机制的重要方法(张强和胡隐樵，2001；胡非等，2003).刘辉志等(2013,2018)从理论研究、观测试验和数值模拟等方面总结了大气边界层物理和大气化学国家重点实验室在2009—2017年期间关于“边界层湍流结构和物质、能量交换规律”的研究成果.近年来，通过采用更精细的探空数据以及新的再分析数据，边界层气象学在研究大尺度空间范围大气边界层高度的气候学统计特征方面取得新进展，特别是在超过4000 m的超高大气边界层高度的特征及其发展机理方面突破了传统认识，为大气边界层参数化发展奠定了基础.下面我们围绕“大气边界层高度及其发展机理”，就大气边界层高度的分析方法、时空特征，极端气候区的特殊边界层结构和发展机理，边界层参数化方案等方面的主要研究进展，以及有待深入研究的问题作简要评述.

图1 大气边界层研究进展的示意图Fig.1 Schematic diagram of the research progress in atmospheric boundary layer (ABL)

1 大气边界层过程机理

由大气边界层气象学(Stull，1988)得知，无论在陆地或海洋之上，大气边界层高度的共同特征是高气压区的边界层高度比低气压区的低，海洋上空的边界层高度时空变化比陆地的小.大气边界层的结构及其演变具有显著的日变化特征.日出后(图2)，太阳辐射对地面加热，近地层产生的热通量增加使湍流混合加强，大气边界层高度升高，边界层内的热力和动力混合均匀，风速、温度和比湿随高度变化较小，此时称为混合层(即对流边界层).日落后，地面长波辐射冷却，湍流输送减弱，近地层形成稳定边界层，其上层混合层从地面向上抬升，大气湍流特性显著减弱，但仍保持白天混合层气象要素分布特征，与上层自由大气明显不同，被称为剩余层.由于大气边界层顶常有逆温层存在，它抑制湍流混合向上发展，从而形成大气边界层与自由大气的界线.

图2 陆地上高压区的大气边界层结构(对“http:∥www.zhihu.com/question/21763748/answer/31231125”图进行了修改)Fig.2 The structure of the atmospheric boundary layer in a high pressure region over land (a modification of http:∥www.zhihu.com/question/21763748/answer/31231125)

地球表面通过混合层与自由大气层进行热量、动量和水汽交换，混合层过程是大气边界层的主要过程(图2).在此过程中，首先，地面层和大气边界层之间的交换在混合层中达到充分混合；其次，充分混合的边界层空气穿越逆温层，进入自由大气层，实现大气边界层内大气要素与自由大气之间的交换.此外，由于混合层顶常有强逆温层和下沉夹卷运动，进入逆温层的湍涡往往会下沉返回混合层，实现逆温层和混合层的交换.

2 大气边界层高度的计算方法

影响大气边界层高度的主要因素有：地表太阳辐射热力扰动过程、地形动力扰动过程、大气边界层顶夹卷热力动力过程以及水平和垂直扩散过程等.因此，大气边界层高度的时空变化是非定常的，无法通过常规气象观测直接获取.大气边界层高度总体上可以通过大气要素廓线观测数据和参数化或模式计算得到(Seibert et al.，2000).依据观测方式的不同，大气要素廓线观测数据方法又分为直接观测(包括无线电探空、系留气球、边界层铁塔和飞机观测等)和遥感观测(包括声雷达、激光雷达、多普勒天气雷达、边界层风廓线雷达、微波辐射计等，可通过地基、机载和星载等多种方式)两种；各种观测平台的数据具有各自的优势与不足(表1).大气边界层高度参数化计算方法主要有：Deardorff (1972)方法，Pielke和Mahrer(1975)方法等；大气边界层高度模式计算方法主要有：MYJ(Mellor Yamada Janjic)和YSU(Yonsei University)参数化方案.

表1 计算大气边界层高度的观测平台及性能(Seibert et al.，2000)Table 1 Measuring platforms and their performances for the ABL height determination (Seibert et al.，2000)

2.1 无线电探空

无线电探空是计算边界层高度最为常用的数据源，其原理为识别探空观测要素或相关物理量自大气边界层转变为自由大气层所具有的特征高度(Seibert et al.，2000；Seidel et al.，2010)，通常由诊断探空观测的温度、湿度和风速廓线得到(表2)，其中基于位温或虚位温廓线判断边界层高度最为常用.如，用位温梯度明显不连续的高度(即位温梯度法)(Liu and Liang，2010；Zhang et al.，2017)、白天对流发展时的混合层顶高度(即气块法)(Holzworth，1964)以及高抬的逆温层顶高度(Seidel et al.，2010)均可定义对流边界层高度；Bradley等(1993)以近地面逆温层顶作为稳定边界层高度；Liu和Liang(2010)通过近地层位温层结状态将边界层细分为稳定、中性和不稳定三类.此外，Richardson数(Ri)算法也常常被用来估算边界层高度(Troen and Mahrt，1986；Vogelezang and Holtslag，1996)，这种方法是将Ri达到指定阈值所对应的最低高度定义为边界层顶，适用于稳定和对流边界层，该算法对探空观测数据的垂直分辨率要求不高，适合于大量数据统计

表2 基于廓线数据的大气边界层高度计算方法(Seidel et al.，2010)Table 2 Calculation methods of the ABL height with vertical profile data (Seidel et al.，2010)

(Seidel et al.，2012；Guo et al.，2016).采用无线电探空数据计算边界层高度方法的优点是垂直分辨率较高，可以识别边界层的精细结构和类型，其不足是数据的水平和时间分辨率较低.

2.2 激光雷达

激光雷达通过观测气溶胶或痕量气体也能够诊断大气边界层高度，其原理是基于大气边界层内气溶胶粒子的输送特征，气溶胶粒子浓度的垂直突变可以识别边界层高度(Torres et al.，2013；Huang et al.，2014).与其他遥感观测相比，激光雷达数据反映了暖气泡自地面产生、上升至破灭的过程，可以识别大气边界层的对流特征，能够更好地测量大气边界层高度(Coulter，1979；Seibert et al.，2000)，是探测和研究边界层高度的一种重要手段.基于激光雷达，邱金恒等(2003)开展了对流层高度的观测试验，王珍珠等(2008)对北京城区大气边界层结构也进行了研究.

2.3 卫星遥感

通过识别卫星信号在大气边界层顶的突变可以确定大气边界层高度.例如，利用CALIPSO卫星传感器接受的激光雷达后向散射信息，诊断气溶胶浓度的垂直分布，从而识别大气边界层的垂直结构(Jordan et al.，2010；McGrath-Spangler and Denning，2012；Korhonen et al.，2014；Zhang et al.，2016).此外，利用卫星反演的大气垂直廓线数据也可以计算边界层高度(Fetzer and Eric，2004；Wood and Bretherton，2004；Martins et al.，2010；Karlsson et al.，2010；Guo et al.，2011；周文等，2018)，例如近年来随着GPS掩星对低层大气探测准确度的明显提高，Ao等(2012)利用梯度法和GPS掩星资料分析了全球大气边界层的时空分布.与无线电探空相比，卫星遥感能够获取同时刻大尺度空间范围的大气边界层信息，尤其对一些无线电探空观测缺少的复杂地形地区(Jordan et al.，2010；Zhang et al.，2016).但是，由于卫星遥感资料在边界层内垂直分辨率较低，因而卫星遥感资料无法应用于研究边界层的细微结构.

2.4 模式参数化

应用模式参数化方法可计算大气边界层高度.典型的有MYJ(Mellor Yamada Janjic)和YSU(Yonsei University)参数化方案.在MYJ方案中，大气边界层高度定义为湍流动能小于临界值0.1 m2·s-2的最低高度；在YSU方案中，定义整体Ri数小于临界值的高度，在张碧辉等(2012)和Miao等(2017)的研究中，整体Ri的临界值取为0.25.在不同数值模式中计算大气边界层高度的方法差异较大(Seibert et al.，2000).Delage(1974)、Wyngaard(1975)及Estournel和Guedalia(1990)采用湍流能量法计算稳定大气边界层高度，Tennekes(1970)基于相似理论计算对流边界层高度.这些方法的优点是可以分析更大范围的边界层物理特征，但是其缺点是不能排除模式带来的误差.

3 大气边界层高度的时空分布特征

3.1 大气边界层日变化特征

各类大气边界层的发生频率具有明显的日变化特征.Seidel等(2010)指出，就全球而言，近地逆温层(稳定边界层)在夜间出现频率高，热带地区为50%，高纬度极区为80%左右；与夜间相比，稳定边界层白天发生的频率显著偏低，在中、低纬地区为15%，极地为50%左右.美国的稳定边界层夜间出现频率超过70%，而白天(09∶00—18∶00 LST)低于5%；与之相反，对流和中性边界层白天峰值时刻(15∶00 LST)的发生频率高达90%，夜间则低于30%(Liu and Liang，2010).对中国而言，14∶00 BJT对流、中性和稳定边界层出现频率分别为70%、26%和4%(Zhang et al.，2017).很显然，东亚的对流边界层在白天占据主导地位，这与美国的相似.

大气边界层高度的频率分布也表现出明显的日变化特征.图3a给出了中国夏季一个日循环内不同时刻大气边界层高度的频率分布形态(Guo et al.，2016)，可以看到：中午时刻(14∶00 BJT)边界层充分发展，边界层高度最高；日落后(20∶00 BJT和02∶00 BJT)边界层由白天的对流型向夜间稳定型过渡，并伴随着边界层高度降低；日出后(08∶00 BJT)，边界层则由夜间稳定型向白天对流型转换，边界层高度开始升高.此外，不同类型的大气边界层高度频率分布形态也不相同，以美国为例(Liu and Liang，2010)，稳定、中性和对流边界层高度出现频率服从三种不同形态的Gamma分布(图3b)：稳定边界层表现为狭窄的Gamma分布，边界层高度分布区间窄，峰值频率(300 m)为16%；对流边界层呈现宽广的Gamma分布，高度分布区间宽，峰值频率(900 m)仅为4%；中性边界层则介于两者之间.

3.2 大气边界层高度空间分布特征

由于地形和下垫面热力性质的差异，大气边界层高度呈现明显的空间变化特征.在印度西南季风爆发前，新德里地区对流边界层高度可以达到4700 m，而在西南季风发生地的班加罗尔地区，对流边界层高度仅900 m左右(Raman et al.，1990)，这说明在不同天气条件下大气边界层高度有明显差异.Zhang等(2017)的研究指出，虽然中国塔克拉玛干沙漠和珠江三角洲具有相似的边界层发生频率，但是塔克拉玛干沙漠的平均大气边界层高度明显偏高.青藏高原地区边界层高度通常在2000～3000 m(叶笃正和高由禧，1979；Xu et al.，2002；周明煜等，2002；Zhang et al.，2003)，高于平原地区的1000～1500 m(赵鸣和苗曼倩，1992).宋正山等(1984)用第一次青藏高原科学试验资料得到高原西部改则边界层高度可以达到3000 m 以上；李家伦等(2000)、刘红燕和苗曼倩(2001)和吕雅琼等(2008)计算的边界层高度则较低，当雄地区为1400～1800 m，纳木错湖区的对流边界层高度最高在1750 m左右.青藏高原地区边界层高度在干季和雨季有明显差异(左洪超等，2004)，那曲边界层高度在干季为2211～4430 m，在雨季为1006～2212 m(李茂善等，2011)，冬季改则地区有超过5000 m的超高边界层(Chen et al.，2013，2016).因此，在不同地点及天气背景下，边界层高度差异较大，反映了青藏高原复杂地形和下垫面状况对边界层结构的影响.

图3 (a) 中国2011年1月—2015年7月夏季02时(BJT；黑线)、 08时(红线)、14时(蓝线)和20时(绿线)大气边界层高度的发生频率以及每个时刻的探空样本数(N)和边界层高度平均值(Guo et al.，2016)； (b) 美国稳定(stable)、中性(neutral)和不稳定(unstable)三类边界层高度的发生频率，括号内数字为各类边界层观测的样本数，平滑曲线为拟合的各类边界层高度频率Gamma分布曲线，(k，s)表示Gamma分布参数值(Liu and Liang，2010)Fig.3 (a) Frequency distribution of the ABL height at 02∶00 BJT (black), 08∶00 BJT (red), 14∶00 BJT (blue), and 20:00 BJT (green) during the summertime from January 2011 to July 2015 in China, in which the number of soundings (N) and the mean value at each observed time are given (Guo et al.，2016); and (b) frequency distribution of the stable boundary layer (stable), neutral boundary layer (neutral), and unstable boundary layer (unstable) heights in United States, in which the sample number is listed in the legend in parentheses, and a smooth curve is drawn for the fitting Gamma distribution with the specified values of parameters (k,s) (Liu and Liang，2010)

Seidel 等(2012)指出日落时刻美国区域大气边界层高度存在较大空间尺度的东西向高度梯度(图4a)，这与太阳高度角的纬向差异有关，在昼夜过渡时段，东部地区较西部更早进入夜间而转为稳定边界层.在中国区域，大气边界层高度存在三种较大尺度的空间分布型(图4b，c，d)，即在日出时刻的西低-东高分布，日落时刻的东低-西高分布，以及中午时刻的南低-北高分布，其中：前两种特征与美国地区的相似(这与太阳高度角纬向差异有关)，而第三种则可能与局地的地表和水文过程有关(Guo et al.，2016；Zhang et al., 2017).由于土壤水分影响着地面感热通量强度，因而边界层的发展对地面湿度很敏感(McCumber and Pielke，1981；Sanchez-Mejia and Papuga，2014；Rihani et al.，2015).中国华南地区土壤湿度高，近地面感热较低，不利于大气边界层发展，而中国西北地区土壤水分偏低，感热偏大，有利于大气边界层发展(Liu et al.，2004；Zhang et al.，2013；Dirmeyer et al.，2014).

图4 (a)夏季平均的美国地区00时(UTC)(Seidel et al.，2012)和中国地区(b)08时、(c)14时、(d)20时(BJT)(Guo et al.，2016)ERA-Interim再分析(彩色阴影)和实测(彩色圆点)边界层高度的空间分布,对图(b)—(d)进行了修改Fig.4 Spatial distributions of summer mean ABL heights from the ERA-Interim reanalysis (color shaded) and sounding observation (color dots) at (a) 00∶00 UTC in United States (Seidel et al., 2012), and at (b) 08∶00 BJT, (c) 14∶00 BJT, and (d) 20∶00 BJT in China (Guo et al., 2016)， a modification of figure (b)—(d)

4 超高对流边界层特征及其成因分析

4.1 超高对流边界层特征

一般认为对流边界层高度应该低于3000 m(Garratt，1992)，但是近年来随着各类野外观测试验的开展，在一些特殊的地理环境地区(如印度新德里、非洲撒哈拉沙漠、中国西北荒漠和青藏高原等地区)(表3)均发现有超过4000 m的特殊边界层(称为超高对流边界层)存在，其中撒哈拉沙漠和青藏高原地区超高对流边界层高度可以发展到5000 m以上.研究指出(张强等，2004；Marsham et al.，2008；张强和王胜，2008；Messager et al., 2010；Chen et al.，2013；杨洋等，2016)，超高对流边界层往往出现在连续多日的烈日晴天，出现前一日的边界层高度超过3000 m，当日早晨近地层有逆温层，其上存在高层的近中性残余层(高度超过2000 m)，当日白天对流发展突破逆温层进入残余层后对流大气边界层高度呈现跳跃式发展，当日午后以干对流运动为主的混合运动使得边界层内位温均匀分布，边界层高度超过4000 m(图5).

图5 青藏高原改则站2008年2月25日不同时刻(a)温度、(b)位温、(c)风速、(d)比湿廓线，廓线采样时间分别为01时(BJT；黑线)、07时(红线)、13时(蓝线)和19时(紫线)，各颜色水平虚线表示相应时刻的对流边界层顶，各颜色水平实线表示相应时刻的对流层顶，SL、RL、ML分别表示稳定层、残余层和混合层(Chen et al.，2013)Fig.5 Profiles of (a) temperature, (b) potential temperature, (c) wind speed, (d) water vapor content at Gaize station of the Tibetan Plateau on 25 Feb 2008. Profiles were recorded at 01∶00 BJT (dark line), 07∶00 BJT (red line), 13∶00 BJT (blue line), and 19∶00 BJT (magenta line). The horizontal dashed lines are for the tops of the convective boundary layer and horizontal solid lines are for the positions of the tropopause. The stable layer (SL), residual layer (RL), and mixed layer (ML) are also marked (Chen et al.，2013)

表3 基于探空数据计算的超高对流边界层高度Table 3 Radiosonde-retrieved heights of super deep convective boundary layer

4.2 超高对流边界层发展成因

地表感热加热是影响超高对流边界层发展的重要因子(张强和王胜，2008； Zhang et al.，2011).在中国西北干旱区，超高对流边界层高度与地面同时次感热通量的相关系数为0.03，而与地面累积感热通量的相关系数高达0.78(图6)，这表明超高对流边界层的发展主要受白天感热通量的累积效应影响(张强等，2019).但是，一些学者研究指出在日际时间尺度上日累积地表感热通量与对流边界层日极大高度的相关系数则没有想象的那样高(Cuesta et al.，2008；张强等，2019).杨洋等(2016)认为地面感热通量能解释30%左右的超高对流边界层发展的总能量；张强等(2019)的计算表明，感热通量能量占西北干旱区(最高大气边界层高度为4170 m)发展所需总能量的40.0%.因此，在超高对流边界层发展过程中地表感热加热并非唯一能量来源.

图6 对流边界层厚度与实时感热通量(a)和累积感热通量(b)的关系(张强等，2019)Fig.6 Relationship between the convective boundary layer thickness and the real-time (a) and cumulative sensible heat fluxes (b)(Zhang et al.，2019)

残余层能量的夹卷效应是影响超高对流边界层发展的另一个重要因子.研究指出，对流边界层上方覆盖大气的层结状况会影响边界层顶附近湍流夹卷速度，进而影响边界层高度的发展(Freire and Dias，2013；Blay-Carreras et al.，2014；Reen et al.，2014).使用边界层单柱发展模式的边界层模拟结果显示，边界层上空大气层结状态对超高对流边界层发展影响最大(Chen et al.，2016).在超高对流边界层出现之前，当地中低层大气存在近中性的残余层，边界层与残余层耦合会呈现跳跃式发展，而在没有近中性残余层出现时，即便地表感热通量较大，当天最大对流边界层高度也偏低(Parker et al.，2005；Zhang et al.，2011；Han et al.，2012；赵采玲等，2016).张强等(2019)进一步给出了残余层夹卷能量的定量计算(图7)，如图所示，残余层的存在使得残余层的近中性位温廓线代替了无残余层时的理想逆位温廓线(AB线段)，残余层夹卷的总能量等于两条位温廓线间的斜线区面积.因此，在超高对流边界层发展的能量机制中必须考虑残余层的夹卷能量效应.

图7 残余层夹卷能量计算方法的示意图，线段AB表示无残余层时的理想逆位温廓线(张强等，2019)Fig.7 Schematic diagram of calculation method for entrainment energy in residual layer. The AB line represents the ideal inverted potential temperature profile without residual layer (Zhang et al.，2019)

通过热力学数值模型可以揭示地表感热加热和浮力夹卷等热力过程在超高对流边界层发展中的作用.赵建华等(2011)针对西北干旱区一次夏季超高对流边界层发展的研究表明，感热对超高对流边界层高度极大值的贡献率为63.4%，而浮力夹卷为18.1%，两者累积贡献超过80%(图8)，是驱动西北干旱区超高对流边界层发展的主要因子.张强等(2019)概括了西北干旱区对流边界层与残余层之间的正反馈循环增长形成超高对流边界层的物理过程.在超高对流边界层形成的持续晴空期内，残余层储存前一天对流边界层能量，成为当天对流边界层发展的重要能量来源，因而残余层的厚度由前一天对流边界层厚度所决定，并且还影响着当日对流边界层高度的增长，边界层与残余层间通过正反馈机制而循环增长，最终形成超高边界层；强烈的地表感热加热作用则是对流边界层与残余层之间正反馈循环增长机制的关键因素.

图8 2006年6月28日西北干旱区(a)感热、潜热日变化，(b)浮力夹卷热通量日变化，(c)热力数值模型模拟的对流边界层高度日变化，(d)浮力夹卷引入后热力数值模型模拟的对流边界层高度日变化(赵建华等，2011)Fig.8 Diurnal variations of (a) sensible and latent heat fluxes, (b) buoyancy entrainment, and convective boundary layer height simulated by the thermal numerical model (c) without and (d) with the buoyancy entrainment in the drought region of northwestern China on 28 June 2006 (Zhao et al.，2011)

5 大气边界层参数化方案

鉴于大气边界层在天气和气候系统中的重要作用，一个成熟的数值模式必须考虑合理的边界层过程描述方案.通常数值模式通过雷诺平均的方法刻画边界层过程，而湍流通量项的出现造成预报方程组无法闭合，因此对于湍流通量项的有效闭合成为边界层参数化方案需要解决的首要问题.此外，由于大气边界层发展受下垫面湍流输送、边界层顶夹卷以及层积云长波辐射冷却等物理过程影响，因而这些物理过程的描述是边界层参数化方案需要解决的另一关键问题.

为了解决垂直湍流通量项的闭合问题，一阶闭合方案基于K理论将湍流通量表达为湍流输送系数与网格物理量梯度的乘积.早期K方案的湍流输送系数K与湍流混合长有关(Louis，1979)，之后的K方案预先给定湍流输送系数K的垂直廓线形式(Troen and Mahrt，1986；Holtslag and Boville，1993).为了便于在中尺度数值模式的边界层方案中刻画湍流通量和其他高阶矩量，高阶湍流闭合技术得以发展，如Rotta(1951)提出了2阶闭合方案，Yamada和Mellor(1975)通过将湍流动能引入湍流输送系数计算，建立了1.5阶K闭合方案，该方案在闭合方程不太复杂的条件下兼顾了对高阶矩的要求，应用比较广泛.然而，由于对流边界层内的湍流输送具有很强的非局地性特征，大尺度湍涡可以引起非局地输送，在对流边界层的中上部会出现逆梯度输送的情况.K方案是局地闭合理论，湍流输送仅与局地垂直梯度有关，未能考虑对大气边界层发展有重要影响的逆梯度输送以及边界层顶夹卷等过程(Wyngaard and Brost，1984；Holtslag and Moeng，1991)，使得计算的边界层高度系统性偏低和水汽在低层堆积，从而造成模式模拟大气边界层低层比实况的湿冷(图9a，b).

图9 模拟的波多黎各圣胡安(San Juan, Puerto Rico；18.3°N, 66°W)7月平均(a)温度、(b)比湿、(c)湍流输送系数、(d)比湿通量和(e)比湿倾向廓线.虚线表示局地边界层方案(Local)模拟结果，实线为非局地边界层方案(Nonlocal)模拟结果.(a)和(b)给出了探空观测结果(带点水平杆)(Holtslag and Boville，1993)Fig.9 The simulated profiles of July mean temperature (a), specific humidity (b), turbulence transfer coefficient (c), specific humidity flux (d), and specific humidity tendency (e) in a grid point near San Juan, Puerto Rico (18.3°N, 66°W). Solid lines are for the local diffusion scheme, and dashed lines are for the nonlocal diffusion scheme. In (a) and (b) radiosonde observational results are also given (dots with horizontal bars) (Holtslag and Boville，1993)

为了弥补K方案对负梯度输送等过程的缺陷，Deardorff(1972b)、Holtslag和Boville(1993)发展了非局地模型，其做法是通过在K理论公式中加入逆梯度项的方式增强湍流垂直输送(图9c，d，e).Pleim(2007)采用直接混合的方法来描述对流边界层的非局地特性；Siebesma等(2007)借鉴浅对流方案中的质量通量方案(Siebesma and Cuijpers，1995)，将局地的K理论与非局地的质量通量结合起来.此外，针对对流边界层顶的夹卷过程，更多的大气边界层方案将夹卷通量转化为对应的湍流输送系数，或者在计算湍流通量时直接显式地引入夹卷通量(Noh et al.，2003；Hong et al.，2006).

然而，非局地大气边界层方案未能考虑层积云顶辐射冷却过程的影响，使非局地模式在有云覆盖情况下湍流垂直输送偏弱，容易产生虚假边界层低云现象.为了解决上述问题，边界层方案可以通过两种方式增加层积云覆盖边界层的云顶辐射冷却过程：一种是在1阶闭合框架下采用双K廓线法，即在原有给定的K廓线外，增加由云顶辐射冷却作为湍流动力来源的湍流混合廓线(Lock et al.，2000；Sun et al.，2016)；另一种则是在2阶和更高阶的闭合框架下，将云顶辐射冷却产生的浮力项通过隐式表达的方式写入湍流动能收支方程(Bretherton et al.，2009；Konor et al.，2009).Sun等(2005)将Mellor/Yamada的湍流2阶闭合方案与中尺度模式MM5(The fifth-generation mesoscale model)耦合，模拟了华南地区的一次暴雨过程.他们指出：与MM5中原有的1.5阶边界层方案相比，2阶闭合方案能够比较好地再现大气边界层中较短时间尺度的风速脉动，从而更好地模拟低涡、低空急流特征及降水分布.

6 结语

影响大气边界层高度的物理过程十分复杂，是大气边界层研究的重要内容之一.本文回顾了大气边界层高度的计算方法、大气边界层的时空分布特征、大气边界层发展成因以及大气边界层参数化方案和大气边界层模式等方面的主要研究进展，从多角度给出了大气边界层高度的诊断识别方法，揭示了混合对流、中性和稳定大气边界层形成机理及结构特征，下垫面状况引起的大气边界层高度的空间差异性，特别是解释了地面感热和混合层顶夹卷对超高大气边界层形成的机理，分析了大气边界层参数化方案关于湍流输送、夹卷和层积云辐射冷却等重要过程的描述缺陷以及对模式预报结果造成的影响.