Fuyuki Saito

Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume

The growth and reduction of Northern Hemisphere ice sheets over the past million years is dominated by an approximately 100,000-year periodicity and a sawtooth pattern (gradual growth and fast termination). Milankovitch theory proposes that summer insolation at high northern latitudes drives the glacial cycles, and statistical tests have demonstrated that the glacial cycles are indeed linked to eccentricity, obliquity and precession cycles. Yet insolation alone cannot explain the strong 100,000-year cycle, suggesting that internal climatic feedbacks may also be at work. Earlier conceptual models, for example, showed that glacial terminations are associated with the build-up of Northern Hemisphere ‘excess ice’, but the physical mechanisms underpinning the 100,000-year cycle remain unclear. Here we show, using comprehensive climate and ice-sheet models, that insolation and internal feedbacks between the climate, the ice sheets and the lithosphere–asthenosphere system explain the 100,000-year periodicity. The responses of equilibrium states of ice sheets to summer insolation show hysteresis, with the shape and position of the hysteresis loop playing a key part in determining the periodicities of glacial cycles. The hysteresis loop of the North American ice sheet is such that after inception of the ice sheet, its mass balance remains mostly positive through several precession cycles, whose amplitudes decrease towards an eccentricity minimum. The larger the ice sheet grows and extends towards lower latitudes, the smaller is the insolation required to make the mass balance negative. Therefore, once a large ice sheet is established, a moderate increase in insolation is sufficient to trigger a negative mass balance, leading to an almost complete retreat of the ice sheet within several thousand years. This fast retreat is governed mainly by rapid ablation due to the lowered surface elevation resulting from delayed isostatic rebound, which is the lithosphere–asthenosphere response. Carbon dioxide is involved, but is not determinative, in the evolution of the 100,000-year glacial cycles.

Insights into spatial sensitivities of ice mass response to environmental change from the SeaRISE ice sheet modeling project I: Antarctica

Sophie Nowicki, Robert A. Bindschadler, Ayako Abe-Ouchi, Andy Aschwanden, Ed Bueler, Hyeungu Choi, Jim Fastook, Glen Granzow, Ralf Greve, Gail Gutowski, Ute Herzfeld, Charles Jackson, Jesse Johnson, Constantine Khroulev, Eric Larour, Anders Levermann, William H. Lipscomb, Maria A. Martin, Mathieu Morlighem, Byron R. Parizek, David Pollard, Stephen F. Price, Diandong Ren, Eric Rignot, Fuyuki Saito, Tatsuru Sato, Hakime Seddik, Helene Seroussi, Kunio Takahashi, Ryan Walker, and Wei Li Wang

Initial results of the SeaRISE numerical experiments with the models SICOPOLIS and IcIES for the Greenland ice sheet

Abstarct SeaRISE (Sea-level Response to Ice Sheet Evolution) is a US-led multi-model community effort to predict the likely range of the contribution of the Greenland and Antarctic ice sheets to sea-level rise over the next few hundred years under global warming conditions. The Japanese ice-sheet modelling community is contributing to SeaRISE with two large-scale, dynamic/thermodynamic models: SICOPOLIS and IcIES. Here we discuss results for the Greenland ice sheet, obtained using both models under the forcings (surface temperature and precipitation scenarios) defined by the SeaRISE effort. A crucial point for meaningful simulations into the future is to obtain initial conditions that are close to the observed state of the present-day ice sheet. This is achieved by proper tuning during model spin-up from the last glacial/interglacial cycle to today. Experiments over 500 years indicate that both models are more sensitive (exhibit a larger rate of ice-sheet mass loss) to future climate warming (based on the A1B emission scenario) than to a doubling in the basal sliding speed. Ice-sheet mass loss varies between the two models by a factor of ~2 for sliding experiments and a factor of ~3 for climate-warming experiments, highlighting the importance of improved constraints on the parameterization of basal sliding and surface mass balance in ice-sheet models.

Results From the Ice-Sheet Model Intercomparison Project-Heinrich Event INtercOmparison (ISMIP HEINO)

Results from the Heinrich Event INtercOmparison (HEINO) topic of the Ice-Sheet Model Intercomparison Project (ISMIP) are presented. ISMIP HEINO was designed to explore internal large- scale ice-sheet instabilities in different contemporary ice-sheet models. These instabilities are of interest because they are a possible cause of Heinrich events. A simplified geometry experiment reproduces the main characteristics of the Laurentide ice sheet, including the sedimented region over Hudson Bay and Hudson Strait. The model experiments include a standard run plus seven variations. Nine dynamic/thermodynamic ice-sheet models were investigated; one of these models contains a combination of the shallow-shelf (SSA) and shallow-ice approximation (SIA), while the remaining eight models are of SIA type only. Seven models, including the SIA-SSA model, exhibit oscillatory surges with a period of ∼1000 years for a broad range of parameters, while two models remain in a permanent state of streaming for most parameter settings. In a number of models, the oscillations disappear for high surface temperatures, strong snowfall and small sediment sliding parameters. In turn, low surface temperatures and low snowfall are favourable for the ice-surge cycles. We conclude that further improvement of ice-sheet models is crucial for adequate, robust simulations of cyclic large-scale instabilities.

Effects of first-order stress gradients in an ice sheet evaluated by a three-dimensional thermomechanical coupled model

Abstract A three-dimensional thermomechanically coupled ice-sheet model including calculation of the first-order stress gradients is developed to perform numerical studies on the effects of stress components neglected in the zeroth-order shallow-ice approximation. Steady-state solutions are obtained with both a shallow-ice and a first-order ice-sheet model for several idealistic, radially symmetric ice-sheet configurations. The results show that the effects of the normal deviatoric-stress gradients on ice thickness are generally small, but the influence on basal temperatures and age profiles at an ice divide is significant.

Thermal structure of Dome Fuji and east Dronning Maud Land, Antarctica, simulated by a three-dimensional ice-sheet model

Abstract Three-dimensional structures of temperature focused on Dome Fuji and east Dronning Maud Land, Antarctica, simulated in a three-dimensional shallow ice model, are reported. With a geothermal heat flux of 54.6 mWm–2, as used in several modelling studies of the Antarctic ice sheet, and an enhancement factor of 1.3, which is smaller than in previous studies, the model result taking into account the glacial cycles is in good agreement with the borehole temperature and surface topography at Dome Fuji. The basal temperature at Dome Fuji must be at or very close to the pressure-melting point. The simulated amplitude of basal temperature through glacial/interglacial cycles is <1 K.

State dependence of climatic instability over the past 720,000 years from Antarctic ice cores and climate modeling

Global cooling in intermediate glacial climate with northern ice sheets preconditions climatic instability with bipolar seesaw. Climatic variabilities on millennial and longer time scales with a bipolar seesaw pattern have been documented in paleoclimatic records, but their frequencies, relationships with mean climatic state, and mechanisms remain unclear. Understanding the processes and sensitivities that underlie these changes will underpin better understanding of the climate system and projections of its future change. We investigate the long-term characteristics of climatic variability using a new ice-core record from Dome Fuji, East Antarctica, combined with an existing long record from the Dome C ice core. Antarctic warming events over the past 720,000 years are most frequent when the Antarctic temperature is slightly below average on orbital time scales, equivalent to an intermediate climate during glacial periods, whereas interglacial and fully glaciated climates are unfavourable for a millennial-scale bipolar seesaw. Numerical experiments using a fully coupled atmosphere-ocean general circulation model with freshwater hosing in the northern North Atlantic showed that climate becomes most unstable in intermediate glacial conditions associated with large changes in sea ice and the Atlantic Meridional Overturning Circulation. Model sensitivity experiments suggest that the prerequisite for the most frequent climate instability with bipolar seesaw pattern during the late Pleistocene era is associated with reduced atmospheric CO2 concentration via global cooling and sea ice formation in the North Atlantic, in addition to extended Northern Hemisphere ice sheets.

Sensitivity of Greenland ice sheet simulation to the numerical procedure employed for ice-sheet dynamics

Abstract The response of the Greenland ice sheet to global warming is simulated by two different numerical approaches, in order to evaluate the sensitivity of the analysis to the numerical structure employed. It is found that the thickness near the margin differs appreciably in these two simulations under identical conditions of modest warming, primarily due to a significant increase in the warming effect by an elevation–ablation feedback mechanism in one of the simulations. The change in ice-sheet volume differs by as much as a factor of two under strong climate-change forcing, demonstrating the need for care in interpreting the results of such climate-change analysis.

Modelled response of the volume and thickness of the Antarctic ice sheet to the advance of the grounded area

Abstract Numerical experiments are performed for the Antarctic ice sheet to study the sensitivity of the ice volume to variations in the area of grounded ice and to changes in the climate during the most recent deglaciation. The effect of the variations in the grounded area is found to be the major source of changes in the ice volume, while the effect of climate change was minor. The maximum possible contribution of the ice-volume change to sea-level rise during the deglaciation is estimated to be 36 m, which covers most values estimated in previous studies. The effect of the advance of the ice-sheet margin over those regions not connected to the major ice shelves contributes one-third of the total ice-volume change, which is comparable to the effect of the grounding of the Filchner–Ronne Ice Shelf and the contribution of the Ross and Amery Ice Shelves together.

Data exchange algorithm and software design of KAKUSHIN coupler Jcup

Abstract A coupler development program is on going as a part of Innovative Program of Climate Change projection for the 21st century (KAKUSHIN Program). The coupler is called Jcup. Basic functions of coupler are data exchange between models that have various grid systems and data conversion (interpolation) calculation. Among these functions, this article focuses on data exchange. Especially, multi-component coupling case that two or more models exist in single execution unit is regarded as the most important topic. What kind of coupling patterns are possible and how such a coupling patterns are realized on multi-component coupling is discussed in this article.