高布錫

(中國(guó)科學(xué)院精密測(cè)量科學(xué)與技術(shù)創(chuàng)新研究院武漢430077)

1 緒論

地月之間距離的長(zhǎng)期變化, 主要是由于月球引起的潮汐耗散引起的. 月球?qū)Ｑ笊铣毕∑鸬奈⑿∫κ窃嫉囊绷? 它正比于潮汐的高度, 而海潮高度反比于地月之間距離的立方[1]. 因此利用開(kāi)普勒第3定律, 可以導(dǎo)出月球的引潮力位正比于月球平均運(yùn)動(dòng)的平方. 根據(jù)月球的平均運(yùn)動(dòng)以及它的加速度, 月潮的能量耗散可以被計(jì)算出來(lái). 如果采用現(xiàn)代的月球平均運(yùn)動(dòng)加速度值?25.90′′·cy?2, 地月距離的演化就能被反推出來(lái)[2].但是計(jì)算結(jié)果顯示, 在28億年前地月距離迅速減小, 甚至趨向于零. 某些學(xué)者推算地月系的演化過(guò)程也得到了同樣的結(jié)果[3–7]. 這個(gè)事件被稱(chēng)為格斯科恩事件, 顯然這是不可能的. 因此采用現(xiàn)在的月球平均運(yùn)動(dòng)加速度值來(lái)反推地月系統(tǒng)的演化, 顯然是不合適的.現(xiàn)在月球平均運(yùn)動(dòng)加速度值為?25.90′′·cy?2, 它大于用古生物鐘給出的結(jié)果. 如果采用月球平均運(yùn)動(dòng)加速度值為?20.72′′·cy?2, 它等于現(xiàn)在月球平均運(yùn)動(dòng)加速度值的0.8倍, 則推算結(jié)果表明, 在35億年前地月距離是0.18百萬(wàn)公里. 因此直到35億年前格斯科恩事件都可以避免.

地球自轉(zhuǎn)速率變化的長(zhǎng)期項(xiàng)主要受3個(gè)因素影響. 最主要的貢獻(xiàn)是月球潮, 第2個(gè)因素是太陽(yáng)潮. 由于地球位相的變化, 月球潮和太陽(yáng)潮總是互相加強(qiáng)或者抵消, Gao[8–9]曾估計(jì)太陽(yáng)潮貢獻(xiàn)的能量耗散僅僅是月球潮的1/6. 第3個(gè)因素是地球極角動(dòng)量矩的變化.太陽(yáng)潮和地球極角動(dòng)量矩變化對(duì)月球自轉(zhuǎn)長(zhǎng)期變化的貢獻(xiàn)可以認(rèn)為是固定的. 因此可以推導(dǎo)出過(guò)去35億年以來(lái)的回歸年日數(shù)變化. 例如在33億年前, 日長(zhǎng)約為9.6 h, 一個(gè)回歸年中有913.2 d.

在大碰撞之后[10], 吸積過(guò)程曾使得月球達(dá)到融化的溫度[11]. 但是根據(jù)熱模型顯示,月殼在50–100百萬(wàn)年間迅速冷卻[12]. 隨后從40到38億年前發(fā)生了后期大轟炸事件[13].月球自轉(zhuǎn)形成與軌道同步應(yīng)該是在后期大轟炸事件之前. 由于早期月球離地球很近以及地球的強(qiáng)大引力場(chǎng), 絕大多數(shù)隕石落在了月球的背面, 它們形成了厚厚的覆蓋層, 并且導(dǎo)致了月球的形狀中心對(duì)月球質(zhì)心有1117 m的偏離(根據(jù)嫦娥月球模型CE-1). 雖然現(xiàn)在的月球沒(méi)有磁場(chǎng), 但是月球巖石顯示在39–36億年前曾有剩余磁場(chǎng)[14]. 所以可以估計(jì)除月核外月球的凝固時(shí)間大約在35億年前.

假定月球在自轉(zhuǎn)離心力和地球的引潮力作用下形成三軸橢球體, 月球表面形成等位面. 再假定月球冷卻后一直保持它的化石形狀至今, 那么根據(jù)月球的三軸橢球體形狀可以算出在35億年前, 地月距離應(yīng)該等于1.7455×108m, 月球的自轉(zhuǎn)周期等于3.652 d[15–16]. 根據(jù)開(kāi)普勒第3定律, 恒星月的長(zhǎng)度應(yīng)該等于8.34 d. 但是后期大轟炸可能部分地破壞了月球的平衡面形狀, 因此上述關(guān)于地月距離以及月球自轉(zhuǎn)周期的數(shù)值可能并不精確.

2 地月距離的演化

月球引起的潮汐摩擦耗散使得地月距離逐漸變大, 并使得地球自轉(zhuǎn)逐漸變慢. 海洋上的潮汐隆起總是略微超前于月球軌道, 并且產(chǎn)生一個(gè)小的引力拖住月球, 這個(gè)力等于總的潮汐摩擦力. 潮汐與地球表面的摩擦力使得地球自轉(zhuǎn)變慢. 月球?qū)Τ毕∑鸬囊φ扔诔毕∑鸬馁|(zhì)量, 而潮汐隆起的質(zhì)量正比于潮汐的高度. 所以可以簡(jiǎn)單地假定潮汐摩擦力正比于潮汐的高度(見(jiàn)圖1). 潮汐的高度反比于r3, 在這r是地月之間的距離[1]. 應(yīng)用開(kāi)普勒第3定律n2r3=f(M+m), 可以導(dǎo)出月潮引起的潮汐摩擦耗散正比于月球的平均運(yùn)動(dòng)n的平方, 其中,f是萬(wàn)有引力常數(shù),M、m分別是地球質(zhì)量. 因此t為時(shí)間. 在這C和ω分別是地球的極慣量矩和自轉(zhuǎn)角速度. 根據(jù)地月系統(tǒng)的角動(dòng)量守恒定律, 可以給出[8–9]：

i是月球軌道的傾角,K是常數(shù). 不考慮太陽(yáng)潮的影響, 假設(shè)地球的極慣量矩C是常數(shù), 利用(1)式, 月球平均運(yùn)動(dòng)n的演化可以用(2)式計(jì)算：

A是一個(gè)常數(shù), 它可以從求得.n0是現(xiàn)在的月球平均運(yùn)動(dòng),0是現(xiàn)在月球平均運(yùn)動(dòng)的長(zhǎng)期加速度. 將(2)式積分, 可以得到：

其中,t0為起始時(shí)間.

圖1 月潮引起的能量耗散正比于海洋潮汐的高度Fig.1 Energy dissipated by lunar tide is proportional to oceanic tidal height

為了計(jì)算現(xiàn)在的月球平均運(yùn)動(dòng)長(zhǎng)期加速度值0, 1995年Yoder[2]根據(jù)多種不同測(cè)量方法得到了數(shù)值相互符合的結(jié)果. 不妨取

由此可算出現(xiàn)在地月平均距離每年增加3.83 cm. 現(xiàn)在月球n0= 13.17636?·d?1, 根據(jù)n0和0的數(shù)值, 利用(3)式可以得到計(jì)算n的公式：

利用開(kāi)普勒第3定律, 可用下式計(jì)算地月距離,

現(xiàn)在的月球平均運(yùn)動(dòng)加速度絕對(duì)值是比較大的, 它比自?shī)W陶紀(jì)(the Ordovician age)以來(lái)用古生物鐘獲得的平均數(shù)值大1.4倍[17]. 1982年Webb[7]采用半球海洋模型發(fā)現(xiàn),假如地球自轉(zhuǎn)速率比現(xiàn)在的數(shù)值增大1倍, 那么在海洋中的潮汐能量耗散將減少約1/3.大陸漂移大約發(fā)生在兩億年前, 它改變了海洋的形狀, 分裂的大陸可能更強(qiáng)烈地阻止潮汐的流動(dòng), 增加了潮汐能量耗散. 因此現(xiàn)在的潮汐能量耗散可能是有史以來(lái)最大的. 但兩億年相對(duì)于35億年的地球演化歷史是個(gè)很短的時(shí)間, 因此在以下的計(jì)算中將忽略由于大陸漂移造成的影響.

利用(6)式并采用絕對(duì)值較小的月球加速度值?20.72′′·cy?2, 就可以算出35億年以來(lái)地月距離的演化. 在圖2中用實(shí)線(xiàn)表示了這個(gè)變化.

圖2 過(guò)去35億年來(lái)地月距離的變化Fig.2 The change of distance between Earth-Moon during past 3.5 billion years

3 地球自轉(zhuǎn)速率的演化

有3種因素影響地球自轉(zhuǎn)速率的長(zhǎng)期變化, 它們是月球和太陽(yáng)引起的潮汐耗散以及地球的極慣量矩的變化. 設(shè)是地球自轉(zhuǎn)角速度的變化率, 有

其中,mt為月潮引起的地球自轉(zhuǎn)速率變化,st為太陽(yáng)潮引起的地球自轉(zhuǎn)速率變化,nt為非潮汐因素引起的地球自轉(zhuǎn)速率變化.

潮汐能量耗散主要是由于月潮引起的. 暫不考慮日潮的影響以及地球極慣量矩的變化, 利用(1)式, 可給出月潮引起的地球的自轉(zhuǎn)速率變化公式,

將(8)式積分, 就可以導(dǎo)出計(jì)算ωmt的公式,

利用(5)、(9)和(10)式, 可以求得上溯到35億年前ω的數(shù)值. 考慮到日月潮汐摩擦以及地球慣量矩的變化, 可以算出, 現(xiàn)在=?5.405×10?22s?2, 相對(duì)于日長(zhǎng)的長(zhǎng)期變化為2.01 ms/cy.

用N表示回歸年日數(shù), 可以用下式求得[8–9],

在圖3中描述了35億年來(lái)回歸年日數(shù)N的變化, 可以看出在34億年前回歸年日數(shù)等于1000 d左右.

圖3 過(guò)去35億年來(lái)回歸年日數(shù)的變化Fig.3 The days in a tropical year during past 3.5 billion years

4 古生物鐘數(shù)據(jù)與理論數(shù)據(jù)的比較

最重要的古生物鐘是珊瑚、貝類(lèi)以及藍(lán)藻形成的疊層石化石. 某些軟體動(dòng)物的殼保存了一些連續(xù)的層, 它們反映了天文現(xiàn)象的周期變化. 周日結(jié)構(gòu)的生長(zhǎng)與太陽(yáng)的照射以及潮汐高度的變化有關(guān). 古生物鐘化石累積的沉積層密度的變化反映了季節(jié)變化, 從而給出了回歸年日數(shù)的數(shù)據(jù). 沉積層密度的變化也反映了潮汐的逐日變化, 從而給出了朔望月日數(shù)的數(shù)據(jù). 因此測(cè)量古生物鐘的層狀結(jié)構(gòu)可以知道古地質(zhì)時(shí)期時(shí)的回歸年日數(shù)和朔望月日數(shù).

在圖3中, 用不同符號(hào)標(biāo)識(shí)了采用不同的古生物鐘化石求得的回歸年日數(shù)變化. 但是珊瑚和軟體動(dòng)物的結(jié)果只能給到5億年前. 疊層石是最古老的古生物鐘, 至今為止,只有疊層石的資料可以提供在寒武紀(jì)時(shí)反映地球自轉(zhuǎn)歷史的資料. Pannella[19]曾研究最古老的疊層石(28億年前)顯示了由于潮汐變化而產(chǎn)生的周期為10、20和40 d的變化.從Biwabik得到的疊層石化石顯示在20億年前一年有800–900 d[20]. 然而更古老的疊層石的數(shù)據(jù)并沒(méi)有能夠提供季節(jié)變化的數(shù)據(jù). 其原因是這些疊層石化石的層狀圖像, 特別是季節(jié)的圖像十分不清楚.

在太古紀(jì)(28億年前)早期太陽(yáng)的光度是微弱的[21], 但是早期的地球大氣富含CO2、CH4和N2. 由于溫室效應(yīng)地球氣候是溫暖的[22], 使得藍(lán)藻大量增加, 藍(lán)藻的光合作用使得大氣中的二氧化碳含量急劇減少, 氧氣成分增加. 由于藍(lán)藻的沉積和粘合形成了疊層石. 太古紀(jì)的鋯石表明在44億年前地球表面有水[13], 為了解釋疊層石的形成同樣也需要假設(shè)在38億年前地球表面有液態(tài)水的海洋. 疊層石的沉積層提供了朔望月日數(shù), 但是由于有周日潮和半日潮, 在1 d中可形成雙沉積線(xiàn). 因此預(yù)先知道合理的天文韻律, 對(duì)于做出正確判斷是很重要的[11]. 另外一種地球古自轉(zhuǎn)的證據(jù)是細(xì)砂粒沉積形成的韻律.在7.5–6.5億年前, 在澳大利亞Neoproterozoic沙礫的沉淀層, 被解釋為是記錄了古潮汐周期的韻律[23–24], 其結(jié)果顯示在6.5億年前每年有13個(gè)月和450 d.

假定回歸年的長(zhǎng)度是固定的, 那么根據(jù)地月距離數(shù)值, 就可以推算出恒星月的長(zhǎng)度.在本文中已經(jīng)求出了回歸年日數(shù), 因此朔望月的日數(shù)(N2)也能夠計(jì)算出來(lái), 計(jì)算結(jié)果見(jiàn)圖4. 可以看出直到28億年前在N2的變化是很小的, 其最大值是在約18億年前, 一個(gè)朔望月約有32.4 d. 根據(jù)古生物鐘資料得到的朔望月日數(shù)與圖4中的曲線(xiàn)符合很好[19,25].在28億年前N2的數(shù)值很快減少到20 d. Pannella[19]、Cao[26]和Zhu等[27]曾經(jīng)在古老的疊層石化石中發(fā)現(xiàn)了40–45 d的潮汐生長(zhǎng)線(xiàn), 作者認(rèn)為它們很可能是朔望月周期的加倍.

圖4 過(guò)去35億年來(lái)朔望月日數(shù)的變化Fig.4 The days in a synodic month during past 3.5 billion years

利用不同的古生物鐘和沙粒沉淀形成的韻律所求得的回歸年日數(shù)被采用不同的符號(hào)標(biāo)志在圖3上. 可以看出直到8億年前, 古生物鐘的結(jié)果與理論結(jié)果都是符合得很好的.但是8億年前疊層石的數(shù)據(jù)大于理論值. 在20億年前回歸年日數(shù)等于800–900 d太大了,是不可能的. 1972年P(guān)annella[19]曾經(jīng)指出在疊層石的9–10個(gè)沉淀線(xiàn)之后緊隨著4–5個(gè)沉淀線(xiàn)是經(jīng)常發(fā)生的. 由于1 d內(nèi)有兩個(gè)潮峰, 并且在太古紀(jì)疊層石化石的季節(jié)韻律是十分不清楚的, 所以這些疊層石給出的回歸年日數(shù)的變化結(jié)果是值得商榷的, 但是給出的朔望月日數(shù)的變化, 即潮汐周期的變化是很有意義的, 它較好地反映了地球的古自轉(zhuǎn)歷史.

5 總結(jié)

30億年以來(lái)潮汐能量耗散成為影響地月系統(tǒng)演化的重要因素. 地月距離變得越來(lái)越遠(yuǎn), 地球的自轉(zhuǎn)速率也變得越來(lái)越慢. 在本文中給出了地月距離和回歸年日數(shù)的演化曲線(xiàn). 理論數(shù)據(jù)與古生物鐘的結(jié)果直到8億年前都基本上符合.

忽略太陽(yáng)潮引起的能量耗散, 地月系統(tǒng)的角動(dòng)量守恒是正確的. 但是地月系統(tǒng)的能量守恒是不正確的. 由于月球潮汐引起的摩擦使得地球自轉(zhuǎn)角速度變慢, 現(xiàn)在數(shù)值為=?6.145×10?22s?2. 因此地球自轉(zhuǎn)的能量耗散率為?3.6×1019erg/s. 同樣根據(jù)月球平均運(yùn)動(dòng)的加速度0=?25.90′′·cy?2可以算出月球軌道運(yùn)動(dòng)的能量增加值為1.22×1018erg/s. 因此只有約3%的地球自轉(zhuǎn)能轉(zhuǎn)化為月球的軌道運(yùn)動(dòng)能. 絕大部分能量(約97%)耗散在地球表面上.

地月系統(tǒng)的演化過(guò)程是十分復(fù)雜的. 大氣、海洋以及大陸的變化都能夠影響地球自轉(zhuǎn)速率的演化, 例如, 在1991年Young[28]曾經(jīng)指出至少在大陸上有11個(gè)主要的冰期被確認(rèn), 冰期存在的時(shí)間大約占地球歷史的5%–10%. 在冰期時(shí)許多海洋都凍結(jié)了, 所以潮汐能量耗散將變得較小, 因此地月系統(tǒng)的動(dòng)力學(xué)演化過(guò)程將變慢. 在本文中, 只是給出了地月系統(tǒng)的大致演化過(guò)程. 還有許多其他的自然因素也會(huì)影響地月系統(tǒng)的動(dòng)力學(xué)演化,它們都需要更多的古地質(zhì)資料才能估計(jì).致謝 感謝鐘敏教授和張子占博士的幫助.