When applied to a snow-covered surface, aerodynamic roughness length, z0, is typically considered as a static parameter within energy balance equations. However, field observations show that z0 changes spatially and temporally, and thus z0 incorporated as a dynamic parameter may greatly improve models. To evaluate methods for characterizing snow surface roughness, we compared concurrent estimates of z0 based on (1) terrestrial light detection and ranging derived surface geometry of the snowpack surface (geometric, z0G) and (2) vertical wind profile measurements (anemometric, z0A). The value of z0G was computed from Lettau’s equation and underestimated z0A but compared well when scaled by a factor of 2.34. The Counihan method for computing z0G was found to be unsuitable for estimating z0 on a snow surface. During snowpack accumulation in early winter, z0 varied as a function of the snow-covered area (SCA). Our results show that as the SCA increases, z0 decreases, indicating there is a topographic influence on this relation.

Historically, snowpack trends have been assessed using one fixed date to represent peak snow accumulation prior to the onset of melt. Subsequent trend analyses have considered the peak snow water equivalent (SWE), but the date of peak SWE can vary by several months due to inter-annual variability in snow accumulation and melt patterns. A 2018 assessment evaluated monthly SWE trends. However, since the month is a societal construct, this current work examines daily trends in SWE, cumulative precipitation, and temperature. The method was applied to 13 snow telemetry stations in Northern Colorado, USA for the period from 1981 to 2018. Temperature trends were consistent among all the stations; warming trends occurred 63% of the time from 1 October through 24 May, with the trends oscillating from warming to cooling over about a 10-day period. From 25 May to 30 September, a similar oscillation was observed, but warming trends occurred 86% of the time. SWE and precipitation trends illustrate temporal patterns that are scaled based on location. Specifically, lower elevations stations are tending to record more snowfall while higher elevation stations are recording less. The largest SWE, cumulative precipitation, and temperature trends were +30 to −70 mm/decade, +30 to −30 mm/decade, and +4 to −2.8 °C/decade, respectively. Trends were statistically significance an average of 25.8, 4.5, and 29.4% of the days for SWE, cumulative precipitation, and temperature, respectively. The trend in precipitation as snow ranged from +/−2%/decade, but was not significant at any station.

The snow surface is at the interface between the atmosphere and Earth. The surface of the snowpack changes due to its interaction with precipitation, wind, humidity, short- and long-wave radiation, underlying terrain characteristics, and land cover. These connections create a dynamic snow surface that impacts the energy and mass balance of the snowpack, blowing snow potential, and other snowpack processes. Despite this, the snow surface is generally considered a constant parameter in many Earth system models. Data from the National Aeronautics and Space Administration (NASA) Cold Land Processes Experiment (CLPX) collected in 2002 and 2003 across northern Colorado were used to investigate the spatial and temporal variability of snow surface roughness. The random roughness (RR) and fractal dimension (D) metrics used in this investigation are well correlated. However, roughness is not correlated across scales, computed here from snow roughness boards at a millimeter resolution and airborne lidar at a meter resolution. Process scale differences were found based on land cover at each of the two measurement scales, as appraised through measurements in the forest and alpine.