Interactive ocean bathymetry and coastlines for simulating the last deglaciation with the Max Planck Institute Earth System Model (MPI-ESM-v1.2)
As ice sheets grow or decay, the net flux of freshwater into the ocean changes and the bedrock adjusts due to isostatic adjustments, leading to variations in the bottom topography and the oceanic boundaries. This process was particularly intense during the last deglaciation due to the high rates of...
| Published in: | Geoscientific Model Development |
|---|---|
| Main Authors: | , |
| Format: | Article in Journal/Newspaper |
| Language: | English |
| Published: |
Copernicus Publications
2018
|
| Subjects: | |
| Online Access: | https://doi.org/10.5194/gmd-11-4677-2018 https://doaj.org/article/9778615cf54d468a8cdf43d4618fb0a9 |
| Summary: | As ice sheets grow or decay, the net flux of freshwater into the ocean changes and the bedrock adjusts due to isostatic adjustments, leading to variations in the bottom topography and the oceanic boundaries. This process was particularly intense during the last deglaciation due to the high rates of ice-sheet melting. It is, therefore, necessary to consider transient ocean bathymetry and coastlines when attempting to simulate the last deglaciation with Earth system models (ESMs). However, in most standard ESMs the land-sea mask is fixed throughout simulations because the generation of a new ocean model bathymetry implies several levels of manual corrections, a procedure that is hardly doable very often for long runs. This is one of the main technical problems towards simulating a complete glacial cycle with general circulation models. For the first time, we present a tool allowing for an automatic computation of bathymetry and land-sea mask changes in the Max Planck Institute Earth System Model (MPI-ESM). The algorithms developed in this paper can easily be adapted to any free-surface ocean model that uses the Arakawa-C grid in the horizontal and z -grid in the vertical including partial bottom cells. The strategy applied is described in detail and the algorithms are tested in a long-term simulation demonstrating the reliable behaviour. Our approach guarantees the conservation of mass and tracers at global and regional scales; that is, changes in a single grid point are only propagated regionally. The procedures presented here are an important contribution to the development of a fully coupled ice sheet–solid Earth–climate model system with time-varying topography and will allow for transient simulations of the last deglaciation considering interactive bathymetry and land-sea mask. |
|---|