Meteorological Knowledge Useful for the Improvement of Snow Rain Separation in Surface Based Models

An accurate precipitation phase determination—i.e., solid versus liquid—is of paramount importance in a number of hydrological, ecological, safety and climatic applications. Precipitation phase can be determined by hydrological, meteorological or combined approaches. Meteorological approaches require atmospheric data that is not often utilized in the primarily surface based hydrological or ecological models. Many surface based models assign precipitation phase from surface temperature dependent snow fractions, which assume that atmospheric conditions acting on hydrometeors falling through the lower atmosphere are invariant. This ignores differences in phase change probability caused by air mass boundaries which can introduce a warm air layer over cold air leading to more atmospheric melt energy than expected for a given surface temperature, differences in snow grain-size or precipitation rate which increases the magnitude of latent heat exchange between the hydrometers and atmosphere required to melt the snow resulting in snow at warmer temperatures, or earth surface properties near a surface observation point heating or cooling a shallow layer of air allowing rain at cooler temperatures or snow at warmer OPEN ACCESS Hydrology 2015, 2 267 temperatures. These and other conditions can be observed or inferred from surface observations, and should therefore be used to improve precipitation phase determination in surface models.


Introduction
The hydrological, ecological, and atmospheric responses to rain are extremely different from snow responses [1] and correct precipitation phase is, therefore, very important for applications across many scientific fields.Truthful identification of the precipitation phase (rain/snow) is of course crucial for the functioning of meteorological models that forecast the precipitation phase itself [2] but also for accurate correction of gauge measured winter precipitation [3] and for land surface models (LSM) predicting snow accumulation and melt [4], glacier and polar ice water balance models [5], models for lake and sea ice growth [6], and climate change models [7].It is also important for models predicting avalanche hazards [8], sublimation of snow in forests [9], urban snowmelt quality [10], winter road safety [11], infiltration into frozen soils [12], survival of mammals and plants under snow cover [13], flooding from rain on snow events [14] etc. Precipitation phase determination is a modeling challenge for both hydrology and meteorology; therefore, a cross discipline approach combining methods and knowledge from both sciences could lead to new insight for both.
The International Association of Hydrological Sciences Decade on Prediction in Ungauged Basins (PUB) had the goal to improve basin hydrology through understanding processes rather than simply calibrating modeling schemes [15,16].In ungauged basins, knowledge of processes would create more robust hydrological solutions compared to the traditional approach to calibrate parameters in one basin, and apply them in others without adjustments for landscape heterogeneity [15].An example of a calibrated parameter used in many hydrological models is the precipitation phase determination scheme (PPDSs) usually based on surface temperature (the 2 m temperature reported from surface meteorological stations) relationships, ignoring important landscape and atmospheric interactions [17].For example, compared to dry landscapes, open water and unfrozen wetlands add heat and moisture to sub or near freezing air in the lowest part of the atmosphere introducing variability into average atmospheric conditions across landscapes.For this reason, ice-free oceans were found to have a warmer PPDS threshold temperature than continental landscapes [18,19].Therefore, it should be accepted that different PPDS schemes will be needed to account for lower atmospheric differences caused by variations in landscape.For example, dew point temperature based PPDS yielded less error than air temperature based PPDS in a saturated mountain environment [20], but an air temperature PPDS performed better in lower drier landscapes [21].
Within hydrological models, the PPDS are important [22], and many different schemes have been used.The most common approaches use surface air temperature (TA) e.g., [23] to determine the precipitation phase.Often TA and precipitation mass are the only measured variables with either lengthy enough historical records, or available at enough measuring stations for many studies requiring a PPDS [24].Unfortunately, TA methods are not physically based as they do not attempt to incorporate physical processes present in the real world e.g., the effects of humidity on latent heat transfer [25] and other hydrometeor processes above the soil surface [26].There are a few hydrological models that attempt to account for moisture in the surface-based PPDS by using surface dew point temperature (TD) [20], surface wet bulb temperature (TW) [27], and surface relative humidity (RH) [1].
All surface-based PPDS approaches assume a constant decrease in air temperature with height (e.g., the Climate High Resolution Model (CHRM)-7.5°C/km) [16] and do not account for interactions either between warm and cold air masses found in the most basic mid-latitude cyclone theories, or between hydrometeors and the air as precipitation falls through the lower atmosphere [28].Therefore, without atmospheric measurements, the most physically based surface PPDS methods can not accurately reproduce physical processes occurring between hydrometeors and the atmosphere they fall through.The phase of precipitation at the ground can be univocally determined only by knowing the initial hydrometeor phase [29], hydrometeor shape and size [26], and the properties of the atmosphere through which the hydrometeor falls.These properties include above ground air temperatures (TAG), height of the 0 °C isotherm (Z0 °C) [17], vertical thickness of freezing and melting layers [26], atmospheric humidity [1], atmospheric stability [30], and interactions between hydrometeors [2,25,29,31].
The equations used for describing the microphysical processes evaporation, sublimation, freezing, and melting [26,29] acting on the falling hydrometeors are calculation intensive and require detailed atmospheric information often not included in current hydrological models [29].
However, some hydrological modelers have begun to integrate meteorological measurements, or add steps that indirectly account for atmospheric changes into their PPDS.These methods use TA, TAG, weather radars, or satellite images.Fassnacht et al. [32] and Gjertsen and Ødegaard [33] e.g., improved precipitation mass by integrating radar precipitation estimates into their hydrological models.However, they chose to use surface TA methods rather than differences in the properties of the radar return pulse that could distinguish precipitation phase.This paper reviews and analyses studies dealing with precipitation phase processes and methods to determine the phase with an aim to suggest approaches that will improve the separation of rain and snow in cold region hydrological models.It starts with a basic overview of fundamental meteorological knowledge to build a baseline understanding of the physical changes between hydrometeors and air in the lower atmosphere.This first section builds to explain terrain/landscape, precipitation rate, and frontal boundary effects on the atmospheric temperature/moisture with height which has an impact on energy needed to drive phase changes in precipitation.This is followed by a short review of PPDS in atmospheric models and a longer review of current PPDS methods in surface models.The advantages and disadvantages of these PPDS methods are discussed along with a proposed formula to simplify surface air temperature-relative humidity PPDS schemes.

Atmospheric Interactions during Phase Change
In general, clouds and precipitation form when an air parcel rises and cools to saturation at the condensation level (the base of clouds).If snow falls through a warm atmospheric layer (>0.0 °C) melting will begin [25,32].Depending on air temperature, humidity and the thickness of a warm layer, latent and sensible heat exchanges between the snow and atmosphere could cause phase change from dry snow, to wet snow, to slush, to rain [26,32,34].
Hydrological models based on surface weather observations do not have enough information to construct a vertical temperature/moisture profile that accurately reflects atmospheric conditions above the surface.A precipitation phase reported from a point location could be the result of many different atmospheric processes e.g., Stewart [25] explained that melting snow at the surface can be a result of (1) the melting process; (2) heat gained from collision with rain drops; (3) refreezing after partial melt; or (4) collisions with supercooled water.The assumption of any of these processes extrapolated over a model surface could lead to large errors.Therefore, the use of a vertical temperature/moisture profile of the environment is preferred in meteorology because it can first determine which hydrometeors are most likely to form in clouds, and then calculate any precipitation phase changes that might occur between and below the cloud level at which the hydrometeor formed.

Lapse Rates
Sampling instruments carried by weather balloons known as radiosondes provide vertical profiles of above ground air pressure, air temperature, dew point, wind speed/direction, and GPS location.Spatial resolution of the observations depends on the number of balloons launched in a given area.Over 800 stations worldwide submit radiosonde data for model intake.The time resolution is twice a day [35] at most locations.The vertical temperature and moisture profile gives useful information; e.g., freezing level height, actual lapse rate compared to the dry and moist adiabatic lapse rates [17], (Figure 1), and the depth of melting and freezing layers.Another reason to use radiosonde data for this application is that TA varies more than TAG since the Earth's surface absorbs ultraviolet solar radiation and emits infrared radiation at different rates through the day.However, representative vertical temperature profiles may not be available at a necessary time, or for some areas the closest observation is hundreds of kilometers away and may not be representative for terrain or landscape induced variations in the lower atmosphere.

Deviations from Average Lapse Rates
The vertical temperature lapse rate under dry atmospheric conditions can be assumed to decrease steadily with height (Figure 2A).A similar assumption is implied through the use of set TA thresholds in hydrological models.These assumptions are incorrect when one or more of the following occurs in the near surface atmosphere; (a) the lapse rate is more unstable than average, decreasing the expected depth of a near surface warm layer (Figure 1); (b) there is an air mass boundary (Figure 2B,C); (c) falling snow is melting or refreezing (Figure 3); and (d) evaporation, sublimation, and/or condensation is occurring.

Air Mass Boundaries
Air mass boundaries introduce deviation from expected lapse rates due to a sharp increase or decrease in TAG (Figure 2B,C).For a given TA, this changes the probability of snow melting before reaching the ground.8 °C/km) and moist (5.9 °C/km) adiabatic lapse rates, give the expected cooling rates for rising air parcels without moisture or temperature advection.The dry and moist adiabatic lapse rates characterize the actual lapse rate as: (A) unstable when the actual TAG decrease is ≥ the dry adiabatic lapse rate (9.8 °C/km).Here, an air parcel that is warmed or lifted will be lighter than the surrounding air and rise whether it is moist or dry.(B) stable when the actual TAG decrease is ≤ the moist adiabatic lapse rate (5.9 °C/km).Here, an air parcel forced to higher altitudes (by terrain, or frontal lift) will be more dense than the surrounding atmosphere and will have a tendency to sink, even if condensation occurs; (C) conditionally unstable when the actual TAG decrease is between 5.9 and 9.8 °C/km.Here, an air parcel must be forced to higher altitudes (over terrain or a frontal slope) to reach its condensation temperature.From that altitude, the moist adiabatic lapse rate decreases slower than the actual lapse rate.Therefore, an air parcel can become lighter than the surrounding air and rise freely.

Isothermal Layers
Melting and refreezing of falling snow causes a near-freezing (+0 °C) isothermal layer (Dmelt) due to sensible and latent heat exchanges (Figure 3) [36].This isothermal layer provides limited energy for phase changes, reducing the total amount of melting or refreezing caused by atmosphere-hydrometeor energy exchanges.

Precipitation Intensity and Duration
The depth of Dmelt increases with precipitation intensity and duration (Figure 3) [17,33,36].Minder et al. [17] found that Dmelt for very light precipitation was about 60 m and could exceed 300 m for high precipitation rates.Gjertsen and Ødegaard's [33] found that the above ground TW lapse rate decreased with increasing precipitation rates from 0.6 to 0.1 °C/100 m.The maximum Dmelt was 598 m, and the minimum was 244 m.However, Minder et al. [17] claims that Dmelt is not an important factor since strong frontal winds usually restrict the average residency period of air in a melting layer to 10 min.
Snow fractions for a given TA can be expected to increase with precipitation intensity and storm duration (Figure 3). Figure 3 illustrates differences between typical low intensity and short duration storms (Figure 3A) and high intensity or long duration storms (Figure 3B).

Evaporation, Sublimation and Condensation
Melting causes a Dmelt to form in the vertical profile that is usually a few hundred meters deep [37,38].The thickness of Dmelt is dependent on snowflake diameter, relative humidity (RH), precipitation rate, and the atmospheric lapse rate [17].Mitra et al. [39] demonstrated that a 10 mm diameter snowflake would have a Dmelt 100 m thicker than a 5 mm snowflake during otherwise identical conditions.The sensible flux is dependent on the temperature difference between the hydrometeor and the ambient air while the latent heat flux is reliant upon the difference in vapor density (i.e., on RH) between the hydrometeor and the ambient air [20,29].Sublimation (2594 J/g) and evaporation (2260 J/g) consume heat energy from the atmosphere reducing available melt energy (334 J/g) while condensation (2260 J/g), ice condensation (2594 J/g), and freezing (334 J/g) release heat energy back into the atmosphere.Melting and refreezing may be combined with evaporation/sublimation, so that hydrometeor latent and sensible heat exchange is balanced by heat flux from the atmosphere.For moist air, i.e., air with high relative humidity, little evaporation/sublimation can be expected due to high vapor density surrounding the hydrometeor while more can be expected for dry air due to low vapor density around the hydrometeor [29].3A,B where initial storm conditions are the same but precipitation rate and duration differ.From bottom to top in the overview, surface temperatures increase to correspond with a precipitation phase change from snow to mixed precipitation to rain.The vertical profiles to the left in the subfigures illustrate the situation in location 2; (A) illustrates the situation for a low intensity or the beginning of a precipitation event, with a thin isothermal layer in the vertical temperature profile; (B) illustrates the situation for a prolonged or high precipitation intensity event.It has increased depth of the isothermal layer but minor change in surface temperature, a more unstable near surface lapse rate, decreasing melt and increasing surface snow areas.

Hydrometeor Interactions
If droplets freeze or ice crystals melt on a hydrometeor, they would release or absorb 334 J/g of latent heat [25].The melting rate of snowflakes thus increases when they collide with cloud or rain droplets [26].This is because a hydrometeor composed of both liquid and solid fractions has a fixed temperature of 0 °C [25].Snow colliding with rain or cloud droplets will undergo melting to chill the droplets to 0 °C [25].
The atmospheric physics of hydrometeor formation and interactions including phase change are very complicated.An understanding of the basic five equation microphysical schemes used by Stewart [25] is suggested.His equations account for: (1) Heat exchange resulting from gas, liquid, or solid phase change; (2) Change in particle mass by condensation; (3) The resulting latent heat flux from condensation; (4) Heat exchange from collision coalescence or through accretion of liquid and ice crystals; and (5) Latent heat of fusion with particles melting or freezing due to contact with each other.Detailed examples of more complicated microphysical schemes are found in e.g., Thériault and Stewart [26] or Lundquist [28].

Terrain Influence
Topography plays an important role in the surface location of rain/snow phase transitions.It is common for phase transition to be associated with a shallow cold layer beneath a warm layer on the lee side of mountains.This is known as cold air damming.This phenomenon can lower Z0 °C [17,25], and may add an extra freezing level in the vertical temperature profile (Figure 4).A phase transition region can sometimes be linked to coastlines due to warm sea surface temperatures.The warm water transfers heat to the near surface atmosphere increasing the available melt energy over the open water, also extending into the near coastal region.Strong winds can push this warm air further inland or offshore, therefore, also moving the location of a phase transition region [25].In most cases, a phase transition region will progress towards adjacent rain areas (Figure 3) [25].

Figure 4.
Phase changes of snow in a melting layer, followed by phase changes in a refreezing layer modified from [26].
Air forced to rise over terrain obstacles may also yield orographically enhanced precipitation.Increased precipitation rates cause a thicker Dmelt, lower Z0 °C and snow elevation (ZS) than upwind [17].This lowering of Z0 °C and ZS will introduce phase errors in a model that simply projects temperature heights from a sounding (or snow elevations from radar returns) without adjusting for variations in terrain.The magnitude of change in Z0 °C and ZS depends on the geometry of the terrain and the properties of the melting layer (including both RH and wind speed) [17].

Atmospheric Precipitation Phase Determination
Most meteorological PPDS approaches are based on measured atmospheric conditions.There are many approaches that range from complex bulk microphysical schemes to simplified empirical approaches using thickness values or freezing levels to determine phase.Snow, wet snow, and/or rain might form when TAG decreases with height, while, a TAG profile consisting of a warm layer (>0 °C) aloft and a cold layer below (<0 °C) usually is needed to form ice pellets, supercooled rain, and freezing rain (Figure 4) [26].

Bulk Microphysical Schemes
Meteorological models commonly use bulk microphysical schemes to determine the hydrometeor phase as precipitation falls to the surface.These schemes assume a size distribution for hydrometeors [26].A hydrometeor's size determines the amount of energy needed for phase change.Small hydrometeors changing phase faster than large hydrometeors is often the cause of mixed phase precipitation [26,29].
There are both sophisticated and simplified bulk microphysical schemes (Table 1).The sophisticated Purdue-Lin scheme incorporates six hydrometeor classes (water vapor, cloud water, cloud ice, rain, snow, and graupel) and mixed-phase processes [40].This scheme applies complex equations to account for evaporation, sublimation, deposition, condensation, aggregation, accretion, Bergeron processes, freezing, and melting.Other bulk microphysical schemes are simpler, e.g., the WSM3 (WFR [Weather Research and Forecasting] Single Moment-3) with three precipitation phases [41] that does not allow mixed phase processes, assuming an instantaneous phase changes to occur at freezing levels.Of note, most microphysical schemes do not account for drizzle formed through the warm rain process in low clouds, usually under inversions.These clouds contain little ice condensation nuclei causing supercooled drizzle without the melting and refreezing layers in Figure 4 [1,33].

Thickness Values
A thickness value is the difference in height (DZ) between two standard pressure levels (e.g., between 500 and 1000 mb) identified by Radiosonde data.DZ is proportional to the changes in pressure (DP) and density.Given a constant DP, warm air is less dense resulting in a larger DZ, while, cold air has the opposite properties (Figure 5).Critical thickness values e.g., 1300 m for 1000-850 mb, 2840 m for 1000-700 mb and 5400 m for 1000-500 mb [43] can be used as rain/snow thresholds.Snow is expected when a DZ is lower than the critical DZ value, and rain if a DZ is higher.Meteorology traditionally used critical DZ values [25,43] when less computational PPDS methods were necessary (Figure 5).The advantage of this method is its simplicity.It is less sensitive to diurnal change than TA but still has many of the same disadvantages.Specifically this method assumes a constant or average change in temperature with height (Figure 2A, Figure 5A,C) but cannot identify cold, warm, (Figure 2B,C and Figure 5B) or isothermal layers (Figure 3B) in the vertical temperature profile affecting precipitation phase [25].

Freezing Levels
The height of Z0 °C is available from radiosonde data and is also used to determine the probability of snow falling through a near surface warm layer (Table 2).Similar approaches exist for above ground TW and TD.On average, a warm layer must be at least 1200 ft (365 m) thick to completely melt all snow [43].This approach does not account for the following factors affecting Dmelt: (a) Actual temperature and wind speed of the warm layer; (b) a possible second warm layer; or (c) snowfall intensity and flake size.

Surface Based Precipitation Phase Determination
Many hydrological models rely solely on automated observing sensors (AOS) for weather data.Most PPDS apply the same empirical formula to all precipitation events reviewed in [21].Of note, one weakness for all surface based PPDS methods is that they are not physically sound, unable to account for microphysical processes occurring above the surface.

Time Step Dependence
Many PPDS schemes use average daily TA for a critical rain/snow threshold (TRS) or slight modifications of it.Ruddell et al. [45] and Schreider et al. [46] e.g., use daily minimum TA.In 1956, the U.S. Army Corps of Engineers [23] determined, from hourly observations, that TRS was normally between 1.1 °C and 1.7 °C while Saelthun [47] found a larger TRS range (−1 °C to 4 °C) using daily average TA.Harder and Pomeroy [29] showed that a daily TA had double the temperature range with 5% or more misclassified precipitation than an hourly time-step.A shorter time step decreases the uncertainty of actual temperature during precipitation caused by e.g., diurnal temperature change [30], and air mass boundaries [48].However, reducing the time step further to 15 min gave similar misclassified precipitation as one hour [29].

Regional Variance
Some studies suggest that TRS may be dependent on local terrestrial properties.Yang et al. [49] and Dai [18] noted changes in TRS due to regional variations in RH, air salinity, and elevation.Using a global 3 hourly synoptic network, Dai [18] estimated the average TRS to be 1.2 °C over land and 1.9 °C over oceans.Similarly, Feiccabrino [19] found a TRS of 2.0 °C over the North Sea, 1.5 °C over the Scandinavian Mountains, and 1.0 °C over most land.Motoyama [50] and Ólafsson and Haraldsdóttir [30] found greater TRS variance in their island studies (0 °C to 3 °C) than found in both broader, and inland studies.Also of note, Ólafsson and Haraldsdóttir, [30] and Kienzle [31] found similar seasonal oscillations for coastal, inland, and lee of mountain TRS Values.

Variations Caused by Threshold Determination Method
The differences in reported TRS values (Table 3) may partly be attributed to the way they have been determined.Feiccabrino et al. [21] and Ye et al. [1] compare measured (visual observation) and modeled precipitation phase, while many studies e.g., Wen et al. [2] compare modeled and measured snow depths.For the latter type, the representation of (a) wind correction of gauge catch; (b) spatial resolution of the watershed; (c) lapse rate; and (d) snow melt etc. in the model has an influence on the determined TRS value.A TRS is defined as the TA at which the likelihood of snow equals the likelihood of rain [24,30].Precipitation is assumed to fall as snow (rain) if the temperature is cooler (warmer) than TRS [2].
Despite knowledge of regional and seasonal variance, many models use a fixed TRS.However, some models e.g., the CHRM and the United States National Weather Service (NWS) snow accumulation and ablation models, have a default TRS that can be changed either for single events or be permanently adjusted [56,57].
The main advantage with TRS is that it is a quick simple method using minimal processing power, and the required parameters are widely available while other options may not be e.g., [24].The main disadvantage is that it does not take into account mixed phase precipitation that often occurs when TA approach TRS.Daly et al. [58] concluded that using a TRS was the main weakness in their study resulting in a lower than observed snow water equivalent by a systematic underestimation of the snow fraction (SF).

Dual Air Temperature Threshold Schemes
To account for the gradual decrease in SF occurring as TA approaches TRS, some models use a dual threshold method whereby all precipitation colder than a snow temperature threshold (TS) is modeled as snow and warmer than a rain temperature threshold (TR) is modeled as rain.The precipitation phase transition zone between TS and TR is considered mixed and proportional to a linear or curvelinear relationship between the defined thresholds [29].Some examples of models using linear dual threshold methods are: Community Atmospheric Model (CAM) Version 3.0 with TS −5 °C and TR 0 °C, The University of British Columbia (UBC) Watershed Model [31] with TS 0.6 °C and TR 3.6 °C.Curvelinear examples are: WATCLASS 2.7 [34], and CLASS 3.1 [4] with TS 0.45 °C and TR 5.97 °C based on the Auer [59] polynomial.
The dual threshold methods have been found to result in less misclassified precipitation than TRS in many comparison studies e.g., [2,21].Of the dual threshold methods, the approach applying a linear decrease in SF requires the least processing power.However, a linear decrease in SF should theoretically have more misclassified precipitation than a curvelinear TA-SF relationship [21].
All methods below will require additional processing power and might require variables not already found in a hydrological model.

Air Mass Boundary Scheme
Air mass boundaries which are known to change lower atmospheric conditions (Figure 2B,C) cannot be measured by AOS, but can be identified at the surface through temperature and wind changes when past, current, and future surface weather observations are compared [60][61][62][63][64]. Cold air mass boundaries were found to have TS and TR values 1 °C warmer than all other precipitation observations [48].When treated separately, total misclassified precipitation was reduced from 7.0% to 5.4% for the temperature range −1°C to 5 °C [48].

Temperature Schemes Including Humidity
PPDS can be based on different types of surface temperature.TA methods do not include humidity, the second-most significant factor in determining precipitation type [27,29,33].Some researchers have found TD [20], TW [27], and a combination of RH and TA [1] to be better predictors of the precipitation phase.However, the assumption that surface humidity is fully representative of the atmospheric humidity is incorrect.It is crucial to remember that the atmosphere above the surface through which a hydrometeor falls is a major factor in melting and refreezing.

Dew Point Temperature Schemes
A TD is the surface air temperature an air parcel would have to be cooled to in order to reach saturation at a constant pressure without moisture exchange [20].Saturation temperatures differ between ice and water, therefore, requiring different formulas for TD [29].
In near-saturated mountain environments, TA ≈ TD ≈ TW and TD was found to be more reliable for phase determination than TA [20].Marks et al. [20] hypothesized that in this moist environment, TD was a proxy for hydrometeor temperature and would therefore identify the liquid or solid phase.Concerning this environment, a TRS > 0 °C was found to misclassify more rain as snow than a threshold of TD = 0 °C.Notably, this method did not work well when there was warm air aloft or when cloud bases were well above the surface [20].
TD thresholds were found to have less spatial variations than TRS [1,21].No overall seasonal changes in TD thresholds were found [1].The TD method can be easily applied when data is supplied by an AOS.However, knowledge of which relationship used by the AOS to derive TD may be necessary.

Wet-Bulb Temperature Schemes
The TW is the temperature an air parcel would be cooled to by evaporation until saturation without exchange of sensible heat [20].A TW scheme approximates the effect of sublimation and evaporation on snowflakes in air warmer and drier than the cloud in which it was formed [20,33].TW is derived from TA, TD, air pressure, and saturation vapor pressure.Snowfall frequently occurs when TA > 0 °C and the near surface RH is < 100%, such that TA > TW > TD [20].Under these unsaturated (drier) conditions, TW may be more effective than TD in a PPDS [20].

Relative Humidity Schemes
There are several studies showing that transition from rain to snow seems to take place at warmer temperatures for dry than for moist air [1,25,27,29,33,50].Matsuo et al. [27] compared observed phase with surface TA and RH in Japan to find snow and rain RH thresholds.All precipitation occurring for conditions plotted below the snow RH threshold were assumed snow, and above the rain RH threshold were assumed rain.Precipitation between the two RH thresholds was considered mixed.They found slightly varying relationships for different locations, but for each 10% drop in RH the TRS would increase by about 1 °C (Figure 6A,B).Similar trends were observed in Gjertsen and Ødegaard's [33] study, where the probability of snow for a given TA is increased with decreasing RH (Figure 7).
Harder and Pomeroy [29] utilized knowledge from research on sublimation of blowing snow to estimate a hydrometeor temperature (Ti) including the effects of RH.This Ti, achieved by an iterative process, was shown to be a better indicator of precipitation phase than TRS.They used energy and mass balances of a sublimating ice sphere and assumed that the inner temperature of the ice sphere was the same as its surface temperature.They also assumed that the Nusselt and Sherwood numbers, indexes for the degree of turbulent transfer from a particle surface to the atmosphere of energy and water vapor respectively, were equal.Hydrometeor temperatures, calculated for different RH-values, using their relationships for the most interesting TA range (−5 °C to + 5 °C) are shown in Figure 8, and confirm that Ti increases with RH.An 80% increase in RH gives a Ti increase of ≈ 4.5 °C for TA ≈ −5 °C and ≈ 6.5 °C for TA ≈ +5 °C (Figure 8).The average Ti thus increases about 0.7 °C per 10% rise in RH in the relevant TA range −5 to +5 °C.The hydrometeor rain/snow threshold temperature Ti = 0 °C, from now on denoted TRH, thus occurs at 0.7 °C warmer TA for each 10% drop in RH (Figures 8 and 9).This agrees well with other reported TRH relationships with about 0.85 °C warmer TA for each 10% RH drop (Figure 9).The TRH values from Gjertsen and Ødegaard [33], were determined by setting TRH as the temperature having 50% SF and for Matsuo et al. [27] as the average of the reported rain and snow humidity thresholds (Figure 6A,B).Since the TRH relationships for TA, and RH from all three studies show similar almost linear relationships (Figure 9) having about 0.85 °C warmer TA per 10% RH decrease, the following simple relationship (for TRH between 0 °C and 5 °C) is suggested; This very simple equation gives TRS (40%) = 4.9 °C and TRS (90%) = 1.65 °C, thus similar results as Matsuo et al. [27], Gjertsen and Ødegaard [33], and Harder and Pomeroy [29] (Figure 9).However, a weakness with RH-based methods is that RH is very sensitive to small errors in TA, e.g., a 0.2 °C error leads to an RH error of 4% at TA = 0 °C [65].Accurate measurements of TA are thus needed.

Summary
The goal here was to increase basic meteorological knowledge of precipitation phase change processes and to point out or identify techniques to improve surface based PPDS.Most hydrological models use daily average surface temperature based empirical formulas to imply precipitation phase without regard for atmospheric processes.These methods like all other surface based methods do not account for variations in atmospheric conditions that cause actual vertical temperature profiles to be uniquely modified by e.g., melting and refreezing.
However, methods developed or modified through knowledge of lower atmospheric interactions between air and hydrometers e.g., relative humidity, air mass boundaries, or precipitation intensity, should allow for statistically better chances of correct precipitation phase determination in surface based models.As a bonus, some of the below methods can be implemented without much difficulty since they only need AOS data already available in hydrological modeling.
Surface based precipitation phase determination methods can be improved by: • Using a shorter a time-step: Harder and Pomeroy [29], showed that by decreasing the time step from one day to one hour resulted in more than a 50% decrease in the temperature range for which mixed-phase (rain and/or snow) precipitation occurred.
This was most likely caused by many precipitation events occurring on a sub-daily timescale (most common exception warm frontal light stratiform precipitation) which often ranges from tens of minutes (intense cumuliform precipitation) to several hours in duration [48].Since temperatures change diurnally with a maximum temperature in the afternoon and a minimum temperature in the morning, a decrease in the temperature sampling time step closer to the time scale to which precipitation events occur would result in a more representative temperature for precipitation phase determination.
• Introducing single (or dual rain/snow) thresholds which vary with relative humidity RH (%): Matsuo et al. [27], Gjertsen and Ødegaard [33], and Harder and Pomeroy [29] all show that TRS increases with decreasing RH.Based on these studies, it is shown that, for the surface air temperature range of 0 °C to 5 °C, TRH could be approximated by 0.75 + 0.085 × (100 − RH).
• Identification of air mass boundaries: Feiccabrino et al. [48] found TS and TR values for a group of cold air mass boundaries to be 1 °C warmer than that of all other observations.When these observation groups were analyzed separately, total misclassified precipitation was reduced by 23% in the temperature range −1 to 5 °C.
This was thought to be the cause of warm air rising over a lower denser cold air mass for most precipitation events with the main exception being cold fronts having precipitation occur ahead of the front in a single air mass.This theory worked in a Central United States continental climate [48] but did not improve precipitation phase identification in the maritime mountain environment of the Scandinavian Peninsula [19].In a maritime or island environment surface air in the winter is modified by the warm ocean temperatures making all but the strongest frontal passages advect warm air over the land.This ocean affect could be the reason for such different results between studies.

• Inclusion of precipitation intensity and/or duration:
No direct study available for confirmation, but Harder and Pomeroy [29] noted that the temperature range between TS and TR narrowed for heavy precipitation events and Minder et al. [17] showed an increase in the isothermal layer depth associated with melting precipitation.The presence of unstable air under the isothermal layer should allow for snow reaching the ground at warmer air temperatures similar to ocean effects from Ólafsson and Haraldsdóttir [30].
The lack of a study in this area would suggest that either negative results were found and therefore not published, or it has not been looked into.Considering that most studies are carried out using a daily or hourly precipitation and temperature time-step, it is possible that cooling of the temperature due to melting and sublimation occurred rapidly enough to keep warm high intensity snow events from being observed often enough to be a notable issue.The isothermal layer can thicken for long duration events (Section 2.2.3), but the surface temperature can also lower with sublimation, evaporation, or, to a lesser extent, modify to being closer to the precipitation temperature.
• Seasonal and/or terrain dependent lapse rates: Dai [18] reported a warm bias (+0.7 °C) for TRS over oceans compared to over land.This phenomena was confirmed in two other studies [19,30], and is thought to be caused by the ocean heating the near surface atmosphere resulting in higher lapse rates over the coast and water.
It is therefore suggested and shown [19] that different snow/rain temperature threshold values should be used for different geographic landscapes.This is due to ocean/water heating the near surface air causing the snow/rain threshold to be warmer over open water than land.The rain/snow threshold differences on land could be caused by terrain effects on cloud formation and precipitation intensity.
This paper was directed mostly towards techniques to improve point precipitation phase identification for possible quality assurance checks, reanalysis of historical data, or the assigning of a precipitation phase when there is no manual observation or automated device present to report a phase.Many of the automated methods/instruments that report precipitation phase use; (1) Doppler radar; (2) optical sensors [66]; (3) hotplates [67] or (4) two dimensional video distrometers (2DVD) [68].However, they have often been put into operation without an in-depth intercomparison which brings into question the quality, consistency, compatibility, and representativeness of their data [66].In response to this, the WMO (World Meteorological Organization) Solid Precipitation InterComparison (WMO-SPICE) study was created.The initial WMO-SPICE findings are just now becoming available e.g., [69,70].
It is also beyond the scope of this paper, however of great value to understand how landscapes such as forest or open fields affect the spatial distribution of snowfall, see e.g., [71].

Conclusions
Fixed rain/snow thresholds in surface based models tend to integrate average atmospheric conditions that control phase change when distinction could be useful.Distinct Atmospheric conditions and terrain features which cause warmer than normal rain/snow temperature thresholds are: • Cold frontal passage • Locations over or near large ice free water bodies • Low relative humidity • Long duration or high intensity precipitation events • Locations on the windward side of mountains (caused by increased precipitation intensity) The above surface and atmospheric weather conditions could all be identified using surface observations while the above terrain features that interact with lower atmospheric conditions can be easily identified on topographic maps.This information is readily available and should therefore be used to improve precipitation phase determination in surface based models.

Figure 1 .
Figure 1.Dry (9.8 °C/km) and moist (5.9 °C/km) adiabatic lapse rates, give the expected cooling rates for rising air parcels without moisture or temperature advection.The dry and moist adiabatic lapse rates characterize the actual lapse rate as: (A) unstable when the actual TAG decrease is ≥ the dry adiabatic lapse rate (9.8 °C/km).Here, an air parcel that is warmed or lifted will be lighter than the surrounding air and rise whether it is moist or dry.(B) stable when the actual TAG decrease is ≤ the moist adiabatic lapse rate (5.9 °C/km).Here, an air parcel forced to higher altitudes (by terrain, or frontal lift) will be more dense than the surrounding atmosphere and will have a tendency to sink, even if condensation occurs; (C) conditionally unstable when the actual TAG decrease is between 5.9 and 9.8 °C/km.Here, an air parcel must be forced to higher altitudes (over terrain or a frontal slope) to reach its condensation temperature.From that altitude, the moist adiabatic lapse rate decreases slower than the actual lapse rate.Therefore, an air parcel can become lighter than the surrounding air and rise freely.

Figure 2 .
Figure 2. Examples of vertical air temperature profiles: (A) a single air mass with a constant lapse rate; (B) an air mass boundary separating a warm over cold air mass; and (C) an air mass boundary separating a cold over warm air mass.

Figure 3 .
Figure 3. Overviews of flat areas are shown to the right in the two Figure3A,B where initial storm conditions are the same but precipitation rate and duration differ.From bottom to top in the overview, surface temperatures increase to correspond with a precipitation phase change from snow to mixed precipitation to rain.The vertical profiles to the left in the subfigures illustrate the situation in location 2; (A) illustrates the situation for a low intensity or the beginning of a precipitation event, with a thin isothermal layer in the vertical temperature profile; (B) illustrates the situation for a prolonged or high precipitation intensity event.It has increased depth of the isothermal layer but minor change in surface temperature, a more unstable near surface lapse rate, decreasing melt and increasing surface snow areas.

Figure 5 .
Figure 5. Example of thickness values (DZ) for 1000 to 850 mb (solid line), and for 1000 to 500 mb (hatched line) for air columns with (A) a cold air mass; (B) an air mass boundary and (C) a warm air mass.

Figure 6 .
Figure 6.Rain (A) and snow (B) thresholds as a function of RH determined from equations in [27].

Figure 7 .
Figure 7. Probability of rain versus surface air temperature for different RH (Figure 7 reprinted from Gjertsen and Ødegaard [33], with permission from Elsevier).

Table 1 .
[42]Model Version 3 Microphysical Schemes from most to least complex table modified from[42].

Table 2 .
Probability of snow based on the height of the freezing layer (TAG = 0 °C).

Table 3 .
Examples of TRS values along with applied region and method used to determine the value.(M) Determined from modeled snow depths, SWE, or snow duration, (DO) determined from direct observation of precipitation phase, and (IO) determined indirectly from e.g., albedo variation.