板塊構造學說（進階）：尤拉極、半空間冷卻與地函對流的力學

板塊在球面上其實「不會直線移動」——你想過嗎？

入門篇裡，我們把板塊想像成在一張平面地圖上滑來滑去的拼圖。這個圖像對於建立直覺很好用，但它藏著一個會誤導人的簡化：地球是球，板塊不是在平面上平移，而是繞著某個軸在球面上旋轉。一旦你接受這件事，許多原本只能死背的現象——為什麼同一條轉形斷層上不同位置的錯動速率不一樣、為什麼太平洋板塊在不同緯度跑的速度差這麼多、為什麼擴張速率沿著中洋脊會系統性變化——立刻就有了一個統一的數學解釋。

這篇進階篇要做的，是把入門篇的「定性故事」升級成「定量機制」。我們會用尤拉極(Euler pole) 描述板塊運動學，用半空間冷卻模型(half-space cooling model) 解釋海底為什麼越老越深，用驅動力的力平衡重新檢視 slab pull 與 ridge push 誰才是老大，最後用瑞利數(Rayleigh number) 把地函對流放回物理脈絡。已經讀過入門的你，準備好把板塊構造從「學說」看成「一套可以計算的力學系統」。

尤拉定理：板塊運動的數學骨架

剛體在球面上的任何運動，都可以表示成繞著一個通過球心的軸的旋轉——這是 18 世紀數學家尤拉(Leonhard Euler) 證明的幾何定理，1960 年代被引入板塊構造，成為描述板塊運動的標準語言。

每一個板塊（相對於另一個板塊）的運動，由兩個量完全決定：

• 尤拉極(Euler pole)：旋轉軸與地表的交點，也叫旋轉極(pole of rotation)。它通常不在地理北極，而是某個特定的經緯度。
• 角速度(angular velocity, $\omega$)：板塊繞此軸旋轉的速率，單位常用「度／百萬年」。

關鍵在於：板塊上任一點 P 的線速度，不是常數，而是隨它到尤拉極的角距離 $\theta$ 變化：

$$v = \omega \, R \sin\theta$$

其中 $R$ 是地球半徑（約 6371 公里），$\theta$ 是 P 點與尤拉極之間的角距離。

這個 $\sin\theta$ 因子是整篇文章最重要的洞見之一。它告訴我們：

• 在尤拉極附近（$\theta \to 0$），$\sin\theta \to 0$，板塊幾乎不動——這裡只是純粹的「轉」。
• 在離尤拉極 90° 的「旋轉赤道」上（$\theta = 90°$），$\sin\theta = 1$，線速度達到最大。

所以同一個板塊在不同位置移動速率不同，並不是因為板塊在「變形」，而純粹是球面旋轉的幾何結果。轉形斷層之所以一定沿著「以尤拉極為中心的小圓」延伸，也是因為相鄰兩板塊只能繞共同尤拉極相對轉動——斷層必須與相對運動方向平行，而那個方向處處是小圓的切線。從衛星測得的轉形斷層走向反推尤拉極位置，正是現代板塊運動學的標準操作。

海底越老越深：半空間冷卻模型

入門篇說「離中洋脊越遠、海底越老」。進階篇要追問：為什麼越老的海底也越深？ 大西洋中洋脊軸部水深約 2.5 公里，但兩側老海盆卻深達 5–6 公里。這不是巧合，而是一條漂亮的物理定律在背後運作。

新生的海洋岩石圈在中洋脊又熱又輕。當它隨擴張往兩側移動、逐漸冷卻，會發生兩件事：變重（熱收縮使密度增加）和變厚（冷卻鋒面往下傳，把更多地函「凍結」成岩石圈）。根據地殼均衡(isostasy)，變重的柱體會下沉得更深。把熱傳導方程式（一維冷卻）解出來，可得到著名的結論：海底深度隨年齡的平方根增加。

$$d(t) = d_0 + C\sqrt{t}$$

其中 $d_0$ 是中洋脊軸部水深，$t$ 是海洋岩石圈的年齡，$C$ 是一個由地函溫度、熱膨脹係數與熱擴散率決定的常數，實測約為 $350 \text{ m}/\sqrt{\text{Myr}}$（百萬年）。

同樣地，岩石圈的厚度也隨 $\sqrt{t}$ 成長：

$$L(t) \approx 2.3\sqrt{\kappa\, t}$$

$\kappa$ 是熱擴散率（約 $10^{-6}\ \mathrm{m^2/s}$）。這就是為什麼老的海洋板塊又厚又重——而「又厚又重」正是它最終會在隱沒帶往下沉的根本原因。冷卻、增厚、變重、隱沒，是同一條物理鏈條。

動手算一下：6 千萬年的海底會比中洋脊深多少？

太平洋上一塊年齡 $t = 64$ Myr 的海洋岩石圈，比中洋脊軸部深多少？用平方根冷卻關係：

$$\Delta d = C\sqrt{t} = 350 \text{ m}/\sqrt{\text{Myr}} \times \sqrt{64\ \text{Myr}}$$

$$\Delta d = 350 \times 8 = 2800 \text{ m}$$

也就是說，這塊老海底比中洋脊軸部深了約 2.8 公里。若軸部水深 2.5 公里，這裡就接近 5.3 公里——和大洋實測深度高度吻合。順帶估一下這塊岩石圈的厚度：

$$L \approx 2.3\sqrt{\kappa t} = 2.3\sqrt{10^{-6}\ \mathrm{m^2/s} \times 64 \times 3.15\times10^{13}\ \mathrm{s}} \approx 2.3 \times 4.5\times10^{4}\ \mathrm{m} \approx 100\ \mathrm{km}$$

剛好就是教科書上「岩石圈厚約 100 公里」那個數字的來源——它不是憑空規定的，而是冷卻時間自然長出來的。注意這個 $\sqrt{t}$ 律在年齡超過約 8 千萬年後會逐漸偏離（海底沒有模型預測的那麼深），暗示老的岩石圈底部可能被地函的小尺度對流「重新加熱」，這正是研究前沿之一。

重新檢視驅動力：slab pull 為什麼是老大

入門篇提到 slab pull（隱沒板塊拉力）與 ridge push（中洋脊推力）。進階篇要把它們量化比較，並說明為什麼現代研究認為 slab pull 主導。

關鍵物理量是負浮力(negative buoyancy)。隱沒的板片又冷又緻密，與周圍較熱的地函有溫差 $\Delta T$，因此產生密度差：

$$\Delta\rho = \rho\,\alpha\,\Delta T$$

$\alpha$ 是熱膨脹係數（約 $3\times10^{-5}\ \mathrm{K^{-1}}$）。單位長度板片受到向下的拉力正比於這個密度差、板片厚度與隱沒深度。對一塊冷卻成熟（厚約 100 km、$\Delta T$ 可達數百度）的隱沒板片，估算出的 slab pull 量級約為 $3\text{–}5\times10^{13}\ \mathrm{N/m}$。

相對地，ridge push 來自中洋脊地形高差造成的側向壓力梯度（本質上是岩石圈內部的重力位能釋放），量級約 $2\text{–}3\times10^{12}\ \mathrm{N/m}$——比 slab pull 小了大約一個數量級。

這個量級差解釋了一個觀測事實：有隱沒帶連接的板塊跑得快，沒有隱沒帶的板塊跑得慢。

• 太平洋板塊、納斯卡板塊、可可板塊：邊緣大量隱沒，速率高達 7–10 cm/yr。
• 北美板塊、歐亞板塊、非洲板塊：以大陸為主、隱沒邊界少，速率僅 1–3 cm/yr。

如果 ridge push 才是主力，這個對比就不該存在。當然，真實的力平衡還要扣掉地函對板片運動的黏滯阻力(slab resistance) 與板塊底部的拖曳(basal drag)。地函對流到底是「推手」還是「阻力」，取決於對流胞與板塊運動方向是否一致——這也是為什麼現代圖像不再說「板塊被對流拖著走」，而說「板塊本身就是地函對流的上邊界冷卻層」：板塊不是被動的乘客，它就是對流系統的一部分。

地函對流的物理量尺：瑞利數

「地函會不會對流」不是靠想像，而是有明確的判準。流體在底部加熱、頂部冷卻時，會不會啟動對流，取決於一個無因次量——瑞利數(Rayleigh number, $Ra$)：

$$Ra = \frac{\rho\, g\, \alpha\, \Delta T\, D^3}{\kappa\, \eta}$$

其中 $g$ 是重力加速度，$\Delta T$ 是上下溫差，$D$ 是地函厚度（約 2900 km），$\kappa$ 是熱擴散率，$\eta$ 是黏滯度。當 $Ra$ 超過一個臨界值（對地函這種上下邊界條件約為 $10^3$ 量級）時，浮力勝過黏滯與熱擴散的阻尼，對流就會發生。

把地函的參數代進去：$\Delta T$ 約數千度、$D \sim 3\times10^6$ m、$\eta \sim 10^{21}\ \mathrm{Pa\cdot s}$，算出來的 $Ra$ 高達 $10^6$–$10^7$，遠遠超過臨界值。所以地函必然在對流——這是一個可以「算」出來的結論，不是猜的。$Ra$ 越大、對流越劇烈，這也解釋了為什麼年輕、更熱（$\Delta T$ 更大、$\eta$ 更小）的地球，板塊運動可能比今天更快。

這裡要破除一個常見迷思：地函「對流」不代表它是液體。地函固態岩石的黏滯度高達 $10^{21}\ \mathrm{Pa\cdot s}$，是水的 $10^{24}$ 倍。它以固態潛變(solid-state creep) 的方式流動——原子在晶格中緩慢擴散、晶格錯位逐步移動，使整體在百萬年尺度下表現得像極黏稠的流體。每年幾公分的板塊速率，本質上就是這種超慢固態流動的地表表現。

看一個例子：用尤拉極算花蓮外海的聚合速率

讓我們把尤拉極的公式用到台灣。菲律賓海板塊相對歐亞板塊的尤拉極，根據 GPS 與地質資料約位於北緯 50° 附近、角速度 $\omega \approx 1.0°/\mathrm{Myr}$（不同模型略有差異，這裡取代表值做數量級估算）。

先把角速度換成弧度／年：

$$\omega = 1.0°/\mathrm{Myr} = \frac{1.0 \times \pi/180}{10^6}\ \mathrm{rad/yr} \approx 1.75\times10^{-8}\ \mathrm{rad/yr}$$

花蓮外海約位於北緯 24°，與北緯 50° 的尤拉極角距離 $\theta \approx 26°$，$\sin\theta \approx 0.44$。代入 $v = \omega R \sin\theta$：

$$v = 1.75\times10^{-8}\ \mathrm{rad/yr} \times 6.371\times10^{6}\ \mathrm{m} \times 0.44$$

$$v \approx 0.049\ \mathrm{m/yr} \approx 4.9\ \mathrm{cm/yr}$$

這個數量級和入門篇提到的「每年約 8 公分」、以及 GPS 實測的數十 mm/yr 是相符的（真實聚合率還疊加了島弧本身變形與多塊微板塊運動，所以實測更高）。重點是：同一個尤拉極，套到不同緯度的位置，就能解釋為什麼台灣不同段落的聚合速率不一樣——北部琉球海溝與南部馬尼拉海溝的速率差異，部分正是 $\sin\theta$ 隨位置變化的幾何效應。這就是把「板塊在球面上旋轉」這個觀念落實到台灣的具體威力。

重點回顧

• 板塊在球面上的運動是繞尤拉極(Euler pole) 的剛體旋轉，任一點線速度 $v = \omega R\sin\theta$ 隨到極的角距離變化——這解釋了沿邊界速率為何系統性改變，以及轉形斷層為何沿小圓延伸。
• 海洋岩石圈冷卻時深度與厚度都隨年齡的平方根成長（$d \propto \sqrt{t}$、$L \propto \sqrt{t}$）；冷卻→增厚→變重，是它最終隱沒的物理根源。
• 力平衡上，slab pull 比 ridge push 大約一個數量級，這由「有隱沒帶的板塊跑得快」的觀測直接印證；現代圖像視板塊為地函對流系統的冷卻上邊界，而非被動乘客。
• 瑞利數 $Ra$ 是地函是否對流的判準；地函的 $Ra \sim 10^6$–$10^7$ 遠超臨界值，必然對流，但它以固態潛變流動，黏滯度高達 $10^{21}\ \mathrm{Pa\cdot s}$，並非熔融液體。
• 把尤拉極公式套到台灣，可從幾何上估算並理解不同段落聚合速率的差異，是球面板塊運動學的在地應用。

深入探討（研究所視角）

流變學分層：為什麼板塊是「板塊」

板塊之所以能當成剛體處理，根源在於地函的流變學(rheology) 是強烈溫度相依的。固態潛變的應變率近似遵循 Arrhenius 型關係：

$$\dot\varepsilon \propto \sigma^n \exp\!\left(-\frac{Q}{RT}\right)$$

其中 $\sigma$ 是應力、$n$ 是應力指數（位錯潛變約 $n\approx3$）、$Q$ 是活化能、$T$ 是絕對溫度。指數項 $\exp(-Q/RT)$ 對溫度極度敏感：冷的岩石圈（低 $T$）黏滯度高出好幾個數量級，於是表現得「硬」；熱的軟流圈（高 $T$）黏滯度驟降，可以流動。所以岩石圈／軟流圈的力學分界，本質上是等溫面而非化學成分界線——這也是為什麼「板塊厚度」會隨冷卻時間（年齡）成長。把這套流變律放進地函對流的數值模擬，正是當代地球動力學的核心工作。

板塊邊界的應變不是一條線：縫合帶的有限變形

入門篇把邊界畫成一條線，但在大陸碰撞與弧陸碰撞區（如台灣），變形其實分布在數十到上百公里寬的帶狀區域。研究這種寬變形帶常用臨界楔(critical taper / critical wedge) 理論：把造山帶想成推土機前方堆起的土楔，楔形的坡度由內部摩擦與底部滑脫面摩擦的平衡決定。台灣的西部麓山帶就是一個近乎教科書級的臨界楔，它的地形坡度、逆衝斷層序列與滑脫面深度，都能用這套力學框架定量預測。配合 GPS 與大地測量的應變率場，可以區分哪些斷層正在累積彈性應變、哪些正在潛移，是地震潛勢評估的前沿。

板片在地函深處的命運：660 公里不連續面與超級地函柱

隱沒板片往下沉時，會在約 660 公里深的相變界面（橄欖石相變為布氏岩 bridgmanite + 鎂方鐵礦的密度躍升處）遇到阻力。地震波層析成像(seismic tomography)顯示，有些板片直接穿透 660 km 界面墜入下部地函、堆積在核幔邊界(core–mantle boundary)上方，形成所謂的「板片墳場(slab graveyard)」；有些則暫時滯留、堆疊在界面上方形成「平躺板片(stagnant slab)」——西太平洋（包含日本、台灣周邊）底下就觀測到顯著的滯留板片。

核幔邊界附近還存在兩個橫跨數千公里的大型低剪波速區(Large Low Shear Velocity Provinces, LLSVPs)，位於非洲與太平洋之下，可能是冷板片長期堆積與化學分異的產物，並被認為與地函柱(mantle plume) 的根源、超大陸的聚合與裂解節律有關。這把板塊構造從「地表現象」延伸成一個貫穿核幔邊界到地表的全地函尺度系統(whole-mantle convection)。地表每年幾公分的板塊運動，與 2900 公里深處冷板片的沉降、熱地函柱的上湧，其實是同一台巨大熱機的不同切面——這正是當代固體地球科學試圖整合的終極圖像。