０ 引言
陆地表面动量输送不仅对全球气候和大气环流变化贡献显著，也是各类尺度大气数值模式物理框架中极其关键的环节。总体输送系数不仅有明显的地区差异，而且动态范围也很大。由于方法差异和观测技术问题，有时不同的人在同一地区得到或使用的总体输送系数也会有数倍之差，甚至量级之别。
总体输送系数的误差是以等同的相对值传递到计算的地表通量上，而大气数值模式所能容忍的地表面通量误差程度极其有限，所以，较准确的总体输送系数是目前模式所必需的。
组合法是由风、温、湿梯度资料、地表辐射分量资料和地热流量资料联合计算地表湍流通量的间接方法，它首先由Thom提出，然而它只能确定感热和潜热通量，无法确定MO通用函数。而胡隐桥等人改进的组合法，不仅能确定感热和潜热通量，同时，还可以确定动量和MO通用函数。这一方法的优点是价格低廉，并且减小了感应头及其支架对气流扰动所造成的湍流通量观测误差。本文就是根据胡隐樵[3]等人改进的组合法，结合江南肥西站的具体资料进行详细的计算。最后，取得了较好的效果。
1 观测站周围环境、观测仪器
肥西边界层观测设在肥西县气象站局（117°.08′E，31°.41′S）观测场内，下垫面为平坦均匀稀疏矮草地，周围较为开阔，西、南面基本为旱地农田，北面为气象局两层办公楼，东面为平房，建筑物距观测场50米外。
超声风温仪、红外水汽和二氧化碳探头均安装在4米高度，观测系统为连续自动观测，由美国生产的CR23X数据采集器采样频率为10Hz，包括U、V、W风速分量，温度、水汽、二氧化碳、气压等，每10分钟生成一个数据文件，经初步处理，将每小时资料合并为一个文件。
肥西站在铁塔的1米、4米、10米高度各伸出2米长的水平臂，两端分别安装温、湿度传感器和风速传感器，此外还安装有四层（0cm,10cm,20cm,40cm深）地温表，装有两个土壤热流量板。
2 计算方法
按照总体输送理论，地面、感热和蒸发潜热的通量可以表示为:
D=ρCDU2 （1）
H=ρCPCHU（Ts0-T）（2）
Q=βρLECEU（qgs-q）（3）
式中ρ是地面空气密度，可表示为海拔高度的降指数。Tso是地面温度（0cm土壤温度），qgs地面饱和比湿。为使计算结果有较好的代表性，地面风速和气温都取为自动站最高观测高度（10m）的值。β是土壤有效因子或称为蒸发系数，它是土壤湿度的升函数。本文由于观测资料有限，不便于计算β的值，这里取β=0.5。LE为水汽的蒸发潜热，取2.500×105为使计算结果有较好的代表性，地面风速和气温都取气象站最高观测高度（10m）值。
又由胡隐樵等人改进的组合法[1]得出的计算通量的表达式可得出动量、感热、潜热的计算式：
D=Ｄ０（Ｆ／α）（4）
Ｈ＝Ｈ０Ｆ（5）
Ｑ＝Ｑ０Ｆ（6）
Ｈ■＝－ρＣ■κ■z■z■■■■■（7）
Ｄ０=ρκ■z■z■（■）2（8）
Ｑ０=ρＣ■κ■z■z■■■■■（9）
F是层结影响函数，α是动量和热量通量层结影响函数的比例，
■≤0时，α=F1/3（10）
■>0时，α=1.0（11）
F=■（12）
求F时，需用到地表面上的能量平衡公式：
A=R-G-P-S （13）[1]
方程中，A为地表面可利用能量，而R，G，P和S分别为净辐射、地流热量、植物光合作用消耗的能量和地表层贮藏的能量。在实际计算过程中，地表面可利用能量一般忽略植物光合作用消耗的能量P，定常条件下也可忽略地表层贮藏的能量S，故
A≈R-G （14）
3 计算及结果分析
实际参加计算的资料是1，4，10三层的风速、温度、湿度、水汽压、地表温度、降水量、气压等7组资料，每十分钟取平均记录，全天144个记录.经过资料的整理和筛选，共有45天的记录可用，从6月6日-7月20日，观测期间，观测站降水很少，出现了大旱的天气，6月9日、10、11、12、13、17、18、22日，7月1、11、13、14、15日有少量降水。其中以7月11日最多。
3.1 整体变化
江南地区，大气状况较好，压强均在1000hp左右，常年有降水，下垫面多为湿润的土壤。以下的图片可表现肥西地区的地表湍流通量变化规律。
图1 6月6日至7月20日的感热随地面温度变化图
从图1可看出，感热的值一般在20-50W/m2之间。除了个别突变的点之外，能发现感热变化的一般规律：地面温度高，则感热高。但是，利用组合法得出的感热值和温度值有突变值。比如7月上旬、中旬的那几天，虽有零星降水，但地面温度还较高，而感热却为负值，不太符合实际，应剔除。
从图2中可看出，潜热变化的大致规律：地面含水率高的，其潜热也大,潜热值的大小在20-220W/m2之间变化。同时也可看出，潜热值的日变化比较明显，这与有无降水有一定的关系，如图中所示，7月11日达到整个观测期的峰值，这与当天有一定量的降水实际情况是相符的。
图2 6月6日至7月20日潜热和
地面0-15厘米含水率的关系图
图3 风速和动量的变化图
从图3可知，对于动量的变化，除了几处（如6月6日至6月9日之间几天，6月21日至6月24日，7月9日前后几天）符合之外，其他各处符合得不是很好，不稳定的点较多。从理论上分析，因为动量和风速的平方成正比的，因此，风速脉动量对动量的影响是比较明显的。从此，也可知，组合法得出的动量数据与风速对应稳定性不太好。当然，还是可从图3中得出：肥西地区6、7月份的动量值在3.8-25.0J之间变化。
3.2 个例分析
这里选取6月18日感热、6月18日潜热的和6月28日的动量变化图说明组合法确定的地面湍流通量的日变化特征以及存在的问题。
注：这里使用气象时间，即：实际上是从6月17日晚20：00至6月18日晚20点为6月18日一天。
图4 6月18日感热随时间变化图
整体观察图4，知基本上符合实际的变化规律。通过分析上图，该地区的感热最大值约为18左右，出现在上午10点至中午12点之间，这比沙漠出现最大值时刻（8点-9点）稍微滞后了一些。最直接的因素当然是下垫面的因素。当然，通过观察可发现9点至10点左右的值是不合理的。从图4中也可得到，感热通量白天为正值，夜晚为负值，大概在子夜的时候达到最小，由正值转变为负值一般在傍晚的6、7点左右
图5 6月28日潜热随时间的变化图
图6 0-15厘米土壤的含水率变化图
潜热通量的日变化的位相特征与感热是相同的。通过对比图5、图6两图可发现，含水量高的时候，水汽蒸发潜热也高，这与理论上是一致的。而且，还可发现，肥西站的蒸发潜热在白天上午10点左右达到最大值，大约是270W/m2。在夜间21点左右达到最小为-20 W/m2。潜热通量的值在白天一般为正值，晚上一般为负值。剔除图中的不合理值，可以看出，从正值变为负值的时刻在晚上的7、8点。
3.3 因子的相关性分析
通过相关性分析表明，对于潜热Q，考察与地面湿度（0-15厘米土壤含水率h）时，其相关系数为ρQh=0.678, 可知相关性较好。
对于感热H，考察感热H与地面温度TS0的相关性时，发现H和TS0的相关系数是ρHTs0=0.55，虽然数值有些小，但这与理论上是基本相符的，即：感热代表的是地面想大气的加热量。从直观上讲，地面温度越高，地面向大气输送的热量的就多。
对于动量D，考察动量D与风速V的相关性，为了使计算有代表性，尽量剔除建筑物的影响，V的速度选取10米高度的风速。通过计算，动量D和风速V的相关系数仅为0.412，达不到基本的概率要求，相关性较差。
3.4 组合法确定的各通量平均值与青藏高原西部各值的比较
这里引用李国平等人、张强等人分别对青藏高原西部狮泉河地区计算结果[2]和黑河流域计算结果[3]，并与肥西地区计算结果进行对比，见下表1。
表1 肥西地区、狮泉河地区的通量平均值对比
从表1可知，肥西地区的感热平均值H（28W/m2）比狮泉河的感热值（170.2W/m2）要小得多，下垫面最重要的因素：狮泉河地处青藏高原西部，多为沙漠地带，白天太阳对地面加热，地面由于其自身的沙质特征，比热容比较小，升温特别快，接着便是地面对空气的输送能量，而且全年降水稀少，全年的情况都没有大的变化，故其感热比较大，最大值可达到380W/m2；而肥西地区的地面特征较湿润，升温较慢，故对空气的加热量也不多。在加上有雨水的调节作用，使得该地区的感热值比较小了。
而对于潜热值Q，肥西地区的蒸发潜热平均值130W/m2，而狮泉河的平均值仅为2.9W/m2，这是因为前者无论是地面上层土壤（0-15cm）还是空气中的含水率都比青藏高原西部的狮泉河沙漠地区的含水率要高得多，无论前者的降水量多与少，前者的水汽蒸发潜热自然要比后者高得多。
至于动量值D，肥西地区的动量值（10N/m2）远比狮泉河地区的（180N/m2）小，从理论上说，这种差别主要取决于风速的差别。但是，事实上江南地区的风速和青藏高原西部的风速差异不会相差那么大的，所以在这里，也有可能是组合法的一个缺陷之一：对于自变量的脉动量依赖性大的函数，组合法很好地保证函数的稳定性。
3.5 总体输送系数的平均值
表2 组合法确定肥西的总体输送系数平均值与
狮泉河和黑河流域值对比表
从表2中可以看出，用组合法计算的总体输送系数CH、CE数量级是10-3，而且感热输送系数CH和潜热输送系数CE非常接近。但是，动量输送系数CD量级较小。对狮泉河地区的计算，李国平等人的算法假设CH=CE（假定水汽廓线与温度廓线函数积分形式相同）[2]。张强等人确定的总体输送系数也是由组合法计算出来的[3]。这三者总体输送系数CH、CE的计算结果没有量级上的差别，只是动量输送系数CD数值相对小。
3.6 问题和讨论
3.6.1 在计算三层各个气象要素的偏微分时，采用何种方法更为稳定，有两种方案，一是：利用牛顿差分公式，得出步长不一致的三点公式，二是：两层两层相减，得出平均值。但用第一种方案时，如何处理三点公式与公式（6）、（7）中Zi和Zi+1不一致的矛盾；用第二种方案时，又会出现奇点，稳定性不能保证。
3.6.2 建筑物对通量的影响是不能忽略的，由于观测场地周围建筑物存在，很可能使通量及通量系数发生不正常的的变化。剔除建筑物的影响会使这些量趋于合理化。
3.6.3 下垫面的气象特征是影响通量的最重要因子，感热、潜热和动量的变化特征在不同的下垫面是截然不同的（沙漠、海洋、戈壁、平原等）。由于，肥西站2001年6、7月份降水稀少，不能代表江南地区的典型气候特征，故计算值有些缺乏一定的代表性。
4 结论
4.1 本文采用组合法结合肥西站的实际观测资料，计算出了通量总体输送系数（CE、CH、CD）、感热、潜热、动量等量。总体而言，组合法确定感热、潜热、动量以及通量总体输送系数的方法是基本可行的，而且相比于其他方法有优点是价格低廉，并且减小了感应头及其支架对气流扰动所造成的湍流通量观测误差，数值的总体误差在量级上没太大的突变。发现总体输送系数的量级是10-3，只是CD偏小，而且，CE和CH很接近。
4.2 通过对感热、潜热的平均值的计算发现，用组合法计算感热和潜热是比较符合的，虽会出现一些异常的数值，但突变的情况比较少。由于观测时期，干旱少雨，雨该地区的平均情况有所差异，所以，感热、潜热和动量的计算值缺乏一定的代表性，特别是计算动量的稳定性较差，数值的突变较多。
4.3 肥西地区的感热通量的最大值约为50W/m2，狮泉河地区的感热最大值可达380W/m2相当于前者的7倍。出现最大值的时间前者（10-12点）也比后者（8-9点）滞后。
4.4 肥西地区的潜热通量的最大值约为280W/m2，狮泉河地区的潜热最大值为15W/m2，前后差别较大，前者是后者的16倍多。根据以上结果推算，6、7月份肥西地区感热的平均值应该是为正值潜热通量的日变化的位相特征与感热是相同的。
4.5 肥西地区感热通量的日变化特征是：白天大，夜间小，上午最大，凌晨（4-6点）左右最小；白天一般是正值，夜间一般是零值或负值。潜热通量也大体一致。
4.6 对于动量通量D，组合法得出的值与风速对应不是很好，这也提出了进一步改进组合法的必要性。
【参考文献】
［1］胡隐樵，等.组合法确定近地面层的通用函数的一般表达式[J].气象学报，1991，49（1）：46-53.
［2］李国平，段延扬，巩远发.青藏高原西部地区总体输送系数和地面通量[J].科学通报，2000，45（8）：865-869.
［3］张强，卫国安，黄荣辉.西北干旱区荒漠戈壁动量和感热总体输送系数[J].科学通报，2001，31（9）：783-791.
［责任编辑：周娜］
相关热词搜索： 通量 组合 长江 下游 系数
