skip to main content


Title: Analysis of Water Vapor Fluxes Over a Seasonal Snowpack Using the Maximum Entropy Production Model
Snow cover plays a key role in the water and energy budgets over cold regions. Understanding and parameterizing water and heat exchange over snow surfaces in hydrologic models remains a major challenge. An innovative approach based on the theory of maximum entropy production (MEP) was developed for modeling energy budgets for snow-covered surfaces. This study generalizes the MEP model to simulate surface water vapor (latent heat) fluxes over an entire snowpack lifecycle, including snow accumulation and melting during the early growing season. The expanded MEP model combines soil evaporation, canopy transpiration, and snow sublimation to evaluate snow water loss during the lifecycle of the snowpack. Two hypotheses are tested: (1) sublimation becomes negligible during snowmelt when snowpack is isothermal (0°C) and (2) transpiration is progressively activated as a function of the air temperature during vegetation awakening. The proposed approach is shown to be effective for modeling the total surface water vapor fluxes over the snowpack's lifecycle. Both the hypotheses are supported by field observations.  more » « less
Award ID(s):
1724633
NSF-PAR ID:
10287801
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of geophysical research
Volume:
126
ISSN:
2169-8953
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Surface fluxes and atmospheric boundary layer budgets of enthalpy and momentum are quantified using bulk and integral methods based on measurements obtained during a research vessel transect across the Terra Nova Bay polynya during an intense late autumnal katabatic wind event. The surface sensible, latent, and net radiation heat fluxes had maximum upward values of 2,500 ± 600, 400 ± 100, and 100 ± 7 Wm−2, respectively, which occurred over open regions that had been cleared of sea ice by the wind stress. In these open areas, sea spray enhanced the sensible and latent heat fluxes (LHFs) by an estimated 106% and 18%, respectively. As sea ice formed on the surface and became thicker in the downwind direction, the sensible and LHFs decreased to 50% of their maximum values over pancake ice and to 5% at the downwind end of the transect, where snow‐covered young ice flows were present. Snow growth removed over 50% of the water vapor that came from the surface. The surface wind stress ranged from mean values of 2.9 Nm−2in the most upwind regions to 0.9 Nm−2at the end of the transect, with smaller scale variations of almost a factor of three due to different surface types (roughnesses) and gustiness. The downwind slowing and turning of the wind vector can be almost entirely explained by frictional and inertial forces, indicating the horizontal pressure gradient was weak. This case could serve as a case study for future modeling studies of coastal polynyas.

     
    more » « less
  2. Abstract

    Several recent studies from both Greenland and Antarctica have reported significant changes in the water isotopic composition of near‐surface snow between precipitation events. These changes have been linked to isotopic exchange with atmospheric water vapor and sublimation‐induced fractionation, but the processes are poorly constrained by observations. Understanding and quantifying these processes are crucial to both the interpretation of ice core climate proxies and the formulation of isotope‐enabled general circulation models. Here, we present continuous measurements of the water isotopic composition in surface snow and atmospheric vapor together with near‐surface atmospheric turbulence and snow‐air latent and sensible heat fluxes, obtained at the East Greenland Ice‐Core Project drilling site in summer 2016. For two 4‐day‐long time periods, significant diurnal variations in atmospheric water isotopologues are observed. A model is developed to explore the impact of this variability on the surface snow isotopic composition. Our model suggests that the snow isotopic composition in the upper subcentimeter of the snow exhibits a diurnal variation with amplitudes in δ18O and δD of ~2.5‰ and ~13‰, respectively. As comparison, such changes correspond to 10–20% of the magnitude of seasonal changes in interior Greenland snow pack isotopes and of the change across a glacial‐interglacial transition. Importantly, our observation and model results suggest, that sublimation‐induced fractionation needs to be included in simulations of exchanges between the vapor and the snow surface on diurnal timescales during summer cloud‐free conditions in northeast Greenland.

     
    more » « less
  3. Abstract. Cold content is a measure of a snowpack's energy deficit and is a linear function of snowpack mass and temperature. Positive energy fluxes into a snowpack must first satisfy the remaining energy deficit before snowmelt runoff begins, making cold content a key component of the snowpack energy budget. Nevertheless, uncertainty surrounds cold content development and its relationship to snowmelt, likely because of a lack of direct observations. This work clarifies the controls exerted by air temperature, precipitation, and negative energy fluxes on cold content development and quantifies the relationship between cold content and snowmelt timing and rate at daily to seasonal timescales. The analysis presented herein leverages a unique long-term snow pit record along with validated output from the SNOWPACK model forced with 23 water years (1991–2013) of quality controlled, infilled hourly meteorological data from an alpine and subalpine site in the Colorado Rocky Mountains. The results indicated that precipitation exerted the primary control on cold content development at our two sites with snowfall responsible for 84.4 and 73.0% of simulated daily gains in the alpine and subalpine, respectively. A negative surface energy balance – primarily driven by sublimation and longwave radiation emission from the snowpack – during days without snowfall provided a secondary pathway for cold content development, and was responsible for the remaining 15.6 and 27.0% of cold content additions. Non-zero cold content values were associated with reduced snowmelt rates and delayed snowmelt onset at daily to sub-seasonal timescales, while peak cold content magnitude had no significant relationship to seasonal snowmelt timing. These results suggest that the information provided by cold content observations and/or simulations is most relevant to snowmelt processes at shorter timescales, and may help water resource managers to better predict melt onset and rate.

     
    more » « less
  4. Abstract

    The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

     
    more » « less
  5. Abstract

    A Lagrangian snow‐evolution model (SnowModel‐LG) was used to produce daily, pan‐Arctic, snow‐on‐sea‐ice, snow property distributions on a 25 × 25‐km grid, from 1 August 1980 through 31 July 2018 (38 years). The model was forced with NASA's Modern Era Retrospective‐Analysis for Research and Applications‐Version 2 (MERRA‐2) and European Centre for Medium‐Range Weather Forecasts (ECMWF) ReAnalysis‐5th Generation (ERA5) atmospheric reanalyses, and National Snow and Ice Data Center (NSIDC) sea ice parcel concentration and trajectory data sets (approximately 61,000, 14 × 14‐km parcels). The simulations performed full surface and internal energy and mass balances within a multilayer snowpack evolution system. Processes and features accounted for included rainfall, snowfall, sublimation from static‐surfaces and blowing‐snow, snow melt, snow density evolution, snow temperature profiles, energy and mass transfers within the snowpack, superimposed ice, and ice dynamics. The simulations produced horizontal snow spatial structures that likely exist in the natural system but have not been revealed in previous studies spanning these spatial and temporal domains. Blowing‐snow sublimation made a significant contribution to the snowpack mass budget. The superimposed ice layer was minimal and decreased over the last four decades. Snow carryover to the next accumulation season was minimal and sensitive to the melt‐season atmospheric forcing (e.g., the average summer melt period was 3 weeks or 50% longer with ERA5 forcing than MERRA‐2 forcing). Observed ice dynamics controlled the ice parcel age (in days), and ice age exerted a first‐order control on snow property evolution.

     
    more » « less