風化、侵蝕與地形（進階）：把峽谷寫成方程式

一條河每年能搬走一座山嗎？

入門篇我們站在太魯閣的步道上，知道了「抬升快、侵蝕也快」這場拉鋸戰。但若有人追問：立霧溪「到底」每年從中央山脈搬走多少噸的岩石？這個數字憑什麼能跟造山運動的抬升「打平」？河水又是依據什麼物理法則，決定哪些石頭搬得動、哪些搬不動？

這些問題把我們從「敘述性地形學(descriptive geomorphology)」推向「定量地形學(quantitative geomorphology)」的世界。在這裡，峽谷不再只是「長期切割的結果」，而是一組可以寫成微分方程、可以用宇宙射線同位素實測、可以放進質量守恆方程式去檢核的動力系統。這篇進階篇要做的，就是把入門篇那句「抬升快侵蝕也快」拆開來，看清楚它背後的搬運門檻、侵蝕定律與造山—剝蝕耦合機制。我們會用到一些公式，但每一條都對應一個你在野外看得到的現象。

顆粒搬得動嗎？臨界剪應力與 Shields 參數

入門篇說「顆粒越大需要的能量越高」，這句話可以寫成精確的物理門檻。河水流過河床時，會對底部顆粒施加一個 底床剪應力(bed shear stress) $\tau_b$。對等深度 $h$、坡度 $S$ 的明渠均勻流，這個剪應力是

$$\tau_b = \rho g h S$$

其中 $\rho$ 是水的密度、$g$ 是重力加速度。注意這裡 $h S$ 的組合：深水或陡坡都會放大剪應力，這解釋了為什麼颱風暴漲（$h$ 暴增）的山區陡流（$S$ 大）搬運力會突然失控。

一顆躺在河床上的顆粒要被「啟動」，剪應力必須超過某個 臨界剪應力(critical shear stress) $\tau_c$。19 世紀末到 20 世紀初的水利學家把這件事無因次化，得到著名的 Shields 參數(Shields parameter)：

$$\tau_* = \frac{\tau_b}{(\rho_s - \rho)\,g\,D}$$

$\rho_s$ 是顆粒密度（石英約 $2650\ \text{kg/m}^3$），$D$ 是顆粒直徑。當 $\tau_*$ 超過臨界值 $\tau_{*c}$（粗顆粒約 $0.045\sim0.06$），顆粒開始移動。把這個條件反解，可以得到「某流量下能搬動的最大粒徑」：

$$D_{\max} = \frac{\tau_b}{(\rho_s-\rho)\,g\,\tau_{*c}} = \frac{\rho\,h\,S}{(\rho_s-\rho)\,\tau_{*c}}$$

這條式子是理解河床粒徑分布的鑰匙。上游 $S$ 大，$D_{\max}$ 大，所以河床鋪滿大礫石；下游 $S$ 趨近於零，$D_{\max}$ 驟減，只剩細砂與泥——這就是入門篇講的「粗顆粒先沉、細顆粒後沉」背後真正的力學原因，而不只是「流速變慢」這種定性說法。

值得一提的是，入門篇提到的 Hjulström 曲線 之所以呈 U 形（極細的黏土反而需要更高流速才能啟動），原因正是黏土顆粒之間有 凝聚力(cohesion)，這部分 Shields 的純重力分析並未涵蓋——這也是定量模型的邊界：模型乾淨，但真實泥砂總有它不肯被簡化的地方。

動手算一下：立霧溪能搬多大的石頭？

假設立霧溪某段在豪雨時水深 $h = 3\ \text{m}$、坡度 $S = 0.01$（1%），取 $\tau_{*c} = 0.05$、$\rho = 1000$、$\rho_s = 2650\ \text{kg/m}^3$、$g = 9.8\ \text{m/s}^2$。先算底床剪應力：

$$\tau_b = \rho g h S = 1000 \times 9.8 \times 3 \times 0.01 = 294\ \text{Pa}$$

再算最大可動粒徑：

$$D_{\max} = \frac{\tau_b}{(\rho_s-\rho)\,g\,\tau_{*c}} = \frac{294}{(2650-1000)\times 9.8 \times 0.05} \approx 0.36\ \text{m}$$

也就是說，這樣一場豪雨足以滾動直徑約 36 公分 的大礫石。若水深翻倍到 $6\ \text{m}$，$\tau_b$ 也翻倍，$D_{\max}$ 直接到約 $0.73\ \text{m}$——這就是為什麼颱風後我們會在河床看到平時動也不動的巨石被搬到下游。搬運力對水深的線性放大，是台灣河流「平時溫吞、暴雨時暴烈」的數學寫照。

河流為何能下切？兩種限制與河川功率

入門篇把「下切」當成既定事實，但下切的速率到底受什麼控制？地形學把基岩河道的侵蝕分成兩種機制：

• 搬運限制(transport-limited)：河床已經堆滿可動的沉積物，侵蝕速率取決於水流「搬得走多少」。這類河道（多見於沖積平原）的演育由輸砂能力主導。
• 分離限制(detachment-limited)：河床是裸露基岩，物質要先被「鑿離」岩體才談得上搬走，侵蝕速率取決於水流「啃得動多少」。台灣這種陡峭、沉積物常被瞬間沖光的基岩峽谷（太魯閣即是），多屬此類。

對分離限制的基岩河道，最廣為使用的是 河川功率侵蝕模型(stream power incision model, SPIM)。其核心量是 單位河川功率(unit stream power)——單位河床面積上水流耗散的位能率。經過推導，下切速率 $E$ 可寫成簡潔的冪律：

$$E = K\,A^{m}\,S^{n}$$

這裡 $A$ 是上游 集水面積(drainage area)（用來代理流量，因為流量大致與集水面積成正比）、$S$ 是河道坡度、$K$ 是綜合了岩性、氣候、水力幾何的 侵蝕係數，$m$、$n$ 是正指數（野外常見 $m/n \approx 0.4\sim0.6$）。這條式子告訴我們：集水面積越大、坡度越陡，下切越快。

看一個例子：坡度—面積關係與裂點

在地形達到「均夷(graded)」的動態平衡時，沿河任一點的下切速率都等於抬升速率 $U$。令 $E = U$ 代入 SPIM 並解出坡度：

$$S = \left(\frac{U}{K}\right)^{1/n} A^{-m/n}$$

兩邊取對數：

$$\log S = \frac{1}{n}\log\!\frac{U}{K} - \frac{m}{n}\log A$$

這是一條直線！把一條河流各點的 $\log S$ 對 $\log A$ 畫出來（稱為 坡度—面積圖），平衡河段會落在斜率為 $-m/n$ 的直線上。更妙的是截距項含有 $U/K$：抬升速率越高，整條直線往上平移。地形學家正是用這個關係從 DEM（數值地形模型）反推一個流域的抬升或岩性差異。

當這條直線出現「斷階」——某一點以下坡度突然變陡——往往就是入門篇提過的 裂點(knickpoint)，代表上游與下游處在不同的均衡狀態，是構造抬升加速或基準面下降的訊號正沿河往上游傳遞。台灣許多河流的裂點分布，被研究者拿來重建近百萬年的抬升史。

山坡怎麼變平滑？地貌的「擴散」

河流負責「切」，但兩條河之間那片廣闊的 山坡(hillslope) 怎麼演育？對緩坡而言，土壤潛移、雨濺、生物擾動等過程的淨效果，可以驚人地用一條 擴散方程式 描述。令地表高程為 $z$，緩坡上的物質通量 $q_s$ 正比於坡度：

$$q_s = -D\,\nabla z$$

$D$ 是 地貌擴散係數(landscape diffusivity)。再套用質量守恆（高程變化 = 抬升 − 通量散度），得到山坡演育方程：

$$\frac{\partial z}{\partial t} = U + D\,\nabla^2 z$$

這跟物理裡的熱傳導方程、化學裡的擴散方程是同一個數學形式！它的意涵很直觀：山坡的凸起處（$\nabla^2 z < 0$）會被削低，凹陷處會被填高，整片地形隨時間趨於平滑——這正是入門篇「削高填低」在山坡尺度的精確版本。

但這個線性擴散只在緩坡成立。一旦坡度逼近 臨界坡度(critical slope angle，約莫摩擦角 $30^\circ\sim40^\circ$），通量會急遽發散、山坡以崩塌(landslide)的方式整批下滑。台灣山區坡度普遍逼近甚至超過臨界值，所以這裡的山坡演育不是溫和的擴散，而是由 地震與颱風觸發的崩塌 主導——這也是台灣崩塌災害如此頻繁的地形學根源。

量出來的侵蝕：宇宙射線同位素

講了這麼多公式，有沒有辦法「直接量」一個流域實際被剝蝕了多快？有，而且是過去三十年地形學最重要的工具之一：宇宙生成核種(cosmogenic nuclide) 定年。

高能宇宙射線打進地表岩石，會在石英的氧、矽原子核內撞出稀有的同位素，最常用的是 $^{10}\text{Be}$（鈹-10）。關鍵在於：宇宙射線只能穿透地表約 數十公分到 1 公尺，所以一顆石英砂粒中累積的 $^{10}\text{Be} $ 濃度，反映了它「待在近地表」的時間長短。

• 若一地剝蝕得 慢，岩石長時間停留在近地表，$^{10}\text{Be}$ 累積得 多。
• 若剝蝕得 快，深處岩石很快被翻到表面又被剝走，曝曬時間短，$^{10}\text{Be}$ 濃度 低。

在穩定狀態下，剝蝕速率 $\varepsilon$ 與表面核種濃度 $N$ 近似成反比：

$$N \approx \frac{P_0}{\dfrac{\lambda}{ } + \dfrac{\rho\,\varepsilon}{\Lambda}} \quad\xrightarrow{\ \text{對 } ^{10}\text{Be} \text{ 半衰期長、忽略衰變}\ }\quad \varepsilon \approx \frac{P_0\,\Lambda}{\rho\,N}$$

其中 $P_0$ 是地表產率、$\Lambda$ 是吸收深度尺度（約 $600\ \text{kg/m}^2$）、$\rho$ 是岩石密度。妙處在於：只要採一把 河砂，量其中石英的 $^{10}\text{Be}$，就能得到整個上游流域 千年到萬年尺度的平均剝蝕速率——這把空間與時間都平均掉了，避開了單點測量的偶然性。台灣的研究用這套方法測出部分流域剝蝕速率高達每年數公釐，是全球少見的「快速剝蝕造山帶」，與入門篇估算的抬升速率同一量級——這就回答了開篇的問題：抬升與剝蝕「打平」，不是巧合，是被野外數據證實的動態平衡。

山為什麼長不高？造山—剝蝕—均衡的耦合

最後我們把尺度拉到整座山脈。一個直覺的迷思是：板塊一直擠，山就會一直長高。實際上山脈高度會趨於一個 穩定態(steady state)，原因是三件事彼此鎖死：

1. 構造抬升 把物質往上推。
2. 剝蝕 隨高度（與坡度、降雨）增強，把物質往下削。
3. 地殼均衡(isostasy)：地殼像漂在地函上的木塊，山根浮力支撐著山體。當頂部被剝走重量，山根會 反彈抬升(isostatic rebound) 補回大部分高度。

這裡有個常被忽略的量化事實：地殼密度約 $2.7\ \text{g/cm}^3$、地函約 $3.3\ \text{g/cm}^3$，由阿基米德原理，剝蝕掉 $1\ \text{km}$ 厚的岩石，山根反彈會補回約

$$\Delta h_{\text{rebound}} = 1\ \text{km} \times \frac{\rho_{\text{crust}}}{\rho_{\text{mantle}}} \approx 1 \times \frac{2.7}{3.3} \approx 0.82\ \text{km}$$

也就是說，山頂淨降低只有約 $0.18\ \text{km}$。這解釋了一個驚人的後果：侵蝕削掉的大多是被均衡反彈補回的高度，淨削低很少，所以山峰反而能比「不考慮均衡」時 存在更久、甚至局部更高。當抬升、剝蝕、均衡三者達成平衡，山脈進入穩定態：高度不再增加，但物質仍源源不絕地「過境」——岩石從深處被抬升、剝蝕、搬運入海，整座山像一條輸送帶，形狀不變而物質長流。台灣正處在這個耦合系統最劇烈的一端，是全球研究造山—剝蝕耦合的天然實驗室。

重點回顧

• 顆粒能否被搬動由 Shields 參數 決定臨界門檻；可動最大粒徑 $D_{\max} \propto hS$，這是河床上游粗、下游細的力學根源，也解釋了暴雨時搬運力為何暴增。
• 基岩河道的下切多屬 分離限制，可用 河川功率模型 $E = K A^m S^n$ 描述；在均衡態下 $\log S$ 對 $\log A$ 呈斜率 $-m/n$ 的直線，裂點 是抬升或基準面變化沿河上溯的訊號。
• 緩坡演育服從 擴散方程 $\partial z/\partial t = U + D\nabla^2 z$（削凸填凹）；但坡度逼近臨界角時改由崩塌主導，這正是台灣的常態。
• 宇宙生成核種（$^{10}\text{Be}$） 能直接實測流域千年尺度的平均剝蝕速率，證實台灣是全球剝蝕最快的造山帶之一。
• 山脈高度受 抬升—剝蝕—地殼均衡 三者耦合而趨於穩定態；均衡反彈會補回約 8 成被剝蝕掉的高度，使山「淨削低很慢、存在很久」。

深入探討（研究所視角）

構造氣候耦合與「降雨能否改變造山帶形狀」。 過去二十年地形學最熱的辯論之一，是 降雨（氣候）能否回過頭來控制構造。SPIM 中的侵蝕係數 $K$ 強烈依賴降雨，於是出現一個誘人假說：迎風坡降雨多、剝蝕強，會像「地形吸塵器(tectonic aneurysm)」把深部物質快速抽到地表，造成局部變質作用與抬升集中——喜馬拉雅 Nanga Parbat、台灣的部分研究都被援引討論。但反方指出，剝蝕速率在快速造山帶常被 構造抬升速率 設定上限，氣候訊號可能被均衡與物質供應「緩衝」掉。這場爭論至今未定，核心難點在於把 $K$ 的氣候依賴從岩性、植被、地震觸發的崩塌中乾淨地分離出來。

隨機—閾值侵蝕(stochastic-threshold incision)。 經典 SPIM 用「平均」流量，但台灣的真相是侵蝕高度集中在少數颱風事件。新一代模型把 SPIM 改寫成對 流量機率分布 的積分，並加入侵蝕門檻 $\tau_c$：只有超過門檻的洪水才做功。其深刻後果是——在門檻效應下，降雨變率(variability) 可能比降雨「總量」更能決定長期侵蝕速率。一個雨量總和相同、但偶有極端暴雨的氣候，會比細水長流的氣候侵蝕得更猛。這對在暖化背景下評估台灣坡地災害與河川輸砂，有直接的政策意涵。

沉積物收支與「源—匯系統(source-to-sink)」。 把單一河道放大到整個系統，現代研究關心的是 質量守恆下的沉積物收支(sediment budget)：山脈剝蝕產出的物質，經河流搬運、暫存於沖積平原（停留時間可達千年），最終匯入海盆。台灣以極小的陸地面積（約佔全球陸地萬分之三）貢獻了全球入海泥砂的可觀比例，是研究 source-to-sink 動力學的關鍵案例。把入門篇的「風化→侵蝕→搬運→堆積」放進這個收支框架，地形學就從描述地貌，升級為追蹤一顆礦物從中央山脈頂峰到台灣海峽海床的全程質量旅程——這也是當代地表動力學(surface dynamics)的研究前沿。