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.
- Award ID(s):
- 1724633
- NSF-PAR ID:
- 10287801
- 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
-
Abstract -
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.
-
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.
-
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.
-
A Lagrangian Snow‐Evolution System for Sea‐Ice Applications (SnowModel‐LG): Part I—Model Description
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.