地球的圈層構造（進階）

為什麼地震波在 660 公里深處會「卡一下」？

讀過入門篇之後，你已經知道地球可以依成分分成地殼、地函、地核，也知道莫氏不連續面(Moho)、古氏不連續面、雷氏不連續面這三道大界線。但如果把一張全球地震波速—深度剖面圖攤開來細看，你會發現事情遠比「三層蛋糕」複雜：在地函內部，波速並不是隨深度平滑增加的，而是在大約 410 公里與 660 公里這兩個深度，各自出現一次明顯的階梯式跳升。地函明明從上到下成分都差不多（都是富鐵鎂的矽酸鹽），為什麼會憑空多出兩道界線？

這個問題的答案，把地球科學從「描述地球有幾層」推進到「理解地球為什麼長這樣」。本篇我們不再重述成分分層，而是聚焦三件入門篇沒細談的事：地函內部的礦物相變(phase transition)如何製造出隱藏的界面、如何用地溫線(geotherm)判斷一個深度是固態還是熔融、以及現代地震學如何用層析成像(tomography)與地球自由振盪(free oscillation)把整顆地球當成一座會發聲的鐘來「聽診」。

地函過渡帶：同一種礦物，換一張臉

先破除一個延伸自入門的迷思。410 與 660 公里的波速跳升，不是成分改變造成的（上下都是橄欖石系統的礦物），而是同一種化學成分的礦物，在高壓下重新排列原子、變成更緊密的晶體結構。化學式沒變，密度卻變了——這就是相變。

地函的主要礦物是橄欖石 $(\mathrm{Mg,Fe})_2\mathrm{SiO}_4$。隨深度增加、壓力升高，它會依序轉變：

• 約 410 公里（壓力約 $14\ \mathrm{GPa}$）：橄欖石(olivine, α 相)轉為瓦茲利石(wadsleyite, β 相)，再到更深處轉為林伍德石(ringwoodite, γ 相)。晶體堆積更緊密，密度與波速階梯式上升。
• 約 660 公里（壓力約 $24\ \mathrm{GPa}$）：林伍德石分解為布里基曼石(bridgmanite，矽酸鎂鈣鈦礦結構) 加 方鎂鐵礦(ferropericlase)。這是一次更劇烈的結構重組，造成 660 公里這道最顯著的不連續面。

夾在 410 與 660 公里之間的這層，地球科學家稱為地函過渡帶(mantle transition zone, MTZ)。它在板塊運動裡扮演關鍵角色：下沉的冷板塊(slab)有時會在 660 公里界面「卡住」暫時攤平，有時則能穿透直墜下部地函——這正是全球地函對流型態的核心爭論之一。

2014 年一項著名研究在巴西鑽石包裹體中找到天然的林伍德石晶體，分析發現其中含有約 1% 重量的水（以氫氧根形式存在於晶格中）。由於過渡帶體積巨大，這暗示地球深部可能鎖著相當於數個海洋質量的「水」——這讓我們對地球水的總量與循環有了全新想像。一顆指甲大的礦物包裹體，竟改寫了整顆行星的水帳。

克勞修斯—克拉佩龍斜率：相變界面為什麼會上下移動

相變界面的深度並非處處固定，而會隨溫度漂移，這對解讀地震影像極其重要。控制相變壓力如何隨溫度變化的，是熱力學裡的克勞修斯—克拉佩龍方程(Clausius–Clapeyron relation)：

$$\frac{dP}{dT} = \frac{\Delta S}{\Delta V}$$

其中 $\Delta S$ 是相變的熵變、$\Delta V$ 是體積變化。這個斜率 $\gamma = dP/dT$ 的正負號，決定了相變界面在「冷區」與「熱區」的反應方向：

• 410 公里相變（橄欖石→瓦茲利石）的克拉佩龍斜率為正（$\gamma \approx +2.9\ \mathrm{MPa/K}$）。在冷的下沉板塊處，較低的溫度讓相變在較淺處就發生，界面被頂高（深度變淺）。
• 660 公里相變（林伍德石分解）的斜率為負（$\gamma \approx -2\ \text{到}\ -3\ \mathrm{MPa/K}$）。在冷板塊處，相變被壓深（界面下凹）。

這個一正一負的組合會產生一個美妙的力學效應：660 公里的負斜率相當於在冷板塊插入處製造了一個浮力屏障，傾向阻擋板塊穿透；而 410 的正斜率則傾向加速下沉。地震層析成像看到的「過渡帶局部增厚或變薄」，其實就是底下溫度異常透過克拉佩龍斜率寫下的指紋——溫度低的地方過渡帶變厚，溫度高的地方變薄。我們因此能用界面的起伏，反推地函深處的冷熱分布。

地溫線 vs. 熔點：判斷一個深度是固是液

入門篇提到「地函是固態、外核是液態」，但為什麼同樣高溫，地函固態、外核卻熔融？關鍵在於比較兩條隨深度變化的曲線：實際溫度分布的地溫線(geotherm)，與物質開始熔化的熔點曲線(solidus / melting curve)。

判準很簡單：

$$ \begin{cases} T_{\text{實際}} < T_{\text{熔點}} & \Rightarrow \text{固態} \\[4pt] T_{\text{實際}} > T_{\text{熔點}} & \Rightarrow \text{液態} \end{cases} $$

地函裡，雖然溫度高達數千度，但熔點隨壓力上升得比地溫線更快，所以地溫線始終壓在熔點曲線之下——固態。到了外核，組成從矽酸鹽換成鐵鎳合金，鐵的熔點曲線位置較低，地溫線此時超過了熔點——液態。再往內核走，壓力極大把鐵的熔點再次抬高到地溫線之上——又變回固態。一條地溫線穿過三條不同的熔點曲線，就解釋了「固—液—固」的順序，這比單純背誦層次深刻得多。

地函對流之所以能進行，還倚賴另一個概念：絕熱地溫線(adiabat)。在對流良好的地函主體，溫度梯度趨近絕熱梯度，可寫成：

$$\left(\frac{dT}{dz}\right)_{\text{ad}} = \frac{\alpha g T}{c_p}$$

其中 $\alpha$ 是熱膨脹係數、$g$ 重力加速度、$c_p$ 定壓比熱。地函絕熱梯度只有約 $0.3\text{–}0.5\ \mathrm{K/km}$，遠小於熱邊界層內陡峭的梯度。換言之，地函中段溫度變化平緩，真正陡的溫度跳變集中在頂部岩石圈與底部 D″ 層這兩個熱邊界層——這正是冷板塊與熱地函柱(mantle plume)誕生的地方。

動手算一下：用波速跳變估算 660 公里界面的密度落差

地震學能測得界面兩側的波速，我們可以反過來估算密度差。已知 660 公里界面附近，P 波速度從約 $v_1 = 10.2\ \mathrm{km/s}$ 跳到 $v_2 = 10.8\ \mathrm{km/s}$。為簡化，假設體積模數 $K$ 與剪切模數大致同比例變化，使得阻抗近似正比於 $\rho v$。界面反射的強弱由地震阻抗(acoustic impedance) $Z = \rho v$ 的對比決定，反射係數為：

$$R = \frac{Z_2 - Z_1}{Z_2 + Z_1} = \frac{\rho_2 v_2 - \rho_1 v_1}{\rho_2 v_2 + \rho_1 v_1}$$

實測 660 界面的 P 波反射係數約 $R \approx 0.05$。設密度由 $\rho_1 = 3.99\ \mathrm{g/cm^3}$ 變為 $\rho_2$，代入：

$$0.05 = \frac{\rho_2 (10.8) - (3.99)(10.2)}{\rho_2 (10.8) + (3.99)(10.2)}$$

令 $a = 10.8\,\rho_2$、$b = 40.70$，則 $0.05(a+b) = a - b$，整理得 $0.95\,a = 1.05\,b$，即 $a = 1.105\,b$：

$$10.8\,\rho_2 = 1.105 \times 40.70 = 44.97 \quad\Rightarrow\quad \rho_2 \approx 4.16\ \mathrm{g/cm^3}$$

密度從 $3.99$ 跳到約 $4.16\ \mathrm{g/cm^3}$，增加約 4%。這個量級與礦物物理實驗測得的林伍德石分解密度躍升一致——一道地震反射的強弱，竟能轉譯成深部礦物的緻密化程度。這就是把地震觀測、熱力學與礦物物理三者扣在一起的力量。

D″ 層與後鈣鈦礦：核—函邊界上的另一層皮

入門篇把核—函邊界(CMB, 約 2900 公里)當成一條乾淨的線，但現代研究發現緊貼 CMB 上方約 200–300 公里厚的地函底部，存在一層性質異常、極不均勻的D″ 層(D-double-prime layer)。這裡是地核（約 $3500\ \mathrm{K}$ 以上）與地函之間最劇烈的熱邊界層，溫度在數百公里內陡升上千度。

2004 年的高壓實驗發現了關鍵拼圖：布里基曼石在 CMB 附近的極端壓力（約 $125\ \mathrm{GPa}$）下，會再次相變為更緻密的後鈣鈦礦(post-perovskite) 結構。這個相變的克拉佩龍斜率為正且陡，因此在較冷的區域（如古老隱沒板塊堆積處）後鈣鈦礦穩定，在較熱的區域又轉回布里基曼石——這解釋了 D″ 層為何有些地方出現地震波速跳升、有些地方沒有。D″ 並不是一道整齊的界面，而是隨溫度斑駁出現的「補丁」。

D″ 之所以重要，是因為它是地函對流的「下熱邊界層」，被認為是地函柱(mantle plume) 的根源。夏威夷、冰島這類「熱點(hotspot)」火山，其岩漿來源可能就一路扎根到 CMB。地球最深的地函，與地表最壯觀的火山，竟由同一根熱柱牽連。

重點回顧

• 410 與 660 公里的波速不連續面源於礦物相變（橄欖石→瓦茲利石→林伍德石→分解為布里基曼石＋方鎂鐵礦），化學成分不變、晶體結構變緻密，與成分分層是兩回事。
• 克勞修斯—克拉佩龍方程 $dP/dT=\Delta S/\Delta V$ 的斜率正負，決定相變界面在冷熱區的起伏：410 為正、660 為負，後者形成阻擋板塊穿透的浮力屏障。
• 一個深度是固是液，取決於地溫線與熔點曲線的相對位置；地函主體近絕熱梯度，陡峭溫變集中在頂部與底部兩個熱邊界層。
• D″ 層與後鈣鈦礦相變揭示核—函邊界並非乾淨界線，而是地函柱（熱點火山源頭）誕生的斑駁熱邊界層。
• 界面的地震反射係數 $R=(Z_2-Z_1)/(Z_2+Z_1)$ 可把波速—密度跳變量化，連結地震觀測與礦物物理。

深入探討（研究所視角）

地震層析成像：把地球當成醫學斷層掃描

入門篇用「陰影區」定性地確認了地核存在，研究所層級則用地震層析成像(seismic tomography) 把地函內部畫成三維溫度／速度地圖。其數學骨架是一個大型線性反問題。觀測量是大量地震—測站對的走時殘差 $\delta t_i$（實測走時與一維參考模型 PREM 預測的差），未知量是各網格的慢度擾動 $\delta s_j = \delta(1/v)$，沿射線路徑 $L_i$ 積分：

$$\delta t_i = \int_{L_i} \delta s(\mathbf{r})\, dl \;\approx\; \sum_j G_{ij}\,\delta s_j$$

寫成矩陣形式 $\mathbf{d} = \mathbf{G}\,\mathbf{m}$，但 $\mathbf{G}$ 通常病態(ill-posed)、射線覆蓋不均，需以正則化(regularization) 求解，例如阻尼最小平方：

$$\mathbf{m} = (\mathbf{G}^\top\mathbf{G} + \lambda^2 \mathbf{I})^{-1}\mathbf{G}^\top\mathbf{d}$$

$\lambda$ 是阻尼參數，權衡資料擬合與模型平滑度。現代「全波形反演(full-waveform inversion)」則進一步用整段地震波形、配合伴隨法(adjoint method)計算敏感度核(sensitivity kernel)，解析度遠勝傳統走時層析。這套方法揭露了 CMB 上方兩塊巨大的低速異常區——位於非洲與太平洋之下的大型低剪切波速省(Large Low-Shear-Velocity Provinces, LLSVPs)，可能是化學成分異於周圍地函的古老堆積體，是當前深部地球動力學最熱門的謎團。

地球自由振盪：把整顆行星當成一座鐘

巨大地震（如 $M_w \gtrsim 8.5$）會激發整顆地球的自由振盪(free oscillation / normal modes)，使地球像被敲響的鐘一樣，以一組離散的本徵頻率持續振盪數天。這些簡正模分為球型模(spheroidal, $_nS_l$)與環型模(toroidal, $_nT_l$)，由地球的徑向密度與彈性結構決定本徵頻率。求解的本徵值問題形式上類似量子力學的束縛態：

$$\mathcal{H}[\,u\,] = \omega^2 \rho\, u$$

其中 $\mathcal{H}$ 是含彈性模數與重力的積分微分算符，$\omega$ 為本徵角頻率，$u$ 為位移本徵函數。最低階的「呼吸模」$_0S_0$ 週期約 20.5 分鐘，整顆地球同步脹縮。透過比對觀測到的數千個本徵頻率與理論值，地球科學家反演出徑向參考模型PREM(Preliminary Reference Earth Model, 1981)——它對密度的約束是純走時資料給不出的，因為簡正模對全球質量分布敏感。值得一提的是，自轉與橢圓度會使本徵頻率分裂(splitting)，類似原子譜線在磁場中的塞曼效應(Zeeman effect)，這個分裂正是探測內核差異自轉與各向異性的精細工具。

連到台灣：密集地震網作為深部探針

台灣位處板塊縫合帶，菲律賓海板塊向西北隱沒，提供了密集且多角度的震源。結合台灣與全球寬頻地震網，研究者得以對東亞地函進行高解析層析，追蹤隱沒板塊在 660 公里界面究竟是停滯攤平、還是穿透墜入下部地函。下次讀到一則地震速報時，不妨想想：那組波不只定位了震央，也正在替全人類更新對地球這座「會發聲的鐘」最深處的認識。