地球科学进展 ›› 2010, Vol. 25 ›› Issue (11): 1148 -1165. doi: 10.11867/j.issn.1001-8166.2010.11.1148

观测数据处理与分析 上一篇    下一篇

大孔径闪烁仪观测数据的处理方法研究
白洁 1,刘绍民 1*,丁晓萍 1,卢俐 1,2   
  1. 1.北京师范大学遥感科学国家重点实验室,地理学与遥感科学学院,北京 100875;2.北京市气象局气象信息中心,北京 100089
  • 收稿日期:2010-03-16 修回日期:2010-09-15 出版日期:2010-11-10
  • 通讯作者: 刘绍民(1967-),男,浙江绍兴人,教授,主要从事陆面过程观测与遥感应用研究.  E-mail:smliu@bnu.edu.cn
  • 基金资助:

    公益性行业(气象)科研专项“大尺度水热通量观测系统的研制与应用研究”(编号:GYHY200706046);国家自然科学基金项目“地表水热通量的时空尺度扩展研究”(编号:40971194)和“基于遥感和数据同化方法的海河流域水文通量预测研究—SP2:不同尺度蒸散量和土壤水分的观测研究”(编号:30911130504)资助.

A Study of  the Processing Method of Large Aperture Scintillometer Observation Data

Bai Jie 1, Liu Shaomin 1, Ding Xiaoping 1, Lu Li 1,2   

  1. 1.State Key Laboratory of Remote Sensing Science, School of Geography, Beijing Normal University, Beijing 100875, China;
    2. Meteorological Information Center of Beijng, Beijng Meteorological Bureau, Beijing 100089, China
  • Received:2010-03-16 Revised:2010-09-15 Online:2010-11-10 Published:2010-11-10

大孔径闪烁仪可以获取公里级尺度的显热/潜热通量,其观测数据在农林业、水文、气象等领域的研究中具有越来越重要的作用。以海河流域2008年密云站和馆陶站大孔径闪烁仪的观测数据为例,讨论了数据筛选与处理方法,以及大气不稳定状态下不同处理方法对显热通量造成的影响。结果表明:空气折射指数结构参数应结合电压信号的方差进行计算、采用逐日的日平均波文比系数进行湿度订正、结合空间权重函数计算有效高度、选取Andreas(1988)稳定度函数计算显热通量等是可信的。通过对不同插补方法的比较可知,采用非线性回归方法和以零值替代的方法对不稳定状态和稳定状态下30 min缺失的显热通量进行插补,采用动态线性回归方法对日显热通量的缺失进行插补是可行的。为了解决由于冬天净辐射和显热通量较小,造成能量平衡方程余项法计算的蒸散量产生较大误差的问题,可采用根据大孔径闪烁仪与涡动相关仪观测的日蒸散量建立的关系式进行估算。基于上述方法的研究,形成了一套较为完整的大孔径闪烁仪观测数据的处理流程,保证了在不同下垫面、不同的天气状况条件下都能获取质量可靠、数据连续的大孔径闪烁仪观测的显热/潜热通量。

Sensible/latent heat fluxes can be obtained by Large Aperture Scintillometer (LAS) over several kilometers, which play an significant role in the analysis  of and application to agricultural and forestry, hydrology and meteorology research. Take LAS observations in Miyun and Guantao stations over Hai River Basin in 2008 as an example. The way of data screening and quality control under unstable conditions as well as the effect of different calculation method on sensible heat flux have been discussed. The results showed: the structure parameter of the refractive index (C2n) should be calculated with the variance of the voltage of structure parameter of the refractive index; Humidity correction can be done with daily Bowen ratio; The effective height of LAS can be calculated with spatial averaged function;  The Andreas(1988) function is  used to get reliable sensible heat flux. Meanwhile, the nonlinear regression method and dynamic linear  regression method have been used to fill the 30min and daily missing data  while the 30  min missing data under stable conditions can be set to zero. The correlation relationship between daily ET measured by LAS and EC has been constructed to estimate daily ET when daily Rn<50 W/m2. Based on the above analysis, a set of LAS data processing scheme has been set up, which ensures continuous and high quality sensible/latent data can be obtained over various surfaces and weather conditions.

中图分类号: 

[1] Chebouni A, Watts C, Lagouarde J P. Estimation of heat and momentum fluxes over complex terrain using a large aperture scintillometer[J].Agricultural and Forest Meteorology,2000, 105: 215-226.
[2] Hoedjes J C B, Zuurbier R M, Watt C J. Large aperture scintillometer used over a homogeneous irrigated area, partly affected by regional advection[J].Boundary-Layer Meteorology, 2002, 105: 99-117.
[3] Hartogensis O K, De Bruin H A R, Van De Wiel B J H. Displaced-beam small aperture scintillometer test. Part II: CASES-99 stable boundary-layer experiment[J].Boundary-Layer Meteorology,2002, 105: 149-176.
[4] Beyrich F, Richter S H, Weisensee U, et al.The LITFASS-98 experiment: Fluxes over a heterogeneous land surface′[C]//Proceedings of the 14th Symposium on Boundary Layer and Turbulence, Aspen, CO, American Meteorological Society, 45 Beacon St., Boston, MA,2000:9-10.
[5] Special Issue of “Bourdary-Layer Meteorology”[M].Springer,2006,121:1-220.
[6] Beyrich F, De Bruin H A R, Meijninger W M L, et al. Results from one year continuous operation of a large aperture scintillometer over a heterogeneous land surface[J].Boundary Layer Meteorology, 2002,105: 85-97.
[7] Meijninger W M L, Beyrich F, Lüdi A, et al.Scintillometer based turbulent fluxes of sensible and latent heat over a heterogeneous land surface—A contribution to LITFASS-2003[J].Boundary-Layer Meteorology, 2006,121: 89-110.
[8] Kleissl J, Gomez J, Hong S H, et al.  Large aperture scintillometer intercomparison study[J].Boundary-Layer Meteorology,2008,128: 133-150.
[9] Savage M J.Estimation of evaporation using a dual-deam surface layer scintillometer and component energy balance measurements[J].Agricultural and Forest Meteorology,2009, 149: 501-517.
[10] Solignac P A, Brut A, Selves J L, et al.Uncertainty analysis of computational methods for deriving sensible heat flux values from scintillometer measurements[J].Atmospheric Measurement Techniques Discussions,2009, 2:1 383-1 417.
[11] Ochs G R, Wilson J J.A second-generation large aperture scintillometer[Z]. NOAA Tech. Memor. ERL ETL-232, NOAA Environmental Research Laboratories, Boulder, CO, USA,1993:24.
[12] Ochs G R, Hill R J.Optical-scintillation method of measuring turbulence inner scale[J].Applied Optics,1985,24:2 430-2 432.[13] Beyrich F. Operational aspects of scintillometry[Z]. 2nd Scintillometer Workshop, Wageningen, 2007.
[14] Ezzahar J, Chehbouni A, Hoedjes J C B, et al. The use of the scintillation technique for monitoring seasonal water consumption of olive orchards in a semi-arid region[J].Agricultural Water Management,2007,89: 173-184.
[15] Von Randow C, Kruijt B, Holtslag A A M, et al. Exploring eddy-covariance and large-aperture scintillometer measurements in an Amazonian rain forest[J].Agricultural and Forest Meteorology,2008,148: 680-690.
[16] Ezzahar J, Chehbouni A, Hoedjes J, et al. Combining scintillometer measurements and an aggregation scheme to estimate area-averaged latent heat flux during the AMMA experiment[J].Journal of Hydrology,  2009,375(1/2):217-226.
[17] Scintec boundary layer scintillometer[Z]. User manual. Scintec, 2007:70.
[18] De Bruin H A R, Kohsiek W, Van den Hurk B J J M. A verification of some methods to determine the fluxes of momentum, sensible heat and water vapor using standard deviation and structure parameter of scalar meteorological quantities[J].BoundaryLayer Meteorology,1993,63: 231-257.
[19] Hoedjes J C B, Chehbouni A, Ezzahar J, et al. Comparison of large aperture scintillometer and eddy covariance measurements: Can thermal infrared data be used to capture footprint-induced differences?[J].Journal of Hydrometeorology,2007,8: 144-159.
[20] Hartogensis O K, Watts C J, Rodriguez J C, et al.Derivation of an effective height for scintillometers: La Poza experiment in northwest Mexico?[J].Journal of Hydrometeorology,2003, 4:915-928.
[21] Marx A, Kunstmann H, Schuttemeyer D,et al. Uncertainty analysis for satellite derived sensible heat fluxes and scintillometer measurements over savannah environment and comparison to mesoscale meteorological simulation results[J].Agricultural and Forest Meteorology,2008,148: 656-667.
[22] Lu Li, Liu Shaomin, Xu Ziwei, et al. Results from measurements of large aperture scintillometer over different surfaces[J].Journal of Applied Meteorological Science,2009,20(2): 171-178.[卢俐, 刘绍民, 徐自为,等. 不同下垫面大孔径闪烁仪观测数据处理与分析[J]. 应用气象学报,2009, 20(2): 171-178.]
[23] Kipp&Zonen. Large aperture scintillometer instruction manual[Z]. Kipp&Zonen, 2007:74.
[24] Green A E, Hayashi Y.Use of the scintillometer technique over a rice paddy[J].Journal of Agricultural Meteorology,1998, 54(3):225-234.
[25] Hemakumara H M, Chandrapala L, Moene A F. Evapotranspiration fluxes over mixed vegetation areas measured from large aperture scintillometer[J].Agricultural Water Management,2003,58: 109-122.
[26] Meijninger W M L, De Bruin H A R. The sensible heat fluxes over irrigated areas in western Turkey determined with a large aperture scintillometer[J].Journal of Hydrology,2000,229: 42-49.
[27] Moene A F. Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation[J].Boundary-Layer Meteorology,2003,107: 635-653.
[28] Asanuma J, Iemoto K. Measurements of regional sensible heat flux over Mongolian grassland using large aperture scintillometer[J].Journal of Hydrology,2007,333: 58-67.
[29] Odhiambo G O, Savage M J. Sensible heat flux by surface layer scintillometry and eddy covariance over a mixed grassland community as affected by bowen ratio and MOST formulations for unstable conditions[J].Journal of Hydrometorology,2009,10:479-492.
[30] Su Z, Timmermans W, Gieske A, et al.'Quantification of landatmosphere exchanges of water, energy and carbon dioxide in space and time over the heterogeneous Barrax site'[J].International Journal of Remote Sensing,2008,29(17):5 215-5 235.
[31] Lu Li. A Study on the Observation, Influencing Factors and Scale Relationships of the Sensible Heat Flux[D].Beijing: Beijing Normal University,2008:163.[卢俐.地表感热通量的观测、影响因子和尺度关系的研究[D].北京:北京师范大学,2008:163.]
[32] Moene A F, De Bruin H A R. Sensible heat flux data derived from the scintillometers[C]//Su Z, Jacobs C,eds. BCRS Report, Advanced Earth Observation Land Surface Climate. Final Report, 01-02,2001:85-90.
[33] De Bruin H A R, Van Den Hurk B J J M, Koshiek W. The scintillation method tested over a dry vineyard area[J].Boundary-Layer Meteorology,1995,76: 25-40.
[34] Xu Ziwei, Liu Shaomin, Xu Tongren, et al. Comparison of the gap filling methods of evapotranspiration measured by eddy covariance system[J].Advances in Earth Science,2009, 24(4):372-382.[徐自为, 刘绍民, 徐同仁,等.涡动相关仪观测蒸散量的插补方法比较[J].地球科学进展,2009, 24(4): 372-382.]
[35] Falge E, Baldocchi D D, Olson R, et al. Gap filling strategies for long term energy flux data sets[J].Agricultural and Forest Meteorology, 2001,107: 71-77.
[36] Alavi N, Warland J S, Berg A A. Filling gaps in evapotranspiration measurements for water budget studies: Evaluation of a Kalman filtering approach[J].Agricultural and Forest Meteorology,2006, 141:57-66.
[37] Ren Lihong, Wang Chenghai, Qiu Chongjian, et al. A study of computing the surface flux in the typical arid region of northwest China by a variational method[J].Chinese Journal of Atmospheric Sciences,2004, 28(2):269-277.[任立宏, 王澄海, 邱崇践,等.利用变分法计算西北典型干旱区地表通量的研究[J].大气科学,2004,28(2):269-277.]
[38] Yang Juan, Zhou Guangsheng, Wang Yunlong, et al. Estimation of sensible and latent heat fluxes of typical steppe in inner Mongolia based on variational method[J].Chinese Journal of Applied Ecology,2006,17(11): 2 046-2 051.[杨娟, 周广胜, 王云龙,等.基于变分方法的内蒙古典型草原水热通量估算[J].应用生态学报,2006,17(11): 2 046-2 051.]
[39] Schmidt A, Wrzesinsky T, Klemm O. Gap filling and quality assessment of CO2 and water vapour fluxes above an urban area with radial basis function neural networks[J].Boundary-Layer Meteorology,2008, 126:389-413.
[40] Dou Zhaoyi, Liu Jianjun. Application of artificial neural networks to interpolation and extrapolation of flux data[J].Journal of Northwest Forestry University,2009, 24(3): 58-62.[窦兆一,刘建军.人工神经网络在通量观测资料插补中的应用[J].西北林学院学报,2009, 24(3): 58-62.]
[41] Verhoef W. Application of Harmonic Analysis of NDVI Time Series (HANTS)[C]//Azzali,  Menenti, eds. Fourier analysis of temporal NDVI in the Southern African and American continents. Report of DLO Winand Staring Centre, Wageningen (The Netherlands), 1996.
[42] Wang T, Ochs G R, Clifford S F. A saturation-resistant optical scintillometer to measure C2n [J].Journal of Optical Society of America,1978,68: 334-338.
[43] Hill R J, Clifford S F, Lawrence R S. Refractive-index and absorption fluctuations in the infrared caused by temperature, humidity, and pressure fluctuations[J].Journal of Optical Society of America,1980, 70(10): 1 192-1 205.
[44] Wesely M L. The combined effect of temperature and humidity fluctuations on refractive index[J].Journal of Applied Meteorology,1976, 15:43-49.
[45] Andreas E L. Two  wavelength method of measuring path-averaged turbulent surface heat fluxes[J].Journal of Atmospheric & Oceanic Technology,1989,6: 280-292.
[46] Andreas E L. Estimating C2n over snow and sea ice from meteorological data[J].Journal of the Optical Society of America,1988,5: 481-495.
[47] De Bruin H A R, Kohsiek W, Van den Hurk B J J M. A verification of some methods to determine the fluxes of momentum, sensible heat and water vapor using standard deviation and structure parameter of scalar meteorological quantities[J].Boundary-Layer Meteorology, 1993,63: 231-257.
[48] Thiermann V, Grassl H. The measurement of turbulent surface layer fluxes by use of bichromatic scintillation[J].Boundary-Layer Meteorology,1992, 58:367-389.
[49] Wyngaard J C, Izumi Y, Collins Jr S A. Behaviour of the refractive index structure parameter near the ground[J].Journal of the Optical Society of America,1971, 61: 1 646-1 650.
[50] Wyngaard J C. On surface-layer turbulence[C]//Proceedings of the Workshop on Micrometeorology. American Meteorological Society, Denver, CO,1973:101-149.
[51] Yang K,Wang J M.A temperature prediction-correction method for estimating surface soil heat flux from soil temperature and moisture data[J].Science in China (series D),2008,51:721-729.

[1] 罗维均, 王世杰, 刘秀明. 喀斯特洞穴系统碳循环的烟囱效应研究现状及展望 *[J]. 地球科学进展, 2014, 29(12): 1333-1340.
[2] 徐自为,黄勇彬,刘绍民. 大孔径闪烁仪观测方法的研究[J]. 地球科学进展, 2010, 25(11): 1139-1147.
[3] 蔡旭晖,朱明佳,刘绍民,徐自为. 大孔径闪烁仪的通量印痕分析与应用[J]. 地球科学进展, 2010, 25(11): 1166-1174.
[4] 郑宁,张劲松,孟平,黄辉,高峻,贾长荣,任迎丰. 基于闪烁仪观测低丘山地人工混交林通量印痕与源区分布[J]. 地球科学进展, 2010, 25(11): 1175-1186.
[5] 朱治林, 孙晓敏, 贾媛媛, 温学发, 张仁华, 袁国富, 唐新斋. 基于大孔径闪烁仪(LAS)测定农田显热通量的不确定性分析[J]. 地球科学进展, 2010, 25(11): 1199-1207.
[6] 王维真,徐自为,李新,王介民,张智慧. 大孔径闪烁仪在黑河流域的应用分析研究[J]. 地球科学进展, 2010, 25(11): 1208-1216.
[7] 艾力·买买提明,何清,霍文,刘新春. 塔克拉玛干沙漠腹地LAS和EC观测感热通量对比分析[J]. 地球科学进展, 2010, 25(11): 1228-1236.
[8] 李远, 孙睿, 刘绍民, 徐自为, 白洁. 大孔径闪烁仪观测数据在陆面模式验证中的应用初探[J]. 地球科学进展, 2010, 25(11): 1237-1247.
[9] 刘雅妮,辛晓洲,柳钦火,周春艳. 基于多尺度遥感数据估算地表通量的方法及其验证分析[J]. 地球科学进展, 2010, 25(11): 1261-1272.
[10] 白洁,刘绍民,丁晓萍. 海河流域不同下垫面上大孔径闪烁仪观测显热通量的时空特征分析[J]. 地球科学进展, 2010, 25(11): 1187-1198.
[11] 卢俐,刘绍民,徐自为,白洁,王介民. 大孔径闪烁仪和涡动相关仪观测显热通量之间的尺度关系[J]. 地球科学进展, 2010, 25(11): 1273-1282.
[12] 张劲松,孟平,郑宁,黄辉,高峻. 大孔径闪烁仪法测算低丘山地人工混交林显热通量的可行性分析[J]. 地球科学进展, 2010, 25(11): 1283-1290.
[13] 涂钢,支克广,支平,付光极,王建平. 大气折射指数的结构参数与北方5~9月降雨的强相关现象[J]. 地球科学进展, 2010, 25(11): 1291-1295.
[14] 双喜,刘绍民,徐自为,王维真. 黑河流域观测通量的空间代表性研究[J]. 地球科学进展, 2009, 24(7): 724-733.
[15] 徐自为,刘绍民,徐同仁,王介民. 涡动相关涡动相关仪观测蒸散量的插补方法比较[J]. 地球科学进展, 2009, 24(4): 372-382.
阅读次数
全文


摘要