Sensitivity of last glacial maximum climate to sea ice conditions in the Nordic Seas

Published reconstructions of last glacial maximum (LGM) sea surface temperatures and sea ice extent differ significantly. We here test the sensitivity of simulated North Atlantic climates to two different reconstructions by using these reconstructions as boundary conditions for model experiments. An atmospheric general circulation model has been used to perform two simulations of the (LGM) and a modern-day control simulation. Standard (CLIMAP) reconstructions of sea ice and sea surface temperatures have been used for the first simulation, and a set of new reconstructions in the Nordic Seas/Northern Atlantic have been used for the second experiment. The new reconstruction is based on 158 core samples, and represents ice-free conditions during summer in the Nordic Seas, with accordingly warmer sea surface temperatures and less extensive sea ice during winter as well. The simulated glacial climate is globally 5.7 K colder than modern day, with the largest changes at mid and high latitudes. Due to more intense Hadley circulation, the precipitation at lower latitudes has increased in the simulations of the LGM. Relative to the simulation with the standard CLIMAP reconstructions, reduction of the sea ice in the North Atlantic gives positive local responses in temperature, precipitation and reduction of the sea level pressure. Only very weak signatures of the wintertime Icelandic Low occur when the standard CLIMAP sea surface temperature reconstruction is used as the lower boundary condition in LGM. With reduced sea ice conditions in the Nordic Seas, the Icelandic Low becomes more intense and closer to its present structure. This indicates that thermal forcing is an important factor in determining the strength and position of the Icelandic Low. The Arctic Oscillation is the most dominant large scale variability feature on the Northern Hemisphere in modern day winter climate. In the simulation of the LGM with extensive sea ice this pattern is significantly changed and represents no systematic large scale variability over the North Atlantic. Reduction of the North Atlantic sea ice extent leads to stronger variability in monthly mean sea level pressure in winter. The synoptic variability appears at a lower level in the simulation when standard reconstructions of the sea surface in the LGM are used. A closer inspection of storm tracks in this model experiment shows that that the synoptic lows follow a narrow band along the ice edge during winter. The trajectories of synoptic lows are not constrained to the sea ice edge to the same degree when the sea ice extent is reduced. Seasonally open waters in the Nordic Seas in the new reconstruction apparently act as a moisture source, consistent with the current understanding of the rapid growth of the Fennoscandian and Barents Ice Sheets, during the LGM. The signal from the intensified thermal forcing in the North Atlantic in Boreal winter is carried zonally by upper tropospheric waves, and thus generates non-local responses to the changed sea ice cover.     

Transient simulation of the last glacial inception. Part I: glacial inception as a bifurcation in the climate system

We study the mechanisms of glacial inception by using the Earth system model of intermediate complexity, CLIMBER-2, which encompasses dynamic modules of the atmosphere, ocean, biosphere and ice sheets. Ice-sheet dynamics are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. We have performed transient experiments starting at the Eemiam interglacial, at 126 ky BP (126,000 years before present). The model runs for 26 kyr with time-dependent orbital and CO₂ forcings. The model simulates a rapid expansion of the area covered by inland ice in the Northern Hemisphere, predominantly over Northern America, starting at about 117 kyr BP. During the next 7 kyr, the ice volume grows gradually in the model at a rate which corresponds to a change in sea level of 10 m per millennium. We have shown that the simulated glacial inception represents a bifurcation transition in the climate system from an interglacial to a glacial state caused by the strong snow-albedo feedback. This transition occurs when summer insolation at high latitudes of the Northern Hemisphere drops below a threshold value, which is only slightly lower than modern summer insolation. By performing long-term equilibrium runs, we find that for the present-day orbital parameters at least two different equilibrium states of the climate system exist—the glacial and the interglacial; however, for the low summer insolation corresponding to 115 kyr BP, we find only one, glacial, equilibrium state, while for the high summer insolation corresponding to 126 kyr BP only an interglacial state exists in the model.

---

Transient simulation of the last glacial inception. Part II: sensitivity and feedback analysis

The sensitivity of the last glacial-inception (around 115 kyr BP, 115,000 years before present) to different feedback mechanisms has been analysed by using the Earth system model of intermediate complexity CLIMBER-2. CLIMBER-2 includes dynamic modules of the atmosphere, ocean, terrestrial biosphere and inland ice, the last of which was added recently by utilising the three-dimensonal polythermal ice-sheet model SICOPOLIS. We performed a set of transient experiments starting at the middle of the Eemiam interglacial and ran the model for 26,000 years with time-dependent orbital forcing and observed changes in atmospheric CO₂ concentration (CO₂ forcing). The role of vegetation and ocean feedback, CO₂ forcing, mineral dust, thermohaline circulation and orbital insolation were closely investigated. In our model, glacial inception, as a bifurcation in the climate system, appears in nearly all sensitivity runs including a run with constant atmospheric CO₂ concentration of 280 ppmv, a typical interglacial value, and simulations with prescribed present-day sea-surface temperatures or vegetation cover—although the rate of the growth of ice-sheets growth is smaller than in the case of the fully interactive model. Only if we run the fully interactive model with constant present-day insolation and apply present-day CO₂ forcing does no glacial inception appear at all. This implies that, within our model, the orbital forcing alone is sufficient to trigger the interglacial–glacial transition, while vegetation, ocean and atmospheric CO₂ concentration only provide additional, although important, positive feedbacks. In addition, we found that possible reorganisations of the thermohaline circulation influence the distribution of inland ice.     

On the causes of glacial inception at 116 kaBP

To explore processes involved in glacial inception at 116 kaBP, the response of an atmospheric general circulation model (AGCM) to changes in lower boundary conditions is investigated. Two 116 kaBP experiments are conducted to examine the importance of sea surface conditions (sea surface temperature and sea ice distribution): one with the present-day sea surface conditions, and the other with 116 kaBP sea surface conditions. These two different sea surface conditions are obtained from simulations using an earth system climate model of intermediate complexity. Perennial snow cover occurred over the Canadian Archipelago under 116 kaBP orbital and CO₂ forcing with present-day "warm" sea surface conditions, and further expanded over northeastern Canada when 116 kaBP "cool" sea surface conditions were applied. The net positive accumulation in northeastern Canada, with little in Alaska, is in good agreement with geological records. Two additional 116 kaBP experiments are conducted to examine the combined importance of sea surface conditions and land surface conditions (vegetation): one with the present-day sea surface and modified land surface conditions, and the other with 116 kaBP sea surface and modified land surface conditions. Modifying vegetation, based on cooling during summer induced by 116 kaBP sea surface conditions, leads to much larger areas of perennial snow cover. Only when 116 kaBP sea surface conditions are applied, is a realistic global net snow accumulation rate obtained. Contrary to the earlier ice age hypothesis, our results suggest that the capturing of glacial inception at 116 kaBP requires the use of "cooler" sea surface conditions than those of the present climate. Also, the large impact of vegetation change on climate suggests that the inclusion of the vegetation feedback is important for model validation, at least, in this particular period of Earth history.     

Carbon cycle variations on the Earth's surface during the last two glacial cycles

Data concerning carbon cycle variations on the earth's surface during the past 200,000 years are reviewed.The variations of the surface temperature (T) and concentration of carbon dioxide (CO₂) in the atmosphere of Antarctica are compared to those of the isotopic ratios of oxygen ¹⁸O/¹⁶O (δ¹⁸O) and of carbon ¹³C/¹²C (°¹³C) of waters in the deep oceans for the two last glacial cycles. This comparison shows that the decrease of the atmospheric CO₂ concentration is accompanied by a carbon transferase from the continental biosphere to the oceanic deep waters. At the glacial maximum this transfer is estimated to be about 500 GtC (1 GtC = 10¹⁵g of carbon) equivalent to 25% of the carbon storage of the biosphere. It occurs mainly in the high latitudes of the Southern Hemisphere by incorporation of CO₂ into particulate matter during photosynthesis. It is shown that the mean oceanic productivity does not increase with a supplementary supply of ions such as phosphate (PO₄³⁻) or nitrate (NO₃⁻) but that the intensity of the thermohaline circulation is certainly reduced. As the warming up of the oceans and the melting of the ice-sheet begin carbon transfer takes place to restore the continental biosphere.Another carbon transfer of a much more important intensity is also at work in the sea shore environment. Its intensity could be sufficient to renew the entire carbon of the continental biospheric, atmospheric and oceanic reservoirs in a length of time comparable to a glacial cycle. This fact shows the importance of studying the deposition of carbon in oceanic zones which are uncovered with drops in sea level. At the present time data on the coastal environment in relation to the global carbon cycle are very scarce and warrants more research in this area.     

The role of land surface dynamics in glacial inception: a study with the UVic Earth System Model

The first results of the UVic Earth System Model coupled to a land surface scheme and a dynamic global vegetation model are presented in this study. In the first part the present day climate simulation is discussed and compared to observations. We then compare a simulation of an ice age inception (forced with 116 ka BP orbital parameters and an atmospheric CO₂ concentration of 240 ppm) with a preindustrial run (present day orbital parameters, atmospheric [CO₂] = 280 ppm). Emphasis is placed on the vegetations response to the combined changes in solar radiation and atmospheric CO₂ level. A southward shift of the northern treeline as well as a global decrease in vegetation carbon is observed in the ice age inception run. In tropical regions, up to 88% of broadleaf trees are replaced by shrubs and C₄ grasses. These changes in vegetation cover have a remarkable effect on the global climate: land related feedbacks double the atmospheric cooling during the ice age inception as well as the reduction of the meridional overturning in the North Atlantic. The introduction of vegetation related feedbacks also increases the surface area with perennial snow significantly.     

High-resolution climate simulation of the last glacial maximum

The climate of the last glacial maximum (LGM) is simulated with a high-resolution atmospheric general circulation model, the NCAR CCM3 at spectral truncation of T170, corresponding to a grid cell size of roughly 75 km. The purpose of the study is to assess whether there are significant benefits from the higher resolution simulation compared to the lower resolution simulation associated with the role of topography. The LGM simulations were forced with modified CLIMAP sea ice distribution and sea surface temperatures (SST) reduced by 1°C, ice sheet topography, reduced CO₂, and 21,000 BP orbital parameters. The high-resolution model captures modern climate reasonably well, in particular the distribution of heavy precipitation in the tropical Pacific. For the ice age case, surface temperature simulated by the high-resolution model agrees better with those of proxy estimates than does the low-resolution model. Despite the fact that tropical SSTs were only 2.1°C less than the control run, there are many lowland tropical land areas 4–6°C colder than present. Comparison of T170 model results with the best constrained proxy temperature estimates (noble gas concentrations in groundwater) now yield no significant differences between model and observations. There are also significant upland temperature changes in the best resolved tropical mountain belt (the Andes). We provisionally attribute this result in part as resulting from decreased lateral mixing between ocean and land in a model with more model grid cells. A longstanding model-data discrepancy therefore appears to be resolved without invoking any unusual model physics. The response of the Asian summer monsoon can also be more clearly linked to local geography in the high-resolution model than in the low-resolution model; this distinction should enable more confident validation of climate proxy data with the high-resolution model. Elsewhere, an inferred salinity increase in the subtropical North Atlantic may have significant implications for ocean circulation changes during the LGM. A large part of the Amazon and Congo Basins are simulated to be substantially drier in the ice age—consistent with many (but not all) paleo data. These results suggest that there are considerable benefits derived from high-resolution model regarding regional climate responses, and that observationalists can now compare their results with models that resolve geography at a resolution comparable to that which the proxy data represent.     

West African palaeosynoptic patterns at the last glacial maximum

Summary The composite method and an atmospheric moisture budget is applied to study the present synoptic summer situation in the West African Sahel. This allows the information of the synoptic scale flow systems at the last glacial maximum to be obtained. Relying on the results of general circulation models and of local geological findings the palaeosynoptic situation was found to consist of the same constituents as today's given pattern. In both time slices the strength of the midtropospheric African Easterly Jet determines the intensity and frequency of the Sahelian synoptic disturbances: the Easterly Waves and the Squall Lines. The glacial rainfall intensity of the Squall Lines, the most effective rain-bearing system of the Sahel, is calculated to be 2/3 of the modern one. Additionally, only half the number of Squall Lines propagated across the West African continent per summer month. Thus the glacial rainfall amount of the Sahel was roughly 30% of the value observed today.With 6 Figures     

Paleoclimate simulations of the mid-Holocene and last glacial maximum by FGOALS

Paleoclimate simulations of the mid-Holocene (MH) and Last Glacial maximum (LGM) by the latest versions of the Flexible Global Ocean-Atmosphere-Land System model, Spectral Version 2 and Grid-point Version 2 (FGOALS-s2 and g2) are evaluated in this study. The MH is characterized by changes of insolation induced by orbital parameters, and the LGM is a glacial period with large changes in greenhouse gases, sea level and ice sheets. For the MH, both versions of FGOALS simulate reasonable responses to the changes of insolation, such as the enhanced summer monsoon in African-Asian regions. Model differences can be identified at regional and seasonal scales. The global annual mean surface air temperature (TAS) shows no significant change in FGOALS-s2, while FGOALS-g2 shows a global cooling of about 0.7 C that is related with a strong cooling during boreal winter. The amplitude of ENSO is weaker in FGOALS-g2, which agrees with proxy data. For the LGM, FGOALS-g2 captures the features of the cold and dry glacial climate, including a global cooling of 4.6 C and a decrease in precipitation by 10%. The ENSO is weaker at the LGM, with a tendency of stronger ENSO cold events. Sensitivity analysis shows that the Equilibrium Climate Sensitivity (ECS) estimated for FGOALS ranges between 4.23 C and 4.59 C. The sensitivity of precipitation to the changes of TAS is～2.3% C-1 , which agrees with previous studies. FGOALS-g2 shows better simulations of the Atlantic Meridional Overturning Circulation (AMOC) and African summer monsoon precipitation in the MH when compared with FGOALS-g1.0; however, it is hard to conclude any improvements for the LGM.     

Sources of continental dust over Antarctica during the last glacial cycle

The soluble and insoluble parts of 4 major components (Al, Ca, K and Mg) of the continental dust input over East Antarctica, as well as size, distribution parameters of the insoluble part of this dust, have been studied along an ice core which spanns the last climatic cycle (160 kyr). These results provide a better understanding of the respective impact of the different potential dust sources. While Al and K were probably entrapped in illite originating from arid areas and in a lesser extent from shallow marine sediments, Ca and Mg inputs were dominated by marine carbonate of exposed continental shelves emissions.     

Forecasts of hydrological parameters over the Mackenzie River Basin: Sensitivity to initial conditions,horizontal resolution and forecast range

Abstract

As part of the Global Energy and Water Cycle Experiment, Canadian global spectral forecast model predictions of surface water and energy fluxes over the Mackenzie River basin are examined. Two nine‐member ensemble forecasts of one month duration are produced with the operational model, for a spring and a summer case, at a horizontal resolution of about 100 km (T95). The sensitivity to initial conditions is measured by the degree to which the individual forecasts in the ensembles vary one from another. The evolution in time of this estimated error (ensemble standard deviation) is determined for the surface energy and water accumulations, averaged over the basin. For comparison the calculations are repeated for the Mississippi basin and over North America. The greatest sensitivity is found for the net accumulation of precipitation minus evaporation. The spring ensemble is redone at a coarser horizontal resolution (T47), and the results are similar. The forecast uncertainty (ensemble standard deviation) of the area‐averages over the basin appear to be unaffected by this change, although the ensemble mean values are sensitive to the change in resolution. The ensemble standard deviation makes a significant, abrupt increase toward the end of the second week into the forecasts. This investigation suggests a need for an improved model, if the forecasts’ useful range is to extend to one month. Available upgrades to the land‐surface, precipitation and evaporation schemes will be used in subsequent work, and the forecasts reported here will serve as a baseline for comparison.     

Atmospheric equilibrium,instability and energy transport at the last glacial maximum

A range of diagnostics from two GCM simulations, one of the present-day climate and one of the last glacial maximum (LGM) is used to gain insight into their different temperature structures and eddy dynamics. There are large local increases in baroclinicity at the LGM, especially in the Atlantic storm track, with large accompanying increases in the low level transient eddy heat flux. However, the differences in the zonal mean are much smaller, and the increases in both baroclinicity and heat flux are confined to low levels. Supplementary experiments with baroclinic wave lifecycles confirm the marked contrast between local and zonal mean behaviour, but do not adequately explain the differences between the zonal mean climates. The total flux of energy across latitude circles in the Northern Hemisphere does not change much during DJF, although its transient component is actually reduced at the LGM (during JJA the transient component is increased). Calculations of total linear eddy diffusivity reveal that changes in the time mean stationary waves are chiefly responsible for the seasonal range of this quantity at the LGM, while they only account for half the seasonal range at the present-day.     

A comparison of PMIP2 model simulations and the MARGO proxy reconstruction for tropical sea surface temperatures at last glacial maximum

Results from multiple model simulations are used to understand the tropical sea surface temperature (SST) response to the reduced greenhouse gas concentrations and large continental ice sheets of the last glacial maximum (LGM). We present LGM simulations from the Paleoclimate Modelling Intercomparison Project, Phase 2 (PMIP2) and compare these simulations to proxy data collated and harmonized within the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface Project (MARGO). Five atmosphere–ocean coupled climate models (AOGCMs) and one coupled model of intermediate complexity have PMIP2 ocean results available for LGM. The models give a range of tropical (defined for this paper as 15°S–15°N) SST cooling of 1.0–2.4°C, comparable to the MARGO estimate of annual cooling of 1.7 ± 1°C. The models simulate greater SST cooling in the tropical Atlantic than tropical Pacific, but interbasin and intrabasin variations of cooling are much smaller than those found in the MARGO reconstruction. The simulated tropical coolings are relatively insensitive to season, a feature also present in the MARGO transferred-based estimates calculated from planktonic foraminiferal assemblages for the Indian and Pacific Oceans. These assemblages indicate seasonality in cooling in the Atlantic basin, with greater cooling in northern summer than northern winter, not captured by the model simulations. Biases in the simulations of the tropical upwelling and thermocline found in the preindustrial control simulations remain for the LGM simulations and are partly responsible for the more homogeneous spatial and temporal LGM tropical cooling simulated by the models. The PMIP2 LGM simulations give estimates for the climate sensitivity parameter of 0.67°–0.83°C per Wm⁻², which translates to equilibrium climate sensitivity for doubling of atmospheric CO₂ of 2.6–3.1°C.     

Simulation of long-term future climate changes with the green McGill paleoclimate model: the next glacial inception

The multi-component “green” McGill Paleoclimate Model (MPM), which includes interactive vegetation, is used to simulate the next glacial inception under orbital and prescribed atmospheric CO₂ forcing. This intermediate complexity model is first run for short-term periods with an increasing atmospheric CO₂ concentration; the model's response is in general agreement with the results of GCMs for CO₂ doubling. The green MPM is then used to derive projections of the climate for the next 100 kyr. Under a constant CO₂ level, the model produces three types of evolution for the ice volume: an imminent glacial inception (low CO₂ levels), a glacial inception in 50 kyr (CO₂ levels of 280 or 290 ppm), or no glacial inception during the next 100 kyr (CO₂ levels of 300 ppm and higher). This high sensitivity to the CO₂ level is due to the exceptionally weak future variations of the summer insolation at high northern latitudes. The changes in vegetation re-inforce the buildup of ice sheets after glacial inception. Finally, if an initial global warming episode of finite duration is included, after which the atmospheric CO₂ level is assumed to stabilize at 280, 290 or 300 ppm, the impact of this warming is seen only in the first 5 kyr of the run; after this time the response is insensitive to the early warming perturbation.     

Surface temperature control in the North and tropical Pacific during the last glacial maximum

Based on coupled modelling evidence we argue that topographically-induced modifications of the large-scale atmospheric circulation during the last glacial maximum may have led to a reduction of the westerlies, and a slowdown of the Pacific subtropical gyre as well as to an intensification of the Pacific subtropical cell. These oceanic circulation changes generate an eastern North Pacific warming, an associated cooling in the Kuroshio area, as well as a cooling of the tropical oceans, respectively. The tropical cooling pattern resembles a permanent La Niña state which in turn forces atmospheric teleconnection patterns that lead to an enhancement of the subtropical warming by reduced latent and sensible cooling of the ocean. In addition, the radiative cooling due to atmospheric CO₂ and water vapor reductions imposes a cooling tendency in the tropics and subtropics, thereby intensifying the permanent La Niña conditions. The remote North Pacific response results in a warming tendency of the eastern North Pacific which may level off the effect of the local radiative cooling. Hence, a delicate balance between oceanic circulation changes, remotely induced atmospheric flux anomalies as well local radiative cooling is established which controls the tropical and North Pacific temperature anomalies during the last glacial maximum. Furthermore, we discuss how the aftermath of a Heinrich event may have affected glacial temperatures in the Pacific Ocean.     

On the relationship between ice volume and sea level over the last glacial cycle

Geological evidence and oxygen-isotope variations in deep-sea cores provide valuable information about the sea-level variations of the past. Ice-volume equivalent is usually computed by using a constant oceanic area. In this paper a relationship is developed between the continental ice-volume variation and the sea-level drop by taking into account the sea-floor topography and, therefore, the variation of the oceanic area. It appears from such calculations that the last glacial maximum ice volume is 7% less than previously estimated, and that the minimal reconstruction of the ice sheets from Hughes et al. (1981) seems the most likely.     

The Southern Westerlies during the last glacial maximum in PMIP2 simulations

The Southern Hemisphere westerly winds are an important component of the climate system at hemispheric and global scales. Variations in their intensity and latitudinal position through an ice-age cycle have been proposed as important drivers of global climate change due to their influence on deep-ocean circulation and changes in atmospheric CO₂. The position, intensity, and associated climatology of the southern westerlies during the last glacial maximum (LGM), however, is still poorly understood from empirical and modelling standpoints. Here we analyse the behaviour of the southern westerlies during the LGM using four coupled ocean-atmosphere simulations carried out by the Palaeoclimate Modelling Intercomparison Project Phase 2 (PMIP2). We analysed the atmospheric circulation by direct inspection of the winds and by using a cyclone tracking software to indicate storm tracks. The models suggest that changes were most significant during winter and over the Pacific ocean. For this season and region, three out four models indicate decreased wind intensities at the near surface as well as in the upper troposphere. Although the LGM atmosphere is colder and the equator to pole surface temperature gradient generally increases, the tropospheric temperature gradients actually decrease, explaining the weaker circulation. We evaluated the atmospheric influence on the Southern Ocean by examining the effect of wind stress on the Ekman pumping. Again, three of the models indicate decreased upwelling in a latitudinal band over the Southern Ocean. All models indicate a drier LGM than at present with a clear decrease in precipitation south of 40°S over the oceans. We identify important differences in precipitation anomalies over the land masses at regional scale, including a drier climate over New Zealand and wetter over NW Patagonia.     

Evolution and modulation of tropical heating from the last glacial maximum through the twenty-first century

Twentieth century observations show that during the last 50?years the sea-surface temperature (SST) of the tropical oceans has increased by ~0.5°C and the area of SST >26.5 and 28°C (arbitrarily referred to as the oceanic warm pool: OWP) by 15 and 50% respectively in association with an increase in green house gas concentrations, with non-understood natural variability or a combination of both. Based on CMIP3 projections the OWP is projected to double during twenty-first century in a moderate CO₂ forcing scenario (IPCC A1B scenario). However, during the observational period the area of positive atmospheric heating (referred to as the dynamic warm pool, DWP), has remained constant. The threshold SST (T _H ), which demarks the region of net heating and cooling, has increased from 26.6°C in the 1950s to 27.1°C in the last decade and it is projected to increase to ~28.5°C by 2100. Based on climate model simulations, the area of the DWP is projected to remain constant during the twenty-first century. Analysis of the paleoclimate model intercomparison project (PMIP I and II) simulations for the Last Glacial maximum and the Mid-Holocene periods show a very similar behaviour, with a larger OWP in periods of elevated tropical SST, and an almost constant DWP associated with a varying T _H . The constancy of the DWP area, despite shifts in the background SST, is shown to be the result of a near exact matching between increases in the integrated convective heating within the DWP and the integrated radiative cooling outside the DWP as SST changes. Although the area of the DWP remains constant, the total tropical atmospheric heating is a strong function of the SST. For example the net heating has increased by about 10% from 1950 to 2000 and it is projected to increase by a further 20% by 2100. Such changes must be compensated by a more vigorous atmospheric circulation, with growth in convective heating within the warm pool, and an increase of subsiding air and stability outside the convective warm pool and an increase of vertical shear at the DWP boundaries. This finding is contrary to some conclusions from other studies but in accord with others. We discuss the similarities and differences at length.