Heterogenous snow cover derived uncertainty in Arctic carbon budget

Author(s):  
Hotaek Park ◽  
Youngwook Kim

<p>The winter of northern Arctic regions is characterized by strong winds that lead to frequent blowing snow and thus heterogeneous snow cover, which critically affects permafrost hydrothermal processes and the associated feedbacks across the northern regions. However until now, observations and models have not documented the blowing snow impacts. The blowing snow process has coupled into a land surface model CHANGE, and the improved model was applied to observational sites in the northeastern Siberia for 1979–2016. The simulated snow depth and soil temperature showed general agreements with the observations. To quantify the impacts of blowing snow on permafrost temperatures and the associated greenhouse gases, two decadal experiments that included or excluded blowing snow, were conducted for the observational sites and over the pan-Arctic scale. The differences between the two experiments represent impacts of the blowing snow on the analytical components. The blowing snow-induced thinner snow depth resulted in cooler permafrost temperature and lower active layer thickness; this lower temperature limited the vegetation photosynthetic activity due to the increased soil moisture stress in terms of larger soil ice portion and hence lower ecosystem productivity. The cooler permafrost temperature is also linked to less decomposition of soil organic matter and lower releases of CO2 and CH4 to the atmosphere. These results suggest that the most land models without a blowing snow component likely overestimate the release of greenhouse gases from the tundra regions. There is a strong need to improve land surface models for better simulations and future projections of the northern environmental changes.</p>

2013 ◽  
Vol 14 (1) ◽  
pp. 203-219 ◽  
Author(s):  
Eric Brun ◽  
Vincent Vionnet ◽  
Aaron Boone ◽  
Bertrand Decharme ◽  
Yannick Peings ◽  
...  

Abstract The Crocus snowpack model within the Interactions between Soil–Biosphere–Atmosphere (ISBA) land surface model was run over northern Eurasia from 1979 to 1993, using forcing data extracted from hydrometeorological datasets and meteorological reanalyses. Simulated snow depth, snow water equivalent, and density over open fields were compared with local observations from over 1000 monitoring sites, available either once a day or three times per month. The best performance is obtained with European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim). Provided blowing snow sublimation is taken into account, the simulations show a small bias and high correlations in terms of snow depth, snow water equivalent, and density. Local snow cover durations as well as the onset and vanishing dates of continuous snow cover are also well reproduced. A major result is that the overall performance of the simulations is very similar to the performance of existing gridded snow products, which, in contrast, assimilate local snow depth observations. Soil temperature at 20-cm depth is reasonably well simulated. The methodology developed in this study is an efficient way to evaluate different meteorological datasets, especially in terms of snow precipitation. It reveals that the temporal disaggregation of monthly precipitation in the hydrometeorological dataset from Princeton University significantly impacts the rain–snow partitioning, deteriorating the simulation of the onset of snow cover as well as snow depth throughout the cold season.


2013 ◽  
Vol 7 (2) ◽  
pp. 631-645 ◽  
Author(s):  
H. Park ◽  
J. Walsh ◽  
A. N. Fedorov ◽  
A. B. Sherstiukov ◽  
Y. Iijima ◽  
...  

Abstract. This study not only examined the spatiotemporal variations of active-layer thickness (ALT) in permafrost regions during 1948–2006 over the terrestrial Arctic regions experiencing climate changes, but also identified the associated drivers based on observational data and a simulation conducted by a land surface model (CHANGE). The focus on the ALT extends previous studies that have emphasized ground temperatures in permafrost regions. The Ob, Yenisey, Lena, Yukon, and Mackenzie watersheds are foci of the study. Time series of ALT in Eurasian watersheds showed generally increasing trends, while the increase in ALT in North American watersheds was not significant. However, ALT in the North American watersheds has been negatively anomalous since 1990 when the Arctic air temperature entered into a warming phase. The warming temperatures were not simply expressed to increases in ALT. Since 1990 when the warming increased, the forcing of the ALT by the higher annual thawing index (ATI) in the Mackenzie and Yukon basins has been offset by the combined effects of less insulation caused by thinner snow depth and drier soil during summer. In contrast, the increasing ATI together with thicker snow depth and higher summer soil moisture in the Lena contributed to the increase in ALT. The results imply that the soil thermal and moisture regimes formed in the pre-thaw season(s) provide memory that manifests itself during the summer. The different ALT anomalies between Eurasian and North American watersheds highlight increased importance of the variability of hydrological variables.


2020 ◽  
Vol 12 (4) ◽  
pp. 645 ◽  
Author(s):  
Sujay Kumar ◽  
David Mocko ◽  
Carrie Vuyovich ◽  
Christa Peters-Lidard

Surface albedo has a significant impact in determining the amount of available net radiation at the surface and the evolution of surface water and energy budget components. The snow accumulation and timing of melt, in particular, are directly impacted by the changes in land surface albedo. This study presents an evaluation of the impact of assimilating Moderate Resolution Imaging Spectroradiometer (MODIS)-based surface albedo estimates in the Noah multi-parameterization (Noah-MP) land surface model, over the continental US during the time period from 2000 to 2017. The evaluation of simulated snow depth and snow cover fields show that significant improvements from data assimilation (DA) are obtained over the High Plains and parts of the Rocky Mountains. Earlier snowmelt and reduced agreements with reference snow depth measurements, primarily over the Northeast US, are also observed due to albedo DA. Most improvements from assimilation are observed over locations with moderate vegetation and lower elevation. The aggregate impact on evapotranspiration and runoff from assimilation is found to be marginal. This study also evaluates the relative and joint utility of assimilating fractional snow cover and surface albedo measurements. Relative to surface albedo assimilation, fractional snow cover assimilation is found to provide smaller improvements in the simulated snow depth fields. The configuration that jointly assimilates surface albedo and fractional snow cover measurements is found to provide the most beneficial improvements compared to the univariate DA configurations for surface albedo or fractional snow cover. Overall, the study also points to the need for improving the albedo formulations in land surface models and the incorporation of observational uncertainties within albedo DA configurations.


2015 ◽  
Vol 16 (4) ◽  
pp. 1736-1741 ◽  
Author(s):  
Sujay V. Kumar ◽  
Christa D. Peters-Lidard ◽  
Kristi R. Arsenault ◽  
Augusto Getirana ◽  
David Mocko ◽  
...  

Abstract Accurate determination of snow conditions is important for several water management applications, partly because of the significant influence of snowmelt on seasonal streamflow prediction. This article examines an approach using snow cover area (SCA) observations as snow detection constraints during the assimilation of snow depth retrievals from passive microwave sensors. Two different SCA products [the Interactive Multisensor Snow and Ice Mapping System (IMS) and the Moderate Resolution Imaging Spectroradiometer (MODIS)] are employed jointly with the snow depth retrievals from a variety of sensors for data assimilation in the Noah land surface model. The results indicate that the use of MODIS data is effective in obtaining added improvements (up to 6% improvement in aggregate RMSE) in snow depth fields compared to assimilating passive microwave data alone, whereas the impact of IMS data is small. The improvements in snow depth fields are also found to translate to small yet systematic improvements in streamflow estimates, especially over the western United States, the upper Missouri River, and parts of the Northeast and upper Mississippi River. This study thus demonstrates a simple approach for exploiting the information from SCA observations in data assimilation.


Author(s):  
Gonzalo Leonardini ◽  
François Anctil ◽  
Vincent Vionnet ◽  
Maria Abrahamowicz ◽  
Daniel F. Nadeau ◽  
...  

AbstractThe Soil, Vegetation, and Snow (SVS) land surface model was recently developed at Environment and Climate Change Canada (ECCC) for operational numerical weather prediction and hydrological forecasting. This study examined the performance of the snow scheme in the SVS model over multiple years at ten well-instrumented sites from the Earth System Model-Snow Model Intercomparison Project (ESM-SnowMIP), which covers alpine, maritime and taiga climates. The SVS snow scheme is a simple single-layer snowpack scheme that uses the force-restore method. Stand-alone, point-scale verification tests showed that the model is able to realistically reproduce the main characteristics of the snow cover at these sites, namely snow water equivalent, density, snow depth, surface temperature, and albedo. SVS accurately simulated snow water equivalent, density and snow depth at open sites, but exhibited lower performance for subcanopy snowpacks (forested sites). The lower performance was attributed mainly to the limitations of the compaction scheme and the absence of a snow interception scheme. At open sites, the SVS snow surface temperatures were well represented but exhibited a cold bias, which was due to poor representation at night. SVS produced a reasonably accurate representation of snow albedo, but there was a tendency to overestimate late winter albedo. Sensitivity tests suggested improvements associated with the snow melting formulation in SVS.


Hydrology ◽  
2020 ◽  
Vol 7 (2) ◽  
pp. 20
Author(s):  
Michael Weber ◽  
Moritz Feigl ◽  
Karsten Schulz ◽  
Matthias Bernhardt

To find the adequate spatial model discretization scheme, which balances the models capabilities and the demand for representing key features in reality, is a challenging task. It becomes even more challenging in high alpine catchments, where the variability of topography and meteorology over short distances strongly influences the distribution of snow cover, the dominant component in the alpine water cycle. For the high alpine Research Catchment Zugspitze (RCZ) a new method for objective delineation of hydrological response units (HRUs) using a time series of high resolution LIDAR derived snow depth maps and the physiographic properties of the RCZ is introduced. Via principle component analysis (PCA) of these maps, a dominant snow depth pattern, that turned out to be largely defined during the (winter) accumulation period was identified. This dominant pattern serves as a reference for HRU delineations on the basis of cluster analyses of the catchment’s physiographic properties. The method guarantees for an appropriate, objective, spatial discretization scheme, which allows for a reliable and meaningful reproduction of snow cover variability with the Cold Regions Hydrological Model — at the same time avoiding significant increase of computational demands. Different HRU schemes were evaluated with measured snow depth and the comparison of their model results identified significant differences in model output and best performance of the scheme which best represents measured snow depth distribution.


2020 ◽  
Author(s):  
Edoardo Raparelli ◽  
Paolo Tuccella ◽  
Rossella Ferretti ◽  
Frank S. Marzano

<p>Italy is a territory characterized by complex orography. Its main mountain chains are the Alps, which identify the northern Italian border, and the Apennines, which cross the entire Italian peninsula ranging from north-west to south-east. The major Apennines peaks reach almost 3000 meters and are located in central Italy, in the Abruzzo region. The near Mediterranean sea is an important source of moisture, which permits to this region to experience a substantial snow cover during winter. Thanks to the orientation of the Apennines chain and the height of its peaks the Abruzzo region is characterized by different climate types. This affects the precipitation patterns and the snowpack evolution, resulting in high regional variability of the snow cover. The goal of this study is to investigate the snow cover evolution in the Abruzzo region, using and comparing different snowpack models. To this end we have used the Weather Research and Forecasting (WRF) model to drive the Noah Land Surface Model (LSM) and the sophisticated three-dimensional snow cover model Alpine3D to simulate the snow cover evolution at regional scale. Noah LSM is already on-line coupled with WRF, but this is not the case for Alpine3D. Thus we have modified and used the interfacing library MeteoIO to force Alpine3D with the meteorological data simulated with WRF, off-line coupling the two models. We have validated the WRF simulation using a dense network of automatic weather stations (AWS), obtaining good agreement between simulated and observed data. We have found that the snow depth simulated with Noah LSM presents a negative bias, caused by the inability of the model to reproduce correctly the snow densification rate. Instead, Alpine3D is capable to better reproduce the observed densification rate, thanks to its more detailed description of the snow metamorphism processes. However, the snow depth simulated with Alpine3D presents a negative bias, caused by an underestimation of the new snow depth, which has a negative impact on the entire simulation.</p>


2007 ◽  
Vol 8 (6) ◽  
pp. 1243-1263 ◽  
Author(s):  
J. E. Cherry ◽  
L-B. Tremblay ◽  
M. Stieglitz ◽  
G. Gong ◽  
S. J. Déry

Abstract A new product, the Pan-Arctic Snowfall Reconstruction (PASR), is developed to address the problem of cold season precipitation gauge biases for the 1940–99 period. The method used to create the PASR is different from methods used in other large-scale precipitation data products and has not previously been employed for estimating pan-arctic snowfall. The NASA Interannual-to-Seasonal Prediction Project Catchment Land Surface Model is used to reconstruct solid precipitation from observed snow depth and surface air temperatures. The method is tested at four stations in the United States and Canada where results are examined in depth. Reconstructed snowfall at Dease Lake, British Columbia, and Barrow, Alaska, is higher than gauge observations. Reconstructed snowfall at Regina, Saskatchewan, and Minot, North Dakota, is lower than gauge observations, probably because snow is transported by wind out of the Prairie region and enters the hydrometeorological cycle elsewhere. These results are similar to gauge biases estimated by a water budget approach. Reconstructed snowfall is consistently higher than snowfall from the 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) but does not have a consistent relationship with snowfall derived from the WMO Solid Precipitation Intercomparison Project correction algorithms. Advantages of the PASR approach include that 1) the assimilation of snow depth observations captures blowing snow where it is deposited and 2) the modeling approach takes into account physical snowpack evolution. These advantages suggest that the PASR product could be a valuable alternative to statistical gauge corrections and that arctic ground-based solid precipitation observing networks might emphasize snow depth measurements over gauges.


2021 ◽  
Author(s):  
Won Young Lee ◽  
Hyeon-Ju Gim ◽  
Seon Ki Park

Abstract. Snow on land surface plays a vital role in the interaction between land and atmosphere in the state-of-the-art land surface models (LSMs) and the real world. Since the snow cover affects the snow albedo and the ground and soil heat fluxes, it is crucial to detect snow cover changes accurately. It is challenging to acquire observation data for snow cover, snow albedo, and snow depth; thus, an excellent alternative is to use the simulation data produced by the LSMs that calculate the snow-related physical processes. The LSMs show significant differences in the complexities of the snow parameterizations in terms of variables and processes considered. Thus, the synthetic intercomparisons of the snow physics in the LSMs will help the improvement of each LSM. This study revealed and discussed the differences in the parameterizations among LSMs related to snow cover fraction, snow albedo, and snow density. We selected the most popular and well-documented LSMs embedded in the Earth System Model or operational forecasting systems. We examined single layer schemes, including the Unified Noah Land Surface Model (Noah LSM), the Hydrology Tiled ECMWF Scheme of Surface Exchanges over Land (HTESSEL), the Biosphere-Atmosphere Transfer Scheme (BATS), the Canadian Land Surface Scheme (CLASS), and multilayer schemes of intermediate complexity including the Community Noah Land Surface Model with Multi-Parameterization Options (Noah-MP), the Community Land Model version 5 (CLM 5), the Joint UK Land Environment Simulator (JULES), and the Interaction Soil-Biosphere-Atmosphere (ISBA). First, we identified that BATS, Noah-MP, JULES, and ISBA reflect the snow depth and roughness length to parameterize snow cover fraction, and CLM 5 accounts for the standard deviation of the elevation value for the snow cover decay function. Second, CLM 5 and BATS are relatively complex, so that they explicitly take into account the solar zenith angle, black carbon, mineral dust, organic carbon, and ice grain size for the determinations of snow albedo. Besides, JULES and ISBA are also complicated model which concerns ice grain size, solar zenith angle, new snow depth, fresh snowfall rate, and surface temperature for the albedo scheme. Third, HTESSEL, CLM 5, and ISBA considered the effects of both wind and temperature in the determinations of the new snow density. Especially, ISBA and JULES considered internal snow characteristics such as snow viscosity, snow temperature, and vertical stress for parameterizing new snow density. The future outlook discussed geomorphic and vegetation-related variables for the further improvement of the LSMs. Previous studies clearly show that spatio-temporal variation of snow is due to the influence of altitude, slope, and vegetation condition. Therefore, we recommended applying geomorphic and vegetation factors such as elevation, slope, time-varying roughness length, vegetation indexes, or optimized parameters according to the land surface type to parameterize snow-related physical processes.


Geosciences ◽  
2018 ◽  
Vol 8 (12) ◽  
pp. 484 ◽  
Author(s):  
Kristi Arsenault ◽  
Paul Houser

Snow depletion curves (SDC) are functions that are used to show the relationship between snow covered area and snow depth or water equivalent. Previous snow cover data assimilation (DA) studies have used theoretical SDC models as observation operators to map snow depth to snow cover fraction (SCF). In this study, a new approach is introduced that uses snow water equivalent (SWE) observations and satellite-based SCF retrievals to derive SDC relationships for use in an Ensemble Kalman filter (EnKF) to assimilate snow cover estimates. A histogram analysis is used to bin the SWE observations, which the corresponding SCF observations are then averaged within, helping to constrain the amount of data dispersion across different temporal and regional conditions. Logarithmic functions are linearly regressed with the binned average values, for two U.S. mountainous states: Colorado and Washington. The SDC-based logarithmic functions are used as EnKF observation operators, and the satellite-based SCF estimates are assimilated into a land surface model. Assimilating satellite-based SCF estimates with the observation-based SDC shows a reduction in SWE-related RMSE values compared to the model-based SDC functions. In addition, observation-based SDC functions were derived for different intra-annual and physiographic conditions, and landcover and elevation bands. Lower SWE-based RMSE values are also found with many of these categorical observation-based SDC EnKF experiments. All assimilation experiments perform better than the open-loop runs, except for the Washington region’s 2004–2005 snow season, which was a major drought year that was difficult to capture with the ensembles and observations.


Sign in / Sign up

Export Citation Format

Share Document