1、书书书龚玺，陈起英，倪长健，等，基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析气象，（）：，（）：（）基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析龚玺，陈起英，倪长健龙海川成都信息工程大学，成都 湖南省株洲市气象局，株洲 气象防灾减灾湖南省重点实验室，长沙 中国气象局地球系统数值预报中心，北京 灾害天气国家重点实验室，北京 提要：由于湍流的时空尺度小、变化率大，缺乏大规模的现场观测，导致对全国范围内的大气湍流特征研究甚少。自 年起，我国开始存储秒级高垂直分辨率探空数据，为深入研究中国各地的湍流特征提供了可能。在用探空资料诊断分析湍流特征时，首先分析了中国地区仪器噪音对分析
2、结果的影响，证实噪音会使对流层平均的湍流耗散率偏小，证明了去除噪音的必要性。在此基础上，采用 分析法详细分析了湍流的强度（湍流耗散率）与下垫面类型的关系。结果显示，在对流层，全年平均草地的湍流最强，然后依次是旱地作物、水田作物、灌木、荒地。在平流层下层依然受到下垫面的影响，全年平均草地的湍流依旧最强，旱地作物次之。湍流在不同下垫面均表现出明显的季节性，且夏季对流层上层的强度和频率最高。研究结果将有助于理解不同下垫面对流层的湍流特征，为航空器飞行安全保障、大气污染扩散模式、天气预报模式的湍流参数化方案改进提供参考。关键词：秒级探空资料，分析法，湍流动能耗散率中图分类号：文献标志码：犇 犗 犐：，
3、犃 犫 狊 狋 狉 犪 犮 狋：，国家重点研发计划（、）和国家自然科学基金项目（）共同资助 年月日收稿；年 月 日收修定稿第一作者：龚玺，主要从事边界层研究 ：通讯作者：陈起英，主要从事数值模式研究 ：第 卷 第期 年 月 气象 ，（），犓 犲 狔狑 狅 狉 犱 狊：，引言大气湍流是大气中重要的中小尺度过程，它影响着大气中的能量平衡和物质传输过程，使能量、动量、物质在不同尺度的运动之间转换，从而调整大气结构，使大气成分在一定的高度范围内均匀分布（，）。在大气边界层，湍流运动相当强，其对地面强迫的响应时间在小时以内（刘辉志等，）。对流层下部的湍流对地面附近动量和大气污染物的向上扩散有重要影响（，
4、），平流层顶的湍流对对流层和平流层之间的物质和能量交换有显著的影响（，），自由大气中的湍流对航空的飞行安全会产生重要影响（，）。因此，充分认识湍流的时空分布特征对于大气污染预报、天气预报和气候预测都相当重要。考虑到大气热力结构、动力结构和湍流动能对湍流混合过程的共同影响（陈德辉等，），湍流运动具有相当强的间歇性和随机性，一般难以对湍流的性质进行准确、定量的描述，因此在数值模式中，常需要根据理论研究结果和观测资料对其效应进行参数化描述。湍流扩散系数是湍流参数化方案需要预报的关键物理量，但是其缺乏实时观测对模式结果来进行验证。湍流扩散系数可以由湍流耗散率进行诊断，因此本文首先利用高垂直分辨率探空资
5、料对湍流动能耗散率进行了研究。湍流动能耗散率是描述湍流运动强度的一个重要物理量，它既可以影响对流层和平流层之间的质量、动量和能量以及化学物质的交换，也能影响各种波（重力波、行星波等）从热带对流层向高层大气的传播过程，进而影响低层与高层大气的耦合（，）。因此，研究湍流动能耗散率的时空分布特征是十分迫切和必要的。从高垂直分辨率探空资料获得湍流动能耗散率的信息，需要采用 分析法。分析法最初由 （）提出，通过对实况海洋位势密度垂直廓线按照稳定层结条件进行垂直方向的重新排序，再将重排获得的稳定位势密度分布和实况位势密度分布比较，用位势密度所在位置的变化信息反演湍流翻转产生的位置和范围，这个方法被广泛应用
6、于海洋湍流特征的分析。（）提出将其用来估计自由大气中的湍流参数，把过去几十年一直用于海洋稳定层结的方法创新性地应用于具有稳定层结特点的自由大气中，利用大气中的保守量位温按照稳定层结要求的分布调整，产生的位置翻转（翻转指对原始位温与按照稳定层结分布要求重新排列产生的位温位置的差异）反演大气的湍流特征。在利用 分析法研究和分析大气湍流耗散率方面，国外已经进行了不少研究工作。（；）对 分析法进行了进一步改进，考虑了仪器噪音的影响，使对湍流识别的准确率得到提高。（）又考 虑 了 湿 度 对 频 率 的 影 响，把 方法的研究范围从晴空扩展到有云大气，至此 分析法估算大气湍流特征的方法趋向成熟。近年，（
7、；）采用美国的个中纬度站点的 年月至 年 月的高 气象 第 卷分辨率探空数据，验证了 分析在大气湍流统计研究中的可信度。年还通过对美国高分辨率探空仪数据集进行 分析，研究了湍流的纬度依赖和地形相关性；同年 （）基于美国的高垂直分辨率探空数据（年月至 年月），在 个业务运行台站用 分析法分析了美国区域自由大气中的 尺度和涡流耗散率，这是少数几个在比较广的时空范围内，利用高分辨率探空观测资料研究湍流参数分布特征的工作。在国 内 的 研 究 中，（）利 用 分析法，基于南海考察船上观测到的热带海洋对流层上部的无线电探测资料，研究了这个考察站对流层上层湍流参数的日变化，及其与中尺度对流系统日变化的密切
8、关系，之后又分别采用 分析法和利用气球垂直上升速度的扰动计算湍流参数的方法，计算南海同样这个考察站对流层和低平流层的湍流动能耗散率和湍流扩散系数，并比较分析了两种方法的优劣。近期，（）使用中国 个站点的高垂直分辨率探空仪数据（年月至 年 月），计算了自由大气中大气湍流的 尺度和湍流动能耗散率，弥补了我国此前没有在较广的时空范围内研究湍流参数分布的不足，但是他们的研究没有分析不同下垫面上湍流参数分布的差异；同时，他们的研究结果虽然考虑了 （）建议的仪器噪音的消除，但是没有详细分析中国的仪器噪音对结果的影响，因此本论文将对这个影响进行讨论。（）最早发现在利用秒级探空资料研究湍流时，仪器噪音的影响会
9、对分析结果产生比较大的影响。仪器噪声会产生虚假的位置翻转和虚假的湍流识别。若将虚假湍流层计入真实的湍流层中，会导致对湍流发生频率的严重高估，从而导致对大气湍流参数和湍流能量的高估，湍流的时空不均匀性会导致扩散系数估计值偏差超过几个数量级。这些虚假翻转很难与湍流产生的真实翻转区分开来，因此对观测资料进行去除仪器噪音过程非常重要（，；，；，；，）。（）基于统计方法提出了去除仪器噪声影响的方法，首先推导了可探测的湍流层大小的阈值与整体趋势噪声比的函数关系，提出了只有反演的垂直尺度超过相同大小正态分布的噪声样本的范围（最大和最小值之差）的规定百分位数，才能将其识别为真实的翻转。本文将把这一方法应用到中
10、国的高垂直分辨率秒级探空资料，研究在利用 分析法对中国区域的湍流耗散率进行分析时仪器噪音对分析结果的影响。已有很多研究证明空气的流动情况受下垫面的影响很大：刘明星等（）利用相同探测设备、相同试验方法在戈壁、草原、郊区和城郊等不同下垫面开展的多次重复湍流观测试验资料，采用相同的湍流数据采集和处理方法，研究了不同下垫面温度和湿度湍流谱特征，发现高粗糙度下垫面的高频段谱曲线高于平坦地形，说明高粗糙度地区的局地湍流耗散较强，戈壁下垫面的低频段湿度谱数值明显高于其他地区；洪雯和王毅勇（）总结了水平非均匀下垫面的大气边界层研究，指出干旱区的绿洲具有明显的“冷湿岛效应”，并可形成边缘逆湿和区域环流，而植被和
11、水域的分布则显著改变了大气边界层结构。同时还阐述了复杂地形下垫面对局地风场和小尺度环流产生的影响、城市热岛效应的产生及特点和 城 市 中 绿 地 和 水 体 的 微 气 候 效 应；龚 玺 等（）利用国家气候中心 年在内蒙古锡林浩特市、河北黄骅市、江西湖口县和星子县开展的 铁塔大气湍流观测试验数据，分析不同下垫面、不同稳定度条件下地面到 的风切变指数变化情况，发现下垫面与风切变指数变化有密切联系。但是之前的研究没有详细分析不同下垫面情况下的湍流参数分布特征。我国的下垫面种类相当丰富，根据国家地球系统科学数据中心的数据，我国的下垫面类型多达 种。空气流动受下垫面影响的方式主要有两种：一是动力作用
12、，如地形起伏、植被或者建筑物等改变地面粗糙度可增加动力湍流；二是下垫面的热力作用，因地形起伏或水陆分布，使得下垫面受热和散热不均匀，从而引起温度场和风场的变化，进而影响湍流的产生，因此下垫面对于湍流的强弱以及频率都会有较强的影响作用。在研究下垫 面 与 湍 流 分 布 特 征 的 关 系 方 面，潘 乃 先 等（）利用三分量风速仪采集风的资料，研究了不同下垫面近地层上的湍流度统计特征，结果表明中性大气条件下湍流度的统计平均值与地表起伏状况存在紧密的关系，并表现出一定的规律性；涂钢等（）分析了 年半干旱区退化草地和农田下垫面近地层湍流通量日、季、年的变化特征，探讨了不同年份的气候背景和下垫面土地
13、覆被对地气交换过程的影响，得出了结论：半干旱区退化草地第期 龚玺等：基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析 和农田近地层湍流通量具有明显的日变化、季变化周期，地气间通量交换年际间的差异主要受当年的气候背景影响，同时还受到降水、下垫面覆被的影响；王少影等（）利用“绿洲系统能量与水分循环过程观测试验”年绿洲、戈壁点的观测资料，分析与讨论了温度、水汽的归一化标准差随稳定度变化的通量方差关系；张璐等（）通过改变鲍恩比的敏感性试验研究不同大气初始状况下，湿润和干旱下垫面湍流特征及其对深对流触发过程的影响，结果表明：干旱下垫面的混合层干而暖，厚度较大，湿润下垫面相反，边界层内的深对流首先
14、在干旱下垫面发生；张蕴帅等（）对比分析了黄河源区两种不同下垫面上（鄂陵湖和湖边草地）对流边界层（）中精细的湍流结构特征，发现两个不同下垫面上模拟结果对模式分辨率的敏感性不同，也模拟出了各种尺度的波对湍流通量的累积贡献。但是这些研究中有的只局限于局部地区和近地层的高度，没有对东亚整个区域以及自由大气高度范围进行研究；有的虽然考虑了一定的范围，但并没有分析湍流动能耗散率的分布情况。因此，本文利用 年月日至 年 月 日我国的高垂直分辨率探空资料，分析在中国区域的仪器噪音对 分析法估算湍流动能耗散率的影响；同时研究东亚不同下垫面特征条件下，湍流动能耗散率的变化特征。资料和方法 高垂直分辨率探空数据本文
15、用于估算湍流动能耗散率的数据来自国家气象信息中心提供的中国无线电探空观测网的高垂直分辨率的探测数据（，；）。如图所示，自 年至今已有 个无线电探空站升级为 波段探空。郝民等（）使用 年月的波段探空资料进行分析，结果发现秒级探空资料在 、等层次上高度分析比传统探空分析更接近 分析，显示出对分析场的改进效果，同时风场在 以上的高层分析中改进显著，表现出高分辨率资料使用在高层分析中对系统模拟更加准确的特征。波段探空仪主要由 （）型高空气象探测雷达和 型数字探空仪组成，气压、温度、风速等气象要素的数据采集由 型数字探空仪完成。型探空仪的测量周期为 ，每分钟采集 次，气球以 的速度上升，因此探测资料的垂
16、直分辨率达到，具有很高的垂直分辨率和实时采集能力（杨湘婧等，；，），其测量的温度、注：红色点为站点，底图为对应的地形高度。图我国 个探空站点的分布 （）气象 第 卷气压、位势高度一致性较好，风向风速与 导航卫星定位测风的结果比较接近（李伟等，）。每天分别在 时和 时（北京时，下同）两次提供温度、相对湿度、压力和风速的垂直廓线。本文使用的高垂直分辨率探测数据的时间段为 年月日至 年 月 日，绝大多数台站都有比较多的观测样本。在本研究中，为了保证足够长的垂直范围，高度范围小于 的探空数据将被舍弃。大多数探空观测的最大高度高于 ，这表明在分析范围内有足够多的观测样本。由于气球的垂直上升速度是变化的，
17、观测的间距是非等距的，为了便于分析研究，本文首先使用三次样条插值把原始数据插值到 等距的垂直格点上。犜 犺 狅 狉 狆 犲分析法空气一般有三种状态，干空气、未饱和和饱和状态。由于水汽对未饱和空气的影响比较小，因此只需要考虑干空气和饱和空气两种情况（，）。本文基于 （）的研究，用探空仪测量的相对湿度（）是否满足两个随高度变化的阈值条件 和 （表），从而确定大气是否饱和（即是否存在云）。判断条件为，如果某一高度范围内的所有层的 （随高度变化），并且此范围内同时存在某一个高度的 （随高度变化），则这一高度范围被确定为有云层。表判别云层存在的相对湿度阈值随高度的变化（犣 犺 犪 狀 犵犲 狋 犪 犾，
18、）犜 犪 犫 犾 犲犞 犪 狉 犻 犪 狋 犻 狅 狀狅 犳 狉 犲 犾 犪 狋 犻 狏 犲犺 狌 犿 犻 犱 犻 狋 狔 狋 犺 狉 犲 狊 犺 狅 犾 犱犳 狅 狉 犼 狌 犱 犵 犻 狀 犵狋 犺 犲 犲 狓 犻 狊 狋 犲 狀 犮 犲狅 犳 犮 犾 狅 狌 犱 狊狑 犻 狋 犺犺 犲 犻 犵 犺 狋（犣 犺 犪 狀 犵犲 狋 犪 犾，）高度 参照 （）、（）的研究方法，针对晴空和有云层分别使用以下公式计算位温：（犻）（犻）犖（犻）狕犵（）（犻）（犻）犖（犻）狕犵（）式（）用于干空气层，式（）用于有云大气层。式中：犵是重力加速度，狕为高度差，犖是干的 频率的平方，可表示为：犖（犻）犵（犻
19、）（）狕犻（）而犖是湿的 频率的平方，可以表示为（，）：犖（犻）犵犜（犻）犜（）狕犻（犻）犔狇（犻）犚 犜（犻）犵狇（犻）狇（）狕犻（）式中：犔是凝结潜热，总混合比狇狇狇，狇是云水冰混合比，狇是水汽饱和混合比。由于无法得到云水冰混合比的信息，且其值相对水汽饱和混合比比较小，因此忽略，近似狇狇。是湿绝热递减率（，）：（犻）犵犮 狇（犻）犮 狇（犻）犮狇（犻）犮 犔狇（犻）犮 犚 犜（犻）狇（犻）（）（）式中：犮 为干空气的比热容，犮为水冰的比热容，犮 为水汽的比热容，狇项同样被忽略。由式（）、式（）计算得出位温后，将其按从小到大的顺序进行重排，得到的重排的位温被称为。如果原廓线上位于高度狕处的值
20、，在廓线中位于高度狕处，则其 位移可以定义为：犔（狕）狕狕（即调整之前的高度减调整之后的高度）。如果在任意高度上发生了位移，则意味着相反方向的位移必然发生在其他某个点。因此，当从大气低层向上积分到累积的 位移满足以下公式时，可确定为湍流层（，）：狀犻犔（犻）（）犽犻犔（犻）（）式中：犽狀（，）。利用式（）和式（）确定湍流层后，计算每个湍流层中 位移的均方根即为 尺度。如 （）所述，方法的基本假设在于 尺度（犔）与 尺度（犔）成线性比例，表示为：犔犮 犔（）式中：犮是一个常数，且在以前的研究中已被广泛讨论（，；，；，）。尺度同时也可以表示为：犔犖（）（）第期 龚玺等：基于高垂直分辨率探空数据对五
21、种典型下垫面湍流分布特征分析 式中：表示湍流耗散率，犖表示 频率。结合式（）和式（），可以得出为：犆犔犖（）式中犆犮，最近的数值试验表明，犆随湍流的时间演变而变化，并随着时间的推移而增加（，）。然而，（）的研究又发现犆呈明显的对数正态分布。因此，犆的平均值通常被视为一个接近常数的值，这使得 分析成为可能，从而从统计角度估计湍流。最近一项基于蒙特卡罗模拟的研究（，）进一步证实了这一点，该研究表明犆保持接近常数，其最佳预测值为 ，因此接下来的分析中都取犆 。为了更好地描述垂直分布的特征，将其分成四个高度间隔进行研究：对流层下部（）、对流层上部（）、平流层下部（）和平流层上部（），它们的 范 围 分
22、 别 为、和 。仪器噪音利用探空观测资料估算湍流时，仪器噪声的影响是一个不可忽视的问题。仪器噪声会产生虚假湍流层，很难将其与真实湍流层区分开来（，；，；，；，），然而在很多早期估算大气湍流的研究中，噪声问题都被忽略了。若将虚假湍流层计入真实的湍流层中，会导致对其发生频率的严重高估，从而导致对大气湍流参数和湍流能量的高估。（；）提出了一种基于统计识别仪器噪声的方法，即使用趋势噪声比（，）来评估测量的噪声水平。可分为两类：一类是局部的趋势噪声比（以下称为局部 ），用于评估各高度处的噪声水平；另一类是用于确定总体质量的整体趋势噪声比（以下称为整体 ）。这里简要介绍去噪声的过程：（）首先将每一个温度的
23、观测廓线均分成长度为 长的区间，由于数据的分辨率为，因此每 个点为一组；（）再将每一个区间分别进行线性拟合，然后从每一段中剔除线性拟合的结果；（）计算剔除拟合结果的温度廓线的一阶差分的方差，该值的一半即是观测温度噪声的方差；（）最后将温度的噪声标准差进行光滑处理。得 到 仪 器 的 温 度 噪 声 标 准 偏 差后，由 （）狆可计算出位温的噪声犜槡狆犜 （）狆（）式中：用重排后的位温，是位温的噪声，狆，犜分别为原始的气压和温度，整体 计算如下：珔（狀）（）（狀）（）式中狀是研究范围内所有样本的总数，整体 是相对于仪器噪声的湍流信号强度的指标。一般来说，整体 越低，说明信号的强度越弱，出现虚假湍
24、流层的概率越高。因此，如果整体 太小（在对流层一般小于 ），则需要进行去噪音处理。最后，在一个湍流层范围内的变化应超过该厚度内同样样本大小的噪声样本范围的。可以识别的湍流层厚度的大小与阈值 之间的统计相关性如 （）研究的图所示。低于标准的湍流层被标记为假，需要被剔除。关于仪器噪音对湍流耗散率的影响将在 节中进行介绍。结果与讨论 仪器噪音对湍流动能耗散率估算的影响由于仪器噪音会使得探测时产生虚假的湍流层，从而使后期计算的湍流耗散率产生偏差。而在使用式（）和式（）进行了整体 的计算后，发现未考虑仪器噪音时，对流层的整体 难以达到规定的（一般为），而平流层则可以达到以上。因此在本节中，我们将主要分析
25、在对流层中中国的探空仪器噪音对于估算的影响。图显示了降噪前后的 年对流层（和）的湍流动能耗散率值的分布情况。需要注意的是，因为值比较小，用对数标度表示，其值是根据式（）计算得出的。由图可知，的范围为。从整体分布看，无论是在对流层低层还是上层，的分布近似对数正态分布，呈正偏斜。在考虑噪声前，在低层和上层的平均值分别是 （）和 （），上层的均值要稍大于低层，这表明上层中的湍流层要 气象 第 卷注：黑线显示了对数正态分布。图 年考虑仪器噪音（）前，（）后的对流层下层（：）、对流层上层（：）的 的百分比分布 （：）（：）（）（）强于低层；而在考虑噪声后，在低层和上层的平均值分别变为了 （）和 （），这
26、意味着噪音产生了一些小尺度的虚假湍流层，使对流层平均的湍流耗散率偏小，证明了去除噪音的必要性。不同下垫面上的湍流动能耗散率水平分布特征图显示了 年 时和 时个高度层的大气湍流耗散率的水平分布。总体而言，对流层上层的湍流最强，对流层下层与平流层下层次之，平流层上层最弱。值得注意的是，以复杂地形为特征的青藏高原在所有高度层都出现了强的峰值，这和在落基山脉地区湍流比较强的结论是一致的（，），这是由于大地形有比较大的粗糙度，对形成湍流的动力作用有较大影响，同时研究表明这个地区的感热通量也比较大，热力作用也对湍流的产生起了很大作用。平流层在这个区域湍流比较强，与青藏高原产生的地形重力波的上传和破碎有关（
27、，）。相反，东北地区在各个高度层都表现出低值，在平流层下层尤为显著，说明该地区的湍流较弱，与这个地区温度偏低，大气常处于稳定状态有关。除平流层上层以外，湍流动能耗散率的年平均值呈现“南高北低”的型式，即在中国南方地区的湍流动能耗散率普遍比北方强，这可能是由于对流云经常占主导地位，频繁的对流作用对湍流有促进作用，尤其是在夏季（，；，），同时由于低纬地区下垫面温度比较高，热力作用比较大。此外，也可以发现除了主要的“南高北低”型式，对流层和平流层下层有比较明显的经向差异，尤其是对流层下层，这与低层大气受下垫面的影响比较大，随高度增加下垫面的影响逐渐减少有关。在平流层上层湍流动能耗散率的分布比较均匀，
28、下垫面的影响几乎看不见，青藏高原的高值区也消失了，说明青藏高原对湍流的影响要低于 。综合而言，对流层（和）大气湍流的强度比平流层（和）强得多。湍流动能耗散率水平分布的南高北低，垂直分布在对流层上层最强，大地形（例如落基山脉）上空为高值区的结论，与在美国和挪威使用相同 分析法基于美国和挪威的高垂直分辨率探空资料的研究结果非常相似（，；，；，）。结合 时和 时不同高度层平均的水平分布可看出，在青藏高原的东西部各存在两个湍流高值区，这与青藏高原西部和东部的主要下垫面类型不同有关，西部主要是灌木，东部主要是草地。位于青藏高原西部的湍流区在 时较弱，在 时则较强，这表明该地区的湍流有较强的日变化；但是这
29、个区域除了在平流层高层以外，各高度层平均强度垂直变化不大，这可能与灌木的影响高度比较高，以及这个区域日出时间较晚有关。而位于青藏高原东部的湍流区 时和 时变化不明显，日变化较弱，但垂直变化较明显，湍流强度在对流层低层最大，对流层高层和平流层下层逐渐减小，在平流层上层基本消失。这表明随着高度的增加，影响其强度的主导因素逐渐减弱，由此推测该地区产生强湍流的主导因素是下垫面。该湍流区及其附近地区下垫面类型均第期 龚玺等：基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析 注：、的垂直范围分别为、和 。图 年（）时，（）时个高度层平均的大气湍流耗散率 的水平分布 （）：，（）：为草地下垫面，但湍
30、流区范围并未与下垫面范围完全一致，可能是因为该地区位于横断山脉地区，随纬度变化高度落差较大，也造成了湍流的差异。而分布在长江中下游地区和华南地区的站点主要为水田和牧场下垫面类型（以下简称水田作物下垫面），该类型下垫面的作物类型多为水生作物，土地含水量大，且该类型站点数量最多，数据具有更高的可信度，虽然个别站点出现了极值，但整体分布情况较为 气象 第 卷均衡，较北方地区而言稍大。西南地区东部的主要下垫面类型为旱地农田和牧场（以下简称旱地作物下垫面），由于旱地农田无灌溉设施，主要靠天然降水种植旱生农作物，因此该下垫面类型的含水量与降水密切相关，中国农田变化对气候影响具有冬季弱、夏季强的季节性变化（
31、曹富强等，），推测旱地下垫面的湍流强度和频率受季节变化影响较大。在内蒙古西部的下垫面类型为灌木下垫面，该区域在各高度层中湍流强度均较均匀，与下垫面条件比较均匀有一定关系。在新疆地区的下垫面类型为荒芜稀疏的植被下垫面（以下简称荒地下垫面），其湍流强度与灌木下垫面类型的站点数据较为接近。湍流的垂直分布与下垫面的关系将在下面的 节中做进一步深入研究。不同下垫面上湍流动能耗散率的垂直分布特征为详细分析下垫面对湍流作用的影响，考虑到中国区域主要的下垫面类型，并根据 节分析的我国湍流水平分布特点，重点选择了个不同类型的下垫面区域：旱地作物、水田作物、草地、灌木、荒地（图），将其中站点的数据样本进行平均来研
32、究讨论。这些区域的选择既保证了下垫面的代表性，又保证了足够的样本数量，并且个区域的纬度差异较小，在一定程度上排除了太阳高度角的季节变化对湍流强度的影响。表给出了选取的种类型下垫面区域的平均海拔、主要植被类型、干湿分区和主要气候类型。可以看出，选择的种下垫面类型中的草地和荒地由于位于我国西部高原地区，平均海拔较其他类型明显更高。种下垫面类型在干湿分区上也有较明显的差异，表明其对应的降水量和土壤含水量也有差别；这也导致了不同类型下垫面主要植被类型的差异。例如草原地区的植物大多是适应半干旱气候条件的草本植物；而荒漠的生态条件极为严酷，夏季炎热干燥，土壤贫瘠，因此植被稀疏，植物种类也很贫乏。本文选取的
33、种下垫面类型对应的气候类型也很多样，几乎涵盖了我国所有的主要气候类型，具有较好的代表性。图给出在个选择的不同下垫面区域，个垂直分层平均的湍流动能耗散率的季节平均和年平均。结合图与表来看，在几乎全部下垫面的所有垂直分层，湍流均表现出了较强的季节性分布，即夏季最强，春秋季次之，冬季最弱。在下垫面影响最注：红圈标出研究选择的下垫面区域，为水田作物下垫面，为旱地作物下垫面，为草地下垫面，为灌木下垫面，为荒地下垫面。图我国 种下垫面类型分布 第期 龚玺等：基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析 表种下垫面类型区域的基本特征犜 犪 犫 犾 犲犅 犪 狊 犻 犮犮 犺 犪 狉 犪 犮 狋
34、犲 狉 犻 狊 狋 犻 犮 狊狅 犳 犳 犻 狏 犲狌 狀 犱 犲 狉 犾 狔 犻 狀 犵 狊 狌 狉 犳 犪 犮 犲 狋 狔 狆 犲 狊下垫面类型平均海拔主要植被类型干湿分区主要气候类型水田作物 水生作物耕地湿润、半湿润温带季风气候旱地作物 旱生作物耕地湿润亚热带季风气候草地 草原半干旱、半湿润高山高原气候灌木 矮小植被半干旱、半湿润温带大陆性气候荒地 多荒漠干旱温带大陆性气候图 年年平均和季节平均的旱地作物（红色）、水田作物（蓝色）、草地（绿色）、灌木（黄色）和荒地（紫色）在个高度层平均的 分布（）对流层下层，（）对流层上层，（）平流层下层，（）平流层上层 （），（），（），（）（）（），
35、（），（），（）大的对流层（和），全年平均草地下垫面的湍流最强，水田作物下垫面、灌木下垫面、荒地下垫面的强度较为接近。旱地作物下垫面的湍流在对流层低层和对流层上层的差异很大，其在对流层低层很强，与草地相当接近，而在对流层高层则比较弱，与水田作物下垫面、灌木下垫面和荒地下垫面比较一致。旱地作物下垫面在对流层下层湍流较强的原因，可能是所选择的区域为南方的旱地，虽然其干湿分区上位于湿润地区，但其主要位于一些山地、丘陵地区，这些地区没有水源保证，且多为坡地，下垫面相对比较干燥，感热通量比较大，而土壤热通量在下垫面能量平衡研究中极为重要，是影响其湍流强度的重要因素（岳平等，）。草地下垫面的湍流最强，原因
36、除了其主要植被类型为草原，能产生较大粗糙度以外，同时该区域的气候类型为高山高原气候，海拔较高，太阳辐射强、地气温差大，是强大的地面感热输送源，更容易产生湍流。除此之外，强湍流也和这个区域同时存在较为复杂的地形有关。从季节平均来看，荒地下垫面在各个季节均表现为最低，与之相反，草地下垫面在除秋季以外的各个季节均为最高。灌木与荒地下垫面除在对流层下层的秋季表现出一定差异以外，分布情况较为重合，这可能是由于两种下垫面在气候环境与地理位置上较为接近，都为温带大陆性气候，灌木虽然较荒地而言降水量更大，但两者的植被类型都为耐旱的低矮植物，对湍流的影响有限。水田作物下垫面稍强，但随四季的分布曲线最为平稳，季节
37、性变化最弱，可能由于水田类型的下垫面降水量大，导致土壤含水量大，热容量大，因此季节变化比较小。在平流层下层，从全年平均来看，草地下垫面的湍流依旧最强，旱地作物下垫面次之，水田作物下垫面稍弱于荒地，其下则是灌木。与对流层一样，灌木和荒地下垫面的湍流分布情况较为重合，但在春夏季出现了一定差异。在平流层上层各类型下垫面的差异较小，也就是说下垫面对湍流的影响只能到平流层下层的高度。草地在平流层上层以下各高度最 气象 第 卷强可能也与这个区域海拔比较高，地形的影响比较大有关。值得注意的是，灌木下垫面的湍流强度相较对流层，降低幅度很大。从季节平均来看，秋冬季各下垫面类型的强弱对比与年均基本一致，而春夏季不
38、同，荒地明显增加，强于除草地之外的下垫面类型。在对流层湍流强度最强的草地下垫面在平流层下层各个季节依然很强，其随季节的变化也更加明显。在平流层上层，由于该层发生的湍流频率和强度都较弱，与下垫面的联系也最弱，因此种下垫面的曲线相较其他高度层较为重合，且季节变化性也最弱，这也是合理的。图给出了在垂直方向上，在不同下垫面的湍流动能耗散率的季节平均和全年平均的垂直廓线。从图上可以看出，种下垫面在 附近出现了较大程度的重合，该高度正好位于对流层顶附近，这说明下垫面对于湍流的影响主要作用范围在对流层，在其之上的高度，下垫面的影响将会在一定程度上减弱。从季节平均来看，各下垫面类型个季节均在 范围内均出现最大
39、值，尤其是旱地作物。除夏季以外，旱地作物的极大值均大于其他下垫面的最大值。因为旱地在四季含水量都比较小，感热通量比较大，且其气候类型为亚热带季风气候，冬季温度相对较高，从而使得其在四季均有湍流高值，而在降水最强的夏季，各下垫面的含水量均会上升，潜热通量变大，湍流加强，因此草地下垫面的湍流在此时强于旱地作物下垫面，夏季还会产生频繁的对流云，而对流云也是湍流的主要来源（，）。水田作物在各个季节的湍流强度均弱于旱地作物和草地，但是其季节性变化最小，因为水田靠人工灌溉，各个季节土壤含水量变化相对较小。灌木和荒地由于其植被类型和气候类型的相似性，湍流强度在各个季节都较为接近，但年平均灌木稍强，这可能与灌
40、木区域的植被覆盖率更大有关。在 范围内，由于下垫面随高度增加产生的影响减弱，各下垫面的湍流均变化缓慢，草地下垫面开始明显强于旱地作物下垫面，其他种下垫面的强度排名则与对流层一致，即水田强于灌木，灌木强于荒地。在 以上，各下垫面的湍流强度均有所上升，但分布线更为重合，说明下垫面的作用难以影响至该高度层，这与图得出的结论一致。图旱地作物（红色）、水田作物（蓝色）、草地（绿色）、灌木（黄色）和荒地（紫色）平均湍流的 分布（）春季，（）夏季，（）秋季，（）冬季，（）年平均 （），（），（），（）（）（），（），（），（），（）第期 龚玺等：基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析 结论
41、本文利用 分析法，基于秒级探空资料，研究了 年期间中国全天空（包括晴天和有云）条件下湍流耗散率的时空分布与下垫面的关系。同时讨论了中国的探空仪器噪音对用 分析法分析湍流的影响。得到以下主要结论：（）仪器噪音会使得探测时产生虚假的湍流层，从而使后期计算的湍流耗散率产生偏差。从整体分布看，无论是在对流层低层还是上层，的分布近似对数正态分布，呈正偏斜。在考虑噪声前，上层中的湍流层要强于低层；而在考虑噪声后，在低层和上层的平均值均变大，这意味着噪音产生了一些小尺度的虚假湍流层，使对流层平均的湍流耗散率偏小，证明了去除噪音的必要性。（）在考虑湿饱和效应的情况下，即在全天空条件下，自由大气的湍流耗散率的年
42、平均分布在空间上表现出明显的“南高北低”和“西高东低”的型式。在中国南方地区的湍流动能耗散率普遍比北方强，这是由于对流云经常占主导地位，尤其是在夏季，频繁的对流作用对湍流有促进作用。而青藏高原上的强湍流可能是由于复杂的下垫面和地形的影响。从垂直分布来看，对流层上层的湍流比其他层次要强，和下垫面的动力和热力作用还可以影响到这一层，以及地形产生的重力波的破碎有关系。（）湍流在不同下垫面均表现出明显的季节性，且夏季对流层上层的强度和频率最高。旱地作物特有的耐旱和蓄水性会使得下垫面的含水量在四季保持稳定，从而使得其在四季均有最大值，而在降水最强的夏季，各下垫面的含水量均会上升，潜热通量变大，湍流加强，
43、因此草地下垫面的湍流在此时强于旱地作物下垫面。夏季还会产生频繁的对流云，而对流云也是湍流的主要来源。由于中国地处东亚季风区，对流具有明显的季节循环，和季风盛行和西太平洋副热带高压的季节性活动有密切关系。对流对湍流的影响还需要更深入的研究。（）通过选定个类型的下垫面具体分析表明，湍流的强度与对应下垫面类型有很大的联系，在受下垫面影响最大的对流层（和），全年平均分布上草地与旱地作物下垫面的湍流都较强，然后依次是水田作物下垫面、灌木下垫面、荒地下垫面。而在平流层下层（），旱地下垫面的湍流强度相较对流层，降低幅度很大，因此其强度开始明显弱于草地下垫面，而其他下垫面的变化与对流层相差不大。在平流层上层（
44、），下垫面难以影响，因此曲线较为重合，季节变化性也最弱。推测湍流的产生这一现象的原因还与下垫面对应的植被覆盖率、降水量、土壤水分含量都有关系，具体原因需要进一步的研究。致谢：感谢郭建平老师提供的程序方法建议。参考文献曹富强，丹利，马柱国，中国农田下垫面变化对气候影响的模拟研究 气象学报，（）：，（）：（）陈德辉，秦成，一个改进的湍流垂直扩散过程参数化方案 气 象 学 报，（）：，（）：（）龚玺，朱蓉，李泽椿，我国不同下垫面的近地层风切变指数研究气象，（）：，（）：（）郝民，田伟红，龚建东，波段秒级探空资料在 同化系统中的应用研究气象，（）：，（）：（）洪雯，王毅勇，非均匀下垫面大气边界层研究进
45、展 南京信息工程大学学报（自然科学版），（）：，（），（）：（）李伟，赵培涛，郭启云，等，中国 型电子探空仪阳江国际比对结果分析气象，（）：，（）：（）刘辉志，王雷，杜群，大气边界层物理研究进展（年）大气科学，（）：，（），（）：（）刘明星，张宏升，宋星灼，等，不同下垫面温度和湿度湍流谱特征研究 北京大学学报（自然科学版），（）：，气象 第 卷 ，（）：（）潘乃先，陈家宜，栾胜基，不同下垫面表面层的湍流度统计特征 北京大学学报（自然科学版），（）：，（）：（）涂钢，刘辉志，董文杰，半干旱区不同下垫面近地层湍流通量特征分析 大气科学，（）：，（）：（）王少影，张宇，吕世华，等，应用通量方差法估算
46、戈壁绿洲下垫面湍流通量的研究大气科学，（）：，（）：（）杨湘婧，徐祥德，陈宏尧，等，波段探空高分辨率廓线中近地层信息分析及相关模型气象，（）：，（）：（）岳平，张强，牛生杰，等，半干旱草原下垫面能量平衡特征及土壤热通量对能量闭合率的影响气象学报，（）：，（）：（）张璐，黄倩，张宏昇，等，干湿地表的湍流特征及其对深对流影响的大涡模拟气象学报，（）：，（）：（）张蕴帅，黄倩，马耀明，等 黄河源区鄂陵湖湖面和湖边草地对流边界层湍流结构特征的大涡模拟研究大气科学，（）：，（）：（），（）：，：，（）：，（）：，（）：，（）：，（）：，：，（）：，：，（）：，（）：，（）：，（）：，（），（）：，：，（）：，（）：，（）：，：，（）：，：，（）：，（），（）：，（）：，第期 龚玺等：基于高垂直分辨率探空数据对五种典型下垫面湍流分布特征分析 ，（）：，（）：，（）：，（）：，：，（）：，：，（）：，？，（）：，（）：，（）：，：，：，（）：，：，（）：，：，（）：（本文责编：俞卫平）气象 第 卷
