Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Results of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958-2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Stable carbon and oxygen isotope records of Cibicidoides wuellerstorfi in DSDP and ODP sites in the Atlantic and Pacific Ocean

Carbon isotopic records from benthic foraminifera are used to map patterns of deep ocean circulation between 3 and 2 million years ago, the interval when significant northern hemisphere glaciation began. The delta18O and delta13C data from four Atlantic sites (552, 607, 610, and 704) and one Pacific site (677) show that global cooling over this interval was associated with increased suppression of North Atlantic Deep Water (NADW) formation. However, the relative strength of NADW production was always greater than is observed during late Pleistocene glaciations when extreme decreases in NADW are observed in the deep North Atlantic. Our data indicate that an increase in the equator-to-pole temperature gradient associated with the onset of northern hemisphere glaciation did not intensify deepwater production in the North Atlantic but rather the opposite occurred. This is not unexpected as it is the "warm high-salinity" characteristic, rather than the "low temperature", of thermocline waters that is critical to the deepwater formation process in this region today.

Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

Geochemistry of planktonic foraminifera from DSDP sites in the Atlantic and Pacific Ocean

The lithium content of planktonic foraminiferal calcite has been determined to evaluate temporal variability of seawater Li concentrations over the past 116 m.y. Mean foraminiferal calcite lithium/calcium in each time interval is no more than 16% greater nor 25% less than the mean Li/Ca of all samples. Li/Ca minima are observed for samples from 50-60 m.y. and 80-90 m.y., with Li/Ca about 25% lower than in adjacent time intervals. At no time during the past 40 m.y does mean Li/Ca appear to be higher than that at present. Subject to the limitations imposed by sample coverage and diagenesis, a similar conclusion holds for the past 116 m.y. Coupled with an oceanic mass balance model for Li, these data suggest that: (1) oceanic Li concentrations and, therefore, high-temperature hydrothermal circulation fluxes during the past 40 m.y. (and perhaps the past 100 m.y.) have not been more than perhaps 30-40% greater than at present for intervals any longer than a million years at most, and (2) these fluxes were not a factor of two higher 100 m.y. ago. By inference, variations in oceanic crustal generation rates over these time periods are similarly limited. Decreases in hydrothermal circulation fluxes and crustal generation rates or fluctuations up to 20% in these rates of a few million years duration are not necessarily ruled out by the Li/Ca data. The lack of variability in Li/Ca over time is not unequivocal evidence that hydrothermal fluxes have not varied because the rates of removal processes may be linked to changes in input fluxes.

Clay mineralogical assemblages, Sm-Nd concentrations, Sm and Nd isotopic ratios and oxygen isotope ratios, OIS number and the age model of ODP Site 105-646

We present 40 Sm-Nd isotope measurements of the clay-size (<2 µm) fractions of sediments from the Southern Greenland rise (ODP-646) that span the last 365 kyr. These data track changes in the relative supply of fine particles carried into the deep Labrador Sea by the Western Boundary Under Current (WBUC) back to the fourth glacial-interglacial cycles. Earlier studies revealed three general sources of particles to the core site: (i) Precambrian crustal material from Canada, Greenland, and/or Scandinavia (North American Shield - NAS), (ii) Palaeozoic or younger crustal material from East Greenland, NW Europe, and/or western Scandinavia (Young Crust - YC) and (iii) volcanic material from Iceland and the Mid-Atlantic Ridge (MAR). Clay-size fractions from glacial sediments have the lowest Nd isotopic ratios. Supplies of young crustal particles were similar during glacial oxygen isotope stages (OIS) 2, 6, and 10. In contrast the mean volcanic contributions decreased relative to old craton material from OIS 10 to OIS 6 and then from OIS 6 to OIS 2. The glacial OIS 8 interval displays a mean Sm/Nd ratio similar to those of interglacials OIS 1, 5, and 9. Compared with other interglacials, OIS 7 was marked by a higher YC contribution but a similar ~30% MAR supply. The overall NAS contribution dropped by a factor of 2 during each glacial/interglacial transition, with the MAR contribution broadly replacing it during interglacials. To decipher between higher supplies and/or dilution, particle fluxes from each end member were estimated. Glacial NAS fluxes were systematically higher than interglacial fluxes. During the time interval examined, fine particle supplies to the Labrador Sea were strongly controlled by proximal ice-margin erosion and thus echoed the glacial stage intensity. In contrast, the WBUC-carried MAR supplies from the eastern basins did not change significantly throughout the last 365 kyr, except for a marked increase in surface-sediments that suggests unique modern conditions. Distal WBUC-controlled inputs from the Northern and NE North Atlantic seem to have been less variable than proximal supplies linked with glacial erosion rate.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Late glacial and deglacial age model and biomarker analyses of sediment core MSM05/5_712-2 from the western continental slope of Svalbard

The transition from last glacial to deglacial and subsequently to modern interglacial climate conditions was accompanied by abrupt shifts in the palaeoceanographic setting in the subpolar North Atlantic. Knowledge about the role that sea ice coverage played during these rapid climate reversals is limited since most marine sediment cores from the higher latitudes provide only a coarse temporal resolution and often poorly preserved microfossils. Here we present a highly resolved reconstruction of sea ice conditions that characterised the eastern Fram Strait - a key area for water mass exchange between the Arctic Ocean and the North Atlantic - for the past 30 ka BP. This reconstruction is based on the distribution of the sea ice biomarker IP25 and phytoplankton derived biomarkers in a sediment core from the continental slope of western Svalbard. During the late glacial (30 ka to 19 ka BP), recurrent advances and retreats of sea ice characterised the study area and point to a hitherto less considered oceanic (and/or atmospheric) variability. A long-lasting perennial sea ice coverage in eastern Fram Strait persisted only at the very end of the Last Glacial Maximum (i.e. from 19.2 to 17.6 ka BP) and was abruptly reduced at the onset of Heinrich Event 1 - coincident with or possibly even inducing the collapse of the Atlantic Meridional Overturning Circulation (AMOC). Further maximum sea ice conditions prevailed during the Younger Dryas cooling event and support the assumption of an AMOC reduction due to increased formation and export of Arctic sea ice through Fram Strait. A significant retreat of sea ice and sea surface warming are observed for the Early Holocene.