近地面层湍流通量观测方法简介.pptVIP

近地面层湍流通量观测方法简介 康 凌 1. 涡旋相关法 2. 空气动力学方法（梯度法） 3. 能量平衡法（Bowen比法） 4. 误差分析 1. 涡旋相关法 1.1 原理 在近地面层大气湍流场中任何一种被动标量属性均满足： 该式是常通量层的条件，是近地面层湍流通量观测的基础。在常通量层内 积分(1.1)式可得： 其中，H为感热通量；LE为潜热通量；?为动量通量；?a为空气密度； Cp为空气的定压比热；Lv为单位质量水汽的汽化热；u’和w’分别为风速在 水平与垂直方向上的脉动；?’为温度脉动；q’为湿度脉动 。 2. 空气动力学方法（梯度法） 2.1 原理 根据Monin-Obukhov相似理论，在粗糙表面上有： 其中，?M=?H=?q=?为Von-Karman常数，其取值范围在0.35 ~ 0.40，一般为 0.39；h为粗糙表面的平均高度；d为零值位移高度； U、T与q分别为平均 风速、温度和比湿；U*、T*与q*分别为特征风速、温度和比湿； L为Monin-Obukhov长度；?M、?H与?q方便为风速、温度和湿度的无因次梯 度（通量廓线关系函数），且是通用函数。 3. 能量平衡法（Bowen比法） 4. 误差分析 * 1.2 测量方法 用超声风温仪可直接测得风速脉动u’、w’与温度脉动?’，用L??湿度计或红外湿度计直接测得湿度脉动q’，我们就可以利用(1.2)?(1.4)式得到H、LE与?。另外，由于超声风温仪在测量?’ 时的局限（ ?的绝对值不准，但不影响相对值?’ ），我们通常在超声风温仪探头下再安装钨铼丝探头来精确测量?’ ，钨铼丝的特点是的精度高（1/1000?C），响应快（200Hz)。测量仪器见Pic1.1 至1.7。 因为?M、?H与?q是稳定度参数z/L或Richardson数Ri的函数，通过对U、T和q的梯度测量，利用(2.1)至(2.3)式就可得到特征速度U*、特征温度T*与特征比湿q* ，然后由下式求得H、LE和?： 由于只有三个方程，但有四个未知数U*、T*、q*和L。所以最常用的解决办法是用观测资料来迭代求解，方法如下： 根据比较常用的Dyer(1974)提出的通量廓线形式： 其中?=z/L， 当稳定条件下时?0；不稳定条件下时?0；中性条件下时?=0； 对三上式积分得如下廓线公式： 其中T0与q0为地表的温度和比湿，z0为粗糙度，?M、?H和?q的形式如下： 先不考虑稳定度修正，将测量的z3高度的风速U3和z1、z2高度的温度T1 、 T2与比湿q1、q2代入下式计算，就可得到初始的U*、T*、q*： 然后计算Monin-Obukhov 长度L： 代入(2.21)至(2.23)式的右边，就得到了新的U*、T*和q*： 再利用(2.20)式计算出新的M-O长度L后，再代入(2.21)至(2.23)式的右边，更新U*、T*和q*，如此迭代，直到精度满意为止。最后由(2.4)至(2.6)式，就得到了H、LE和?。 2.2 测量方法 风速梯度是利用微气象专用的风速表测量，这种风速表异于气象站使用的普通风速表，其特点是精度高（误差0.1m/s）、启动风速小（0.3m/s）。 温度梯度和湿度梯度的测量则是利用Bowen比系统，其特点是利用它自身的翻转系统，在测量过程中每10分种翻转一次来消除温度探头的误差，从 从而得到高精度的干、湿球温度（误差0.01°C）。测量系统见Pic2.1至2.9。 3.1 原理 根据地表能量平衡方程： 其中Rn为净辐射通量，G为地表土壤热通量，规定Rn与G向下为正，向上为负；?Qs为能量储存项（植物光合作用的化学储存，总为正）；H与LE则分别为感热通量和潜热通量（向上为正，向下为负）。 由于?Qs Rn和G，故可忽略不计，(3.1)式则可写成： Bowen比的定义为： 将(3.3)式代入(3.2)式得： 其中： SI为总辐射；SO为地表反射辐射； LI为大气逆辐射；LO为地表长波辐射；G5为5cm处的土壤热通量。Ts则为不同深度的土壤温度（通常为0、5、10、15、20和40cm，有时达80cm），CW为土壤容积热容量（可推算出来）。 因为SI、SO、 LI、LO、G5和Ts都可实测，所以只要知道了B值，用(3.4)和(3.5)式就可得到H与LE。根据(3.3)与(1.2)及(1.3)式有： 根据湍流输送的K理论： 其中Kh和Kq分别为