本文是对近20年新发展的双波段光闪烁方法的一个综述。陆面过程模式常常是基于局地或斑块尺度上的观测发展的， 其与大气模式较大网格尺度的不匹配， 显然会对后者的效能带来很大影响。如何扩展野外实验站点的代表性尺度已成为当前更好地了解陆面过程， 进而改善陆面过程模式与大气模式网格尺度匹配问题的关键。可用的面积平均通量观测方法， 包括以涡动相关方法为主的多点微气象观测、 飞机观测、 卫星和地面遥感等5种。其中， “光闪烁方法”是当前最为可行的、 可以大到10 km尺度的感热通量和潜热通量观测方法， 特别是， 它可以应用于复杂下垫面包括山谷地区和城市等。光闪烁方法的理论涉及电磁波传输和大气湍流。文章从折射指数、 结构参数、 湍流谱等基本概念开始， 对由对数光强方差计算折射指数的结构参数等基本公式的推导， 由光程权重函数、 空间谱权重函数、 光强的时间序列谱等对闪烁仪主要工作尺度的了解， 由折射指数结构参数计算温度、 湿度的结构参数的方法， 以及利用近地层相似理论计算感热通量和潜热通量等光闪烁方法的理论、 公式和计算步骤等做了较系统的阐述。进而， 在介绍双波段闪烁仪通量足迹函数之后， 对光闪烁方法与涡动相关方法从特征、 优势与缺点三方面做了比较； 指出结合使用涡动相关和光闪烁两种方法， 可以更好地进行面积平均通量分析， 进而用于模式的发展和检验， 以及更好的流域尺度的能量和水循环研究。最后， 对光闪烁方法的应用， 从较均匀下垫面、 复杂下垫面、 城市地区、 遥感模式的地面“真值”及在大气模式中的应用等几方面做了介绍； 特别是， 结合黑河流域阿柔和大满两站2020年部分资料的分析， 彰显了双波段闪烁仪对较大尺度感热、 潜热通量观测的明显优势。但是， 有关方法特别是观测水汽通量的微波闪烁仪的研制究竟为时尚短； 相关硬件、 软件、 资料处理方法等， 许多地方都还需要研究改进。相对于涡动相关方法， 光闪烁方法的理论和数据处理都更为复杂； 有关台站和资料分析人员， 需要更好的物理和微气象学基础。国内闪烁仪在城市地区的应用， 至今仍开展较少； 青藏高原和内地一些湖面蒸发研究中的难点， 也有望在应用光闪烁方法的过程中有所突破。基于光闪烁法等多种面积平均通量观测， 对陆面过程研究做尺度扩展， 并借以推动中—大尺度大气模式的发展， 更是我们殷切期盼的。
This is a review of the Optical-Microwave Scintillometer （OMS） system newly developed in last two decades， which can measure area averaged sensible and latent heat fluxes over a scale of 1 -10 km， especially over heterogeneous surfaces such as cross a valley or over urban areas.Among the methods of area averaged flux measurements， such as eddy-covariance based multi-point observations， air craft observations， satellite， and surface remote sensing etc.， scintillometry is probably the most feasible technique in getting areal fluxes up to 10 kilometers.The basic theory of scintillometry includes electromagnetic wave propagation， atmosphere turbulence， and micrometeorology， which are more sophisticated than that of the popular eddy covariance （EC） system.Based on the introduction of concepts such as refractive index， structure parameter and turbulence spectra etc.， basic scintillometer theories and equations are presented briefly， including： （1） calculation of structure parameter of refractive index via the variance of received light log-intensities； （2） the understanding of the working scale of scintillometry via the light-path weighting function， the spatial spectral weighting function， and the temporal spectral characteristics； （3） the derivation of structure parameters of temperature and humidity via the structure parameters of refractivity； （4） the calculation of fluxes by using the typical functions of Monin-Obukhov similarity； （5） the footprint analysis of scintillometry.A comparison between scintillometry and EC are presented in three aspects： ‘Characteristics’， ‘Advantages’ and ‘Weakness’.It is clear that a combined use of EC & Scintillometry can provide better area averaged fluxes， and， refined flux aggregation schemes.Then， the application of scintillometry is introduced for rather homogeneous surfaces， complicated surfaces， urban areas， and， the ‘ground truth’ for remote sensing and the application of areal averaged fluxes in atmospheric models.The example utilizations of the OMS systems in the Arou alpine-meadow station and the Zhangye oasis station， of the Heihe River basin， clearly show the advantages of scintillometry over the EC in the measurements of larger scale evapotranspiration.Nowadays there are hundreds flux stations operating over various climate regions and surface states of the world.Upscaling of the spatial representativeness of these stations becomes a key point to better understanding the land surface processes， and improving the spatial matching between land surface models and meso/large scale atmospheric models.Comparatively， the time of development of scintillometry， particularly the microwave scintillometers in measuring water vapor fluxes， are still short.Further improvements of relevant hardware， data sampling， and data processing software etc.are still needed.
[1]AndreasE L， 1988.Atmospheric stability from scintillation measurements［J］.Applied Optics， 27： 2241-2246.
[2]AndreasE L， 1990.Selected papers on turbulence in a refractive medium［C］//SPIE Milestones Series， 25， 693pp.SPIE Optical Engineering Press， Bellingham， Wash.， USA.
[3]AvissarR， PielkeR A， 1989.A parameterization of heterogeneous land surface for atmospheric numerical models and its impact on regional meteorology［J］， Monthly Weather Review， 117： 2113-2136.
[4]BastiaanssenW， MenentiM， FeddesR， alet， 1998.A remote sensing surface energy balance algorithm for land （SEBAL）.1.Formulation［J］.Journal of Hydrology， 212： 198-212.
[5]BeyrichF， MengelkampH T， 2006.Evaporation over a heterogeneous land surface： EVA_GRIPS and the LITFASS-2003 experiment-an overview［J］.Boundary-Layer Meteorology， 121： 5-32.
[6]BeyrichF， BangeJ， HatogensisO， alet， 2012.Towards a validation of scintillometer measurements： The LITFASS-2009 experiment［J］.Boundary-Layer Meteorology， 144： 83-112.
[7]Beyrich F， Dereszynski P， van Kesteren B， 2013.Some aspects and results of scintillometer long-term operation.4th Workshop on Scintillometers and Applications［Z］.Tübingen， Germany， 7-9 October， 2013.
[8]BooneA， de RosnayP， BasalmoG， alet， 2009.The AMMA Land Surface Model Intercomparison Project［J］.Bulletin of the American Meteorological Society， 90（12）： 1865-1880.
[9]CliffordS F， 1971.Temporal-frequency spectra for a spherical wave propagating through atmospheric turbulence［J］.Journal of the Optical Society of America， 61（10）： 1285-1292.
[10]de BruinH， van den HurkB， KohsiekW， 1995.The Scintillation Method tested over a vineyard area［J］.Boundary-Layer Meteorology， 76： 25-40.
[11]de BruinH， WangJ M， 2017.Scintillometry： A review［M/OL］.［2021-02-15］.https： //www.researchgate.net/project/ Scintillometry.
[12]de VreseP， SchulzJ， HagemannS， 2016.On the representation of heterogeneity in land-surface-atmosphere coupling［J］.Boundary-Layer Meteorology， 160： 157-183.
[13]DickinsonR E， OlesonK W， BonanG， alet， 2006.The Community Land Model and its climate statistics as a component of the Community Climate System Model［J］， Journal of Climate， 19： 2302-2324.DOI： 10.1175/jcli3742.1.
[14]EvansJ G， 2009.Long-path scintillometry over complex terrain to determine areal-averaged sensible and latent heat fluxes［D］.Ph.D.thesis， The University of Reading， 176.
[15]EzzaharJ， ChehbouniA， HoedjesJ C B， alet， 2007.The use of the scintillation technique for monitoring seasonal water consumption of olive orchards in a semi-arid region［J］.Agricultural Water Management， 89（3）： 173-184.
[16]GiorgiF， AvissarR， 1997.Representation of heterogeneity effects in earth system modelling： Experience from land surface modelling［J］.Reviews of Geophysics， 35： 413-437.
[17]Henderson-SellersA， PitmanA J， LoveP K， alet， 1995.The Project for Intercomparison of Land Surface Parameterization Schemes （PILPS）： Phase 2 and Phase 3［J］.Bulletin of the American Meteorological Society， 74（4）： 489-503.
[18]HemakumaraH M， ChandrapalaL， MoeneA F， 2003.Evapotranspiration fluxes over mixed vegetation areas measured from large aperture scintillometer［J］.Agricultural Water Management， 58： 109-122.
[19]HillR J， BohlanderR A， CliffordS F， alet， 1988.Turbulence-induced millimeter-wave scintillation compared with micrometeorological measurements［J］.IEEE Transactions on Geoscience and Remote Sensing， 26（3）： 330-342.
[20]KleisslJ， HongS H， HendrickxJ M H， 2009.New Mexico scintillometer network supporting remote sensing and hydrologic and meteorological models［J］.Bulletin of the American Meteorological Society， 90（2）： 207-218.DOI： 10.1175/2008BAMS2480.1.
[21]KljunN， CalancaP， RotachM W， alet， 2015.The simple two-dimensional parameterization for Flux Footprint Prediction （FFP）［J］.Geoscientific Model Development Discussions， 8（8）： 6757-6808.DOI： 10.5194/gmdd-8-6757-2015.
[22]KolmogorovA N， 1941.The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers［J］.C.R.Acad.Sci.URSS， 30： 301-305.
[23]KooijmansL M J， HartogensisO K， 2016.Surface-layer similarity functions for dissipation rate and structure parameters of temperature and humidity based on eleven field experiments［J］.Boundary-Layer Meteorology， 160： 501-527.
[24]KormannR， MeixnerF X， 2001.An analytical footprint model for non-neutral stratification［J］.Boundary-Layer Meteorology， 99： 207-224.
[25]KustasH， NormanJ， 1999.Evaluation of soil and vegetation heat flux predictions using a simple two-source model with radiometric temperatures for partial canopy cover［J］.Agricultural and Forest Meteorology， 94（1）： 13-29.
[26]LawrenceD M， FisherR A， KovenC D， alet， 2019.The Community Land Model Version 5： Description of new features， benchmarking， and impact of forcing uncertainty［J］.J.Advances in Modeling earth Systems， 11， 4245-4287.
[27]LeclercM Y， FokenT， 2014.Footprints in Micrometeorology and Ecology［M］.Springer， Heidelberg， New York， Dordrecht， London， 239 pp.
[28]LeeS H， LeeJ H， KimB Y， 2015.Estimation of turbulent sensible heat and momentum fluxes over a heterogeneous urban area using a large aperture scintillometer［J］.Advances in Atmospheric Sciences， 32： 1-14.
[29]LiX， ChenG D， LiuS M， alet， 2013.Heihe Watershed Allied Telemetry Experimental Research （HiWATER）： Scientific objectives and experimental design［J］.Bulletin of the American Meteorological Society， 94（8）： 1145-1160.
[30]LiuS M.Xu Z W， Wang W Z， et al， 2011.A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem［J］.Hydrology & Earth System Sciences， 15（4）： 1291-1306.
[31]LiuS M， LiX， XuZ W， alet， 2018.The Heihe Integrated Observatory Network： A basin scale land surface processes observatory in China［J］.Vadose Zone Journal， 17： 1-21.
[32]LuH， ZhengD H， YangK， alet， 2020.Last-decade progress in understanding and modeling the land surface processes on the Tibetan Plateau［J］.Hydrology and Earth System Sciences， 24： 5745-5758.
[33]LudiA， BeyrichF， MatzlerC， 2005， Determination of the turbulent temperature-humidity correlation from scintillometric measurements［J］ Boundary-Layer Meteorology， 117： 525-550.
[34]MahrtL， VickersD， SunJ， alet， 2001.Calculation of area-averaged fluxes： Application to BOREAS［J］.Environmental Science， 40： 915-920.
[35]McJannetD， CookF， McGloinR， alet， 2013.Long-term energy flux measurements over an irrigation water storage using scintillometry［J］.Agricultural and Forest Meteorology， 168： 93-107.
[36]MengelkampH， BeyrichF， HeinemannG， alet， 2006.Evaporation over a heterogeneous land surface-The EVA-GRIPS project［J］.Bulletin of the American Meteorological Society， 87： 775-786.
[37]MoninA， YaglomA， 1975.Statistical fluid mechanics： mechanics of turbulence.Volume 2［M］.The MIT Press， Cambridge， Massachusetts， and London， England， 862.
[38]OwensJ C， 1967.Optical refractive index of air： Dependence on pressure， temperature and composition［J］.Applied Optics.6（1）： 51-59.
[39]SellersP J， DickinsonR E， RandallD A， alet， 1997a.Modeling the exchanges of energy， water， and carbon between continents and the atmosphere［J］.Science， 275（5299）： 502-509.
[40]SellersP J， HeiserhM D， Ha11 F G， et al， 1997b.The impact of using area-averaged land surface properties--topography， vegetation condition， soil wetness-in calculations of intermediate scale （approximately 10 km’） surface-atmosphere heat and moisture fluxes［J］.Journal of Hydrology， 190： 269-361.
[41]SchuttemeyerD， SchillingsC， MoeneA F， alet， 2007.Satellite-based actual evapotranspiration over drying semiarid terrain in west Africa［J］.Journal of Applied Meteorology & Climatology， 46： 97-111.
[42]SchuttemeyerD， MoeneA F， HoltslagA A M， alet， 2008.Evaluation of two land surface schemes used in terrains of increasing aridity in West Africa［J］.Journal of Hydrometeorology.9（2）： 173-193.
[43]SteeneveldG J， TolkL F， MoeneA F， alet， 2011.Confronting the WRF and RAMS mesoscale models with innovative observations in the Netherlands： evaluating the boundary layer heat budget［J］.Journal of Geophysical Research Atmospheres， 116： D23114.
[44]SuB， 2002.The Surface Energy Balance System （SEBS） for estimation of turbulent heat fluxes［J］.Hydrology and Earth System Sciences， 6： 85-100.
[45]TatarskiiV I， 1961.Wave propagation in a turbulent medium［M］.McGraw-Hill Book Company Inc.， New York， 285pp（中译本： 湍流大气中波的传播理论 ［M］.译者： 温景嵩等， 科学出版社， 1978， 232 pp）.
[46]van KesterenB， HartogensisO K， Van DintherD， alet， 2013.Measuring H2O and CO2 fluxes at field scales with scintillometry： Part I-Introduction and validation of four methods［J］.Agricultural and Forest Meteorology， 178/179： 75-87
[47]WardH C， EvansJ G， HartogensisO K， alet， 2013.A critical revision of the estimation of the latent heat flux from two-wavelength scintillometry［J］.Quarterly Journal of the Royal Meteorological Society， 139： 1912-1922.
[48]WardH C， EvansJ G， GrimmondC S B， alet， 2015.Infrared and millimetre-wave scintillometry in the suburban environment-Part 1： Structure parameters［J］.Atmospheric Measurement Techniques， 8： 1385-1405.
[49]WardH C， 2017.Scintillometry in urban and complex environments： A review［J］.Measurement Science & Technology， 28（6）： 064005.DOI： 10.1088/1361-6501/aa5e85.
[50]Wang J M， 2018.Area-averaged Flux Measurements and Scintillometry［Z］.3rd Intl.Potsdam GHG Flux Workshop， Nanjing， China， 22-25 October， 2018.
[51]WangT， OchsG R， CliffordS F， 1978.A saturation-resistant optical scintillometer to measure Cn2［J］.Journal of the Optical Society of America 68（3）： 334-338.
[52]WeselyM L， 1976.The combined effect of temperature and humidity fluctuations on refractive index［J］.Journal of Applied Meteorology， 15（1）： 43-49.
[53]WyngaardJ C， IzumiY， CollinsS A， 1971.Behavior of the refractive-index-structure parameter near the ground［J］.Journal of the Optical Society of America， 61（12）： 1646-1650.
[54]ZhangH， ZhangH S， CaiX H， alet， 2016.Contribution of Low-Frequency Motions to Sensible Heat Fluxes over Urban and Suburban Areas［J］.Boundary-Layer Meteorology， 161（1）： 183-201.
[55]ZhaoP， XuX D， ChenF， alet， 2018.The third atmospheric scientific experiment for understanding the earth-atmosphere coupled system over the Tibetan Plateau and its effects［J］.Bulletin of the American Meteorological Society， 99（4）： 757-776.
[56]蔡旭晖， 2008.湍流微气象观测的印痕分析方法及其应用拓展［J］.大气科学， 32（1）： 123-132.
[57]卢俐， 刘绍民， 徐自为， 等， 2009.不同下垫面大孔径闪烁仪观测数据处理与分析［J］.应用气象学报， 20（2）： 171-178.
[58]王介民， 1999.陆面过程实验和地气相互作用研究—从HEIFE到IMGRASS 和GAME-Tibet/ TIPEX［J］.高原气象， 18（3）： 280-294.
