Crop–water interactions define productivity in water-limited dryland agricultural production systems in cold regions. Despite the agronomic and economic importance of this relationship there are challenges in quantifying crop water use efficiency (WUE). To understand dynamics driving crop water use and agricultural productivity in these environments, observations of evapotranspiration, carbon assimilation, meteorology, and crop growth were collected over 17 site-years at 5 agricultural sites in the sub-humid continental Canadian Prairies. Eddy-covariance (EC) derived water and carbon fluxes provided a means to comprehensively assess the WUE of current agricultural practices by both physiological (WUEP: g C kg−1 H2O) and agronomic (WUEY): kg yield mm H2O−1 hectare−1) approaches. Mean field scale WUEY for grain yields were 10.4 (Barley), 10.2 (Wheat), 6.0 (Canola), 19.3 (Peas), 12.2 (Lentils) and for silage/forage crops were 23.0 (Barley), 11.9 (Forage), and 20.7 (Corn) (kg yield mm H2O−1 hectare−1). An assessment of environmental factors and their covariance with WUE, utilising a conditional inference tree approach, demonstrated that WUE decreased when crops were under greater evapotranspiration demands. EC-based areal WUE approaches, measuring fluxes over footprints of hundreds of square metres, were compared with more commonly reported point-scale water balance residual approaches (WUEWB) and demonstrated consistently smaller magnitudes. WUEWB was greater than EC-estimated WUEY by an average of 52% and 65% for grain and forage/silage crops respectively. WUEWB also had greater variability than EC estimates, with standard deviations 188% and 128% greater than Barley and Wheat crops, respectively. This comparison highlights the scale dependency of WUE estimation methods, demonstrates considerable uncertainty in point scale water balance approaches due to spatial variability in crop–water interactions, and shows how this variability can be accounted for by EC observations. This improves the understanding of WUE and quantifies its variability in cold continental water-limited climates and provides a means to diagnose improved agricultural water management.

Abstract. Vegetation has a tremendous influence on snow processes and snowpack dynamics, yet remote sensing techniques to resolve the spatial variability of sub-canopy snow depth are not always available and are difficult from space-based platforms. Unmanned aerial vehicles (UAVs) have had recent widespread application to capture high-resolution information on snow processes and are herein applied to the sub-canopy snow depth challenge. Previous demonstrations of snow depth mapping with UAV structure from motion (SfM) and airborne lidar have focussed on non-vegetated surfaces or reported large errors in the presence of vegetation. In contrast, UAV-lidar systems have high-density point clouds and measure returns from a wide range of scan angles, increasing the likelihood of successfully sensing the sub-canopy snow depth. The effectiveness of UAV lidar and UAV SfM in mapping snow depth in both open and forested terrain was tested in a 2019 field campaign at the Canadian Rockies Hydrological Observatory, Alberta, and at Canadian prairie sites near Saskatoon, Saskatchewan, Canada. Only UAV lidar could successfully measure the sub-canopy snow surface with reliable sub-canopy point coverage and consistent error metrics (root mean square error (RMSE) <0.17 m and bias −0.03 to −0.13 m). Relative to UAV lidar, UAV SfM did not consistently sense the sub-canopy snow surface, the interpolation needed to account for point cloud gaps introduced interpolation artefacts, and error metrics demonstrated relatively large variability (RMSE<0.33 m and bias 0.08 to −0.14 m). With the demonstration of sub-canopy snow depth mapping capabilities, a number of early applications are presented to showcase the ability of UAV lidar to effectively quantify the many multiscale snow processes defining snowpack dynamics in mountain and prairie environments.

Abstract. Vegetation has a tremendous influence on snow processes and snowpack dynamics, yet remote sensing techniques to resolve the spatial variability of sub-canopy snow depth are not always available and are difficult from space-based platforms. Unmanned aerial vehicles (UAVs) have had recent widespread application to capture high-resolution information on snow processes and are herein applied to the sub-canopy snow depth challenge. Previous demonstrations of snow depth mapping with UAV structure from motion (SfM) and airborne lidar have focussed on non-vegetated surfaces or reported large errors in the presence of vegetation. In contrast, UAV-lidar systems have high-density point clouds and measure returns from a wide range of scan angles, increasing the likelihood of successfully sensing the sub-canopy snow depth. The effectiveness of UAV lidar and UAV SfM in mapping snow depth in both open and forested terrain was tested in a 2019 field campaign at the Canadian Rockies Hydrological Observatory, Alberta, and at Canadian prairie sites near Saskatoon, Saskatchewan, Canada. Only UAV lidar could successfully measure the sub-canopy snow surface with reliable sub-canopy point coverage and consistent error metrics (root mean square error (RMSE) <0.17 m and bias −0.03 to −0.13 m). Relative to UAV lidar, UAV SfM did not consistently sense the sub-canopy snow surface, the interpolation needed to account for point cloud gaps introduced interpolation artefacts, and error metrics demonstrated relatively large variability (RMSE<0.33 m and bias 0.08 to −0.14 m). With the demonstration of sub-canopy snow depth mapping capabilities, a number of early applications are presented to showcase the ability of UAV lidar to effectively quantify the many multiscale snow processes defining snowpack dynamics in mountain and prairie environments.

Field-based assessment of transpiration phenology in boreal tree species is a significant challenge. Here we develop an objective approach that uses stem radius change and its correlation with sapwood temperature to determine the timing of phenological changes in transpiration in mixed evergreen species. We test the stem-temp approach using a five year stem-radius dataset from black spruce (Picea mariana) and jack pine (Pinus banksiana) trees in Saskatchewan (2016–2020). We further compare transpiration phenological transition dates from this approach with tower-based phenological assessment from green chromatic coordinate derived from phenocam images, eddy-covariance-derived evapotranspiration and carbon uptake, tower-based measurements of solar-induced chlorophyll fluorescence and snowmelt timing. The stem-temp approach identified the start and end of four key transpiration phenological phases: (i) the end of temperature-driven cycles indicating the start of biological activity, (ii) the onset of stem rehydration, (iii) the onset of transpiration, and (iv) the end of transpiration-driven cycles. The proposed method is thus useful for characterizing the timing of changes in transpiration phenology and provides information about distinct processes that cannot be assessed with canopy-level phenological measurements alone.

Abstract. The US Northern Great Plains and the Canadian Prairies are known as the world’s breadbaskets for its large spring wheat production and exports to the world. It is essential to accurately represent spring wheat growing dynamics and final yield and improve our ability to predict food production under climate change. This study attempts to incorporate spring wheat growth dynamics into the Noah-MP crop model, for a long time period (13-year) and fine spatial scale (4-km). The study focuses on three aspects: (1) developing and calibrating the spring wheat model at point-scale, (2) applying a dynamic planting/harvest date to facilitate large-scale simulations, and (3) applying a temperature stress function to assess crop responses to heat stress amid extreme heat. Model results are evaluated using field observations, satellite leaf area index (LAI), and census data from Statistics Canada and the US Department of Agriculture (USDA). Results suggest that incorporating a dynamic planting/harvest threshold can better constrain the growing season, especially the peak timing and magnitude of wheat LAI, as well as obtain realistic yield compared to prescribing a static province/state-level map. Results also demonstrate an evident control of heat stress upon wheat yield in three Canadian Prairies Provinces, which are reasonably captured in the new temperature stress function. This study has important implications for estimating crop production, simulating the land-atmosphere interactions in croplands, and crop growth’s responses to the raising temperatures amid climate change.

Estimating soil moisture (SM) over the circumpolar boreal forest would have numerous applications including wildfire risk detection, and weather prediction. Evaluation of satellite derived SM retrievals in boreal ecoregions is hindered by available in situ SM observation networks. To address this, an SM monitoring network was established in a boreal forest region in Saskatchewan, Canada. The network is unique as there are no other SM network of similar size in the boreal forest. The network consisted of 17 SM stations within a single Soil Moisture Active Passive (SMAP) satellite observation pixel ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$33\times 33$ </tex-math></inline-formula> km). We present an analysis of the sensitivity and accuracy of SMAP SM products in a boreal forest environment over a two-year period in 2018 and 2019. Results show current SMAP radiometer-based L2 SM products have higher correlation with the in situ lower mineral layer SM than with the top organic layer, although the overall correlation is low. Correlations between in situ mineral layer SM and SMAP brightness-temperature (TB) products are higher than those observed with the SMAP SM product, suggesting current SMAP SM retrieval from the TB using the <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\tau $ </tex-math></inline-formula> – <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\omega $ </tex-math></inline-formula> model introduces large uncertainties in the SM estimation, possibly from uncertain vegetation and surface parameters in the retrieval model. Results show SM can be retrieved using the <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\tau $ </tex-math></inline-formula> – <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\omega $ </tex-math></inline-formula> model with reasonable accuracy over the boreal forest provided the vegetation and soil parameters are optimized. The SM retrieval using a dual channel <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\tau $ </tex-math></inline-formula> – <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\omega $ </tex-math></inline-formula> model, which utilize both horizontally and vertically polarized SMAP TB, performs better than that with a single channel algorithm (SCA), using optimized parameters.

Water use efficiency (WUE) can be calculated using a range of methods differing in carbon uptake and water use variable selection. Consequently, inconsistencies arise between WUE calculations due to complex physical and physiological interactions. The purpose of this study was to quantify and compare WUE estimates (harvest or flux-based) for alfalfa (C3 plant) and maize (C4 plant) and determine effects of input variables, plant physiology and farming practices on estimates. Four WUE calculations were investigated: two “harvest-based” methods, using above ground carbon content and either precipitation or evapotranspiration (ET), and two “flux-based” methods, using gross primary productivity (GPP) and either ET or transpiration. WUE estimates differed based on method used at both half-hourly and seasonal scales. Input variables used in calculations affected WUE estimates, and plant physiology led to different responses in carbon assimilation and water use variables. WUE estimates were also impacted by different plant physiological responses and processing methods, even when the same carbon assimilation and water use variables were considered. This study highlights a need to develop a metric of measuring cropland carbon-water coupling that accounts for all water use components, plant carbon responses, and biomass production.

The boreal forest is a major contributor to the global climate system, therefore, reducing uncertainties in how the forest will respond to a changing climate is critical. One source of uncertainty is the timing and drivers of the spring transition. Remote sensing can provide important information on this transition, but persistent foliage greenness, seasonal snow cover, and a high prevalence of mixed forest stands (both deciduous and evergreen species) complicate interpretation of these signals. We collected tower-based remotely sensed data (reflectance-based vegetation indices and Solar-Induced Chlorophyll Fluorescence [SIF]), stem radius measurements, gross primary productivity, and environmental conditions in a boreal mixed forest stand. Evaluation of this data set shows a two-phased spring transition. The first phase is the reactivation of photosynthesis and transpiration in evergreens, marked by an increase in relative SIF, and is triggered by thawed stems, warm air temperatures, and increased available soil moisture. The second phase is a reduction in bulk photoprotective pigments in evergreens, marked by an increase in the Chlorophyll-Carotenoid Index. Deciduous leaf-out occurs during this phase, marked by an increase in all remotely sensed metrics. The second phase is controlled by soil thaw. Our results demonstrate that remote sensing metrics can be used to detect specific physiological changes in boreal tree species during the spring transition. The two-phased transition explains inconsistencies in remote sensing estimates of the timing and drivers of spring recovery. Our results imply that satellite-based observations will improve by using a combination of vegetation indices and SIF, along with species distribution information.

Atmospheric boundary layer (ABL) dynamics over glaciers mediate the response of glacier mass balance to large‐scale climate forcing. Despite this, very few ABL observations are available over mountain glaciers in complex terrain. An intensive field campaign was conducted in June 2015 at the Athabasca Glacier outlet of Columbia Icefield in the Canadian Rockies. Observations of wind and temperature profiles with novel kite and radio‐acoustic sounding systems showed a well‐defined mesoscale circulation developed between the glacier and snow‐free valley in fair weather. The typical vertical ABL structure above the glacier differed from that expected for “glacier winds”; strong daytime down‐glacier winds extended through the lowest 200 m with no up‐valley return flow aloft. This structure suggests external forcing at mesoscale scales or greater and is provisionally termed an “icefield breeze.” A wind speed maximum near the surface, characteristic of a “glacier wind,” was only observed during night‐time and one afternoon. Lapse rates of air temperature down the glacier centerline show the interaction of down‐glacier cooling driven by sensible heat loss into the ice, entrainment and periodic disruption and warming. Down‐glacier cooling was weaker in “icefield breeze” conditions, while in “glacier wind” conditions, stronger down‐glacier cooling enabled large increases in near‐surface temperature on the lower glacier during periods of surface boundary layer (SBL) disruption. These results raise several questions, including the impact of Columbia Icefield on the ABL and melt of Athabasca Glacier. Future work should use these observations as a testbed for modeling spatio‐temporal variations in the ABL and SBL within complex glaciated terrain.

Abstract. The interior of western Canada, like many similar cold mid- to high-latitude regions worldwide, is undergoing extensive and rapid climate and environmental change, which may accelerate in the coming decades. Understanding and predicting changes in coupled climate–land–hydrological systems are crucial to society yet limited by lack of understanding of changes in cold-region process responses and interactions, along with their representation in most current-generation land-surface and hydrological models. It is essential to consider the underlying processes and base predictive models on the proper physics, especially under conditions of non-stationarity where the past is no longer a reliable guide to the future and system trajectories can be unexpected. These challenges were forefront in the recently completed Changing Cold Regions Network (CCRN), which assembled and focused a wide range of multi-disciplinary expertise to improve the understanding, diagnosis, and prediction of change over the cold interior of western Canada. CCRN advanced knowledge of fundamental cold-region ecological and hydrological processes through observation and experimentation across a network of highly instrumented research basins and other sites. Significant efforts were made to improve the functionality and process representation, based on this improved understanding, within the fine-scale Cold Regions Hydrological Modelling (CRHM) platform and the large-scale Modélisation Environmentale Communautaire (MEC) – Surface and Hydrology (MESH) model. These models were, and continue to be, applied under past and projected future climates and under current and expected future land and vegetation cover configurations to diagnose historical change and predict possible future hydrological responses. This second of two articles synthesizes the nature and understanding of cold-region processes and Earth system responses to future climate, as advanced by CCRN. These include changing precipitation and moisture feedbacks to the atmosphere; altered snow regimes, changing balance of snowfall and rainfall, and glacier loss; vegetation responses to climate and the loss of ecosystem resilience to wildfire and disturbance; thawing permafrost and its influence on landscapes and hydrology; groundwater storage and cycling and its connections to surface water; and stream and river discharge as influenced by the various drivers of hydrological change. Collective insights, expert elicitation, and model application are used to provide a synthesis of this change over the CCRN region for the late 21st century.

Representing climate-crop interactions is critical to earth system modeling. Despite recent progress in modeling dynamic crop growth and irrigation in land surface models (LSMs), transitioning thes...

Abstract. Vegetation has a tremendous influence on snow processes and snowpack dynamics, yet remote sensing techniques to resolve the spatial variability of sub-canopy snow depth are not always available and are difficult from space-based platforms. Unmanned aerial vehicles (UAVs) have had recent widespread application to capture high-resolution information on snow processes and are herein applied to the sub-canopy snow depth challenge. Previous demonstrations of snow depth mapping with UAV structure from motion (SfM) and airborne lidar have focussed on non-vegetated surfaces or reported large errors in the presence of vegetation. In contrast, UAV-lidar systems have high-density point clouds and measure returns from a wide range of scan angles, increasing the likelihood of successfully sensing the sub-canopy snow depth. The effectiveness of UAV lidar and UAV SfM in mapping snow depth in both open and forested terrain was tested in a 2019 field campaign at the Canadian Rockies Hydrological Observatory, Alberta, and at Canadian prairie sites near Saskatoon, Saskatchewan, Canada. Only UAV lidar could successfully measure the sub-canopy snow surface with reliable sub-canopy point coverage and consistent error metrics (root mean square error (RMSE) <0.17 m and bias −0.03 to −0.13 m). Relative to UAV lidar, UAV SfM did not consistently sense the sub-canopy snow surface, the interpolation needed to account for point cloud gaps introduced interpolation artefacts, and error metrics demonstrated relatively large variability (RMSE<0.33 m and bias 0.08 to −0.14 m). With the demonstration of sub-canopy snow depth mapping capabilities, a number of early applications are presented to showcase the ability of UAV lidar to effectively quantify the many multiscale snow processes defining snowpack dynamics in mountain and prairie environments.

Abstract. Vegetation has a tremendous influence on snow processes and snowpack dynamics, yet remote sensing techniques to resolve the spatial variability of sub-canopy snow depth are not always available and are difficult from space-based platforms. Unmanned aerial vehicles (UAVs) have had recent widespread application to capture high-resolution information on snow processes and are herein applied to the sub-canopy snow depth challenge. Previous demonstrations of snow depth mapping with UAV structure from motion (SfM) and airborne lidar have focussed on non-vegetated surfaces or reported large errors in the presence of vegetation. In contrast, UAV-lidar systems have high-density point clouds and measure returns from a wide range of scan angles, increasing the likelihood of successfully sensing the sub-canopy snow depth. The effectiveness of UAV lidar and UAV SfM in mapping snow depth in both open and forested terrain was tested in a 2019 field campaign at the Canadian Rockies Hydrological Observatory, Alberta, and at Canadian prairie sites near Saskatoon, Saskatchewan, Canada. Only UAV lidar could successfully measure the sub-canopy snow surface with reliable sub-canopy point coverage and consistent error metrics (root mean square error (RMSE) <0.17 m and bias −0.03 to −0.13 m). Relative to UAV lidar, UAV SfM did not consistently sense the sub-canopy snow surface, the interpolation needed to account for point cloud gaps introduced interpolation artefacts, and error metrics demonstrated relatively large variability (RMSE<0.33 m and bias 0.08 to −0.14 m). With the demonstration of sub-canopy snow depth mapping capabilities, a number of early applications are presented to showcase the ability of UAV lidar to effectively quantify the many multiscale snow processes defining snowpack dynamics in mountain and prairie environments.

Abstract. Local-scale advection of energy from warm snow-free surfaces to cold snow-covered surfaces is an important component of the energy balance during snow-cover depletion. Unfortunately, this process is difficult to quantify in one-dimensional snowmelt models. This paper proposes a simple sensible and latent heat advection model for snowmelt situations that can be readily coupled to one-dimensional energy balance snowmelt models. An existing advection parameterization was coupled to a conceptual frozen soil infiltration surface water retention model to estimate the areal average sensible and latent heat advection contributions to snowmelt. The proposed model compared well with observations of latent and sensible heat advection, providing confidence in the process parameterizations and the assumptions applied. Snow-covered area observations from unmanned aerial vehicle imagery were used to update and evaluate the scaling properties of snow patch area distribution and lengths. Model dynamics and snowmelt implications were explored within an idealized modelling experiment, by coupling to a one-dimensional energy balance snowmelt model. Dry, snow-free surfaces were associated with advection of dry air that compensated for positive sensible heat advection fluxes and so limited the net influence of advection on snowmelt. Latent and sensible heat advection fluxes both contributed positive fluxes to snow when snow-free surfaces were wet and enhanced net advection contributions to snowmelt. The increased net advection fluxes from wet surfaces typically develop towards the end of snowmelt and offset decreases in the one-dimensional areal average melt energy that declines with snow-covered area. The new model can be readily incorporated into existing one-dimensional snowmelt hydrology and land surface scheme models and will foster improvements in snowmelt understanding and predictions.

Spring snowmelt is the most important hydrological event in agricultural cold regions, recharging soil moisture and generating the majority of annual runoff. Melting agricultural snowcovers are pat...

Abstract Forest clearings are common features of evergreen forests and produce snowpack accumulation and melt differing from that in adjacent forests and open terrain. This study has investigated the challenges in specifying the turbulent fluxes of sensible and latent heat to snowpacks in forest clearings. The snowpack in two forest clearings in the Canadian Rockies was simulated using a one-dimensional (1D) snowpack model. A trade-off was found between optimizing against measured snow surface temperature or snowmelt when choosing how to specify the turbulent fluxes. Schemes using the Monin–Obukhov similarity theory tended to produce negatively biased surface temperature, while schemes that enhanced turbulent fluxes, to reduce the surface temperature bias, resulted in too much melt. Uncertainty estimates from Monte Carlo experiments showed that no realistic parameter set could successfully remove biases in both surface temperature and melt. A simple scheme that excludes atmospheric stability correction was required to successfully simulate surface temperature under low wind speed conditions. Nonturbulent advective fluxes and/or nonlocal sources of turbulence are thought to account for the maintenance of heat exchange in low-wind conditions. The simulation of snowmelt was improved by allowing enhanced latent heat fluxes during low-wind conditions. Caution is warranted when snowpack models are optimized on surface temperature, as model tuning may compensate for deficiencies in conceptual and numerical models of radiative, conductive, and turbulent heat exchange at the snow surface and within the snowpack. Such model tuning could have large impacts on the melt rate and timing of the snow-free transition in simulations of forest clearings within hydrological and meteorological models.

Abstract Greenhouse gas (GHG) emissions from agricultural soils in the Canadian Prairie region are generally low and, due to dry, well aerated soil conditions, can be quite variable. Compared to dryland (rainfed) crop production, irrigated cropping has potential to contribute greater quantities of soil derived nitrous oxide (N2O), carbon dioxide (CO2), and methane (CH4) to the atmosphere as producers target higher yields by minimizing soil moisture limitations and applying greater amounts of nitrogen fertilizers. However, the actual GHG dynamics from irrigated soils in this region are not well understood as there have been few field-based studies in the semi-arid prairies of western Canada. The goal of this study was to identify how emissions of soil derived N2O, CO2, and CH4 are influenced by changes in soil temperature, water status, and nitrogen rates brought about by irrigated crop management. This was achieved through continuous, in-situ monitoring of soil conditions and chamber-based measurements of soil GHG flux. The most notable change in soil conditions brought about by irrigation was elevated moisture levels, which appeared to influence the flux dynamics of all three agricultural greenhouse gases—specifically, a reduction in CH4 uptake and periodic increases in CO2 and N2O emissions. Despite the reduced soil moisture limitation, annual N2O emissions from the irrigated cropping system were much lower than those calculated using the current Canadian National GHG Inventory Reporting. This suggests that annual emissions are limited more by N availability rather than moisture deficits, as the current method for emissions accounting assumes. Consequently, our results indicate that emissions from irrigated cropping systems in the semi-arid Canadian Prairies are overestimated by the current inventory approach. Moreover, because irrigated crop production involves more than just the application of water, our results demonstrate that a more systems-oriented approach to GHG accounting is required to capture the combined effects of water-soil-crop management on GHG emissions from irrigated cropping systems.

Abstract On the Canadian Prairies, agricultural practices result in millions of hectares of standing crop stubble that gradually emerges during snowmelt. The importance of stubble in trapping wind-blown snow and retaining winter snowfall has been well demonstrated. However, stubble is not explicitly accounted for in hydrological or energy balance snowmelt models. This paper relates measurable stubble parameters (height, width, areal density, and albedo) to the snowpack energy balance and snowmelt with the new, physically based Stubble–Snow–Atmosphere Model (SSAM). Novel process representations of SSAM quantify the attenuation of shortwave radiation by exposed stubble, the sky and vegetation view factors needed to solve longwave radiation terms, and a resistance scheme for stubble–snow–atmosphere fluxes to solve for surface temperatures and turbulent fluxes. SSAM results were compared to observations of radiometric snow-surface temperature, stubble temperature, snow-surface solar irradiance, areal-average turbulent fluxes, and snow water equivalent from two intensive field campaigns during snowmelt in 2015 and 2016 over wheat and canola stubble in Saskatchewan, Canada. Uncalibrated SSAM simulations compared well with these observations, providing confidence in the model structure and parameterization. A sensitivity analysis conducted using SSAM revealed compensatory relationships in energy balance terms that result in a small increase in net snowpack energy as stubble exposure increases.

The breakup of snow cover into patches during snowmelt leads to a dynamic, heterogeneous land surface composed of melting snow, and wet and dry soil and plant surfaces. Energy exchange with the atmosphere is therefore complicated by horizontal gradients in surface temperature and humidity as snow surface temperature and humidity are regulated by the phase change of melting snow unlike snow-free areas. Airflow across these surface transitions results in local-scale advection of energy that has been documented as sensible heat during snowmelt, while latent heat advection has received scant attention. Herein, results are presented from an experiment measuring near-surface profiles of air temperature and humidity across snow-free to snow-covered transitions that demonstrates that latent heat advection can be the same order of magnitude as sensible heat advection and is therefore an important source of snowmelt energy. Latent heat advection is conditional on an upwind source of water vapor from a wetted snow-free surface.

Abstract The snow surface temperature (SST) is essential for estimating longwave radiation fluxes from snow. SST can be diagnosed using finescale multilayer snow physics models that track changes in snow properties and internal energy; however, these models are heavily parameterized, have high predictive uncertainty, and require continuous simulation to estimate prognostic state variables. Here, a relatively simple model to estimate SST that is not reliant on prognostic state variables is proposed. The model assumes that the snow surface is poorly connected thermally to the underlying snowpack and largely transparent for most of the shortwave radiation spectrum, such that a snow surface energy balance among only sensible heat, latent heat, longwave radiation, and near-infrared radiation is possible and is called the radiative psychrometric model (RPM). The RPM SST is sensitive to air temperature, humidity, ventilation, and longwave irradiance and is secondarily affected by absorption of near-infrared radiation at the snow surface that was higher where atmospheric deposition of particulates was more likely to be higher. The model was implemented with neutral stability, an implicit windless exchange coefficient, and constant shortwave absorption factors and aerodynamic roughness lengths. It was evaluated against radiative SST measurements from the Canadian Prairies and Rocky Mountains, French Alps, and Bolivian Andes. With optimized and global shortwave absorption and aerodynamic roughness length parameters, the model is shown to accurately predict SST under a wide range of conditions, providing superior predictions when compared to air temperature, dewpoint, or ice bulb calculation approaches.