Brilli M, Lerche L, Ciaranfi N, Turi B
Earth Science Department, University of Rome La Sapienza, Italy.
Appl Radiat Isot. 2000 Apr;52(4):957-64. doi: 10.1016/s0969-8043(99)00079-2.
Quantitative signal processing methods have been applied to a delta(18)O profile for a land-based stratigraphic section, extending from the upper part of lower Pleistocene to the lower part of middle Pleistocene. The section is well exposed with a continuous succession of muds and muddy silts, about 400 m thick, located in the southernmost part of Bradano Trough, near Montalbano Jonico in Basilicata (south Italy). The sampled part of the section is about 240 m thick, in which a foram benthic species (Cassidulina carinata) is continuously available for oxygen isotope ratio measurements. The aim of the data treatment is to discover how much of the Earth's orbital periodic movements, precession and obliquity, which represent the dominant periodicities in paleoclimatic variations from the base of the Pleistocene until 0.735 Myr BP, are responsible for the oscillations observed in the oxygen-18 record of the Montalbano Jonico section. A time framework of the section was constructed on the basis of calcareous nannofossil biostratigraphic analyses, preliminary magnetostratigraphic results and oxygen isotope correlation with the record from DSDP s607 (isotope data collected in the NOAA World Data Center). The resulting time-scale extends from 1.15 to 0.74 Myr. Power spectrum analysis was performed on the isotope data to illuminate the most important periodicity components of the Montalbano Jonico record. The periodic components of 41,000 and 21,000 yr are present in this record; the former associated with periodic changes in the tilt of Earth's axis and the latter with periodic changes with the precession of the equinoxes, as predicted by the astronomical theory of ice ages. They are, however, not the most important components of the power spectrum, in which a lower frequency component contains most of the variance. This low-frequency component is centered at a period around 208,000 yr. This periodicity seems not to be attributable to any known astronomical or paleoclimatic phenomenon. An attempt was made to verify if this periodicity was due to the composite effect of precession and obliquity signals together at different frequencies from their forcing frequencies. In order to investigate this effect, isotope data have been parameterized in terms of a sum of simple functions of precession and obliquity signals with unknown coefficients. The coefficients are estimated from the time series with the assumption that the best coefficients are those which minimize the 'noise' i.e. the difference between the data function and the precession and obliquity functions. Cross-spectra analyses were also performed on the data and the precession signal and on the data and the obliquity signal. The power spectrum of the residual 'noise' functions and the cross-spectra demonstrate that precession and obliquity signals are not in phase with the data at their forcing frequencies and so damp. The precession and obliquity signals were then shifted towards lower frequencies at equally spaced lags, the resultant 'noise' power spectra were plotted for every combination of lags of precession and obliquity. The results of this data processing demonstrate that it is possible to have a combination of precession and obliquity cyclicities that could be responsible for the signal with 208 kyr periodicity.
定量信号处理方法已应用于一个陆基地层剖面的δ¹⁸O曲线,该剖面从更新世早期上部延伸至更新世中期下部。该剖面出露良好,有连续的泥质和泥质粉砂层序,厚约400米,位于布拉达诺海槽最南端,靠近巴西利卡塔(意大利南部)的蒙塔尔巴诺乔尼科。该剖面的采样部分厚约240米,其中一种有孔底栖生物(卡氏卡西杜林娜)可连续用于氧同位素比值测量。数据处理的目的是确定地球轨道的周期性运动,即岁差和倾角,在更新世开始至距今0.735百万年前的古气候变化中代表主要周期,它们在蒙塔尔巴诺乔尼科剖面的氧-18记录中观测到的振荡中占多大比例。该剖面的时间框架是基于钙质超微化石生物地层分析、初步的磁性地层结果以及与深海钻探计划607站记录的氧同位素对比(在美国国家海洋和大气管理局世界数据中心收集的同位素数据)构建的。所得时间尺度从1.15百万年延伸至0.74百万年。对同位素数据进行了功率谱分析,以阐明蒙塔尔巴诺乔尼科记录中最重要的周期成分。该记录中存在41000年和21000年的周期成分;前者与地轴倾斜的周期性变化有关,后者与二分点岁差的周期性变化有关,正如冰期天文理论所预测的那样。然而,它们并不是功率谱中最重要的成分,其中低频成分包含了大部分方差。这个低频成分集中在约208000年的周期附近。这个周期性似乎不归因于任何已知的天文或古气候现象。有人试图验证这个周期性是否是由于岁差和倾角信号在与它们的强迫频率不同的频率下共同作用的复合效应。为了研究这种效应,同位素数据已根据岁差和倾角信号的简单函数之和进行参数化,系数未知。系数是从时间序列中估计出来的,假设最佳系数是那些使“噪声”最小化的系数,即数据函数与岁差和倾角函数之间的差异。还对数据与岁差信号以及数据与倾角信号进行了互谱分析。残余“噪声”函数的功率谱和互谱表明,岁差和倾角信号在其强迫频率下与数据不同步,因此会衰减。然后将岁差和倾角信号以等间距滞后向低频方向移动,针对岁差和倾角滞后的每种组合绘制所得的“噪声”功率谱。这种数据处理的结果表明,有可能存在岁差和倾角周期的组合,它们可能是208千年周期信号的成因。