On the influence of microphysics parametrization on the rainfall rates in numerical models of clouds

Accuracy of formulas for growth by accretion and evaporation of rain in bulk parametrization of these processes for the case of light and moderate precipitation is investigated. It is done by comparison of results from two simple models: with bulk approach and with exact calculations of growth or evaporation of drops in each size bin separately. Growth by accretion is accurately represented in bulk parametrization but rain evaporation is overpredicted. Corrected formula for rain evaporation is suggested.     

Estimating the snow water equivalent from snow depth measurements in the Swiss Alps

The snow water equivalent (SWE) characterizes the hydrological significance of snow cover. However, measuring SWE is time-consuming, thus alternative methods of determining SWE may be useful. SWE can be calculated from snow depth if the bulk snow density is known. Thus, a reliable estimation method of snow densities could (a) potentially save a lot of effort by, at least partly, sampling snow depth instead of SWE, and would (b) allow snow hydrological evaluations, when only snow depth data are available. To generate a useful parameterization of the bulk density a large dataset was analyzed covering snow densities and depths measured biweekly over five decades at 37 sites throughout the Swiss Alps. Four factors were identified to affect the bulk snow density: season, snow depth, site altitude, and site location. These factors constitute a convenient set of input variables for a snow density model developed in this study. The accuracy of estimating SWE using our model is shown to be equivalent to the variability of repeated SWE measurements at one site. The technique may therefore allow a more efficient but indirect sampling of the SWE without necessarily affecting the data quality.     

Modeling the Phase Transition Associated with Melting Snow in a 1D Kinematic Framework: Sensitivity to the Microphysics

A simple 1D kinematic cloud model coupled to a two-moment bulk microphysics scheme is used to perform quasi-idealized simulations of snow, with a prescribed upper boundary snow field based on observed radar reflectivity and temperature, falling into a low-level melting layer. The model realistically simulates the formation of a nearly isothermal layer below the melting level, the surface precipitation rate, and the phase transition from liquid to solid, consistent with observations for this case. A series of test runs is performed to examine the sensitivity of modeling the timing and duration of the phase transition period to details of specific parameterization aspects related to snow in the microphysics scheme. The sensitivity tests include varying the number of prognostic moments, the mass–diameter relation, the fall velocity–diameter relation, the treatment of aggregation, and the lower limit for the slope of the size distribution. It is shown that the simulated transition period, for such a case with the initial melting level being close to the surface, can be quite sensitive to model parameters specified within realistic ranges and/or ranges within our physical understanding.     

Modelling the poroelasticity of rocks and ice

For many geophysical and astrophysical applications the relationship between physical bulk properties and porosity is needed. As part of the preparation for the 'Rosetta' comet rendezvous mission, a simple model for the porosity dependence of the elastic properties of granular media, i.e. the elastic moduli and the propagation velocities of elastic waves, has been developed based on textural properties and the contact stiffness of the constituent particles. It is shown that the derived relationships fit very well with sandstone data. The model is also consistent with data for snow and ice and is in agreement with the transformation mechanisms from snow to ice. A short review shows the relevance to another physical bulk property, the thermal conductivity.     

SOME ASPECTS OF THE INTERACTION OF BLOWING SNOW WITH THE ATMOSPHERIC BOUNDARY LAYER

Several possible effects of blowing snow on the atmospheric boundary layer are investigated, mostly within the general framework of the Prairie Blowing Snow Model (PBSM). The processes of snow saltation and suspension are first described. Variations to the drift density profile are tested and the effects of stratification and density variation calculations are evaluated. Despite high density gradients of blowing snow, stratification effects on turbulence and the velocity profiles can generally be neglected. However, with saltating or suspended snow in a constant shear stress layer, part of the shear stress is carried by the particles. A highly simplified, single-phase approach, based on the density variation of the air–snow mixture coupled to a simple turbulent stress–strain relationship, is used to illustrate this. Sublimation rates in a column of blowing snow are calculated using the PBSM and results are compared with those obtained with a modified formulation which incorporates a spectrum of sublimating particles of varying sizes at each height in a steady-state surface boundary layer and different specifications of the ventilation velocity.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

Charred forests accelerate snow albedo decay: parameterizing the post‐fire radiative forcing on snow for three years following fire

As large, high‐severity forest fires increase and snowpacks become more vulnerable to climate change across the western USA, it is important to understand post‐fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post‐fire radiative forcing effects on snow to improve hydrologic modelling of snow‐dominated watersheds having experienced severe forest fires. Following a 2011 high‐severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow, and micrometeorological conditions, sampled snow surface debris, and modelled snowpack energy and mass balance in adjacent burned forest (BF) and unburned forest sites. For three winters following the fire, charred debris in the BF reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared with the unburned forest. We demonstrate a new parameterization of post‐fire snow albedo as a function of days‐since‐snowfall and net snowpack energy balance using an empirically based exponential decay function. Incorporating our new post‐fire snow albedo decay parameterization in a spatially distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the BF, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. Copyright © 2016 John Wiley & Sons, Ltd.     

The nutrient cycle through snow and ice, a review

This paper reviews the merging of the nutrient cycle with the water cycle in the seasonal alpine snow cover, emphasizing physical processes at the snowpack and snow grain scale. Nutrients are incorporated into snowflakes growing in the atmosphere, they are part of the dry deposition from the atmosphere to the snowpack and they reach the snow as plant litter. The physical processes of the accumulation of nutrients and their redistribution in and on the snow grains and in the pore space of the snow matrix are described.¶The first flush of meltwater that reaches the soil carries a solution of nutrients and acids in a concentration several times higher than bulk values, an effect that increases with the age of the snow and the number of melt/freeze cycles and is more pronounced for sulfate than for chloride. Species that are attached to insoluble particles will be concentrated near the snow surface and will display peak concentrations in the final fraction of meltwater.     

Measurements of powder snow avalanche- nature -

Snow cloud growth rates of powder snow avalanche were obtained with analysing the pictures recorded in Ryggfonn, Norway. Although results showed wide scattering, as far as the data obtained in runout zone are concerned, they roughly agreed with the water tank experiments by Beghin and Olagne (1991).Air movement in snow avalanche cloud was measured with an ultra-sonic anemometer in Kurobe Canyon, Japan. It showed rising current existed near the front and downward at the trail. Comparing with drifting snow threshold and particle suspension criterion, entrainment, suspension and deposition of snow particles in the snow cloud were discussed.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

A COMPARISON OF EVAPORATION FROM SNOW AND SOIL SURFACES

ABSTRACT

During the spring of 1961, evaporation from snow and soil surfaces was measured in the central Rocky Mountains near Fraser, Colorado. Measurements were made in natural forest openings at 9,000 feet elevation. Evaporation from wet soil surfaces greatly exceeded evaporation from nearby snow. There was little evidence of transfer of vapor from soil to nearby patches of snow, but as areas of bare, wet soil increased and evaporation amounts from such surfaces increased, evaporation from snow decreased. It was concluded that, as greater amounts of water evaporated from soils, the vapor pressure of the air was raised sufficiently to reduce evaporation from snow. Since transfer of vapor from soil to snow appeared small at best, evaporation losses from snow and soil surfaces essentially constituted a total moisture loss from the area.     

Estimation of snow ablation under a dust layer covering a wide range of albedo

Ablation processes of snow under a thin dust cover are complicated compared with those under a thick cover, mainly owing to the effects of aggregation (redistribution) of dust particles on the conditions of surface melting. Aggregation of dust particles causes the snow surface to brighten after the initial dust configuration, thus affecting the relationship between initial dust concentration and surface albedo. In order to estimate snow ablation rate under a thin dust cover, we used a composite energy balance model in which the surface albedo is taken as a measured input variable. The estimated results of snow ablation agreed reasonably well with the observation, considering the measurement errors inherited in the snow depressions. Comparison of the two cases, that is, one considering the aggregation of dust particles (observation: albedo variable) and the other without aggregation (assumption: albedo constant), showed that the ablation rates were noticeably lower on the former case. This suggests that the aggregation of dust particles induces a reduction of snow ablation. Copyright © 2002 John Wiley & Sons, Ltd.     

Observation and estimation of evaporation from the ground surface of the cryosphere in eastern Asia

The characteristics of evaporation from the ground surface of Asian cryosphere sites are presented, as estimated by the lysimeter method, a profile method, and a heat budget method. The observation sites were located on the eastern Tibetan Plateau, in the Qilian and Tianshan Mountains of China, and in eastern Siberia. The lysimeter method has been demonstrated to be a reliable observation technique for estimating daily evaporation from the land surface, given suitable experiment design and operation. Daily mean evaporation varied within the range of 0·3 to 3·5 mm on the permafrost surface, and regional differences in evaporation were strongly related to surface soil moisture. Locally, topography, by way of its influence on surface soil moisture, was found to control evaporation systematically. Seasonality of ground evaporation in permafrost regions is dominated by thaw–freeze cycles at the surface; evaporation from the melting permafrost surface is up to four to seven times greater than that from frozen ground. In forested terrain, the interception of precipitation can reduce daily evaporation by 60 to 70%. Sublimation from the snow surface was observed at some sites in the range of 0·2 to 1·0 mm daily; atmospheric conditions, such as wind speed and saturation deficit, were dominant factors in determining snow sublimation. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of spatial variability in snow water equivalent for a high mountain catchment

Multivariate statistical analysis was used to explore relationships between catchment topography and spatial variability in snow accumulation and melt processes in a small headwater catchment in the Spanish Pyrenees. Manual surveys of snow depth and density provided information on the spatial distribution of snow water equivalent (SWE) and its depletion over the course of the 1997 and 1998 melt seasons. A number of indices expressing the topographic control on snow processes were extracted from a detailed digital elevation model of the catchment. Bivariate screening was used to assess the relative importance of these topographic indices in controlling snow accumulation at the start of the melt season, average melt rates and the timing of snow disappearance. This suggested that topographic controls on the redistribution of snow by wind are the most important influence on snow distribution at the start of the melt season. Furthermore, it appeared that spatial patterns of snow disappearance were largely determined by the distribution of snow water equivalent (SWE) at the start of the melt season, rather than by spatial variability in melt rates during the melt season. Binary regression tree models relating snow depth and disappearance date to terrain indices were then constructed. These explained 70–80% of the variance in the observed data. As well as providing insights into the influence of topography on snow processes, it is suggested that the techniques presented herein could be used in the parameterization of distributed snowmelt models, or in the design of efficient stratified snow surveys. Copyright © 2003 John Wiley & Sons, Ltd.     

Intercomparison of measurements of bulk snow density and water equivalent of snow cover with snow core samplers: Instrumental bias and variability induced by observers

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra-taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.     

Surface-active substances on rainwater and atmospheric particles

The surface tension lowering by surface-active substances has been measured on rainwater, melted snow, and dispersions of atmospheric particles in water, with a film balance and a tensiometer. The precipitation water was sampled during 1979, 1980, and 1981 in the city of Frankfurt/Main. From measurements with the film balance technique, normalized concentrations of insoluble and weakly soluble surface-active substances have been estimated. Soluble surface-active substances were determined from measurements with a tensiometer. It was found that the normalized concentration of theinsoluble and weakly soluble surface-active material on rainwater or melted snow shows a maximum during late spring of about 2.5 · 10⁻⁷ moles/l and a minimum during wintertime of about 5 · 10⁻⁸ moles/l. These concentrations are too low to influence significantly the condensation of water vapour on cloud droplets or the evaporation of water from them. Thesoluble surface-active material on rainwater or melted snow was found to have concentrations of the order of 2 · 10⁻⁶ moles/l. These concentrations are also too small to have a significant influence on cloud physical processes.     

Impact of spatial variations of land surface parameters on regional evaporation: a case study with remote sensing data

Most precipitation in watersheds is consumed by evaporation, thus techniques to appraise regional evaporation are important to assess the availability of water resources. Many algorithms to estimate evaporation from remotely sensed spectral data have been developed in the recent past. In addition to differences in the physical parameterization of surface fluxes, these algorithms have different solutions for describing spatial variations of the parameters in the soil–vegetation–atmosphere–transfer (SVAT) continuum. In this study, the necessity to spatially distinguish SVAT parameters for computing surface heat fluxes is analysed for the Naivasha watershed in the Kenyan Rift Valley. Landsat Thematic Mapper (TM) spectral data have been used to first delineate the watershed into 15 hydrological units using surface temperature, normalized difference vegetation index and surface albedo as attributes. Thereafter, semi‐empirical relationships between these TM‐based parameters and other SVAT parameters have been applied to compute the spatial variation of SVAT parameters and the associated evaporation from the different hydrological units. The impact of using watershed‐constant or watershed‐distributed SVAT parameters on the fluxes is analysed. The determination of watershed averaged evaporation with area‐aggregated SVAT parameters is feasible without significant loss of accuracy. Distributed evaporation in heterogeneous watersheds, however, can be investigated only with remote sensing flux algorithms that can account for spatially variable air temperature, surface roughness, surface albedo and the stability correction of the temperature profile due to buoyancy. Erroneous results can be expected if area‐aggregated SVAT parameters are used to calculate local evaporation. As most of the recently developed remote sensing flux algorithms are based on areal constant SVAT parameters, direct applications in watersheds are still limited. Copyright © 2001 John Wiley & Sons, Ltd.