论文部分内容阅读
水热交换是地气间相互作用最为直接和显著的物理过程，一直是气象学、气候学、水文学、生态学等相关学科共同关注的重点和难点研究内容。显热通量（Sensible HeatFlux）是地表湍流运动和能量交换的驱动因子，准确测算区域显热通量显得尤为重要。但迄今为止，在植被分布非均一和地形起伏的复杂下垫面条件下，即便是公里尺度的通量观测技术及数据代表性的研究工作仍有许多问题需要解决。大孔径闪烁仪法（LargeAperture Scintillometers Methods）在公里尺度水热通量的测算研究中，具有广阔的应用前景。但目前大孔径闪烁仪法只能直接观测到大气折射指数的结构参数（the StructureParameter of the Refractive Index of Air，Cn2），测算显热通量需要进一步计算温度结构参数（the Structure Parameter of Temperature，CT2），并根据大气湍流相似理论（Monin-Obukhov Similarity Theory，MOST），建立CT2的MOST普适函数结合，气象数据才能得到最终结果。由于MOST要求地表相对平坦均一，该方法应用于山地森林等复杂或非均一植被显热通量的测算时，应当首先研究实际观测区域湍流运动特征，但至今为止，尚无详尽的研究文献报道。本研究以华北南部低丘山地栓皮栎-侧柏-刺槐混交林为研究对象，在研究实际观测区域湍流运动特征的基础上，判别湍流关键参数是否适合MOST适用条件，再建立适合于实际观测区域的温度结构参数的普适函数，并进行验证，旨在为该方法有效测算山地森林生态系统水热通量提供理论依据，以期促进公里尺度水热通量研究的发展，最终为深入研究了解森林生态系统在能量平衡和调节气候变化中的作用提供工作基础。主要结果如下：（1）湍流强度和风速的相关性很强。当U<3m/s时，湍流强度最大，并随着风速的增大明显减小。一年生长季中两个典型时段（2011年3月15日至4月15日，叶片萌芽生长期；2011年7月15日至8月15日，叶片生长完全期）的平均湍流强度分别为Iu=0.326，Iv=0.306， Iw=0.189和Iu=0.324， Iv=0.296， Iw=0.173。（2）各主要气象要素在山地人工林生态系统复杂下垫面条件下基本满足MOST适用条件。无因次风速方差在不稳定和稳定条件下均符合1/3次幂相似规律，垂直方向比水平方向上的拟合效果好，近中性层结下由σu/u*=A, σv/u*=B, σw/u*=C式计算得到A、B和C的平均值在春季观测期为2.55、2.06和1.30，在夏季观测期为2.61、2.45和1.21，以上观测结果均与类似经典的陆面实验结果接近；同时无量纲湍流动能（Turbulent KineticEnergy,TKE）在不稳定和稳定条件下均符合1/3次幂相似规律；无因次温度方差在春季观测时期不同稳定度条件下均符合2/3次幂相似规律，在夏季观测时期不同稳定度条件下分别符合2/3次和1/3次幂相似规律。湍流谱结构基本符合-5/3规律，实际观测区森林冠层上方大气湍流是非各向同性的。（3）本研究建立的适合本研究区域实际条件的温度结构参数的普适函数准确性较高。研究MOST适用性发现，当热力因子较小时，大气处于非稳定状态，实测值离散程度最大；当热力作用较大时，大气处于稳定状态，离散程度小，说明了近中性条件下机械湍流对LAS观测中的MOST适用性有一定影响，而在对流不稳定时影响较小或可忽略。热力湍流越强，MOST适用性越好，反之则越差。湍流结构的实际情况可作为MOST适用性的判断方法，计算得到的普适函数为，fTJiYuan=4.91（16.51ζ）2/3z/L <0；fTJiYuan=6.43(1+2.20ζ2/3)z/L>0；LAS实际观测与湍流谱方法计算得到的温度结构参数的一致性较好，大部分情况下相关系数都达到了0.9以上。（4）大气层结稳定时更要注意观测信息的空间代表性分析；在风向近似垂直于光径路线时，LAS观测得到的通量信息最全面完整。在大气层结稳定条件下，因湍流活动较弱，观测数据的源区面积相对较大；非稳定层结时，湍流活动旺盛，源区面积较小。大气层结稳定时最小的通量源区面积比不稳定时最大的源区面积大14.7%。（5）湍流理论分析对于提高LAS观测森林生态系统显热通量的准确性具有重要意义。确定了迭代法为该复杂下垫面条件下LAS测算显热通量的计算方法。分析重叠区域印痕权重比值大小以及各风向源区后，建立了适合本观测区域的基于单点通量观测数据扩展到公里观测尺度上的方程为Hec={0.77f1/（0.77f1+0.57f2）}Hec1+{0.572f2/（0.77f1+0.57f2）}Hec2在重叠源区，对由建立的普适函数（FTJiYuan）计算得到的显热通量（计算值）与涡度相关法观测得到的显热通量（实测值）进行相关性分析，结果表明：二者线性相关系数可达0.90。如采用其他5种广泛应用的普适函数计算显热通量，其计算值与实测值的相关系数为0.27-0.63，说明建立的温度结构参数普适函数具有较高的准确性。