Jeremy G. Fyke

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...

Future climate warming increases Greenland ice sheet surface mass balance variability

The integrated surface mass balance (SMB) of the Greenland ice sheet (GrIS) has large interannual variability. Long-term future changes to this variability will affect GrIS dynamics, freshwater fluxes, regional oceanography, and detection of changes in ice volume trends. Here we analyze a simulated 1850–2100 GrIS SMB time series from the Community Earth System Model, currently the only global climate model that realistically simulates GrIS SMB. We find a significant increase in interannual integrated SMB variability over time, which we attribute primarily to a shift to a high-variability melt-dominated SMB regime due to GrIS ablation area growth. We find temporal increases to characteristic ablation and accumulation area-specific SMB variabilities to be of secondary importance. Since ablation area SMB variability is driven largely by variability in summer surface melt, variability in the climate processes regulating the energy fluxes that control melting will likely increasingly determine future GrIS SMB variability.

Surface Melting over Ice Shelves and Ice Sheets as Assessed from Modeled Surface Air Temperatures

Abstract Summer surface melting plays an important role in the evolution of ice shelves and their progenitor ice sheets. To explore the magnitude of surface melt occurring over modern ice shelves and ice sheets in a climate scenario forced by anthropogenic emissions of carbon dioxide (CO2), a coupled climate model was used to simulate the distribution of summer melt at high latitudes and project the future evolution of high-melt regions in both hemispheres. Forcing of the climate model with CO2 emissions resulting from combustion of the present-day fossil-fuel resource base resulted in expansion of high-melt regions, as defined by the contour marking 200 positive degree-days per year, in the Northern Hemisphere and the Antarctic Peninsula and the introduction of high summer melt over the Ross, Ronne-Filchner, and Amery ice shelves as well as a large portion of the West Antarctic Ice Sheet (WAIS) and most of the Greenland Ice Sheet (GIS) by the year 2500. Capping CO2 concentrations at present-day levels av...

The Effect of Potential Future Climate Change on the Marine Methane Hydrate Stability Zone

Abstract The marine gas hydrate stability zone (GHSZ) is sensitive to temperature changes at the seafloor, which likely affected the GHSZ in the past and may do so in the future in response to anthropogenic greenhouse gas emissions. A series of climate sensitivity and potential future climate change experiments are undertaken using the University of Victoria Earth System Climate Model (UVic ESCM) with resulting seafloor temperature changes applied to a simple time-dependent methane hydrate stability model. The global GHSZ responds significantly to elevated atmospheric CO2 over time scales of 103 yr with initial decreases of the GHSZ occurring after 200 yr in shallow high-latitude seafloor areas that underlie regions of sea ice loss. The magnitude and rate of GHSZ change is dependent primarily upon the thermal diffusivity of the seafloor and the magnitude and duration of the seafloor temperature increase. Using a simple approximation of the amount of carbon stored as hydrate in the GHSZ, estimates of carbo...

Comment on "Climate sensitivity estimated from temperature reconstructions of the Last Glacial Maximum".

Schmittner et al. (Reports, 9 December 2011, p. 1385) report a new, low estimate of equilibrium climate sensitivity based on a comparison of Last Glacial Maximum climate model simulations and paleoproxy data. Here, we show that exclusion of questionable comparison points and constructive changes to model design are both likely capable of altering the most probable value of equilibrium climate sensitivity suggested in Schmittner et al.

The pattern of anthropogenic signal emergence in Greenland Ice Sheet surface mass balance

Surface mass balance (SMB) trends influence observed Greenland Ice Sheet (GrIS) mass loss, but the component of these trends related to anthropogenic forcing is unclear. Here we study the simulated spatial pattern of emergence of an anthropogenically derived GrIS SMB signal between 1850 and 2100 using the Community Earth System Model. We find emergence timing heterogeneity, with a bimodal structure reflecting interior snowfall increases against a background of low SMB variability, and peripheral surface melting increases against a backdrop of high SMB variability. We also find a nonemerging intermediate region. We conclude that (1) a bimodal pattern of GrIS SMB change will unambiguously reflect the impact of anthropogenic forcing; (2) present-day peripheral and interior SMB trends likely have an underlying anthropogenically forced component; (3) local emergence occurs well before emergence of a spatially integrated signal; and (4) the GrIS summit region may be an ideal location for monitoring regional/global climate change.

A probabilistic analysis of cumulative carbon emissions and long-term planetary warming

Efforts to mitigate and adapt to long-term climate change could benefit greatly from probabilistic estimates of cumulative carbon emissions due to fossil fuel burning and resulting CO2-induced planetary warming. Here we demonstrate the use of a reduced-form model to project these variables. We performed simulations using a large-ensemble framework with parametric uncertainty sampled to produce distributions of future cumulative emissions and consequent planetary warming. A hind-cast ensemble of simulations captured 1980–2012 historical CO2 emissions trends and an ensemble of future projection simulations generated a distribution of emission scenarios that qualitatively resembled the suite of Representative and Extended Concentration Pathways. The resulting cumulative carbon emission and temperature change distributions are characterized by 5–95th percentile ranges of 0.96–4.9 teratonnes C (Tt C) and 1.4 °C–8.5 °C, respectively, with 50th percentiles at 3.1 Tt C and 4.7 °C. Within the wide range of policy-related parameter combinations that produced these distributions, we found that low-emission simulations were characterized by both high carbon prices and low costs of non-fossil fuel energy sources, suggesting the importance of these two policy levers in particular for avoiding dangerous levels of climate warming. With this analysis we demonstrate a probabilistic approach to the challenge of identifying strategies for limiting cumulative carbon emissions and assessing likelihoods of surpassing dangerous temperature thresholds.

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

LIVVkit: An extensible, python-based, land ice verification and validation toolkit for ice sheet models: LAND ICE VERIFICATION AND VALIDATION

To address the pressing need to better understand the behavior and complex interaction of ice sheets within the global Earth system, significant development of continental-scale, dynamical ice sheet models is underway. Concurrent to the development of the Community Ice Sheet Model (CISM), the corresponding verification and validation (V&V) process is being coordinated through a new, robust, Python-based extensible software package, the Land Ice Verification and Validation toolkit (LIVVkit). Incorporated into the typical ice sheet model development cycle, it provides robust and automated numerical verification, software verification, performance validation, and physical validation analyses on a variety of platforms, from personal laptops to the largest supercomputers. LIVVkit operates on sets of regression test and reference data sets, and provides comparisons for a suite of community prioritized tests, including configuration and parameter variations, bit-for-bit evaluation, and plots of model variables to indicate where differences occur. LIVVkit also provides an easily extensible framework to incorporate and analyze results of new intercomparison projects, new observation data, and new computing platforms. LIVVkit is designed for quick adaptation to additional ice sheet models via abstraction of model specific code, functions, and configurations into an ice sheet model description bundle outside the main LIVVkit structure. Ultimately, through shareable and accessible analysis output, LIVVkit is intended to help developers build confidence in their models and enhance the credibility of ice sheet models overall.