Berkeley code captures retreating Antarctic ice

April 3, 2013 in Earth / Earth Sciences
A simulation of the retreating Amundsen Sea Embayment with AMR. Green identifies regions that are modeled at 2km resolution, blue represents 1km, yellow is 500 meters and brown is 250 meters. The grounding line is shown as a red line. Credit: Antony Payne and Stephen Cornford, University of Bristol

Satellite observations suggest that the shrinking West Antarctic ice sheet is contributing to global sea level rise. But until recently, scientists could not accurately model the physical processes driving retreat of the ice sheet. Now, a new ice sheet model—called Berkeley-ISICLES (BISICLES)—is shedding light on these details.

Developed in a collaboration between at Lawrence Berkeley National Laboratory (Berkeley Lab), and at the Los Alamos National Laboratory and the University of Bristol, BISICLES allows researchers to model regions of interest—like the retreating edge of an ice-sheet—at sub-kilometer resolution, while employing computationally cheaper coarser resolution in areas that don't need such fine detail. This allows an accurate, high-resolution view of processes like glacier surges, ice streams, and grounding-line migration, at a relatively low computational cost.

This video shows one possible scenario for Antarctica's rapidly changing Amundsen Sea Embayment at two different model resolutions over a period of 218 years (from 1982 to 2200). In all three clips, the grounding line—where ice flowing into the sea separates from the bedrock and begins to float—is shown as a red line, and the speed of the flowing ice is depicted by brown coloring (the darker the brown, the faster the ice is flowing).

"When you watch these simulations side-by-side it is apparent that you need very fine resolution, better than one kilometer, to accurately simulate ice streams and grounding line migration," says Dan Martin, a Berkeley Lab computational scientist who helped develop BISICLES. Martin and his collaborators generated the simulations using the BISICLES code and up to 800 on the National Energy Research Scientific Computing Center's (NERSC's) Hopper system.

Martin worked with Gunther Weber of Berkeley Lab's Visualization Group to create this movie with VisIt from simulations generated by the University of Bristol researchers, Payne and Cornford. VisIt is a free interactive tool for visualizing scientific data. Weber's contribution to this project is supported by the Scalable Data Management, Analysis and Visualization (SDAV) project, which is funded by the Department of Energy's Scientific Discovery Through Advanced Computing (SciDAC) program.

Berkeley Lab contributes to BISICLES code

In its Fourth Assessment Report in 2007, the United Nations Intergovernmental Panel on Climate Change (IPCC) deemed contemporary state-of-the-art ice sheet models inadequate for making realistic projections about rise. In response, the Ice Sheet Initiative for Climate Extremes (ISICLES), funded by the DOE Advanced Scientific Computing Research Program, sought to leverage the expertise of computational scientists to help improve the models.

As a part of this program, computational scientists at Berkeley Lab (Dan Martin, Esmond Ng, Sam Williams, and Woo-Sun Yang) formed a collaboration with climate scientists at the Los Alamos National Laboratory (William Lipscomb) and the University of Bristol (Stephen Cornford and Tony Payne), which led to the development of a code called Berkeley-ISICLES (BISICLES).

Martin developed scalable Adaptive Mesh Refinement (AMR) ice sheet modeling tools, built on the Chombo software framework, and linked them to the existing Community Ice Sheet Model code (CISM). With AMR, researchers can now model dynamic points of interest at extremely high resolution, while maintaining a coarser resolution for areas that don't change as quickly. This enables accurate models at relatively low computational cost. The collaboration was so successful that the DOE continued supporting this work as part of PISCEES (Predicting Ice Sheet and Climate Evolution at Extreme Scales), a multi-institutional SciDAC Applications Partnership, which aims to further improve models.

More information:

Provided by Lawrence Berkeley National Laboratory

"Berkeley code captures retreating Antarctic ice" April 3, 2013