• 气候与全球变化 •

### 几种典型地表粗糙度计算方法的比较研究

1. 1. 中国科学院地理科学与资源研究所 生态系统网络观测与模拟重点试验室,北京 100101;
2. 中国科学院研究生院,北京 100049
• 收稿日期:2006-11-13 修回日期:2007-05-15 出版日期:2007-09-25 发布日期:2007-09-25
• 作者简介:周艳莲(1980-),女,湖北人,博士研究生。GIS专业。 *通讯作者 : 孙晓敏,研究员。主要从事地表通量观测与尺度转换研究。E-mail:sunxm@igsnrr.ac.cn
• 基金资助:

国家自然科学基金资助项目(40471099)

### Comparative research on four typical surface roughness length calculation methods

ZHOU Yan-lian1,2, SUN Xiao-min1, ZHU Zhi-lin1, WEN Xue-fa1, TIAN Jing1,2, ZHANG Ren-hua1

1. 1. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research,CAS, Beijing 100101, China;
2. Graduate University of the CAS, Beijing 100049, China
• Received:2006-11-13 Revised:2007-05-15 Online:2007-09-25 Published:2007-09-25
• Supported by:

国家自然科学基金资助项目(40471099)

Abstract: There are four typical surface roughness length calculation methods, which are iterative method by least-square fit, Newton iterative method, Temperature Variance Method (TVM method) and a method proposed by Martano. The former two methods need wind speed and temperature profiles in several levels, while the latter two methods use three-dimensional sonic anemometer data in a single level. In this paper, the surface roughness length in various atmospheric stratifications over forest underlying in Changbai Mountains Experimental Station is calculated by the four methods respectively, with wind speed and temperature profiles and three-dimensional sonic anemometer data in 2003. Discrepancy among the four methods and uncertainty of each method are analyzed. The results indicate that, except Martano method, there is minor difference between surface roughness lengths calculated by the other three methods, while roughness length calculated by Martano method is obviously much smaller. By uncertainty analysis of Martano method, it is found out that smaller roughness length calculated by Martano method is probably related to data amount, and calculation results change obviously with data amount. Otherwise, by making uncertainty analysis for the other three methods, it is found out that as for iterative method by least-square fit and Newton iterative method, calculation roughness length is influenced by friction velocity u*, and if u* changes 1% between two neighborhood measurement heights, the calculated roughness length would change by 8.3%. As for TVM method, parameter of C1 changes in such a certain range from 0.9 to 1.05 that the calculated roughness length changes by about 29.9%. By calculating surface roughness length in various atmospheric stratifications with four methods, it is verified that surface roughness length changes with atmospheric stratifications. Surface roughness length in unstable stratification is the smallest, the one in stable stratification is the highest, and in the neutral stratification is medium. Also in this paper specific values in various stratifications are calculated, The results are 2.642, 2.103 and 1.616m for stable, neutral and unstable stratification roughness length, respectively.