William J. Sacks

Interactive Crop Management in the Community Earth System Model (CESM1): Seasonal Influences on Land–Atmosphere Fluxes

AbstractThe Community Earth System Model, version 1 (CESM1) is evaluated with two coupled atmosphere–land simulations. The CTRL (control) simulation represents crops as unmanaged grasses, while CROP represents a crop managed simulation that includes special algorithms for midlatitude corn, soybean, and cereal phenology and carbon allocation. CROP has a more realistic leaf area index (LAI) for crops than CTRL. CROP reduces winter LAI and represents the spring planting and fall harvest explicitly. At the peak of the growing season, CROP simulates higher crop LAI. These changes generally reduce the latent heat flux but not around peak LAI (late spring/early summer). In midwestern North America, where corn, soybean, and cereal abundance is high, simulated peak summer precipitation declines and agrees better with observations, particularly when crops emerge late as is found from a late planting sensitivity simulation (LateP). Differences between the CROP and LateP simulations underscore the importance of simul...

Implementation and Initial Evaluation of the Glimmer Community Ice Sheet Model in the Community Earth System Model

AbstractThe Glimmer Community Ice Sheet Model (Glimmer-CISM) has been implemented in the Community Earth System Model (CESM). Glimmer-CISM is forced by a surface mass balance (SMB) computed in multiple elevation classes in the CESM land model and downscaled to the ice sheet grid. Ice sheet evolution is governed by the shallow-ice approximation with thermomechanical coupling and basal sliding. This paper describes and evaluates the initial model implementation for the Greenland Ice Sheet (GIS). The ice sheet model was spun up using the SMB from a coupled CESM simulation with preindustrial forcing. The models sensitivity to three key ice sheet parameters was explored by running an ensemble of 100 GIS simulations to quasi equilibrium and ranking each simulation based on multiple diagnostics. With reasonable parameter choices, the steady-state GIS geometry is broadly consistent with observations. The simulated ice sheet is too thick and extensive, however, in some marginal regions where the SMB is anomalousl...

Integration of process-based soil respiration models with whole-ecosystem CO2 measurements

We integrated soil models with an established ecosystem process model (SIPNET, simplified photosynthesis and evapotranspiration model) to investigate the influence of soil processes on modelled values of soil CO2 fluxes (RSoil). Model parameters were determined from literature values and a data assimilation routine that used a 7-year record of the net ecosystem exchange of CO2 and environmental variables collected at a high-elevation subalpine forest (the Niwot Ridge AmeriFlux site). These soil models were subsequently evaluated in how they estimated the seasonal contribution of RSoil to total ecosystem respiration (TER) and the seasonal contribution of root respiration (RRoot) to RSoil. Additionally, these soil models were compared to data assimilation output of linear models of soil heterotrophic respiration. Explicit modelling of root dynamics led to better agreement with literature values of the contribution of RSoil to TER. Estimates of RSoil/TER when root dynamics were considered ranged from 0.3 to 0.6; without modelling root biomass dynamics these values were 0.1–0.3. Hence, we conclude that modelling of root biomass dynamics is critically important to model the RSoil/TER ratio correctly. When soil heterotrophic respiration was dependent on linear functions of temperature and moisture independent of soil carbon pool size, worse model-data fits were produced. Adding additional complexity to the soil pool marginally improved the model-data fit from the base model, but issues remained. The soil models were not successful in modelling RRoot/RSoil. This is partially attributable to estimated turnover parameters of soil carbon pools not agreeing with expected values from literature and being poorly constrained by the parameter estimation routine. We conclude that net ecosystem exchange of CO2 alone cannot constrain specific rhizospheric and microbial components of soil respiration. Reasons for this include inability of the data assimilation routine to constrain soil parameters using ecosystem CO2 flux measurements and not considering the effect of other resource limitations (for example, nitrogen) on the microbe biomass. Future data assimilation studies with these models should include ecosystem-scale measurements of RSoil in the parameter estimation routine and experimentally determine soil model parameters not constrained by the parameter estimation routine.

Greenland Surface Mass Balance as Simulated by the Community Earth System Model. Part I: Model Evaluation and 1850–2005 Results

AbstractThe modeling of the surface mass balance (SMB) of the Greenland Ice Sheet (GIS) requires high-resolution models in order to capture the observed large gradients in the steep marginal areas. Until now, global climate models have not been considered suitable to model ice sheet SMB owing to model biases and insufficient resolution. This study analyzes the GIS SMB simulated for the period 1850–2005 by the Community Earth System Model (CESM), which includes a new ice sheet component with multiple elevation classes for SMB calculations. The model is evaluated against observational data and output from the regional model Regional Atmospheric Climate Model version 2 (RACMO2). Because of a lack of major climate biases, a sophisticated calculation of snow processes (including surface albedo evolution) and an adequate downscaling technique, CESM is able to realistically simulate GIS surface climate and SMB. CESM SMB agrees reasonably well with in situ data from 475 locations (r = 0.80) and output from RACMO2...

A Multiscale and Multidisciplinary Investigation of Ecosystem- Atmosphere CO2 Exchange over the Rocky Mountains of Colorado

A significant fraction of Earth consists of mountainous terrain. However, the question of how to monitor the surface–atmosphere carbon exchange over complex terrain has not been fully explored. This article reports on studies by a team of investigators from U.S. universities and research institutes who carried out a multiscale and multidisciplinary field and modeling investigation of the CO2 exchange between ecosystems and the atmosphere and of CO2 transport over complex mountainous terrain in the Rocky Mountain region of Colorado. The goals of the field campaign, which included ground and airborne in situ and remote-sensing measurements, were to characterize unique features of the local CO2 exchange and to find effective methods to measure regional ecosystem–atmosphere CO2 exchange over complex terrain. The modeling effort included atmospheric and ecological numerical modeling and data assimilation to investigate regional CO2 transport and biological processes involved in ecosystem–atmosphere carbon exch...

Improving the Representation of Polar Snow and Firn in the Community Earth System Model

In Earth system models, terrestrial snow is usually modeled by the land surface component. In most cases, these snow models have been developed with an emphasis on seasonal snow. Questions about future sea level rise, however, prompt the need for a realistic representation of perennial snow, as snow processes play a key role in the mass balance of glaciers and ice sheets. Here we enhance realism of modeled polar snow in the Community Land Model (CLM), the land component of the Community Earth System Model (CESM), by implementing (1) new parametrizations for fresh snow density, destructive metamorphism, and compaction by overburden pressure, (2) by allowing for deeper snow packs, and (3) by introducing drifting snow compaction, with a focus on the ice sheet interior. Comparison with Greenlandic and Antarctic snow density observations show that the new physics improve model skill in predicting firn and near-surface density in the absence of melt. Moreover, compensating biases are removed and spurious subsurface melt rates at ice sheets are eliminated. The deeper snow pack enhances refreezing and allows for deeper percolation, raising ice temperatures up to 15 oC above the skin temperature.

Description and Evaluation of the Community Ice Sheet Model(CISM) v2.1

We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3D thermomechanical model, written mainly in Fortran 90/95, that solves equations for the momentum balance and thickness and temperature evolution of ice sheets. CISM’s velocity solver incorporates a hierarchy of Stokes-flow approximations, including shallow-shelf, depth-integrated higher-order, and 3D higher-order. CISM also includes a suite of test cases, links to third-party solver libraries, and parameterizations of physical processes such as basal sliding and iceberg calving. The model has been 5 verified for standard test problems, including the ISMIP-HOM experiments for higher-order models, and has participated in the initMIP–Greenland initialization experiment. In multi-millennial simulations with modern climate forcing on a 4-km grid, CISM reaches a steady state that is broadly consistent with observed flow patterns of the Greenland ice sheet. CISM has been integrated into version 2.0 of the Community Earth System Model, where it is being used for Greenland simulations under past, present and future climates. The code is open-source with extensive documentation, and remains under active development. 10