Comparing transient, accelerated, and equilibrium simulations of the last 30 000 years with the GENIE-1 model

• Main Authors: D. J. Lunt, M. S. Williamson, P. J. Valdes, T. M. Lenton, R. Marsh
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2006
• Subjects: Environmental pollution; TD172-193.5; Environmental protection; TD169-171.8; Environmental sciences; GE1-350; Southern Ocean; ice core; Ice Sheet; Sea ice
• Online Access: https://doaj.org/article/e1d25875221f4b99a6f7af77c293cc69

We examine several aspects of the ocean-atmosphere system over the last 30 000 years, by carrying out simulations with prescribed ice sheets, atmospheric CO 2 concentration, and orbital parameters. We use the GENIE-1 model with a frictional geostrophic ocean, dynamic sea ice, an energy balance atmosphere, and a land-surface scheme with fixed vegetation. A transient simulation, with boundary conditions derived from ice-core records and ice sheet reconstructions, is compared with equilibrium snapshot simulations, including the Last Glacial Maximum (21 000 years before present; 21 kyrBP), mid-Holocene (6 kyrBP) and pre-industrial. The equilibrium snapshot simulations are all very similar to their corresponding time period in the transient simulation, indicating that over the last 30 000 years, the model's ocean-atmosphere system is close to equilibrium with its boundary conditions. However, our simulations neglect the transfer of fresh water from and to the ocean, resulting from the growth and decay of ice sheets, which would, in reality, lead to greater disequilibrium. Additionally, the GENIE-1 model exhibits a rather limited response in terms of its Atlantic Meridional Overturning Circulation (AMOC) over the 30 000 years; a more sensitive AMOC would also be likely to lead to greater disequilibrium. We investigate the method of accelerating the boundary conditions of a transient simulation and find that the Southern Ocean is the region most affected by the acceleration. The Northern Hemisphere, even with a factor of 10 acceleration, is relatively unaffected. The results are robust to changes to several tunable parameters in the model. They also hold when a higher vertical resolution is used in the ocean.