The role of snow cover affecting boreal-arctic soil freeze–thaw and carbon dynamics

Northern Hemisphere permafrost affected land areas contain about twice as much carbon as the global atmosphere. This vast carbon pool is vulnerable to accelerated losses through mobilization and decomposition under projected global warming. Satellite data records spanning the past 3 decades indicate widespread reductions (~ 0.8–1.3 days decade−1) in the mean annual snow cover extent and frozen-season duration across the pan-Arctic domain, coincident with regional climate warming trends. How the soil carbon pool responds to these changes will have a large impact on regional and global climate. Here, we developed a coupled terrestrial carbon and hydrology model framework with a detailed 1-D soil heat transfer representation to investigate the sensitivity of soil organic carbon stocks and soil decomposition to climate warming and changes in snow cover conditions in the pan-Arctic region over the past 3 decades (1982–2010). Our results indicate widespread soil active layer deepening across the pan-Arctic, with a mean decadal trend of 6.6 ± 12.0 (SD) cm, corresponding to widespread warming. Warming promotes vegetation growth and soil heterotrophic respiration particularly within surface soil layers (≤ 0.2 m). The model simulations also show that seasonal snow cover has a large impact on soil temperatures, whereby increases in snow cover promote deeper (≥ 0.5 m) soil layer warming and soil respiration, while inhibiting soil decomposition from surface (≤ 0.2 m) soil layers, especially in colder climate zones (mean annual T ≤ −10 °C). Our results demonstrate the important control of snow cover on northern soil freeze–thaw and soil carbon decomposition processes and the necessity of considering both warming and a change in precipitation and snow cover regimes in characterizing permafrost soil carbon dynamics

Snow model verification using ensemble prediction and operational benchmarks

Hydrologic model evaluations have traditionally focused on measuring how closely the model can simulate various characteristics of historical observations. Although advancing hydrologic forecasting is an often-stated goal of numerous modeling studies, testing in a forecasting mode is seldom undertaken, limiting information derived from these analyses. One can overcome this limitation through generation, and subsequent analysis, of ensemble hindcasts. In this study, long-range ensemble hindcasts are generated for the available period of record for a basin in southwestern Idaho for the purpose of evaluating the Snow-Atmosphere-Soil Transfer (SAST) model against the current operational benchmark, the National Weather Service's (NWS) snow accumulation and ablation model SNOW17. Both snow models were coupled with the NWS operational rainfall runoff model and ensembles of seasonal discharge and weekly snow water equivalent (SWE) were evaluated. Ensemble predictions from both the SAST and SNOW17 models were better than climatology forecasts, for the period studied. In most cases, the accuracy of the SAST-generated predictions was similar to the SNOW17-generated predictions, except during periods of significant melting. Differences in model performance are partially attributed to initial condition errors. After updating the SWE state in the snow models with the observed SWE, the forecasts were improved during the first 2-4 weeks of the forecast window and the skills were essentially equal in both forecasting systems for the study watershed. Climate dominated the forecast uncertainty in the latter part of the forecast window while initial conditions controlled the forecast skill in the first 3-4 weeks of the forecast. The use of hindcasting in the snow model analysis revealed that, given the dominance of the initial conditions on forecast skill, streamflow predictions will be most improved through the use of state updating. © 2008 American Meteorological Society

Operational snow modeling: Addressing the challenges of an energy balance model for National Weather Service forecasts

Prediction of snowmelt has become a critical issue in much of the western United States given the increasing demand for water supply, changing snow cover patterns, and the subsequent requirement of optimal reservoir operation. The increasing importance of hydrologic predictions necessitates that traditional forecasting systems be re-evaluated periodically to assure continued evolution of the operational systems given scientific advancements in hydrology. The National Weather Service (NWS) SNOW17, a conceptually based model used for operational prediction of snowmelt, has been relatively unchanged for decades. In this study, the Snow-Atmosphere-Soil Transfer (SAST) model, which employs the energy balance method, is evaluated against the SNOW17 for the simulation of seasonal snowpack (both accumulation and melt) and basin discharge. We investigate model performance over a 13-year period using data from two basins within the Reynolds Creek Experimental Watershed located in southwestern Idaho. Both models are coupled to the NWS runoff model [SACramento Soil Moisture Accounting model (SACSMA)] to simulate basin streamflow. Results indicate that while in many years simulated snowpack and streamflow are similar between the two modeling systems, the SAST more often overestimates SWE during the spring due to a lack of mid-winter melt in the model. The SAST also had more rapid spring melt rates than the SNOW17, leading to larger errors in the timing and amount of discharge on average. In general, the simpler SNOW17 performed consistently well, and in several years, better than, the SAST model. Input requirements and related uncertainties, and to a lesser extent calibration, are likely to be primary factors affecting the implementation of an energy balance model in operational streamflow prediction. © 2008 Elsevier B.V. All rights reserved

Implementation of a soil albedo scheme in the CABLEv1.4b land surface model and evaluation against MODIS estimates over Australia

Land surface albedo, the fraction of incoming solar radiation reflected by the land surface, is a key component of the Earth system. This study evaluates snow-free surface albedo simulations by the Community Atmosphere Biosphere Land Exchange (CABLEv1.4b) model with the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Satellite Pour L'Observation de la Terre (SPOT) albedo. We compare results from offline simulations over the Australian continent. The control simulation has prescribed background snow-free and vegetation-free soil albedo derived from MODIS whilst the experiments use a simple parameterisation based on soil moisture and colour, originally from the Biosphere Atmosphere Transfer Scheme (BATS), and adopted in the Common Land Model (CLM). The control simulation, with prescribed soil albedo, shows that CABLE simulates overall albedo over Australia reasonably well, with differences compared to MODIS and SPOT albedos within ±0.1. Application of the original BATS scheme, which uses an eight-class soil classification, resulted in large differences of up to −0.25 for the near-infrared (NIR) albedo over large parts of the desert regions of central Australia. The use of a recalibrated 20-class soil colour classification from the CLM, which includes a higher range for saturated and VIS (visible) and NIR soil albedos, reduced the underestimation of the NIR albedo. However, this soil colour mapping is tuned to CLM soil moisture, a quantity which is not necessarily transferrable between land surface models. We therefore recalibrated the soil color map using CABLE's climatological soil moisture, which further reduced the underestimation of the NIR albedo to within ±0.15 over most of the continent as compared to MODIS and SPOT albedos. Small areas of larger differences of up to −0.25 remained within the central arid parts of the continent during summer; however, the spatial extent of these large differences is substantially reduced as compared to the simulation using the default eight-class uncalibrated soil colour map. It is now possible to use CABLE coupled to atmospheric models to investigate soil-moisture–albedo feedbacks, an important enhancement of the model

Impacts of model calibration on high-latitude land-surface processes: PILPS 2(e) calibration/validation experiments

In the PILPS 2(e) experiment, the Snow Atmosphere Soil Transfer (SAST) land-surface scheme developed from the Biosphere-Atmosphere Transfer Scheme (BATS) showed difficulty in accurately simulating the patterns and quantities of runoff resulting from heavy snowmelt in the high-latitude Torne-Kalix River basin (shared by Sweden and Finland). This difficulty exposes the model deficiency in runoff formations. After representing subsurface runoff and calibrating the parameters, the accuracy of hydrograph prediction improved substantially. However, even with the accurate precipitation and runoff, the predicted soil moisture and its variation were highly "model-dependent". Knowledge obtained from the experiment is discussed. © 2003 Elsevier Science B.V. All rights reserved

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a 3-Dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of Solar System and exoplanetary terrestrial planets. Its parent model, known as ModelE2 (Schmidt et al. 2014), is used to simulate modern and 21st Century Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (slowly rotating to more rapidly rotating than modern Earth, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the Solar System such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn's moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents we can then expand its capabilities further to those exoplanetary rocky worlds that have been discovered in the past and those to be discovered in the future. We discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.Comment: Revisions since previous draft. Now submitted to Astrophysical Journal Supplement Serie

Evaluation of the Land Surface Water Budget in NCEP/NCAR and NCEP/DOE Reanalyses using an Off-line Hydrologic Model

The ability of the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) reanalysis (NRA1) and the follow-up NCEP/Department of Energy (DOE) reanalysis (NRA2), to reproduce the hydrologic budgets over the Mississippi River basin is evaluated using a macroscale hydrology model. This diagnosis is aided by a relatively unconstrained global climate simulation using the NCEP global spectral model, and a more highly constrained regional climate simulation using the NCEP regional spectral model, both employing the same land surface parameterization (LSP) as the reanalyses. The hydrology model is the variable infiltration capacity (VIC) model, which is forced by gridded observed precipitation and temperature. It reproduces observed streamflow, and by closure is constrained to balance other terms in the surface water and energy budgets. The VIC-simulated surface fluxes therefore provide a benchmark for evaluating the predictions from the reanalyses and the climate models. The comparisons, conducted for the 10-year period 1988–1997, show the well-known overestimation of summer precipitation in the southeastern Mississippi River basin, a consistent overestimation of evapotranspiration, and an underprediction of snow in NRA1. These biases are generally lower in NRA2, though a large overprediction of snow water equivalent exists. NRA1 is subject to errors in the surface water budget due to nudging of modeled soil moisture to an assumed climatology. The nudging and precipitation bias alone do not explain the consistent overprediction of evapotranspiration throughout the basin. Another source of error is the gravitational drainage term in the NCEP LSP, which produces the majority of the model\u27s reported runoff. This may contribute to an overprediction of persistence of surface water anomalies in much of the basin. Residual evapotranspiration inferred from an atmospheric balance of NRA1, which is more directly related to observed atmospheric variables, matches the VIC prediction much more closely than the coupled models. However, the persistence of the residual evapotranspiration is much less than is predicted by the hydrological model or the climate models

Do we (need to) care about canopy radiation schemes in DGVMs? Caveats and potential impacts

Dynamic global vegetation models (DGVMs) are an essential part of current state-of-the-art Earth system models. In recent years, the complexity of DGVMs has increased by incorporating new important processes like, e.g., nutrient cycling and land cover dynamics, while biogeophysical processes like surface radiation have not been developed much further. Canopy radiation models are however very important for the estimation of absorption and reflected fluxes and are essential for a proper estimation of surface carbon, energy and water fluxes. The present study provides an overview of current implementations of canopy radiation schemes in a couple of state-of-the-art DGVMs and assesses their accuracy in simulating canopy absorption and reflection for a variety of different surface conditions. Systematic deviations in surface albedo and fractions of absorbed photosynthetic active radiation (faPAR) are identified and potential impacts are assessed. The results show clear deviations for both, absorbed and reflected, surface solar radiation fluxes. FaPAR is typically underestimated, which results in an underestimation of gross primary productivity (GPP) for the investigated cases. The deviation can be as large as 25% in extreme cases. Deviations in surface albedo range between −0.15 ≤ Δα ≤ 0.36, with a slight positive bias on the order of Δα ≈ 0.04. Potential radiative forcing caused by albedo deviations is estimated at −1.25 ≤ RF ≤ −0.8 (W m−2), caused by neglect of the diurnal cycle of surface albedo. The present study is the first one that provides an assessment of canopy RT schemes in different currently used DGVMs together with an assessment of the potential impact of the identified deviations. The paper illustrates that there is a general need to improve the canopy radiation schemes in DGVMs and provides different perspectives for their improvement