4.1 導論與 4.2 震源機制
內容整合了地震學第四章 4.1 與 4.2 的重點，並加入地震分析報告與震源機制球互動網頁。
4.1 導論 (Introduction)
闡述地震學從「波傳研究」轉向「震源物理」的演進
1. 地震學的核心定義與演進
- 早期焦點（傳播學）：側重於地震波的傳播（propagation）過程，將地震波作為探針，藉以推求地球內部的速度結構與分層。
- 現代焦點（震源物理）：轉向地震波的產生（generation）機制，利用輻射出的波動研究震源本身的物理過程（如破裂滑動、能量釋放）。
- 地震與斷層的強關聯性：地震幾乎毫無例外地發生在斷層（faults）上。斷層是岩石圈內部一側相對於另一側發生相對滑動的幾何不連續面。
2. 彈性回跳理論 (Elastic Rebound Theory)
由美國科學家 H. Reid 在深入調查 1906 年舊金山大地震後提出。核心概念是：板塊運動先在斷層兩側累積彈性應變，等到應力超過強度時，斷層突然破裂並滑移，釋放儲存的能量，形成地震波。
1
板塊運動蓄能遠離斷層的板塊邊緣持續以穩定速率相對運動。
2
摩擦力鎖定斷層面因強大摩擦力而處於 locked 狀態，阻止兩側岩體滑動。
3
彈性應變累積岩體逐漸彎曲變形，應變與能量持續累積。
4
突發破裂釋能剪應力超過抗剪強度極限後，斷層瞬間破裂並發生滑移。
觀念上的重大突破：在這套理論提出之前，常把地表斷層破裂看成地震的副作用；Reid 的理論確認了「斷層的突發破裂與滑移」本身就是地震的本質原因。
地震週期 (Seismic Cycle) 的四個階段
大地震通常發生在板塊邊界，反映了斷層面經歷數百至數千年能量累積與釋放的循環過程。
- 地震間期 (Interseismic stage)：占整個週期的絕大部分時間。遠離斷層處穩定運動，斷層被摩擦力鎖定，應變能緩慢且持續累積；部分斷層可能伴隨微弱的無震 creeping 蠕動。
- 震前期 (Preseismic stage)：主震破裂前夕，可能伴隨微小的前震（foreshocks）、微地動或地殼局部變形加速等潛在震前變形。
- 同震期 (Coseismic phase)：僅持續數秒至數十秒。斷層發生快速破裂滑動，向外輻射強烈地震波。
- 震後期 (Postseismic phase)：地震過後系統尚未恢復穩定，伴隨餘震（aftershocks）與暫態餘滑（afterslip），持續數月至數年，之後進入下一個週期。
三種基本斷層類型與滑移角對應
- 走向滑移斷層：左移 $\lambda = 0^\circ$；右移 $\lambda = 180^\circ$。
- 正斷層：$\lambda = 270^\circ$（或 $-90^\circ$）。
- 逆衝／逆斷層：$\lambda = 90^\circ$；若傾角低於 $45^\circ$，通常稱為 thrust。
真實地震多為斜滑斷層（oblique-slip），滑移角常介於上述特殊值之間。
2. 體波初動 (First Motions) 與節面 (Nodal Planes)
- 初動極性（Polarity）：P 波到達測站時的初始運動方向，是推求震源機制最簡單且經典的方法。
- 壓縮（Compression）：垂直分量地震圖上表現為初始波峰向上（Up / +）。
- 舒張/擴張（Dilatation）：垂直分量地震圖上表現為初始波谷向下（Down / -）。
- 四象限分佈與雙節面：剪切滑動會使 P 波初動在震源周圍呈四象限交替對稱分佈，兩個正交交界面稱為節面（nodal planes）。
多值性（Ambiguity）挑戰：滑動在斷層面與滑動在輔助面的初動場在波動力學上完全相同，因此單憑初動極性無法判定哪個才是真正的斷層面。
判定輔助手段：地質觀測、餘震分佈、大地測量變形場、波形破裂直接效應（directivity effects）。
3. 等效體力：雙力偶模型 (Double-Couple Source)
複雜的斷層面剪切滑動可簡化為無淨力矩的雙力偶等效體力源。這解釋了地震波輻射圖樣對斷層面與輔助面具備對稱性。
零軸（Null axis / B axis）：垂直於 $\mathbf{\hat{n}}$ 與 $\mathbf{\hat{d}}$ 的軸，也就是斷層面與輔助面的交線。
體波輻射振幅公式
在無限均勻彈性介質中，雙力偶源在遠場產生的位移公式如下。
$$u_r(t, \mathbf{x}) = \frac{1}{4\pi\rho\alpha^3 r} F(t - r/\alpha) \sin 2\theta \cos\phi$$
$$u_\theta = \frac{1}{4\pi\rho\beta^3 r} F(t - r/\beta) \cos 2\theta \cos\phi$$
$$u_\phi = \frac{1}{4\pi\rho\beta^3 r} F(t - r/\beta) (-\cos\theta \sin\phi)$$
- $\frac{1}{r}$ 是幾何擴散衰減項。
- $F(t)$ 是地震矩率函數（seismic moment rate function），即地震矩 $M(t)$ 的時間導數。
- $\sin 2\theta \cos\phi$ 描述 P 波四瓣狀輻射圖樣；節面上振幅為零。
- S 波沒有簡單零值節面，但其位移方向與 P 波節面有高度幾何對稱。
- 因分母含波速三次方，且 $\alpha \approx \sqrt{3}\beta$，S 波平均波幅通常比 P 波大約 5 倍。
標量地震矩（Scalar Seismic Moment, $M_0$）
$$M_0 = \mu \bar{D} S$$ 其中 $\mu$ 為剪切模量，$\bar{D}$ 為平均滑移量，$S$ 為破裂面積。
$$M_0 = \mu \bar{D} S$$ 其中 $\mu$ 為剪切模量，$\bar{D}$ 為平均滑移量，$S$ 為破裂面積。
4. 射出角 (Take-off Angle) 與下半球投影 (Lower Focal Hemisphere)
為了將全球測站觀測資料對應回震源球面上，必須考量射線傳播幾何。
- 射出角 $i$：射線在震源處出發時，與垂直向下方向的夾角。
- 司乃爾定律：$$p = \frac{\text{d}T}{\text{d}\Delta} = \frac{r_0 \sin i}{v_0}$$
- 測站距離 $\Delta$ 越遠，射出角 $i$ 越小，越接近垂直向下。
- 通常將全球測站投影於震源下方的下半球焦點投影球（Lower Focal Hemisphere）上。
5. 立體網投影與沙灘球表示法 (Beachballs)
使用立體投影網（stereonet）將三維下半球的壓縮與舒張分佈投影到二維平面，形成常見的震源機制解「沙灘球（beachball）」圖樣。
典型震源機制球特徵對照表
| 斷層類型 | 沙灘球投影特徵 | 物理與應力特徵 |
|---|---|---|
| 純走向滑移斷層 | 完美的「十字棋盤狀」四象限，兩個正交節面投影為通過球心的直線。 | 最大主應力與最小主應力皆位於水平面上，中性軸（Null/B 軸）垂直。 |
| 純逆衝／逆斷層 | 三象限「貓眼／眼睛狀」圖樣，中心大片區域為壓縮區，兩側為舒張區。 | 水平受強烈擠壓，上盤往上爬升，使震源上方受壓，T 軸垂直。 |
| 純張裂正斷層 | 三象限「貓眼狀」圖樣，中心大片區域為舒張區，兩側為壓縮區。 | 水平受拉張，上盤沿斷層面下滑，P 軸垂直（重力為主應力）。 |
| 斜滑／混合斷層 | 四象限邊界呈不對稱的彎曲弧線，黑白比例隨走向與傾向滑動分量而變。 | 結合水平切向滑移與垂直傾向滑動，應力軸傾斜。 |
6. 震源機制解與構造應力場的物理關係
- P 軸：最大主壓應力方向，平分舒張象限。
- T 軸：最小主壓應力（或最大張應力）方向，平分壓縮象限。
- B 軸 / Null 軸：垂直於 P 軸與 T 軸，也就是兩個正交節面的交線。
在立體網上的作圖求法：先找出斷層面與輔助面的節面大圓，再找極點，接著繪製通過兩個極點的垂直大圓，最後在大圓上兩極點正中間的位置得到 P 軸與 T 軸。
應力場垂直關係對照
- P 軸垂直 → 重力主導，形成正斷層。
- T 軸垂直 → 水平強烈擠壓，形成逆衝斷層。
- Null 軸垂直 → 以水平剪切為主導，形成走向滑移斷層。
7. 斷層幾何變換解析數學（必考重點公式）
在地震矩張量分析中，常需要進行走向、傾角、滑移角與法向量、滑移向量之間的解析互換。
（1）向量坐標表示
$$\mathbf{\hat{n}} = \begin{pmatrix} n_1 \\ n_2 \\ n_3 \end{pmatrix} = \begin{pmatrix} -\sin\delta \sin\phi_f \\ -\sin\delta \cos\phi_f \\ \cos\delta \end{pmatrix}$$
$$\mathbf{\hat{d}} = \begin{pmatrix} d_1 \\ d_2 \\ d_3 \end{pmatrix} = \begin{pmatrix} \cos\lambda \cos\phi_f + \sin\lambda \cos\delta \sin\phi_f \\ -\cos\lambda \sin\phi_f + \sin\lambda \cos\delta \cos\phi_f \\ \sin\lambda \sin\delta \end{pmatrix}$$
$$\mathbf{\hat{b}} = \mathbf{\hat{n}} \times \mathbf{\hat{d}} = \begin{pmatrix} \sin\lambda \cos\phi_f \cos\delta - \cos\lambda \sin\phi_f \cos\delta \\ \sin\lambda \sin\phi_f \cos\delta + \cos\lambda \cos\phi_f \cos\delta \\ \cos\lambda \sin\delta \end{pmatrix}$$
$$\mathbf{\hat{t}} = \frac{1}{\sqrt{2}}(\mathbf{\hat{n}} + \mathbf{\hat{d}}), \quad \mathbf{\hat{p}} = \frac{1}{\sqrt{2}}(\mathbf{\hat{n}} - \mathbf{\hat{d}})$$
（2）雙正交節面轉換公式
$$\cos\lambda_1 = \sin\delta_2 \sin(\phi_{f1} - \phi_{f2})$$
$$\cos\lambda_2 = \sin\delta_1 \sin(\phi_{f2} - \phi_{f1})$$
$$\cos\delta_2 = \sin\lambda_1 \sin\delta_1$$
$$\cos\delta_1 = \sin\lambda_2 \sin\delta_2$$
$$\tan\delta_1 \tan\delta_2 \cos(\phi_{f1} - \phi_{f2}) = -1$$
正交節面變換的計算流程
- 先用 $\cos\delta_2 = \sin\lambda_1 \sin\delta_1$ 計算第二節面的傾角。
- 再用 $\sin\lambda_2 = \frac{\cos\delta_1}{\sin\delta_2}$ 與 $\cos\lambda_2$ 確定滑移角象限。
- 利用 $\tan\delta_1 \tan\delta_2 \cos(\phi_{f1} - \phi_{f2}) = -1$ 解出走向差。
- 若 $\delta_2 > 90^\circ$，依幾何慣例做平移轉換。
外部互動網站
這兩個連結會直接放在頁面中；若對方網站允許嵌入，會直接顯示，否則也能從按鈕開啟。