The sea level conundrum: insights from paleo studies

Empirical Constraints on Future Sea Level Rise; Bern, Switzerland, 25–29 August 2008; Eustatic sea level (ESL) rise during the 21st century is perhaps the greatest threat from climate change, but its magnitude is contested. Geological records identify examples of nonlinear ice sheet response to climate forcing, suggesting a strategy for refining estimates of 21st-century sea level change. In August 2008, Past Global Changes (PAGES), International Marine Past Global Change Study (IMAGES), and the University of Bern cosponsored a workshop to address this possibility. The workshop highlighted several ways that paleoceanography studies can place limits on future sea level rise, and these are enlarged upon here.

Evidence for a widespread climatic anomaly at around 9.2 ka before present.

The 8.2 ka event was triggered by a meltwater pulse (MWP) into the North Atlantic and resultant reduction of the thermohaline circulation (THC). This event was preceded by a series of at least 14 MWPs; their impact on early Holocene climate has remained almost unknown. A set of high-quality paleoclimate records from across the Northern Hemisphere shows evidence for a widespread and significant climatic anomaly at ∼9.2 ka B.P. This event has climatic anomaly patterns very similar to the 8.2 ka B.P. event, cooling occurred at high latitudes and midlatitudes and drying took place in the northern tropics, and is concurrent with an MWP of considerable volume (∼8100 km3). As the 9.2 ka MWP occurs at a time of enhanced baseline freshwater flow into the North Atlantic, this MWP may have been, despite its relatively small volume, sufficient to weaken THC and to induce the observed climate anomaly pattern.

Impact of strong deep ocean stratification on the glacial carbon cycle

During the Last Glacial Maximum, the climate was substantially colder and the carbon cycle was clearly different from the late Holocene. According to proxy data deep oceanic δ13C was very low, and the atmospheric CO2 concentration also reduced. Several mechanisms have been proposed to explain these changes, but none can fully explain the data, especially the very low deep ocean δ13C values. Oceanic core data show that the deep ocean was very cold and salty, which would lead to enhanced deep ocean stratification. We show that such an enhanced stratification in the coupled climate model CLIMBER-2 helps get very low deep oceanic δ13C values. Indeed the simulated δ13C reaches values as low as −0.8‰ in line with proxy data evidences. Moreover it increases the oceanic carbon reservoir leading to a small, yet robust, atmospheric CO2 drop of approximately 10 ppm.

Sea-surface temperature records of Termination 1 in the Gulf of California: Challenges for seasonal and interannual analogues of tropical Pacific climate change

Centennial-scale records of sea-surface temperature and opal composition spanning the Last Glacial Maximum and Termination 1 (circa 25–6 ka) are presented here from Guaymas Basin in the Gulf of California. Through the application of two organic geochemistry proxies, the U37K′ index and the TEX86H index, we present evidence for rapid, stepped changes in temperatures during deglaciation. These occur in both temperature proxies at 13 ka (∼3°C increase in 270 years), 10.0 ka (∼2°C decrease over ∼250 years) and at 8.2 ka (3°C increase in <200 years). An additional rapid warming step is also observed in TEX86H at 11.5 ka. In comparing the two temperature proxies and opal content, we consider the potential for upwelling intensity to be recorded and link this millennial-scale variability to shifting Intertropical Convergence Zone position and variations in the strength of the Subtropical High. The onset of the deglacial warming from 17 to 18 ka is comparable to a “southern hemisphere” signal, although the opal record mimics the ice-rafting events of the north Atlantic (Heinrich events). Neither the modern seasonal cycle nor El Niño/Southern Oscillation patterns provide valid analogues for the trends we observe in comparison with other regional records. Fully coupled climate model simulations confirm this result, and in combination we question whether the seasonal or interannual climate variations of the modern climate are valid analogues for the glacial and deglacial tropical Pacific.

Last glacial vegetation of northern Eurasia

In order to investigate the potential role of vegetation changes in megafaunal extinctions during the later part of the last glacial stage and early Holocene (42–10 ka BP), the palaeovegetation of northern Eurasia and Alaska was simulated using the LPJ-GUESS dynamic vegetation model. Palaeoclimatic driving data were derived from simulations made for 22 time slices using the Hadley Centre Unified Model. Modelled annual net primary productivity (aNPP) of a series of plant functional types (PFTs) is mapped for selected time slices and summarised for major geographical regions for all time slices. Strong canonical correlations are demonstrated between model outputs and pollen data compiled for the same period and region. Simulated aNPP values, especially for tree PFTs and for a mesophilous herb PFT, provide evidence of the structure and productivity of last glacial vegetation. The mesophilous herb PFT aNPP is higher in many areas during the glacial than at present or during the early Holocene. Glacial stage vegetation, whilst open and largely treeless in much of Europe, thus had a higher capacity to support large vertebrate herbivore populations than did early Holocene vegetation. A marked and rapid decrease in aNPP of mesophilous herbs began shortly after the Last Glacial Maximum, especially in western Eurasia. This is likely implicated in extinction of several large herbivorous mammals during the latter part of the glacial stage and the transition to the Holocene.

Validating NCAR-CCSM last glacial maximum sea surface temperature in the tropical and South Atlantic with proxy-data

Last Glacial Maximum simulated sea surface temperature from the Paleo-Climate version of the National Center for Atmospheric Research Coupled Climate Model (NCAR-CCSM) are compared with available reconstructions and data-based products in the tropical and south Atlantic region. Model results are compared to data proxies based on the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface product (MARGO). Results show that the model sea surface temperature is not consistent with the proxy-data in all of the region of interest. Discrepancies are found in the eastern, equatorial and in the high-latitude South Atlantic. The model overestimates the cooling in the southern South Atlantic (near 50 degrees S) shown by the proxy-data. Near the equator, model and proxies are in better agreement. In the eastern part of the equatorial basin the model underestimates the cooling shown by all proxies. A northward shift in the position of the subtropical convergence zone in the simulation suggests a compression or/and an equatorward shift of the subtropical gyre at the surface, consistent with what is observed in the proxy reconstruction. (C) 2008 Elsevier B.V. All rights reserved

Changes in the intermediate water mass formation rates in the global ocean for the Last Glacial Maximum, mid-Holocene and pre-industrial climates

The paleoclimate version of the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR-CCSM3) is used to analyze changes in the water formation rates in the Atlantic, Pacific, and Indian Oceans for the Last Glacial Maximum (LGM), mid-Holocene (MH) and pre-industrial (PI) control climate. During the MH, CCSM3 exhibits a north-south asymmetric response of intermediate water subduction changes in the Atlantic Ocean, with a reduction of 2 Sv in the North Atlantic and an increase of 2 Sv in the South Atlantic relative to PI. During the LGM, there is increased formation of intermediate water and a more stagnant deep ocean in the North Pacific. The production of North Atlantic Deep Water (NADW) is significantly weakened. The NADW is replaced in large extent by enhanced Antarctic Intermediate Water (AAIW), Glacial North Atlantic Intermediate Water (GNAIW), and also by an intensified of Antarctic Bottom Water (AABW), with the latter being a response to the enhanced salinity and ice formation around Antarctica. Most of the LGM intermediate/mode water is formed at 27.4 < sigma(theta) < 29.0 kg/m(3), while for the MH and PI most of the subduction transport occurs at 26.5 < sigma(theta) < 27.4 kg/m(3). The simulated LGM Southern Hemisphere winds are more intense by 0.2-0.4 dyne/cm(2). Consequently, increased Ekman transport drives the production of intermediate water (low salinity) at a larger rate and at higher densities when compared to the other climatic periods.

Impact of suborbital climate changes in the North Atlantic on ice sheet dynamics at the Mid-Pleistocene Transition

[1] Early and Mid-Pleistocene climate, ocean hydrography and ice sheet dynamics have been reconstructed using a high-resolution data set (planktonic and benthicδ18O time series, faunal-based sea surface temperature (SST) reconstructions and ice-rafted debris (IRD)) record from a high-deposition-rate sedimentary succession recovered at the Gardar Drift formation in the subpolar North Atlantic (Integrated Ocean Drilling Program Leg 306, Site U1314). Our sedimentary record spans from late in Marine Isotope Stage (MIS) 31 to MIS 19 (1069–779 ka). Different trends of the benthic and planktonic oxygen isotopes, SST and IRD records before and after MIS 25 (∼940 ka) evidence the large increase in Northern Hemisphere ice-volume, linked to the cyclicity change from the 41-kyr to the 100-kyr that occurred during the Mid-Pleistocene Transition (MPT). Beside longer glacial-interglacial (G-IG) variability, millennial-scale fluctuations were a pervasive feature across our study. Negative excursions in the benthicδ18O time series observed at the times of IRD events may be related to glacio-eustatic changes due to ice sheets retreats and/or to changes in deep hydrography. Time series analysis on surface water proxies (IRD, SST and planktonicδ18O) of the interval between MIS 31 to MIS 26 shows that the timing of these millennial-scale climate changes are related to half-precessional (10 kyr) components of the insolation forcing, which are interpreted as cross-equatorial heat transport toward high latitudes during both equinox insolation maxima at the equator.

Toward explaining the Holocene carbon dioxide and carbon isotope records: Results from transient ocean carbon cycle-climate simulations

[1] The Bern3D model was applied to quantify the mechanisms of carbon cycle changes during the Holocene (last 11,000 years). We rely on scenarios from the literature to prescribe the evolution of shallow water carbonate deposition and of land carbon inventory changes over the glacial termination (18,000 to 11,000 years ago) and the Holocene and modify these scenarios within uncertainties. Model results are consistent with Holocene records of atmospheric CO2 and δ13C as well as the spatiotemporal evolution of δ13C and carbonate ion concentration in the deep sea. Deposition of shallow water carbonate, carbonate compensation of land uptake during the glacial termination, land carbon uptake and release during the Holocene, and the response of the ocean-sediment system to marine changes during the termination contribute roughly equally to the reconstructed late Holocene pCO2 rise of 20 ppmv. The 5 ppmv early Holocene pCO2 decrease reflects terrestrial uptake largely compensated by carbonate deposition and ocean sediment responses. Additional small contributions arise from Holocene changes in sea surface temperature, ocean circulation, and export productivity. The Holocene pCO2 variations result from the subtle balance of forcings and processes acting on different timescales and partly in opposite direction as well as from memory effects associated with changes occurring during the termination. Different interglacial periods with different forcing histories are thus expected to yield different pCO2 evolutions as documented by ice cores.

On the relationship between Nd isotopic composition and ocean overturning circulation in idealized freshwater discharge events

Using a cost-efficient climate model, the effect of changes in overturning circulation on neodymium isotopic composition,ϵNd, is systematically examined for the first time. Idealized sequences of abrupt climate changes are induced by the application of periodic freshwater fluxes to the North Atlantic (NA) and the Southern Ocean (SO), thus mainly affecting either the formation of North Atlantic Deep Water (NADW) or Antarctic Bottom Water (AABW). Variations in ϵNd reflect weakening and strengthening of the formation of NADW and AABW, changes in ϵNdof end-members are relatively small. Relationships betweenϵNd and the strength of NADW or AABW are more pronounced for AABW than for NADW. Atlantic patterns of variations in ϵNd systematically differ between NA and SO experiments. Additionally, the signature of changes in ϵNd in the Atlantic and the Pacific is alike in NA but opposite in SO experiments. Discrimination between NA and SO experiments is therefore possible based on the Atlantic pattern of variations in ϵNd and the contrariwise behavior of ϵNd in the Atlantic and the Pacific. In further experiments we examined the effect of variations in magnitudes of particle export fluxes. Within the examined range, and although settling particles represent the only sink of Nd, their effects on ϵNd are relatively small. Our results confirm the large potential of ϵNd as a paleocirculation tracer but also indicate its limitations of quantitative reconstructions of changes in the Atlantic Meridional Ocean Circulation.