中尺度数值大气模式WRF在水文气象领域的研究
开篇：润墨网以专业的文秘视角，为您筛选了一篇中尺度数值大气模式WRF在水文气象领域的研究范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！
摘要：近几年中尺度数值大气模式wrf迅速发展，其应用越来越广。为阐释WRF模式的机理，揭示其发展方向，为相关领域的研究人员提供参考，介绍了WRF模式系统构成，探讨和总结了物理参数化方案、数据同化方法、以及合理的空间尺度对WRF模式在水文气象领域的模拟预报效果的影响。相关研究表明：（1）针对不同区域、不同时段的研究应选取不同的物理参数化方案组合；（2）目前WRF模式常用的数据同化方法是三维变分数据同化，混合数据同化方法是改进WRF模式模拟预报效果的发展方向；（3）并非水平分辨率越高，模拟效果越好，应依据实际情况选用合理的空间尺度；（4）WRF模式的模拟效果好、预报精度高，未来有着更加广阔的应用前景。
关键词：WRF模式；物理参数化方案；数据同化；空间尺度；应用进展
中图分类号：P458.2 文献标志码：A 文章编号：16721683（2015）06102506
Abstract：With the rapid development of mesoscale numerical atmospheric model WRF in recent years， its application is wider and wider.In order to explain the mechanism， reveal the development direction and provide reference for related fields researchers， WRF model system is introduced， it is discussed that the influence of parameterized physical processes， data assimilation and reasonable spatial scale on the simulation effect of WRF model.Related studies have shown that： （1） different parameterized physical processes need to be chosen in different region and different time；（2） the common data assimilation method is 3DVar data assimilation at present， hybrid data assimilation may be a better method to improve the simulation effect of WRF model；（3） reasonable spatial scales need to be chosen because not the higher horizontal resolution， the better simulation effect for all research；（4） WRF model has good simulation effect， high prediction accuracy， there is more wide application.
Key words：WRF model；parameterized physical processes；data assimilation；spatial scale；application progress
为解决区域性气象问题，减少或避免中尺度天气造成的自然灾害，到20世纪80年代，中尺度数值大气模式逐渐成熟。20世纪90年代，部分中尺度数值大气模式已经非常先进，比较有代表性的是美国的Eta模式、MM5模式、RAMS模式、RSM模式、COAMPS模式以及WRF模式等，英国的UKMO模式，加拿大的MC2模式，法国的MESONH模式，日本的JRSM模式。近几年我国也在自主研发中尺度数值大气模式，其中中科院大气物理研究所的REM模式和中国气象科学研究院的GRAPES模式应用较广。与全球模式和区域模式相比，中尺度数值大气模式的研究区域更小，网格距一般为10 km以下，能够更精细地模拟出局部地区动量和能量的湍流输送，提高局部地区气象要素的模拟精度，可用来描述局部强对流天气，如龙卷风、雷雨等。
在众多中尺度数值大气模式中，新一代中尺度数值大气模式WRF以其较高的预报精度和灵活的分辨率，且免费对外开放，在水文气象领域具有广阔的发展与应用前景[1]。与其前身MM5模式相比，WRF模式拥有更多元化的模式动力框架和物理参数化方案，且同时支持单向和双向嵌套。在模拟和预报我国不同地区、不同性质的强降水过程中，WRF模式的较早版本在模拟天气系统和降水落区上已明显优于MM5模式[2]，其对于高度场、风场等物理量的模拟效果也优于MM5[3]。与其它中尺度大气模式相比，WRF模式在预报各种天气中都具有较好的性能，能够比较成功地再现中尺度过程中环流形势的演变和雨带分布特征[45]。
WRF模式结合先进的数值方法和资料同化技术，采用改进的物理过程方案，同时具有多重嵌套及易定位于不同地理位置的能力，已逐渐成为改进从云尺度到天气尺度等不同尺度重要天气特征预报精度的工具。但WRF模式在我国的研究与应用仍处于起步阶段，对其基本理论与应用方向认识不够全面。本文对WRF模式在物理参数化方案的选择、数据同化方法以及空间尺度上结合国内外研究进行了探讨和梳理，为WRF模式在水文气象领域的进一步应用提供参考。
1 WRF模式介绍
WRF模式是由美国国家大气研究中心（NCAR）、国家大气海洋局预报系统实验室（FSL）、国家环境预报中心（NCEP）、俄克拉荷马大学暴雨分析中心等科研单位共同研发的新一代中尺度数值大气模式，于2000年第一版，目前最新版是2014年的V3.6.1。WRF模式分为ARW和NMM两种，前者用于科学研究，后者用于实际业务，分别由NCEP和NCAR管理维持，本文仅探讨WRF ARW。
WRF模式是完全可压缩的非静力中尺度模式，由预处理系统、主程序、后处理三部分构成。通过收集外部数据，利用预处理系统对数据进行处理、控制与引导，最终形成namelist.wps文件，进而导入主程序进行计算，最终通过后处理辅助分析模拟结果，包括数据的可视化操作等。WRF采用高度模块化和分层设计，用户可根据实际情况做出不同选择，使比较模式性能和进行集合预报成为可能。模式系统组成与流程分别见图1、图2。
从云尺度到天气尺度是WRF模式预报的重点，主要考虑1～10 km的水平分辨率。最新版V3.6.1水平方向采用大气模式常用的ArakawaC网格，随着水平分辨率的提高，该网格优势更加凸显，垂直方向采用地形跟随的质量坐标（即从0～1的eta分层），时间积分采用时间分裂积分法，声波项采用时间分裂小步长方案，非声波项多采用三阶RungeKutta算法，主要的控制方程有动量守恒方程、质量守恒方程、能量守恒方程、气体状态方程、水汽混合比守恒方程等，其控制方程以通量形式表达，均提高了模拟的准确性。更多WRF模式的相关机理介绍和适用说明可参考NCAR提供的WRF模式使用者手册[6]。
2 物理参数化方案
WRF模式包含的物理过程较多（主要物理过程见表1），不同物理参数化方案的组合对模拟结果会产生较大影响，直接影响预报的准确性[78]，因此物理参数化方案的选择是WRF模式研究的重点。气象过程在不同地区、不同时段具有唯一性，目前气象学中的机理研究并不能完全准确地描述大气运动过程，目前对WRF模式物理参数化的研究需要进行对比试验，才能确定方案的优劣，缺乏机理分析。由于微物理过程、积云对流参数化过程、陆面过程、边界层对温度、湿度、气压、风速等主要气象因素以及降雨产生过程的影响比较显著，本文主要探讨这几个物理过程的参数化方案的选择。
2.1 微物理过程
微物理过程主要是指云粒子的形成、增长以及产生降水的微观物理过程，涉及相变潜热的释放与吸收，其选取结果影响积云对流发生发展条件，从而影响积云降水预报，在WRF模式的模拟中对降水模拟结果的影响较大。各微物理过程在不同地区、不同场次、不同类型的降水模拟中，表现出不同的适用性。张少济等[9]采用12种微物理方案，对山西省一次强降雨过程进行模拟，表明不同微物理过程模拟的降水情况差别较大，WSM6方案、M2M方案和Thompson方案较好地模拟了此次强降水过程；诸葛丰林等[10]利用9种微物理过程对江苏省里下河地区一次暴雨进行模拟，表明各方案基本模拟出了雨带的走向、降水落区，但降水强度以及降水中心均有不同程度的偏差；Kim等[11]选用了6种微物理方案分析了其对降雨模拟的精确性的影响，结果表明Lin方案和Thompson方案模拟对流雨的精度不高；Efstathiou等[12]选用不同的微物理方案对Chalkidiki半岛进行了暴雨分析，分析选用YSU边界层时，Eta方案的模拟结果优于Lin方案；Rajeevan等[13]选用了4种微物理方案对印度Gadanki地区的一次降水进行了模拟，结果表明Thompson方案的模拟结果略优于Lin、WSM6、Morrison方案。
2.2 积云对流过程
积云对流过程伴随着云团生成、发展、旺盛以及由于温度不同时形成的对流，云团消失则降水结束，因此其与降水过程联系紧密，不同的积云对流方案对不同地区、不同场次的降水模拟效果不同。Jankov等[14]通过对不同物理方案进行大量的模拟比较，认为模拟结果对积云对流方案是最敏感的，且不同积云对流方案中KF方案的模拟结果最优；Rama等[15]对印度的三次暴雨进行了模拟，表明BMJ方案更好地模拟了暴雨的中心和强度； Biswas等[16]选用了3种积云对流方案对大西洋和北太平洋东部的强降水进行了模拟，结果显示SAS方案更好完成预报；廖镜彪等[17]对广东一次降水进行模拟，结果表明积云对流参数化方案KF与微物理方案Lin组合条件下，模拟结果与实测结果最接近；屠妮妮等[18]选用4种积云对流参数化方案对四川省的一次区域性暴雨进行模拟，结果表明KF和GD总体预报效果好于BMJ和SAS。
2.3 陆面过程
陆面过程直接影响地表水热通量，并通过影响模拟的低层大气环流和水汽辐合量，从而影响气温等气象要素和降雨的发生。模拟时应结合实际，选择合适的陆面过程。Jimenez等[19]认为WRF的模拟结果一定程度上依赖于陆面过程方案；Khvorostyanov等[20]选用了4种陆面方案对Paris地区进行了模拟对比，结果表明各方案的适用性不同；Jin等[21]研究了4种不同陆面方案的选取对美国西部一段时期气温模拟的影响，发现结构最复杂的CLM3方案模拟效果优于STD、Noah、RUC3个方案；曾新民等[22]采用3种陆面方案对江南及华南地区的高温天气进行了模拟，结果表明模拟高温对不同陆面方案比较敏感，SLAB方案模拟最好，RUC次之，Noah再次之；李安泰等[23]选取3种陆面方案对舟曲“8.8”暴雨天气过程进行了模拟，结果表明采用陆面方案明显优于不采用陆面方案，其中PX陆面方案比SLAB、Noah方案更接近实况；张瑛等[24]采用4种陆面方案对江西省“6.19”暴雨过程进行了模拟研究，发现强降水的特性在无陆面方案情况下偏差较大，而SLAB、Noah、RUC、PX方案的模拟结果均较好的反映了降水的情况，模拟结果相差不大。
2.4 边界层
边界层在WRF模式中布局影响低层大气要素，而且通过垂直输送也可对高层大气产生影响，进而影响对温度、风速、湿度等主要气象要素的模拟，因而模拟时应根据具体的模拟要素，分析并选择合适的边界层参数化方案。Flaounas等[25]对西非季风个例进行了模拟，研究发现选取不同的边界层参数化方案对温度的影响最大，YSU方案模拟的温度较实测温度高，而MYJ方案模拟的温度相对较低；Cha等[26]指出YSU方案降低了边界层的垂直混合作用，因此比MRF方案模拟效果好；Miglietta等[27]用2种边界层参数化方案模拟地中海东部的一日风场，分析认为YSU方案优于MYJ方案；王子谦等[28]基于4种边界层参数化方案分别对东亚夏季风进行模拟研究，结果表明BL方案和MYJ方案对东亚夏季风及其降水的模拟优于YSU和ACM2方案；张龙等[29]利用3种边界层参数化方案对兰州冬季几个气象要素进行了模拟，结果表明总体上MYJ方案对低层风温、感热等模拟结果优于YSU和ACM2；张小培等[30]对安徽黄山及周边复杂地形地区的气象要素进行了模拟，结果表明YSU、ACM2、MYN、BL方案对气温、露点温度、风速、湿度等的模拟各有优劣且误差均较小。
6 讨论与结论
（1）对于绝大多数气象要素，WRF模式具有较高的模拟精度和运算效率，未来在水文气象领域有更广阔的应用前景。
（2）物理参数化方案和空间尺度的选取是WRF模式在应用过程中的重点和难点，不同地区、不同时段的气象要素模拟所适用的最优化的物理参数化方案不一定相同，水平分辨率也并非越高越好，尤其是与其他模型耦合时，应安排不同尺度组合方案进行试验，研究最优的耦合尺度。
（3）目前WRF模式最常用的数据同化方法是三维数据同化，可以大幅提升模式的模拟和预报精度，为能进一步兼顾计算效率和计算精度的提高，混合数据同化将成为未来WRF模式数据同化方法的发展趋势。
参考文献（References）：
[1]Skamaraock W C，Klemp J B，Dudhia J，et al.A description of the advanced research WRF version[Z].NCAR Technical Notes，2008.
[2]孙健，赵平.用WRF与MM5模拟1998年三次暴雨过程的对比分析[J].气象学报，2003，61（6）： 692701.（SUN Jian，ZHAO Ping.Simulation and analysis of three heavy rainfall processes in 1998 with WRF and MM5[J].Acta Meteorologica Sinica，2003，61（6）：692701.（in Chinese））
[3]刘宁微，王奉安.WRF和MM5模式对辽宁暴雨模拟的对比分析[J].气象科技，2006，34（4）：364369.（LIU Ningwei，WANG parative analysis of heavy rainfall simulations with two mesoscale models in Liaoning province[J].Meteorological Science and Technology，2006，34（4）：364369.（in Chinese））
[4]Welsh P，Wildman A，Shaw B，et al.Implementing the weather research and forecast （WRF） model with local data assimilation in a NWSWFO，Seattle，January 1116，2004[C].USA：Seattle，2004.
[5]张芳华，马旭林，杨克明.2003年6月24日-25日江南特大暴雨数值模拟和诊断分析[J].气象，2004，30（1）：2832.（ZHANG Fanghua，MA Xulin，YANG Keming.Numerical simulation and diagnostic analysis of a heavy rainfall in Jiangnan area during 2425 June 2003[J].Meteorological Monthly，2004，30（1）：2832.（in Chinese））
[6]National Center for Atmospheric Research （NCAR） （2014），User’s Guide for the Advanced Research WRF （ARW） Modelling System Version3.6，http：//www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3.6/ARWUsersGuideV3.6.1.pdf
[7]Gallus J W，Bresch J parison of impacts of WRF dynamic core，physics package，and initial conditions on warm season rainfall forecasts[J].Monthly Weather Review，2006，134：26322641.
[8]Evans J P，Ekstrom M，Ji F.Evaluating the performance of a WRF physics ensemble over SouthEast Australia[J].Climate Dynamics，2012，39：12411258.
[9]张少济，孙科，张华，等.WRF模式微物理方案对山西省一次强降雨过程模拟的影响[J].水利水电技术，2013，44（9）：811.（ZHANG Shaoji，SUN Ke，ZHANG Hua，et al.Impact from microphysical scheme of WRF model on simulation of a heavy rainfall process in Shanxi Province[J].Water Resources and Hydropower Engineering，2013，44（9）：811.（in Chinese））
[10]诸葛丰林，郑有飞，吴荣军，等.江苏省里下河地区一次暴雨过程的模拟研究[J].自然灾害学报，2014，23（5）：164176.（ZHUGE Fenglin，ZHENG Youfei，WU Rongjun，et al.Study on simulation of a rain process in Lixia River area，Jiangsu Province[J].Journal of Natural Disasters，2014，23（5）：164176.（in Chinese））
[11]Kim J H，Shin D B，Kummerow C.Impacts of a priori databases using six WRF microphysics schemes on passive microwave rainfall retrievals[J].Journal of Atmospheric and Oceanic Technology，2013，30：23672381.
[12]Efstathiou G A，Zoumakis N M，Melas D，et al.Sensitivity of WRF to boundary layer parameterizations in simulating a heavy rainfall event using different microphysical schemes[J].Atmospheric Research，2013，132133：125143.
[13]Rajeevan M，Kesarkar A，Thampi S B，et al.Sensitivity of WRF cloud microphysics to simulations of a severe thunderstorm event over southeast India[J].Annales Geophysicae，2010，28（2）：603619.
[14]Jankov，Isidora，et al.An investigation of IHOP convective system predictability using a matric of 19 WRF Members[M].84th AMS Annual Meeting，Seattle，U.S.A Jan.1015，2004.
[15]Rama R Y V，Hatwar H R，Kamal S A，et al.An experiment using the high resolution eta and WRF models to forecast heavy precipitation over India[J].Pure & Applied Geophysics，2007，164（8）：15931615.
[16]Biswas M K，Bernardet L，Dudhia J.Sensitivity of hurricane forecasts to cumulus parameterizations in the HWRF model[J].Geophysical Research Letters，2014，41（24）：91139119.
[17]廖镜彪，王雪梅，夏北成，等.WRF模式中微物理和积云参数化方案的对比试验[J].热带气象学报，2012，28（4）：461470.（LIAO Jingbiao，WANG Xuemei，XIA Beicheng，et al.The effects of different physics and cumulus parameterization schemes in WRF on heavy rainfall simulation in PRD[J].Journal of Tropical Meteorology，2012，28（4）：461470.（in Chinese））
[18]屠妮妮，何光碧，张利红.WRF模式中不同积云对流参数化方案对比试验[J].高原山地气象研究，2011，31（2）：1825.（TU Nini，HE Guangbi，ZHANG Lihong.Simulation test on the effects of various cumulus parameterization schemes in WRF model[J].Plateau and Mountain Meteorology Research，2011，31（2）：1825.（in Chinese））
[19]Jimenez P A，Dudhia J，GonzalezRouco J F，et al.A revised scheme for the WRF surface layer formulation[J].Monthly Weather Review，2012，140（3）：898918.
[20]Khvorostyanov D V，Menut L，Dupont J C，et al.The role of WRF land surface schemes on weather simulations in Paris area[A].Proceedings of ISARS 2010 conference[C].2010.
[21]Jin J，Miller N L，Schlegel N.Sensitivity study of four land surface schemes in the WRF model[J].Advances in Meteorology，2010，2010：111.
[22]曾新民，吴志皇，熊仕焱，等.WRF模式短期高温天气模拟对陆面方案的敏感性[J].中国科学D辑：地球科学，2011，41（9）：13751384.（ZENG Xinmin，WU Zhihuang，XIONG Shiyan，et al.Sensitivity of simulated shortrange hightemperature weather to land surface by WRF[J].Sci China Earth Sci，2011，41（9）：13751384.（in Chinese））
[23]李安泰，王亚明，何宏让，等.耦合不同陆面方案的WRF模式对“8.8”舟曲暴雨过程的模拟[J].气象与减灾研究，2014，37（1）：2128.（LI Antai，WANG Yaming，HE Hongyang，et al.Numeric simulations of rainstorm in ZhouQu county on 8 august 2010 using the WRF model with different land surface schemes[J].Meteorology and Disaster Reduction Research，2014，37（1）：2128.（in Chinese））
[24]张瑛，肖安，马力，等.WRF耦合4个陆面过程对“6.19”暴雨过程的模拟研究[J].气象，2011，37（9）：10601069.（ZHANG Ying，XIAO An，MA Li，et al.Simulation of “19 June 2010” heavy rainfall by using WRF coupled with four land surface processes[J].Meteorological Monthly，2011，37（9）：10601069.（in Chinese））
[25]Flaounas E，Bastin S，Janicot S.Regional climate modelling of the 2006 West African monsoon：sensitivity to convection and planetary boundary layer parameterisation using WRF[J].Climate Dynamics，2011，36：10831105.
[26]Cha D H，Lee D K，Hong S Y.Impact of boundary layer processes on seasonal simulation of the East Asian summer monsoon using a Regional Climate Model[J].Meteorol Atmos Phys，2008，100：5372.
[27]Miglietta M M，Zecchetto S，Biasio F D.WRF model and ASARretrieved 10m wind field comparison in a case study over eastern Mediterranean Sea[J].Adv Sci Res，2010，4：8388.
[28]王子谦，段安民，吴国雄.边界层参数化方案及海气耦合对WRF模拟东亚夏季风的影响[J].中国科学D辑：地球科学，2014，44（3）：548562.（WANG Ziqian，DUAN Anmin，WU Guoxiong.Impacts of boundary layer parameterization schemes and airsea coupling on WRF simulation of the East Asian summer monsoon[J].Sci China Earth Sci，2014，44（3）：548562.（in Chinese））
[29]张龙，张镭，王颖，等.基于WRF不同边界层方案的黄土高原丘陵冬季地面气象要素日变化模拟分析[J].干旱气象，2012，30（2）：158168.（ZHANG Long，ZHANG Lei，WANG Ying，et al.Simulation of daily variations about surface meteorological factor in loess plateau region in winter based on different PBL schemes of WRF[J].Journal of Arid Meteorology，2012，30（2）：158168.（in Chinese））
[30]张小培，银燕.复杂地形地区WRF模式四种边界层参数化方案的评估[J].大气科学学报，2013，36（1）：6876.（ZHANG Xiaopei，YIN Yan.Evaluation of the four PBL schemes in WRF model over complex topographic areas[J].Transactions of Atmospheric Sciences，2013，36（1）：6876.（in Chinese））
[31]Dimet F X L，Talagrand O.Variational algorithms for analysis and assimilation of meteorological observations：Theoretical aspects[J].Tellus Series ADynamic Meteorology and Oceanography，1986，38A（2）：97110.
[32]Liu J，Bray M，Han D.Exploring the effect of data assimilation by WRF3DVar for numerical rainfall prediction with different types of storm events[J].Hydrological Processes，2013，27：36273640.
[33]Maiello I，Ferretti R，Gentile S，et al.Impact of radar data assimilation for the simulation of a heavy rainfall case in central Italy using WRF3DVAR[J].Atmospheric Measurement Techniques，2014，7：29192935.
[45]Done J，Davis C A，Weisman M.The next generation of NWP：explicit forecasts of convection using the weather research and forecasting （WRF） model[J].Atmospheric Science Letters，2004，5：110117.
[46]Schwartz C S，Kain J S，Weiss S J，et al.Nextday convectionallowing WRF model guidance：a second look at 2km versus 4km grid spacing[J].Monthly Weather Review，2009，137：33513372.
[47]Ma P L，Rasch P J，Fast J D，et al.Assessing the CAM5 physics suite in the WRFChem model：implementation，resolution sensitivity，and a first evaluation for a regional case study[J].Geoscientific Model Development，2014，7：755778.
[48]郝春沣，贾仰文，王浩.气象水文模型耦合研究及其在渭河流域的应用[J].水利学报，2012，43（9）：10411049.（HAO Chunfeng，JIA Yangwen，WANG Hao.Atmospheric and hydrological models'coupling and application in flood forecasting of the Weihe Basin[J].Journal of Hydraulic Engineering，2012，43（9）：10411049.（in Chinese））
[49]彭艳，周建中，贾梦，等.三峡库区陆气耦合研究与应用[J].水文，2014，34（3）：1116.（PENG Yan，ZHOU Jianzhong，JIA Meng，et al.A coupled atmospherichydrological modeling system and its application for Three Gorges reservoir region[J].Journal of China Hydrology，2014，34（3）：1116.（in Chinese））
[50]Wagner S，Fersch B，Kunstmann H，et al.Coupled atmospherichydrological modeling for feedback investigations in the Poyang lake catchment，China[C].AGU Fall Meeting Abstracts，2012.
[51]Givati A，Lynn B，Liu Y，et al.Using the WRF model in an operational streamflow forecast system for the Jordan river[J].Journal of Applied Meteorology and Climatology，2012，51：285299.
[52]Hamill T M.Performance of operational model precipitation forecast guidance during the 2013 Colorado frontrange floods[J].Monthly Weather Review，2014，142：26092618.
[53]高洋.WRF模式对2008年1月我国南方冻雨极端天气过程的数值模拟研究[D].北京：中国气象科学研究院，2011.（GAO Yang.WRF simulations of the freezing weather in Jan 2008 over the southern China[D].Beijing：Chinese Academy of Meteorological Sciences，2011）
[54]Jimenez P A，Dudhia J.Improving the representation of resolved and unresolved topographic effects on surface wind in the WRF model[J].Journal of Applied Meteorology and Climatology，2012，51（2）：300316.
[55]文小航，廖小罕，袁文平，等.中国东北半干旱区能量水分循环的同化模拟[J].中国科学D辑：地球科学，2014，44（12）：27682784.（WEN Xiaohang，LIAO Xiaohan，YUAN Wenping，et al.Numerical simulation and data assimilation of the waterenergy cycle over semiarid northeastern China[J].Sci China Earth Sci，2014，44（12）：27682784.（in Chinese））
[56]崔琳琳，胡松，杨红，等.绿潮早期聚集期间天气过程分析[J].海洋环境科学，2014，33（6）：941946.（CHUI Linlin，HU Song，YANG Hong，et al.Weather process analysis for the period during the early aggregation of green tide[J].Marine Environmental Science，2014，33（6）：941946.（in Chinese））
[57]侯爱中，倪广恒，丛振涛，等.耦合城市冠层模型对暴雨模拟精度的影响[J].水力发电学报，2013，32（6）：3135.（HOU Aizhong，NI Guangheng，CONG Zhengtao，et al.Impacts of coupling urban canopy model on storm simulation[J].Journal of Hydroelectric Engineering，2013，32（6）：3135.（in Chinese））
[58]朱坚，杨涛，黄丹青，等.基于WRFNoahLSM耦合模型的江苏省淮北地区土壤旱情评估[J].水电能源科学，2014，32（10）：1013.（ZHU Jian，YANG Tao，HUANG Danqing，et al.Soil drought assessment of Huaibei regions in Jiangsu province based on WRFNoahLSM coupled model[J].Water Resources and Power，2014，32（10）：1013.（in Chinese））
[59]赵求东.WRF+DHSVM融雪径流预报模式研究[D].新疆：新疆大学，2008.（ZHAO Qiudong.The study of WRF+DHSVM snowmelt runoff forecasting model[D].Xinjiang：Xinjiang University，2008.（in Chinese））
[60]邱辉，訾丽.WRF模型在山洪灾害预警预报中的试验应用[J].人民长江，2013，44（13）：59.（Qiu Hui，Zi Li.Trial application of WRF model in early warning and forecast of mountain torrent disaster[J].Yangtze River，2013，44（13）：59.（in Chinese））
[61]郑子彦，张万昌，徐精文.山区流域暴雨洪水的数值模拟[J].山地学报，2012，30（2）：222229.（ZHENG Ziyan，ZHANG Wanchang，XU Jingwen.Numeric simulation of torrential rainfall and flashflood in a mountainous basin[J].Journal of Mountain Science，2012，30（2）：222229.（in Chinese））
[62]Warner J C，Armstrong B，He R，et al.Development of a coupled oceanatmospherewavesediment transport （COAWST） modeling system[J].Ocean Modeling，2010，35：230244.