Michael Haugeneder

In absence of the high-frequency measurements of wind components, sonic temperature and water vapour required by the eddy covariance (EC) method, Monin–Obukhov similarity theory (MOST) is often used to calculate heat fluxes. However, MOST requires assumptions of stability corrections and roughness lengths. In most environments and weather situations, roughness length and stability corrections have high uncertainty. Here, we revisit the modified Bowen-ratio method, which we call C-method, to calculate the latent heat flux over snow. In the absence of high-frequency water vapour measurements, we use sonic anemometer data, which have become much more standard. This method uses the exchange coefficient for sensible heat flux to estimate latent-heat flux. Theory predicts the two exchange coefficients to be equal and the method avoids assuming roughness lengths and stability corrections. We apply this method to two datasets from high mountain (Alps) and polar (Antarctica) environments and compare it with MOST and the three-layer model (3LM). We show that roughness length has a great impact on heat fluxes calculated using MOST and that different calculation methods over snow lead to very different results. Instead, the 3LM leads to good results, in part due to the fact that it avoids roughness length assumptions to calculate heat fluxes. The C-method presented performs overall better or comparable to established MOST with different stability corrections and provides results comparable to the direct EC method. An application of this method is provided for a new station installed in the Pamir mountains.

The near-surface boundary layer over patchy snow is highly heterogeneous and dynamic. Layers of opposing stability coexist within only a few horizontal meters. Conventional experimental methods to investigate this layer suffer from limitations related to the fixed positions of eddy covariance sensors. To overcome these difficulties, we set up a centimeter-resolution large eddy simulation of flow across an idealised transition from bare ground to snow. We force the simulation with high-frequency eddy covariance data recorded during a field campaign. We show that the model can represent the real flow by comparing it to independent eddy covariance data. However, the simulation underestimates vertical wind speed fluctuations, especially at high frequencies. Sensitivity analyses show that this is influenced by grid resolution and surface roughness representation but not much by subgrid-scale parameterization. Nevertheless, the model can reproduce the experimentally observed plumes of warm air intermittently detaching from bare ground and being advected over snow. This process is highly dynamic, with time scales of only a few seconds. We can show that the growth of a stable internal boundary layer adjacent to the snow surface can be approximated by a power law. With low wind speeds, deeper stable layers develop, while strong wind speeds limit the growth. Even close to the surface, the buoyancy fluxes are heterogeneous and driven by terrain variations, which also induce the frequent decoupling of a thin layer adjacent to the snow surface. Our simulations point the path towards generalizing point-based and aerial measurements to three dimensions.

The near-surface boundary layer above patchy snow cover in mountainous terrain is characterized by a highly complex interplay of various flows on multiple scales. In this study, we present data from a comprehensive field campaign that cover a period of 21 days of the ablation season in an alpine valley, from continuous snow cover until complete melt out. We recorded near-surface eddy covariance data at different heights and investigated spectral decompositions. The topographic setting led to the categorisation of flows into up and down valley flows, with a down valley Föhn event in the middle of the observation period. Our findings reveal that the snow cover fraction is a major driver for the structure and dynamics of the atmospheric layer adjacent to the snow surface. With bare ground emerging, stable internal boundary layers (SIBL) developed over the snow. As the snow coverage decreased, the depth of the SIBL decreased below 1 m and spectra of air temperature variance showed a transition towards turbulent time scales, which were caused by the intermittent advection of shallow plumes of warm air over the snow surface. The intermittent advection could also be observed visually with high spatio-temporal resolution measurements using a thermal infrared camera. While the shallow advection only affected the lowest measurement level at 0.3 m, the measurements above at 1 m, 2 m, and 3 m indicate that the distribution of eddy size and, thus, the turbulence structure, did not distinctly change with height.

Snow is a vital part of water resources, and sublimation may remove 10%–90% of snowfall from the system. To improve our understanding of the physics that govern sublimation rates, as well as how those rates might change with the climate, we deployed an array of four towers with over 100 instruments from NCAR’s Integrated Surface Flux System from November 2022 to June 2023 in the East River watershed, Colorado, in conjunction with the U.S. Department of Energy’s Surface Atmosphere Integrated Field Laboratory (SAIL) and the National Oceanic and Atmospheric Administration (NOAA)’s Study of Precipitation, the Lower Atmosphere and Surface for Hydrometeorology (SPLASH) campaigns. Mass balance observations, snow pits, particle flux sensors, and terrestrial lidar scans of the evolving snowfield demonstrated how blowing snow influences sublimation rates, which we quantified with latent heat fluxes measured by eddy-covariance systems at heights 1–20 m above the snow surface. Detailed temperature profiles at finer resolutions highlighted the role of the stable boundary layer. Four-stream radiometers indicated the important role of changing albedo in the energy balance and its relationship to water vapor losses. Collectively, these observations span scales from seconds to seasons, from boundary layer turbulence to valley circulation to mesoscale meteorology. We describe the field campaign, highlights in the observations, and outreach and education products we are creating to facilitate cross-disciplinary dialogue and convey relevant findings to those seeking to better understand Colorado River snow and streamflow.

SIGNIFICANCE STATEMENT: Snow provides over 80% of water for the overallocated Colorado River, and in recent years, less runoff has occurred per unit snowfall. Sublimation, the conversion of ice to water vapor, results in less water for runoff, but due to a historic lack of observations, this process is hard to constrain. Variations in how sublimation is represented in models have led to a large divergence of projected water resource availability for the Colorado River basin over both current and future climates. The field campaign described here provides the first comprehensive examination of how snow accumulates, blows around, evolves, and sublimates over 8 months in the Colorado Rocky Mountains, providing a critical benchmark for process understanding and model development.

Dynamic downscaling of atmospheric forcing data to the hectometer resolution has shown increases in accuracy for landsurface models, but at great computational cost. Here we present a validation of a novel intermediate complexity atmospheric model, HICAR, developed for hectometer scale applications. HICAR can run more than 500x faster than conventional atmospheric models, while containing many of the same physics parameterizations. Station measurements of air temperature, wind speed, and radiation, in combination with data from a scanning Doppler wind LiDAR, are compared to 50 m resolution HICAR output during late spring. We examine the model’s performance over bare ground and melting snow. The model shows a smaller root mean squared error in 2 m air temperature than the driving model, and approximates the 3D flow features present around ridges and along slopes. Timing and magnitude of changes in shortwave and longwave radiation also show agreement with measurements. Nocturnal cooling during clear nights is overestimated at the snow covered site. Additionally, the thermal wind parameterization employed by the model typically produces excessively strong surface winds, driven in part by this excessive nocturnal cooling over snow. These findings highlight the utility of HICAR as a tool for dynamically downscaling forcing datasets, and expose the need for improvements to the snow model used in HICAR.

Mountain snowpack forecasting relies on accurate mass and energy input information in relation to the snowpack. For this reason, coupled snow–atmosphere models, which downscale input fields to the snow model using atmospheric physics, have been developed. These coupled models are often limited in the spatial and temporal extents of their use by computational constraints. In addressing this challenge, we introduce HICARsnow, an intermediate-complexity coupled snow–atmosphere model. HICARsnow couples two physics-based models of intermediate complexity to enable basin-scale snow and atmospheric modeling at seasonal timescales. To showcase the efficacy and capability of HICARsnow, we present results from its application to a high-elevation basin in the Swiss Alps. The simulated snow depth is compared throughout the snow season to aerial lidar data. The model shows reasonable agreement with observations from peak accumulation through late-season melt-out, representing areas of high snow accumulation due to redistribution processes, as well as melt patterns caused by interactions between radiation and topography. HICARsnow is also found to resolve preferential deposition, with model outputs suggesting that parameterizations of the process using surface wind fields may only be inappropriate under certain atmospheric conditions. The two-way coupled model also improves surface air temperatures over late-season snow, demonstrating added value for the atmospheric model as well. Differences between observations and model outputs during the accumulation season indicate a poor representation of redistribution processes away from exposed ridges and steep terrain and a low bias in albedo at high elevations during the ablation season. Overall, HICARsnow shows great promise for applications in operational snow forecasting and in studying the representation of snow accumulation and ablation processes.

The lateral transport of heat above abrupt (sub-)metre-scale steps in land surface temperature influences the local surface energy balance. We present a novel experimental method to investigate the stratification and dynamics of the near-surface atmospheric layer over a heterogeneous land surface. Using a high-resolution thermal infrared camera pointing at synthetic screens, a 30Hz sequence of frames is recorded. The screens are deployed upright and horizontally aligned with the prevailing wind direction. The screen’s surface temperature serves as a proxy for the local air temperature. We developed a method to estimate near-surface two-dimensional wind fields at centimetre resolution from tracking the air temperature pattern on the screens. Wind field estimations are validated with near-surface three-dimensional short-path ultrasonic data. To demonstrate the capabilities of the screen method, we present results from a comprehensive field campaign at an alpine research site during patchy snow cover conditions. The measurements reveal an extremely heterogeneous near-surface atmospheric layer. Vertical profiles of horizontal and vertical wind reflect multiple layers of different static stability within 2m above the surface. A dynamic, thin stable internal boundary layer (SIBL) develops above the leading edge of snow patches protecting the snow surface from warmer air above. During pronounced gusts, the warm air from aloft entrains into the SIBL and reaches down to the snow surface adding energy to the snow pack. Measured vertical turbulent sensible heat fluxes are shown to be consistent with air temperature and wind profiles obtained using the screen method and confirm its capabilities to investigate complex in situ near-surface heat exchange processes.

High-resolution (< 1 km) atmospheric modeling is increasingly used to study precipitation distributions in complex terrain and cryosphere–atmospheric processes. While this approach has yielded insightful results, studies over annual timescales or at the spatial extents of watersheds remain unrealistic due to the computational costs of running most atmospheric models. In this paper we introduce a high-resolution variant of the Intermediate Complexity Atmospheric Research (ICAR) model, HICAR. We detail the model development that enabled HICAR simulations at the hectometer scale, including changes to the advection scheme and the wind solver. The latter uses near-surface terrain parameters which allow HICAR to simulate complex topographic flow features. These model improvements clearly influence precipitation distributions at the ridge scale (50 m), suggesting that HICAR can approximate processes dependent on particle–flow interactions such as preferential deposition. A 250 m HICAR simulation over most of the Swiss Alps also shows monthly precipitation patterns similar to two different gridded precipitation products which assimilate available observations. Benchmarking runs show that HICAR uses 594 times fewer computational resources than the Weather Research and Forecasting (WRF) atmospheric model. This gain in efficiency makes dynamic downscaling accessible to ecohydrological research, where downscaled data are often required at hectometer resolution for whole basins at seasonal timescales. These results motivate further development of HICAR, including refinement of parameterizations used in the wind solver and coupling of the model with an intermediate-complexity snow model.