Spectral analysis of seabed elevation from the Danish Wadden Sea

Bedforms both reflect and influence shallow water hydrodynamics and sediment dynamics. A correct characterization of their spatial distribution and dimensions is required for the understanding, assessment and prediction of numerous coastal processes. A method to parameterize geometrical characteristics using two-dimensional (2D) spectral analysis is presented and tested on seabed elevation data from the Knudedyb tidal inlet in the Danish Wadden Sea, where large compound bedforms are found. The bathymetric data were divided into 20x20 m areas on which a 2D spectral analysis was applied. The most energetic peak of the 2D spectrum was found and its energy, frequency and direction were calculated. A power-law was fitted to the average of slices taken through the 2D spectrum; its slope and y-intercept were calculated. Using these results the test area was morphologically classified into 4 distinct morphological regions. The most energetic peak and the slope and intercept of the power-law showed high values above the crest of the primary bedforms and scour holes, low values in areas without bedforms, and intermediate values in areas with secondary bedforms. The secondary bedform dimensions and orientations were calculated. An area of 700x700 m was used to determine the characteristics of the primary bedforms. However, they were less distinctively characterized compared to the secondary bedforms due to relatively large variations in their orientations and wavelengths. The method is thus appropriate for morphological classification of the seabed and for bedform characterization, being most efficient in areas characterized by bedforms with regular dimensions and directions.

Short-term spectral analysis of ODP Hole 114-704B

Short-term spectral analysis was carried out on geochemical logging data from ODP Site 704. The FFT was used to compute the amplitude spectra of short-term overlapping segments to produce depth-period-amplitude spectrograms of the logging data. The spectrograms provided a means of evaluating the significance of the observed periodic components. The periodic components that were consistently present and prominent across a given record interval were considered to be significant. Changes in the spectrogram characteristics seem to reflect changes in either lithology, sedimentation rates, or hiatuses and may therefore provide useful information to aid in stratigraphic and paleoenvironmental studies. The dominant periodicity during the late Pleistocene and Brunhes Chron (0.97 to 0.47 Ma) was determined to be > 100,000 yr whereas the upper Matuyama Chron was dominated by the 41,000-yr periodicity. These periodicities suggest that the sedimentation patterns within the upper Matuyama Chron (0.98-1.78 Ma) were influenced by the Milankovitch obliquity cycle and those within the latest Matuyama-Brunhes Chron (<0.98 Ma) by the eccentricity cycle. The Brunhes/Matuyama boundary therefore represents a major discontinuity. Periodicities observed within the lower Matuyama and the upper Gauss Chron did not correlate with any of the periodicities within the Milankovitch frequency bands.

Grain size analysis on CRP sediment cores from the Ross Sea, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-l.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

Stable isotope analysis of benthic foraminifera and sedimentology on ODP Site 162-984

Sediments recovered at ODP Site 984 on the Reykjanes Ridge provided multicentennial-scale records (SST, planktic and benthic delta18O, IRD and magnetic susceptibility) of Late Pliocene climate change over the onset of Northern Hemisphere glaciation (NHG), 2.95-2.82 Ma. Short-term climate variations prior and after the onset of continent-wide glaciation were compared to test the hypothesis, whether Dansgaard-Oeschger (DO) cycles may have been triggered by continental ice breakouts. By means of spectral analyses for two selected interglacial stages prior to and after NHG (G15 and G1), we found that climate variability resembled that of the Holocene and the mid-Pliocene warm period. In contrast, DO-like periodicities near 1470, 2900, and 4400 yr indeed only occurred in glacial stages after the onset of NHG (G14, G6, and 104), but hardly in stage G20 prior to the onset. These results suggest a causal link between DO cycles and the Late Pliocene onset of major NHG and ice breakouts in the North Atlantic. This data set provides all primary data and spectral analysis related to this scientific work.

Dating analysis from DSDP Hole 90-594

Date-32 is a fast and easily used computer program developed to date Quaternary deep-sea cores by associating variations in the earth's orbit with recurring oscillations in core properties, such as carbonate content or isotope composition. Starting with known top and bottom dates, distortions in the periodicities of the core properties due to varying sedimentation rates are realigned by fast Fourier analysis so as to maximise the spectral energy density at the orbital frequencies. This allows age interpolation to all parts of the core to an accuracy of 10 kyrs, or about 1.5% of the record duration for a typical Brunhes sequence. The influence of astronomical forcing is examined and the method is applied to provide preliminary dates in a high-resolution Brunhes record from DSDP Site 594 off southeastern New Zealand.

Microfossil cyclicity in Late Albian sediments of Kirchrode I borehole

We studied the biological response to orbital forcing in marine Upper Albian sediments recovered from the 245 m-long Kirchrode I borehole in the Lower Saxony basin in northwestern Germany. Results from quantitative analysis of planktonic and benthic foraminifera, of calcareous nannofossils, and radiolaria were used for this study. Spectral analysis in the depth domain indicates for the high sedimentation rate part of the Upper Albian dominant periods with wavelengths of 10±13 m, 5±6 m, and 2±3 m, which we interpret to represent the biological response to orbital forcing in the Milankovitch frequency bands eccentricity, obliquity, and precession, respectively. In addition, a low amplitude 40±50 m cycle was found, which would represent the long-term eccentricity variation of roughly 400 ka. Microfossil cyclicity does not change significantly within the whole core indicating sedimentation rates of 11±12 cm/ka on an average, with variations between 3.5 and 13 ka. Microfossils show greater variability in their abundance changes than the physical and chemical parameters and also greater power in the higher-frequency bands (obliquity and precession). While most of the planktonic foraminifer species studied are dominated by variations in the obliquity, most benthic foraminifer species show an additional strong influence of precession. These differences in the cyclicity of the abundance changes are interpreted as reflecting a stronger influence of low latitude water in the deep waters of the Late Albian Lower Saxony basin than in the shallow waters. This basin was part of a wide, 'Boreal' epicontinental sea, which was connected to the Tethys to the south via the Polish strait and via the Paris basin, and which was connected with the North Atlantic and Arctic to the north. In analogy to results from analysis of data from the Late Neogene, strong effects of precession interpreted as being more characteristic for changes/influences triggered in the low latitudes and those of obliquity to be more characteristic for influences from the high latitudes. The presence of a relatively strong eccentricity cycle, not only in the compound parameters, but also in the abundance changes of single species during the Late Albian means that there must have been a non-linear response to orbital forcing and internal feedbacks.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Silt grain-size analysis of ODP holes from the western equatorial Atlantic

Grain-size records of the terrigenous and calcareous silt fraction, preservation of planktic foraminifera, and benthic foraminiferal stable-isotope data (delta13C, delta18O values of C. wuellerstorfi) at ODP Site 927 on the Ceará Rise (5°27.7'N, 44°28.8'W), are used to reconstruct variations in the history of bottom current strength, ventilation, and carbonate corrosiveness of deep waters during the time interval from 0.8 to 0.3 Ma. Glacial periods are characterized by generally smaller mean sizes of the terrigenous sortable silt fraction (mean(SS)), lower delta13C values, and poorer preservation of planktic foraminifera compared to interglacials. This indicates lower bottom current speeds, larger nutrient contents and more corrosive deep water. By contrast, larger mean(SS) sizes, higher delta13C values, and well preserved planktic foraminifera indicate strong circulation and a well ventilated deep-water mass during interglacials. The observed changes are most likely related to the weakening and strengthening of circulation of Lower North Atlantic Deep Water (LNADW). Cross-spectral analysis between the mean(SS) and benthic delta18O records reveals that minima in mean(SS) occur about 7.6 k.y. after the maximum in ice volume. This indicates a considerable lag time between ice-shield induced changes in LNADW production and subsequent changes in the velocity of LNADW flow in the western equatorial Atlantic. Striking changes in bottom current speed occur regularly during glacial to interglacial transitions. Extremely fine mean(SS) minima point to an almost complete shutdown of bottom current vigor in response to a cessation of LNADW production caused by an enhanced melt water release during the initial phases of deglaciation. However, each of the fine minima extremes is followed by a rapid shift to very high mean(SS) values that indicate strong bottom currents, and hence, vigorous LNADW flow during the early interglacials. After the onset of glacial Stage 12, generally poorer carbonate preservation and higher variability is registered. This coincides with a global decrease in carbonate preservation during the mid-Brunhes (mid-Brunhes dissolution event). Detailed grain-size analysis of the calcareous fine fraction (<63 µm) revealed a considerable reduction of particles in the fraction from 7 to 63 µm during periods of enhanced dissolution. This indicates a preferential dissolution of larger planktic foraminiferal fragments which leads to an enrichment of coccoliths in the calcareous fine fraction.

Upper Cretaceous sediments of ODP Hole 122-762C

Well-developed Campanian to Maestrichtian pelagic cyclic sediments were recovered from Hole 762C on the Exmouth Plateau, off northwest Australia, during Ocean Drilling Program Leg 122. The cycles consist of nannofossil chalk (light beds) and clayey nannofossil chalk (dark beds). Both light and dark beds are strongly to moderately bioturbated, alternate on a decimeter scale, and exhibit gradual boundaries. Bioturbation introduces materials from a bed of one color into an underlying bed of another color, indicating that diagenesis is not responsible for the cyclicity. Differences in composition between the light and dark beds, revealed by calcium carbonate measurement and X-ray diffraction analysis, together with trace fossil evidence, indicate that the cycles in the sediments are a depositional feature. Diagenetic processes may have intensified the appearance of the cycles. Spectral analysis was applied to the upper Campanian to lower Maestrichtian cyclic sediments to examine the regularity of the cycles. Power spectra were calculated from time series using Walsh spectral analysis. The most predominant wavelengths of the color cycles are 34-41 cm and 71-84 cm. With an average sedimentation rate of 1.82 cm/k.y. in this interval, we found the time durations of the cycles to be around 41 k.y. and 21 k.y., respectively, comparable to the obliquity and precession periods of the Earth's rotation, which strongly suggests an orbital origin for the cycles. On the basis of sedimentological evidence and plate tectonic reconstruction, we propose the following mechanism for the formation of the cyclic sediments from Hole 762C. During the Late Cretaceous, when there was no large-scale continental glaciation, the cyclic variations in insolation, in response to cyclic orbital changes, controlled the alternation of two prevailing climates in the area. During the wetter, equable, and warmer climatic phases under high insolation, more clay minerals and other terrestrial materials were produced on land and supplied by higher runoff to a low bioproductivity ocean, and the dark clayey beds were deposited. During the drier and colder climatic phases under low insolation, fewer clay minerals were produced and put into the ocean, where bioproductivity was increased and the light beds were deposited.

Stable isotopes of benthic foraminifera from Site ODP 184-1146

We generated a high resolution (~8 ky) benthic record from a West Pacific marginal basin to investigate the detailed structure and spectral characteristics of deep water isotope fluctuations during the middle Miocene. The benthic record from ODP Site 1146 allows unprecedented resolution of the structure of the middle Miocene delta13C excursion, as well as tighter control on the chronology of climatic events. Spectral analysis of the variance in the delta18O and delta13C records from ODP Site 1146 reveals spectral power concentrated in the eccentricity band (400-, ~100-ky) over the time interval between 13 and 17 Ma. The amplitude evolution in the 400-ky band is strikingly similar to that of the long eccentricity in Laskar's solution. There is an abrupt switch to the obliquity band in the delta18O record at -14.9 Ma, suggesting a shift in the ocean/climate response to orbital forcing (from low latitude eccentricity to high latitude obliquity forcing). The obliquity signal is pervasive in the delta18O record until -13.9 Ma, when a sharp increase in delta18O values indicates a major climatic transition. Comparison of delta18O and delta13C profiles from DSDP Site 588 (SW Pacific Ocean), ODP Site 761 (E Indian Ocean) and ODP Site 1146 (South China Sea) reveals significantly cooler deep water in the NE Indian Ocean throughout the middle Miocene and a restricted deep water exchange between the Pacific Ocean and Indian Ocean.