Integrated Data Sets of the EU FP5 Research Project ORFOIS: Origin and fate of biogenic particle fluxes in the ocean and their interactions with atmospheric CO2 concentrations as well as the marine sediment

For a reliable simulation of the time and space dependent CO2 redistribution between ocean and atmosphere an appropriate time dependent simulation of particle dynamics processes is essential but has not been carried out so far. The major difficulties were the lack of suitable modules for particle dynamics and early diagenesis (in order to close the carbon and nutrient budget) in ocean general circulation models, and the lack of an understanding of biogeochemical processes, such as the partial dissolution of calcareous particles in oversaturated water. The main target of ORFOIS was to fill in this gap in our knowledge and prediction capability infrastructure. This goal has been achieved step by step. At first comprehensive data bases (already existing data) of observations of relevance for the three major types of biogenic particles, organic carbon (POC), calcium carbonate (CaCO3), and biogenic silica (BSi or opal), as well as for refractory particles of terrestrial origin were collated and made publicly available.

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.

Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

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.

(Table 1) Diatom abundance and preservation in ODP Hole 130-806B samples

We report on diatom abundance and preservation stratigraphy in the uppermost four cores of Ocean Drilling Program Hole 806B, which span the Quaternary period. Changes in diatom abundance and preservation show a rather complicated pattern, with much noise at high frequencies. However, in the cycles corresponding to eccentricity and obliquity variations, the picture is quite clear. Abundance and preservation follow glacial-interglacial cycles, with lowest abundances and poorest preservation observed in sediments that correspond to glacial stages. Seventy taxa compose the diatom assemblage of Hole 806B (from Samples 130-806B-1H-1, 8-9 cm, to -4H-7, 73-74 cm) with Azpeitia nodulifera as the dominant member. This species exhibits significant size variations related to glacial and interglacial stages during the Pleistocene. The distribution of power in the Fourier spectrum of the diatom signal (in the time domain) displays the expected Milankovitch frequencies (at 100,41, and 24-18 k.y.). It also shows concentration at various "odd" frequencies, especially at 62 k.y., suggesting a complicated response of productivity (and silicate chemistry) to climatic forcing.

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.