Advanced Search
Article Contents

Impact of Rain Snow Threshold Temperature on Snow Depth Simulation in Land Surface and Regional Atmospheric Models

Fund Project:

doi: 10.1007/s00376-012-2192-7

  • This study investigates the impact of rain snow threshold (RST) temperatures on snow depth simulation using the Community Land Model (CLM) and the Weather Research and Forecasting model (WRF-coupled with the CLM and hereafter referred to as WRF_CLM), and the difference in impacts. Simulations were performed from 17 December 1994 to 30 May 1995 in the French Alps. Results showed that both the CLM and the WRF_CLM were able to represent a fair simulation of snow depth with actual terrain height and 2.5℃ RST temperature. When six RST methods were applied to the simulation using WRF_CLM, the simulated snow depth was the closest to observations using 2.5℃ RST temperature, followed by that with Pipes', USACE, Kienzle's, Dai's, and 0℃ RST temperature methods. In the case of using CLM, simulated snow depth was the closest to the observation with Dai's method, followed by with USACE, Pipes', 2.5℃ RST temperature, Kienzle's, and 0℃ RST temperature method. The snow depth simulation using the WRF_CLM was comparatively sensitive to changes in RST temperatures, because the RST temperature was not only the factor to partition snow and rainfall. In addition, the simulated snow related to RST temperature could induce a significant feedback by influencing the meteorological variables forcing the land surface model in WRF_CLM. In comparison, the above variables did not change with changes in RST in CLM. Impacts of RST temperatures on snow depth simulation could also be influenced by the patterns of temperature and precipitation, spatial resolution, and input terrain heights.
    摘要: This study investigates the impact of rain snow threshold (RST) temperatures on snow depth simulation using the Community Land Model (CLM) and the Weather Research and Forecasting model (WRF-coupled with the CLM and hereafter referred to as WRF_CLM), and the difference in impacts. Simulations were performed from 17 December 1994 to 30 May 1995 in the French Alps. Results showed that both the CLM and the WRF_CLM were able to represent a fair simulation of snow depth with actual terrain height and 2.5℃ RST temperature. When six RST methods were applied to the simulation using WRF_CLM, the simulated snow depth was the closest to observations using 2.5℃ RST temperature, followed by that with Pipes', USACE, Kienzle's, Dai's, and 0℃ RST temperature methods. In the case of using CLM, simulated snow depth was the closest to the observation with Dai's method, followed by with USACE, Pipes', 2.5℃ RST temperature, Kienzle's, and 0℃ RST temperature method. The snow depth simulation using the WRF_CLM was comparatively sensitive to changes in RST temperatures, because the RST temperature was not only the factor to partition snow and rainfall. In addition, the simulated snow related to RST temperature could induce a significant feedback by influencing the meteorological variables forcing the land surface model in WRF_CLM. In comparison, the above variables did not change with changes in RST in CLM. Impacts of RST temperatures on snow depth simulation could also be influenced by the patterns of temperature and precipitation, spatial resolution, and input terrain heights.
  • 加载中
  • Auer A. H.,1974: The rain versus snow threshold temperatures. Weatherwise, 27(2), 67-67, doi: 10.1080/00431672.1974.9931684.
    Barnett T. P.,L. Dumenil,U. Schlese,E. Roeckner,M. Latif, 1989: The effect of Eurasian snow cover on regional and global climate variations. J. Atmos. Sci., 46(5), 661-685.
    Belair S.,R. Brown,J. Mailhot,B. Bilodeau,L. P. Crevier, 2003: Operational implementation of the ISBA land surface scheme in the Canadian regional weather forecast model. Part II: Cold season results. Journal of Hydrometeorology, 4(2), 371-386.
    Boone A.,P. Etchevers, 2001: An intercomparison of three snow schemes of varying complexity coupled to the same land surface model: Local-scale evaluation at an Alpine site. Journal of Hydrometeorology, 2(4), 374-394.
    Braun L. N.,H. Lang, 1986: Simulation of snowmelt runoff in lowland and lower alpine regions of Switzerland. Modelling Snowmelt-Induced Processes IAHS Publ, No. 155, 125-140.
    Chevallier P.,Y. Caballero,R. Gallaire,R. Pillco, 2004: Flow modelling in a high mountain valley equipped with hydropower plants: Rio zongo valley, Cordillera Real, Bolivia. Hydrological Processes, 18(5), 939-957, doi: 10.1002/hyp.1339.
    Collins W. D.,P. J. Rasch,B. A. Boville,J. J. Hack,J. R. McCaa,D. L. Williamson,J. T. Kiehl,B. Briegleb, 2004: Description of the NCAR community atmosphere model (CAM3). Tech. Rep. NCAR/TN-464+STR, 226pp.
    Collins, W. D.,Coauthors, 2006: The community climate system model version 3 (CCSM3). J. Climate, 19(11), 2122-2143.
    Dai A.,2008: Temperature and pressure dependence of the rain-snow phase transition over land and ocean. Geophys. Res. Lett., 35(12), doi: 10.1029/2008GL033295.
    Dickinson R. E.,A. Henderson-Sellers,P. J. Kennedy, 1993: Biosphere atmosphere transfer scheme (BATS) version 1e as coupled to the NCAR community climate model. NCAR Tech. Note NCAR/TN-387 STR 77-77, 80pp.
    Dudhia J.,1989: Numerical study of convection observed during the winter monsoon experiment using a mesoscale two dimensional model. J. Atmos. Sci., 46(20), 3077-3107.
    Essery R.,P. Etchevers, 2004: Parameter sensitivity in simulations of snowmelt. J. Geophys. Res., 109(D20), doi: 10.1029/2004JD005036.
    Essery R.,E. Martin,H. Douville,A. Fernandez,E. Brun, 1999: A comparison of four snow models using observations from an alpine site. Climate Dyn., 15(8), 583-593.
    Fassnacht S. R.,E. D. Soulis, 2002: Implications during transitional periods of improvements to the snow processes in the land surface scheme—Hydrological model WATCLASS. Atmos.-Ocean, 40(4), 389-403.
    Feiccabrino J.,A. Lundberg, 2008: Precipitation phase discrimination in Sweden. Proc.of the 65th Eastern Snow ConferenceþÿþÿFairlee, Vermont, USA.[Fern dez, 1998] Fernández, A., 1998: An energy balance model of seasonal snow evolution. Physics and Chemistry of the Earth, 23(5-6), 661-666, doi: 10.1016/S0079-1946(98)00107-4.[Gilliesetal.,2012]Gillies,R.R.,S.Y.Wa,andW.R.Huang,2012:ObservationalandsupportivemodellinganalysesofwinterprecipitationchangeinChinaoverthelasthalfcentury.Int.J.Climatol.,32(5),747-758,doi:10.1002/Joc.2303.
    Jin J. M.,L. J. Wen, 2012: Evaluation of snowmelt simulation in the weather research and forecasting model. J. Geophys. Res., 117, doi: 10.1029/2011jd016980.
    Kain J. S.,2004: The Kain-Fritsch convective parameterization: An update. J. Appl. Meteor., 43(1), 170-181.
    Kanamitsu M.,W. Ebisuzaki,J. Woollen,S. K. Yang,J. J. Hnilo,M. Fiorino,G. L. Potter, 2002: NCEP-DOE AMIP-II Reanalysis (R-2). Bull. Amer. Meteor. Soc., 83(11), 1631-1643, doi: 10.1175/Bams-83-11-1631.
    Kienzle S. W.,2008: A new temperature based method to separate rain and snow. Hydrological Processes, 22(26), 5067-5085.
    Koren V.,J. Schaake,K. Mitchell,Q. Y. Duan,F. Chen,J. M. Baker, 1999: A parameterization of snowpack and frozen ground intended for NCEP weather and climate models. J. Geophys. Res., 104(D16), 19569-19585.
    L'Hote, Y., P. Chevallier, A. Coudrain, Y. Lejeune,P. Etchevers, 2005: Relationship between precipitation phase and air temperature: comparison between the Bolivian Andes and the Swiss Alps. Hydrological Sciences Journal, 50(6), 989-997.
    Loth B.,H. F. Graf,J. M. Oberhuber, 1993: Snow cover model for global climate simulations. J. Geophys. Res., 98(D6), 10451-10464.
    Lundquist J. D.,P. J. Neiman,B. Martner,A. B. White,D. J. Gottas,F. M. Ralph, 2008: Rain versus snow in the Sierra Nevada, California: Comparing doppler profiling radar and surface observations of melting level. Journal of Hydrometeorology, 9(2), 194-211.
    Marks D. G.,A. H. Winstral, 2007: Finding the rain/snow transition elevation during storm events in mountain basins. Abstract in Joint Symposium JHW001: Interactions between snow, vegetation, and the atmosphere, the 24th General Assembly of the IUGG, Perugia, Italy, July 2-13.
    Mlawer E. J.,S. J. Taubman,P. D. Brown,M. J. Iacono,S. A. Clough, 1997: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res., 102(D14), 16663-16682.
    Morrison H.,J. A. Curry,V. I. Khvorostyanov, 2005: A new double-moment microphysics parameterization for application in cloud and climate models. Part I: Description. J. Atmos. Sci., 62(6), 1665-1677.
    Motoyama H.,1990: Simulation of seasonal snowcover based on air-temperature and precipitation. J. Appl. Meteor., 29(11), 1104-1110.
    Noh Y.,W. G. Cheon,S. Y. Hong,S. Raasch, 2003: Improvement of the K-profile model for the planetary boundary layer based on large eddy simulation data. Bound.-Layer Meteor., 107(2), 401-427.
    Oleson, K. W.,Coauthors, 2008: Improvements to the Community Land Model and their impact on the hydrological cycle. J. Geophys. Res., 113, G01021, doi: 01010.01029/02007JG000563.
    Pipes A.,M. C. Quick, 1977: UBC watershed model users guide. Department of civil engineering, University of British Columbia: Vancouver, British Columbia, Canada.
    Rohrer M.,1989: Determination of the transition air temperature from snow to rain and intensity of precipitation. Presented at Workshop on Precipitation Measurement, St. Moritz, Switzerland, 475-482.
    Schlosser C. A.,A. Robock,K. Y. Vinnikov,N. A. SperanskayaþÿþÿY. K. Xue, 1997: 18-year land-surface hydrology model simulations for a midlatitude grassland catchmenin Valdai, Russia. Mon . Wea. Rev., 125(12), 3279-3296.
    Schreider S. Y.,P. H. Whetton,A. J. Jakeman,A. B. Pittock, 1997: Runoff modelling for snow-affected catchments in the Australian alpine region, eastern Victoria. J. Hydrol., 200(1-4), 1-23, doi: 10.1016/S0022-1694(97)00006-1.
    Shewchuk S. R.,1997: Surface mesonet for BOREAS. J. Geophys. Res., 102(D24), 29077-29082, doi: 10.1029/96jd03875.
    Skamarock W. C.,J. B. Klemp, 2008: A time-split nonhydrostatic atmospheric model for weather research and forecasting applications. J. Comput. Phys., 227(7), 3465-3485.
    Strasser U.,P. Etchevers,Y. Lejeune, 2002: Inter-comparison of two snow models with different complexity using data from an alpine site. Nordic Hydrology, 33(1), 15-26.
    Subin Z. M.,W. J. Riley,J. Jin,D. S. Christianson,M. S. Torn,L. M. Kueppers, 2011: Ecosystem feedbacks to climate change in California: Development, testing, and analysis using a coupled regional atmosphere and land surface model (WRF3-CLM3.5). Earth Interactions, 15, 1-38., doi:http://dx.doi.org/10.1175/2010EI331.1.
    Sun S. F.,Y. K. Xue, 2001: Implementing a new snow scheme in simplified simple biosphere model. Adv. Atmos. Sci., 18(3), 335-354, doi: 10.1007/BF02919314.
    U.S. Army Corps of Engineers, 1956: Summary report of the snow investigation—Snow hydrology, North Pacific Division report, Portland, Oregon, 437pp.
    Walsh J. E.,W. H. Jasperson,B. Ross, 1985: Influences of snow cover and soil-Moisture on monthly air-temperature. Mon. Wea. Rev., 113(5), 756-768.
    Yang Z. L.,R. E. Dickinson,A. Robock,K. Y. Vinnikov, 1997: Validation of the snow submodel of the biosphere-atmosphere transfer scheme with Russian snow cover and meteorological observational data. J. Climate, 10(2), 353-373.
    Yeh T.-C.,R. T. Wetherald,S. Manabe, 1983: A model study of the short-term climate and hydrologic effects of sudden snowcover removal. Mon. Wea. Rev., 111, 1013-1024.
  • [1] Xuejie GAO, Ying SHI, Zhenyu HAN, Meili WANG, Jia WU, Dongfeng ZHANG, Ying XU, Filippo GIORGI, 2017: Performance of RegCM4 over Major River Basins in China, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 441-455.  doi: 10.1007/s00376-016-6179-7
    [2] WANG Aihui, ZENG Xubin, 2009: Improving the Treatment of the Vertical Snow Burial Fraction over Short Vegetation in the NCAR CLM3, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 877-886.  doi: 10.1007/s00376-009-8098-3
    [3] SONG Xiang, ZENG Xiaodong, ZHU Jiawen, 2013: Evaluating the Tree Population Density and Its Impacts in CLM-DGVM, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 116-124.  doi: 10.1007/s00376-012-1271-0
    [4] Yuanyuan WANG, Zhenghui XIE, Binghao JIA, 2016: Incorporation of a Dynamic Root Distribution into CLM4.5: Evaluation of Carbon and Water Fluxes over the Amazon, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 1047-1060.  doi: 10.1007/s00376-016-5226-8
    [5] Lian LIU, Massimo MENENTI, Yaoming MA, Weiqiang MA, 2022: Improved Parameterization of Snow Albedo in WRF + Noah: Methodology Based on a Severe Snow Event on the Tibetan Plateau, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1079-1102.  doi: 10.1007/s00376-022-1232-1
    [6] LI Weiping, LUO Yong, XIA Kun, LIU Xin, 2008: Simulation of Snow Processes Beneath a Boreal Scots Pine Canopy, ADVANCES IN ATMOSPHERIC SCIENCES, 25, 348-360.  doi: 10.1007/s00376-008-0348-2
    [7] LI Weiping, SUN Shufen, WANG Biao, LIU Xin, 2009: Numerical Simulation of Sensitivities of Snow Melting to Spectral Composition of the Incoming Solar Radiation, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 403-412.  doi: 10.1007/s00376-009-0403-7
    [8] Huanlian LI, Huijun WANG, Dabang JIANG, 2017: Influence of October Eurasian Snow on Winter Temperature over Northeast China, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 116-126.  doi: 10.1007/s00376-016-5274-0
    [9] Xuke LIU, Xiaojing JIA, Min WANG, Qifeng QIAN, 2022: The Impact of Tibetan Plateau Snow Cover on the Summer Temperature in Central Asia, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1103-1114.  doi: 10.1007/s00376-021-1011-4
    [10] Haoxin ZHANG, Naiming YUAN, Zhuguo MA, Yu HUANG, 2021: Understanding the Soil Temperature Variability at Different Depths: Effects of Surface Air Temperature, Snow Cover, and the Soil Memory, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 493-503.  doi: 10.1007/s00376-020-0074-y
    [11] GU Jianfeng, Qingnong XIAO, Ying-Hwa KUO, Dale M. BARKER, XUE Jishan, MA Xiaoxing, 2005: Assimilation and Simulation of Typhoon Rusa (2002) Using the WRF System, ADVANCES IN ATMOSPHERIC SCIENCES, 22, 415-427.  doi: 10.1007/BF02918755
    [12] Xianghui KONG, Aihui WANG, Xunqiang BI, Dan WANG, 2019: Assessment of Temperature Extremes in China Using RegCM4 and WRF, ADVANCES IN ATMOSPHERIC SCIENCES, 36, 363-377.  doi: 10.1007/s00376-018-8144-0
    [13] BI Yun, CHEN Yuejuan, ZHOU Renjun, YI Mingjian, DENG Shumei, 2011: Simulation of the Effect of an Increase in Methane on Air Temperature, ADVANCES IN ATMOSPHERIC SCIENCES, 28, 129-138.  doi: 10.1007/s00376-010-9197-x
    [14] LIAO Zhijie, ZHANG Yaocun, 2013: Simulation of a Persistent Snow Storm over Southern China with a Regional Atmosphere-Ocean Coupled Model, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 425-447.  doi: 10.1007/s00376-012-2098-4
    [15] BI Yun, CHEN Yuejuan, ZHOU Renjun, YI Mingjian, DENG Shumei, 2011: Simulation of the Effect of Water-vapor Increase on Temperature in the Stratosphere, ADVANCES IN ATMOSPHERIC SCIENCES, 28, 832-842.  doi: 10.1007/s00376-010-0047-7
    [16] XU Ying, GAO Xuejie, SHEN Yan, XU Chonghai, SHI Ying, F. GIORGI, 2009: A Daily Temperature Dataset over China and Its Application in Validating a RCM Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 763-772.  doi: 10.1007/s00376-009-9029-z
    [17] Jianfeng WANG, Ricardo M. FONSECA, Kendall RUTLEDGE, Javier MARTÍN-TORRES, Jun YU, 2020: A Hybrid Statistical-Dynamical Downscaling of Air Temperature over Scandinavia Using the WRF Model, ADVANCES IN ATMOSPHERIC SCIENCES, 37, 57-74.  doi: 10.1007/s00376-019-9091-0
    [18] Ji-Hyun HA, Dong-Kyou LEE, 2012: Effect of Length Scale Tuning of Background Error in WRF-3DVAR System on Assimilation of High-Resolution Surface Data for Heavy Rainfall Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 29, 1142-1158.  doi: 10.1007/s00376-012-1183-z
    [19] MA Xiaoyan, GUO Yufu, SHI Guangyu, YU Yongqiang, 2004: Numerical Simulation of Global Temperature Change during the 20th Century with the IAP/LASG GOALS Model, ADVANCES IN ATMOSPHERIC SCIENCES, 21, 227-235.  doi: 10.1007/BF02915709
    [20] ZHANG Jie, Laurent LI, ZHOU Tianjun, and XIN Xiaoge, 2013: Variation of Surface Temperature during the Last Millennium in a Simulation with the FGOALS-gl Climate System Model, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 699-712.  doi: 10.1007/s00376-013-2178-0

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 19 September 2012
Manuscript revised: 10 December 2012
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Impact of Rain Snow Threshold Temperature on Snow Depth Simulation in Land Surface and Regional Atmospheric Models

    Corresponding author: WEN Lijuan, gamevictory@gmail.com
  • 1. Key Laboratory of Land Surface Process and Climate Change in Cold and Arid Regions, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000; 
  • 2. Laboratory of Arid Climatic Changing and Reducing Disaster of Gansu Province, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000; 
  • 3. Institut für Umweltplanung, Gottfried Wilhelm Leibniz Universität, Hannover, 30419, Germany; 
  • 4. Department of Plants, Soils, and Climate, Utah State University, Logan, UT, 84322, USA
Fund Project:  The study is under the auspices of the National Natural Science Foundation of China (Grant No. 41275014) and the National Basic Research Program of China (973 Program) Project (Grant No. 2010CB950503). We acknowledge the support for using the computing resources at Supercomputing Center of Cold and Arid Region Environment and Engineering Research Institute of Chinese Academy of Sciences.

Abstract: This study investigates the impact of rain snow threshold (RST) temperatures on snow depth simulation using the Community Land Model (CLM) and the Weather Research and Forecasting model (WRF-coupled with the CLM and hereafter referred to as WRF_CLM), and the difference in impacts. Simulations were performed from 17 December 1994 to 30 May 1995 in the French Alps. Results showed that both the CLM and the WRF_CLM were able to represent a fair simulation of snow depth with actual terrain height and 2.5℃ RST temperature. When six RST methods were applied to the simulation using WRF_CLM, the simulated snow depth was the closest to observations using 2.5℃ RST temperature, followed by that with Pipes', USACE, Kienzle's, Dai's, and 0℃ RST temperature methods. In the case of using CLM, simulated snow depth was the closest to the observation with Dai's method, followed by with USACE, Pipes', 2.5℃ RST temperature, Kienzle's, and 0℃ RST temperature method. The snow depth simulation using the WRF_CLM was comparatively sensitive to changes in RST temperatures, because the RST temperature was not only the factor to partition snow and rainfall. In addition, the simulated snow related to RST temperature could induce a significant feedback by influencing the meteorological variables forcing the land surface model in WRF_CLM. In comparison, the above variables did not change with changes in RST in CLM. Impacts of RST temperatures on snow depth simulation could also be influenced by the patterns of temperature and precipitation, spatial resolution, and input terrain heights.

摘要: This study investigates the impact of rain snow threshold (RST) temperatures on snow depth simulation using the Community Land Model (CLM) and the Weather Research and Forecasting model (WRF-coupled with the CLM and hereafter referred to as WRF_CLM), and the difference in impacts. Simulations were performed from 17 December 1994 to 30 May 1995 in the French Alps. Results showed that both the CLM and the WRF_CLM were able to represent a fair simulation of snow depth with actual terrain height and 2.5℃ RST temperature. When six RST methods were applied to the simulation using WRF_CLM, the simulated snow depth was the closest to observations using 2.5℃ RST temperature, followed by that with Pipes', USACE, Kienzle's, Dai's, and 0℃ RST temperature methods. In the case of using CLM, simulated snow depth was the closest to the observation with Dai's method, followed by with USACE, Pipes', 2.5℃ RST temperature, Kienzle's, and 0℃ RST temperature method. The snow depth simulation using the WRF_CLM was comparatively sensitive to changes in RST temperatures, because the RST temperature was not only the factor to partition snow and rainfall. In addition, the simulated snow related to RST temperature could induce a significant feedback by influencing the meteorological variables forcing the land surface model in WRF_CLM. In comparison, the above variables did not change with changes in RST in CLM. Impacts of RST temperatures on snow depth simulation could also be influenced by the patterns of temperature and precipitation, spatial resolution, and input terrain heights.

1 Introduction
  • Snow can modify regional and, possibly, remote hydroclimatic environments by changing the surface energy and water balance (Yeh et al., 1983; Walsh et al., 1985; Barnett et al., 1989; Yang et al., 1997; Essery et al., 1999). Accurate numerical modeling of snow is a potential way to not only understand local climate but also facilitate seasonal prediction. The separation of precipitation into snow and rain remains a key challenge in the snow simulation process. Many factors, such as thickness and temperature of the atmospheric boundary layer, position of the 0℃ isotherm, cloud types, air mass, and humidity, determine the precipitation morphology (Kienzle, 2008). However, including all these factors for partitioning rain and snow in models will render complicated model physics, undermining the overall simulation efficiency of the model. In order to simplify the modeling process and improve the simulation efficiency, empirical studies have adopted climatologically parameterized values to determine the phase of precipitation. Rain snow threshold (RST) temperature is one such parameter, which is based on minimum air temperature, dew point temperature, or air temperature. The minimum air temperature was used in a rainfall-runoff model to determine the precipitation type in Australian Alpine region (Schreider et al., 1997). (Marks and Winstral, 2007) found that dew point temperature was more reliable than air temperature as a predictor of the precipitation phase in a mountain in Idaho, USA. However a study from Sweden (Feiccabrino and Lundberg, 2008) showed that air temperature is a better indicator than dew point temperature. And the RST air temperature has been reported to be the commonly used parameter (Motoyama, 1990; Yang et al., 1997; Gillies et al., 2012). Taking reference from the above examples, we employed RST air temperature for the current analysis. In our study, RST temperature referred to the air temperature. The RST temperature of 2.5℃ was first determined by (Auer, 1974) based on nearly 1000 weather observations, which separated solid and liquid precipitations in equal probabilities. (Auer, 1974) also illustrated that 0℃ (6.1℃) was the lowest (highest) temperature for rain (snow) to exist. Regional variations in RST temperature were also observed. For example, RST temperature was observed to be about 0℃ at Hokkaido (northern region) and 2℃-3℃ at Honshu (southern region) in central Japan (Motoyama, 1990). In other mountainous regions, RST temperatures varied with elevation (Lundquist et al., 2008).

    Based on observations related to RST temperatures, climate models have adopted different parameterizations to categorize the precipitation phase. General circulation models typically employ a constant RST temperature; for instance, 2.2℃ was applied as the RST temperature for the Biosphere-Atmosphere Transfer Scheme (BATS) (Dickinson et al., 1993; Yang et al., 1997). For lowland and lower Alpine regions of Switzerland, 0.5℃ was considered as the optimal transitional temperature for partitioning snowfall and rainfall in the Hydrologiska Byråns Vattenbalansavdelning (HBV) runoff model, developed by the Swedish Meteorological and Hydrological Institute (Braun and Lang, 1986). In the tropical Andes Cordillera, a distinct RST temperature of -1.5℃ was employed to distinguish snowfall and rainfall events in the Interactions between Soil, Biosphere, and Atmosphere (ISBA) model (Boone and Etchevers, 2001; Chevallier et al., 2004). The Community Atmospheric Model version 3.0 determined the percent of snowfall as a linear function of air temperature between 0℃ and -5℃ (Collins et al., 2004) and assumed that all precipitations occurred as snow (rain) when the temperature was less than -5℃(higher than 0℃). The Canadian University of British Columbia (UBC) Watershed Model (Pipes and Quick, 1977; Kienzle, 2008) applied a linear approach to depict a mixed precipitation of snow and rain, when the RST temperatures were regulated between 0.6℃ and 3.6℃.

    Apparently, there has been a wide ranging difference between RST temperatures used in different models. Simulated snow and related energy budget were sensitive to the changes in RST temperatures (Loth et al., 1993; Fassnacht and Soulis, 2002) (Loth et al., 1993) found significant differences between the simulated snow depths and snow water equivalents with three different RST temperature methods. (Fassnacht and Soulis, 2002) show that warmer RST temperatures produced more snow, larger latent heat flux and ground heat flux during the melt period, along with increased surface heat flux. However, up to now only few studies have examined the impacts of different RST temperature methods on snow simulation. Even fewer have delved into studying the impacts with offline land surface models (forced with observed meteorological data, without considering the interaction between atmosphere and land surface) and regional atmospheric models coupled with land surface models (considering the interaction between atmosphere and land surface). Our study intended to address this research gaps with an objective to study the impact and sensitivity of RST temperature methods to snow simulation using an offline land surface model (Community Land Model version 3.5, CLM) and a regional atmospheric model (Weather Research and Forecasting model version 3.2, WRF, which was recently coupled with the CLM—WRF_CLM). Section 2 of this article describes six RST methods applied to the models for comparing their impacts on snow depth simulation. Section 3 introduces the observation data and models, simultaneously providing a detailed insight into the model settings and experiments. Section 4 illustrates the snow depth simulation from 17 December 1994 (when snow started to appear) to 30 May 1995 (when snow melted out). Six popular RST temperature methods were applied to the offline land surface model CLM and regional atmospheric model WRF_CLM to demonstrate their impacts on snow depth simulation and study the feedback on the meteorological fields. In addition, sensitivity of RST temperatures (with same intervals) to snow depth simulation is discussed. Other factors that possibly affect how RST temperatures influence snow simulation are also addressed. The last section includes discussion of results and conclusions.

2 RST temperature methods
  • Six RST temperature methods were applied to the models to compare the impacts of RST methods on snow simulation in the study. The first method was based on about 1000 weather observations (Auer, 1974) and set the RST temperature as 2.5℃. Precipitation occurred as snow (rain) when the temperature was cooler (warmer) than the threshold temperature. The second method employed 0℃ RST temperature, as was used in case of the Noah Model (Koren et al., 1999). In the third method, a linear function of air temperature was used in the Canadian UBC Watershed Model (Pipes et al.,1977; Kienzle et al.,2008):

    where rp is the ratio of falling rain to precipitation and T is the air temperature (℃).

    The fourth method was similar to function (1), except the upper (lower) limit temperature and the slope, which is as follows (U.S. Army Corps of Engineers, 1956) :

    where Tt (2.6℃) is the threshold temperature at which rainfall and snowfall were equal; TR is the temperature range at which rainfall and snowfall could exist together, and it was set as 13.3℃. The sixth method was based on the three-hourly synoptic weather reports from over 15 000 land stations and many ships globally (Dai, 2008). Here, the percentage of rainfall was a hyperbolic tangent function of temperature:

    where a, b, c, and d are the annual parameters for land and equal -48.2292, 0.7205, 1.1662, and 1.0223, respectively.

3 Data and models
  • The data (Essery et al., 1999) used to force the offline land surface model CLM and validate the WRF_CLM model included hourly records of air temperature, humidity, wind speed, shortwave radiation, long wave radiation, and precipitation from 17 December 1994 to 30 May 1995, which were retrieved from the Centre d'Etudes de la Neige situated at Col de Porte (45°N, 6°E, 1320 m) over a short grassland (Fig. 1). The observed data for temperature and humidity were collected at 2 m above ground and the wind speed was measured at 2.5 m. Snow depth was monitored on an hourly basis using an ultrasonic sensor. Surface temperature and albedo were calculated using long and shortwave radiation measurements. This data set has also been reported to be employed in case of various other models, such as Chameleon Surface Model (CHASM), Simplified Simple Biosphere Model (SSIB), BATS, ISBA, Energy balance Snow Cover Integrated Model (ESCIMO), and others (Fernández, 1998; Sun and Xue, 2001; Strasser et al., 2002; Belair et al., 2003; Essery and Etchevers, 2004).

    Figure 1.  Topography of the simulated domain and the observation site (red circle).

    Figure 2.  Observed and simulated snow depths with different RST temperature methods by CLM and WRF_CLM.

    The aforementioned observed hourly interval atmospheric data were used to drive CLM and evaluate the effects of different RST temperatures on the snow depth simulation in offline land surface and regional atmospheric models. The land surface model CLM was developed by the National Center for Atmospheric Research (NCAR) and was used extensively for offline and coupled model simulations in varied landscapes globally. (Oleson and Coauthors, 2008) and (Collins et al, 2006) described the model in detail. In the present study, the offline simulation was performed during the period from 17 December 1994 (when snow started to appeared) to the end of May in 1995 (when all snow melted) using the CLM model.

    Further, the newly coupled regional atmospheric model WRF_CLM was applied to study the effects of RST temperatures on snow simulation and on the interaction between land surface and atmosphere. Details on the model can be obtained from the study of (Subin et al., 2011). The WRF model is a limited-area, nonhydrostatic, primitive-equation model with multiple options for various physical parameterization schemes (Skamarock and Klemp, 2008). The options selected for atmospheric physics in our work are listed in Table 1. The period of simulations using the WRF_CLM was in accordance with the period of study using the CLM. The simulated domain was centered at 46°N, 9°W, with 20 km horizontal grid spacing (Fig. 1). The grid point dimension was 60×60. The model had 30 vertical layers. The initial and lateral boundary conditions were provided by National Centers for Environmental Prediction (NCEP) reanalysis data version II (Kanamitsu et al., 2002), which are updated every six hours.

    To compare the impacts of different RST temperatures using offline and regional atmospheric models, we employed two models (CLM and WRF_CLM) to perform a couple of simulations. However, simulation by the WRF_CLM with default settings was not so satisfactory. As explained in the next section, replacement of model defaulted terrain height with the real value can improve the snow simulation ability in WRF_CLM. So the comparison of RST impact with WRF_CLM is based on the replacement. To evaluate the impacts of RST temperatures on snow simulation under different terrain heights, we performed sensitivity runs with different RST temperatures and different terrain heights. The details of CLM and WRF_CLM experiments are listed in Table 2.

4 Results
  • With default model settings (with 2.5℃ RST temperature), the CLM could capture the observed snow depth variation fairly well (Fig. 2a). However, it tended to underestimate the snow depth. In general, it can be used for the snow study. Hereafter, the experiment will be referred to as CLM-2.5 (Table 2).

    Figure 3.  Observed and simulated snow depths and 2-m-height air temperature with 2.5℃ RST temperature and default/changed terrain height.

    The simulation, using the WRF_CLM model with default model settings (with the input terrain height at about 20-km horizontal spatial resolution and 2.5℃ RST temperature), will be referred to as OH-2.5 in the rest of the script (Table 2). The maximum simulated snow depth was nearly 4 m in OH-2.5 (Fig. 3a), while the observation value was only 1.9 m. The simulated temperature was always lower than the observation value (Fig. 3c) and the simulated precipitation was slightly high in general (figure not shown). The inaccuracy in the results of experiment OH-2.5 may be attributed to incorrect input value of the terrain height for the station (1860 m was used in the model, whereas the actual value was 1320 m). Use of higher terrain heights led to a cooling of the surface atmosphere and lowering of the saturated vapor pressure, inducing more precipitation and increasing snow depth. This snow accumulation made the atmosphere even cooler, causing more precipitation and further snow accumulation. Thus, a positive feedback loop was formed. Owing to a low temperature, the snow melted quite slowly, thereby preventing a rise in temperature. As a consequence, a simulated snow depth of 3 m still existed in mid-May (Fig. 3a) and the simulated 2-m height temperature was around 0℃ (Fig. 3c). However, there was no observation value for snow and the observed maximum daily temperature was about 10℃.

    Application of the accurate topography could possibly improve the snow simulation (Jin and Wen, 2012). Hence, the simulation called CH-2.5 was performed, that is almost same with the OH-2.5 experiment, but we replaced the terrain height of the observation site in the model with the real value. From the root mean square deviation (RMSD) values and the bias in Table 3, it can be seen that the simulation of CH-2.5 was improved significantly. The CH-2.5 experiment captured snow depth variation fairly well, except that snow was a little underestimated and melted a little bit earlier (Figs. 2b and 3b). The simulated temperature was more realistic in the CH-2.5 experiment (Fig. 3d). In summary, the CH-2.5 experiment showed better simulation than the OH-2.5 one. The model was found to be suitable for adoption in similar research assignments if the real terrain height value could be adjusted well across the whole domain. Therefore, henceforth the study with WRF_CLM was performed with the real terrain height value, except for some sensitivity analyses.

  • The separation of precipitation into snow and rain in models was found to affect the simulation of water and energy balance notably. Presently, lots of RST temperature methods are used in numerical models. To find the differences in the impacts of these popularly used RST temperature methods on snow simulation, six different RST temperature methods (described in section 2) were applied to the offline land surface model CLM and the regional atmospheric model WRF_CLM for comparison. The corresponding experiments are listed in Table 2. The offline land surface model CLM with all these different RST temperature methods could simulate the snow depth fairly well (Fig. 2), as it employed observed forcing. The difference caused by RST temperature with the CLM model reflects primarily on the magnitude of the simulated snow depth. All methods in general overestimated the snow depth, except Dai's method and 0℃ threshold temperature when using the CLM model. The average maximum snow depth was simulated using Kienzle's method, followed by CLM-2.5. The finest simulation was achieved with Dai's method (Tables 4 and 5) and the worst was in case of CLM-0.

    Figure 4.  Observed and simulated snow depths with 3℃, 2℃, 1℃, and 0℃ RST temperatures by CLM and WRF_CLM.

    In the experiments using the regional atmospheric model WRF_CLM with the actual terrain height, the simulated snow depth with all different RST temperatures (Fig. 2) were underestimated and snow melted earlier. These results were different from those obtained using the offline CLM. The finest simulation was obtained in CH-2.5 experiment (Tables 4 and 5); however, snow disappeared 20 days earlier. This was followed by simulations with Pipes', USACE, and Kienzle's methods (Tables 4 and 5), although the values were quite close to each other; simulation of less snow caused its melting 33, 34, 34 days earlier, respectively.

    The accuracy of simulation with the offline land surface model CLM and regional atmospheric model WRF_CLM varied with the change of the RST temperature. The simulated snow had different responses to different RST temperatures with CLM and WRF_CLM models. This is possibly because RST temperatures partitioned precipitation in different forcing heights in offline land surface model CLM and regional atmospheric model WRF_CLM. The forcing heights in the land surface model were 2 m and 54 m in CLM and WRF_CLM, respectively. In addition, certain other factors influenced the simulation in regional atmospheric model WRF (such as spatial resolution, input terrain height, system error for simulating temperature, precipitation, radiation, and assumingly more).

    The difference in simulated snow depths between CH-2.5 and CH-0 was higher than that between CLM-2.5 and CLM-0 (Fig. 2). Additionally, the simulated snow depth in regional atmospheric model WRF_CLM was more sensitive to RST temperatures than that in offline land surface model CLM. This was because not only RST temperature was a key factor for partitioning the snowfall and rainfall in regional atmospheric model WRF_CLM, but also simulated snow related to RST temperature could induce the large feedback by changing the forcing meteorological variables, that is reversely always the same for different RST temperatures in offline land surface model CLM. In order to explain how sensitive the simulated snow was to the change in RST temperatures, sensitivity experiments were conducted with intervals of 1℃ RST temperature, using the regional atmospheric model WRF_CLM and offline land surface model CLM. Such an analysis would help reflect the effect of RST temperatures on snow simulation, in addition to the energy and water interaction between land and the atmosphere.

  • From the beginning of snow accumulation toward the end of its melting, the average observed snow depth was 0.87 m, whereas the average simulated snow depth in CLM-3, CLM-2, CLM-1, and CLM-0 were 1.22, 1.01, 0.89, and 0.48 m, respectively (Fig. 4a). Snow depths simulations of CLM-3, CLM-2, and CLM-1 were high as compared with the observation. The simulation in CLM-1 was the closest to the observed values perhaps because the actual RST temperature was around 1℃. Long-term observations of a meteorological station in Davos (1590 m, 45 km north of Morteratsch) showed mixed precipitation most appeared between 0.75℃ and 1.5℃, with a standard deviation of approximately 0.3℃-0.5℃ (Rohrer, 1989). On average, 50% rain (snow) fall occurred between 0.5℃ and 1℃ (L'Hote et al., 2005). The difference in snow depths simulated with offline land surface model CLM was relatively small between experiments with a 1℃ interval in the RST temperature range from 3℃ to 1℃. The difference in average snow depths was only 0.21 (0.12) m between RST temperatures 3℃ and 2℃ (2℃ and 1℃). However, the difference between experiments with 1℃ and 0℃ RST temperatures could be as large as 0.41 m. Discrepancies in the results of experiments using the same RST temperature intervals were possibly induced by the different distribution patterns of the observed temperature and precipitation forcing the offline land surface model CLM (Fig. 5). Nearly one-sixth of the accumulated precipitation in the snow accumulation period occurred between 0℃ and 1℃, which almost equaled the sum of precipitation falling between 1℃ and 2℃ and between 2℃ and 3℃. It is important to note that the simulated snow depths were not the same in the above experiments. However, their simulated time of accumulation and end of ablation data were closely related, as the forcing meteorological data were same.

    Figure 5.  Distribution of observed precipitation and temperature.

    The highest value of simulated snow depth with the offline land surface model CLM and regional atmospheric model WRF_CLM was obtained with the maximum value of RST temperature when the major portion of precipitation consisted of snow (Fig. 4). In general, the simulated snow depths in all experiments with WRF_CLM model were underestimated (Fig. 4b). The average simulated snow depths were approximately 0.68, 0.53, 0.14, and 0.05 m at RST temperature 3℃, 2℃, 1℃, and 0℃. The average simulated snow decreased by about 93% from CH-0 to CH-3, and by nearly 61% from CLM-3 to CLM-0. The simulated snow melting out time for the experiments employing the WRF_CLM also varied. Using the WRF_CLM, the earliest ablation was observed at the beginning of April and the latest toward the beginning of May. On the other hand, all the simulations by the CLM (using different RST temperatures) reflected melting out of snow during mid-May.

    Figure 6.  Simulated results by WRF_CLM with 3℃, 2℃, 1℃, and 0℃ RST temperatures.

  • With original input terrain height from model, the difference in average snow depths between OH-3 and OH-0 was about 0.31 m (Fig. 7a), which was low compared with the difference (0.63 m) between CH-3 and CH-0 (Fig. 4b). This reason for this immense difference could be attributed to the relatively low temperature simulated by WRF_CLM model with the original input terrain height (Fig. 7b), which was on average about -4.1℃ in both OH-3 and OH-0. The false original terrain height made the temperature fall below 0℃ during the simulated period. Therefore, with 0℃ and 3℃ RST temperatures, practically comparable separation of the precipitation into rain or snow took place. Conversely, after changing the input terrain height to the actual value, the simulated average 2-m height temperatures were 1.7℃ and -0.2℃ in CH-0 and CH-3, respectively (Fig. 6a). Referring to the above experiments, one can possibly conclude that the impact of RST temperatures on snow depth simulation was perhaps affected by other model parameters in addition to the terrain height.

5 Discussions and conclusions
  • In this study, two major objectives were addressed. The first was to validate the snow simulation ability of the recently coupled regional atmospheric model WRF_CLM and the second was to illustrate the impact of different RST temperature methods on snow depth simulation based on the offline land surface model CLM (forced by observation, not considering the interaction between atmosphere and land surface) and the regional atmospheric model WRF_CLM (considering the interaction between atmosphere and land surface). In addition, the difference in impacts estimated by the above two models was also deliberated. As the results reflected, both the CLM and the WRF_CLM demonstrated good simulation capacity for snow depth with the actual terrain height and the default 2.5℃ RST temperature. Snow depth simulated by the offline land surface model CLM using Dai's method was the closest to the observation, followed by that with USACE, Pipes', 2.5℃ RST temperature, Kienzle's, and 0℃ RST temperature methods. The snow depth simulated by the WRF_CLM with 2.5℃ RST temperature was the closest to the observation, followed by that with Pipes', USACE, Kienzle's, Dai's, and 0℃ RST temperature methods. The difference in performance of the above two models could be explained based on factors such as spatial resolution, input terrain height, system error for the simulation of temperature, precipitation, radiation, and more, which affected simulation of the regional atmospheric model WRF. In case of the offline land surface model CLM, the forcing data were fed with observed values. The difference in snow depth simulations with the same intervals of RST temperature by the regional atmospheric model WRF_CLM was much higher than that by the offline land surface model CLM, not only because the RST temperature decided the partition of snowfall and rainfall, but also because the simulated snow related to RST temperature could induce a significant feedback by changing the forcing meteorological variables in the regional atmospheric model WRF_CLM. Reversely, the forcing data were always the same for different RST temperatures in the offline land surface model CLM. The pattern of temperature and precipitation, along with other factors influencing the temperature simulation, could influence the effect of RST temperatures on snow simulation in the regional atmospheric model WRF_CLM. Therefore, a realistic treatment of model parameters, such as terrain height and spatial resolution, is highly recommended. In addition, representation of precipitation phase with suitable RST temperatures is vital for snow simulation in the models (such as the regional atmospheric model WRF_CLM), considering the interaction between atmosphere and land surface.

    In addition, the impact of different RST temperature methods on the snow depth simulation and difference between the impacts estimated by land surface and regional atmospheric models were also tested with data from the Boreal Ecosystem Atmosphere Study (Shewchuk, 1997) and Valdai, Russia (Schlosser et al., 1997). It was noted that the accuracy rate of simulated results was not consistent for different sites with different RST temperature methods. However, a broad conclusion can be drawn stating that simulated results with the regional atmospheric model were more sensitive to the changes in RST temperatures than those with the offline land surface model. Other factors (pattern of temperature and precipitation, spatial resolution, and input terrain heights) were shown to have considerable effects on simulation sensitivity of the regional atmospheric model, in addition to the changes of RST temperatures.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint