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.

Mineralogy of ODP Hole 117-722B (Appendix)

Bulk mineralogy of the terrigenous fraction of 99 samples from ODP Site 722 on the Owen Ridge, western Arabian Sea, has been determined by x-ray diffraction, using an internal standard method. The sampling interval, approximately 4.3 k.y., provides a detailed mineralogic record for the past 500 k.y. Previous studies have identified important modern continental sediment sources and the mineral assemblages presently derived from each. These studies have also demonstrated that most of this material is supplied by southwest and northwest winds during the summer monsoon. A variety of marine and terrestrial records and general circulation model (GCM) simulations have indicated the importance of monsoonal circulation during the Pleistocene and Holocene and have demonstrated increased aridity during glacial times and increased humidity during inter glacials. The mineralogic data generated here were used to investigate variations in source area weathering conditions during these environmental changes. Terrigenous minerals present include smectite, illite, palygorskite, kaolinite, chlorite, quartz, plagioclase feldspar, and dolomite. This mineralogy is consistent with the compositions of source areas presently supplying sediment to the Arabian Sea. An R-mode factor analysis has identified four mineral assemblages present throughout the past 500 k.y.: quartz/chlorite/dolomite (Factor 1), kaolinite/plagioclase/illite (Factor 2), smectite (Factor 3), and palygorskite/dolomite (Factor 4). Chlorite, illite, and palygorskite are extremely susceptible to chemical weathering, and a spectral comparison of these factors with the eolian mass accumulation rate (MAR) record from Hole 722B (an index of dust source area aridity) indicates that Factors 1, 2, and 4 are directly related to changes in aridity. Because of these characteristics, Factors 1,2, and 4 are interpreted to originate from arid source regions. Factor 3 is interpreted to record more humid source conditions. Time-series of scores for the four factors are dominated by short-term (10-100 k.y.) variability, and do not correlate well to glacial/interglacial fluctuations in the time domain. These characteristics suggest that local climatic shifts were complex, and that equilibrium weathering assemblages did not develop immediately after climatic change. Spectral analysis of factor scores identifies peaks at or near the primary Milankovitch frequencies for all factors. Factor 1 (quartz/chlorite/dolomite), Factor 2 (kaolinite/plagioclase/illite), and Factor 4 (illite/palygorskite) are coherent and in phase with the MAR record over the 23, 41, and 100 k.y. bands, respectively. The reasons for coherency at single Milankovitch frequencies are not known, but may include differences in the susceptibilities of minerals to varying time scales of weathering and/or preferential development of suitable continental source environments by climatic changes at the various Milankovitch frequencies.

Geochemistry, physical properties, diatom abundance and datums of ODP Site 186-1150 and Hole 186-1151A

Detection of climate response to orbital forcing during Cenozoic long-term global cooling is a key to understanding the behavior of Earth's icehouse climate. Sedimentary rhythm, which is a rhythmic or cyclic variation in the sequence of sediments and sedimentary rocks, is useful for quantitative reconstruction of Earth's evolution during geological time. In this study, we attempt to (1) identify sources of natural gamma ray (NGR) emissions of core recovered during Ocean Drilling Program (ODP) Leg 186 by analyses of physical properties, major element concentrations, diatom abundances, and total organic carbon contents, (2) integrate whole-core NGR intensity of recovered core with wireline logging NGR measurements in order to construct a continuous sedimentary sequence, and (3) discuss changes in the NGR signal in the time domain. This attempt gives us preliminary information to discuss climate stability in relation to orbital forcing thorough geologic time. NGR values are obtained mainly by indirectly measuring the amount of terrigenous minerals including potassium and related elements in the sediments. NGR intensity is also affected by high porosity, which in these sediments was related to the amount of diatom valves. NGR signals might be a proxy of the intensity of the East Asian monsoon off Sanriku. A continuous sedimentary record was constructed by integration of the whole-core NGR intensity measured in sediments obtained from the drilled holes with that measured directly in the borehole by wireline logging, then using a stratigraphic age model to convert to a time series covering 1.3-9.7 Ma with a short break at ~5 Ma. High sedimentation rate (H) stages were identified in the sequence, related to intervals of low-amplitude precession and eccentricity variations. The transition of the dominant periodicities through the four H stages may correlate to major shifts in the climate system, including the onset of major Northern Hemisphere glaciation, the initial stage of the East Asian monsoon intensification, and the onset of the East Asian monsoon with uplift of the Himalayas and the Tibetan Plateau.

Modelled global ocean temperature from 6 ka to present

We compare the ocean temperature evolution of the Holocene as simulated by climate models and reconstructed from marine temperature proxies. This site provides informations about the Holocene temperature trends as simulated by the models. We use transient simulations from a coupled atmosphere-ocean general circulation model, as well as an ensemble of time slice simulations from the Paleoclimate Modelling Intercomparison Project. The general pattern of sea surface temperature (SST) in the models shows a high latitude cooling and a low latitude warming. The proxy dataset comprises a global compilation of marine alkenone- and Mg/Ca-derived SST estimates. Independently of the choice of the climate model, we observe significant mismatches between modelled and estimated SST amplitudes in the trends for the last 6000 years. Alkenone-based SST records show a similar pattern as the simulated annual mean SSTs, but the simulated SST trends underestimate the alkenone-based SST trends by a factor of two to five. For Mg/Ca, no significant relationship between model simulations and proxy reconstructions can be detected. We tested if such discrepancies can be caused by too simplistic interpretations of the proxy data. We tested different seasons and depths in the model to compare the proxy data trends, and can reconcile only part of the mismatches on a regional scale. We therefore considered the additional environmental factor changes in the planktonic organisms' habitat depth and a time-shift in the recording season to diagnose whether invoking those environmental factors can help reconciling the proxy records and the model simulations. We find that invoking shifts in the living season and habitat depth can remove some of the model-data discrepancies in SST trends. Regardless whether such adjustments in the environmental parameters during the Holocene are realistic, they indicate that when modeled temperature trends are set up to allow drastic shifts in the ecological behavior of planktonic organisms, they do not capture the full range of reconstructed SST trends. Our findings indicate that climate model and reconstructed temperature trends are to a large degree only qualitatively comparable, thus providing a challenge for the interpretation of proxy data as well as the models' sensitivity to orbital forcing.