地球的運動與四季（進階）：均時差、日射積分與被月球穩定的轉軸

為什麼太陽最早升起的日子，不是白晝最長的那一天？

如果你已經讀過入門篇，你大概很篤定：夏至(summer solstice)白晝最長、冬至(winter solstice)白晝最短，太陽的高度與晝夜長短全由黃赤交角 $23.5°$ 一手安排。這聽起來乾淨俐落。但真實的天空藏著一些不肯乖乖聽話的小細節。

打開任何一份精確的日出日落表你會發現：北半球日出最早的日子並不在夏至，而是落在六月初；日落最晚的日子也不在夏至，而是在六月底之後。同樣地，冬至前後，「最晚日出」與「最早日落」也各自錯開幾天。白晝確實在夏至最長、冬至最短，但日出與日落的極值卻偏離了——這個錯位看似微小，背後卻牽動著地球公轉速度的不均勻、軌道的橢圓、以及一條叫做「均時差(equation of time)」的隱形曲線。

入門篇把地球的運動講成一台精準的時鐘。這篇進階文章，我們要拆開這台時鐘的擒縱機構，看看為什麼「太陽時間」其實天天都在和我們的手錶吵架。

恆星日與太陽日：那消失的四分鐘

入門篇提過一句話：恆星日(sidereal day)約 $23$ 小時 $56$ 分，太陽日(solar day)才是 $24$ 小時。但這多出來的約 $4$ 分鐘從哪裡來？我們把它推導清楚。

地球一邊自轉，一邊公轉。所謂恆星日，是地球相對於遙遠恆星轉滿 $360°$ 所需的時間；所謂太陽日，是太陽連續兩次通過同一條子午線（例如連續兩個正午）的時間。問題在於：當地球自轉滿一整圈時，它同時也沿著軌道前進了一小段，太陽相對於地球的方向因此移動了一點點，地球必須多轉一點角度才能讓太陽再次回到正南方。

一年約 $365.25$ 天，地球公轉一圈 $360°$，所以平均每天公轉前進：

$$\Delta\theta = \frac{360°}{365.25} \approx 0.9856°/\text{天}$$

地球自轉角速度約 $360°/(23.934\ \text{h})$，要多補上這 $0.9856°$ 所需的時間：

$$\Delta t = \frac{0.9856°}{360°} \times 23.934\ \text{h} \times 60 \approx 3.93\ \text{分鐘}$$

這就是太陽日比恆星日長約 $4$ 分鐘的由來。換句話說，恆星每天提早約 $4$ 分鐘升起，一年累積下來剛好多繞天空一圈——這正是為什麼夏夜與冬夜看到的星座完全不同，也是為什麼天文觀測要用恆星時(sidereal time)而非民用時來對準目標。

值得注意：這裡我說的是「平均」太陽日。真實的太陽日長度其實每天都不一樣，會在 $24$ 小時上下擺盪達數十秒。要理解這個擺盪，得先認識太陽在天球上的兩種「不老實」。

均時差：太陽為什麼不守時

我們手錶上的時間，是一個假想的「平均太陽(mean sun)」——它以完全均勻的速度繞天球赤道運行，每天精準走 $24$ 小時。但天上真正的太陽（視太陽，apparent sun）跑得時快時慢。兩者的差距，就是均時差：

$$\text{EoT} = (\text{視太陽時}) - (\text{平均太陽時})$$

均時差在一年中於約 $-14$ 分鐘到 $+16$ 分鐘之間擺動，由兩個獨立效應疊加而成。

效應一：軌道偏心率(eccentricity)。 地球軌道是橢圓，依克卜勒第二定律(Kepler's second law)，地球在近日點(perihelion，一月初)附近跑得快、在遠日點(aphelion，七月初)附近跑得慢。公轉快的時候，太陽在天球黃道上「東移」得快，太陽日就被拉長；公轉慢時則縮短。這個效應的週期是一年，振幅約 $\pm 7.7$ 分鐘。

效應二：黃赤交角(obliquity)。 即使地球以等速沿黃道前進，由於黃道相對天球赤道傾斜 $23.5°$，太陽沿黃道的均勻移動「投影」到赤道上（我們量時間是沿赤道量的）並不均勻。在二分點附近，黃道與赤道斜交，赤道方向的分量較小；在二至點附近，黃道平行於赤道，投影分量反而較大。這個效應的週期是半年，振幅約 $\pm 9.9$ 分鐘。

把兩者相加，均時差呈現一條一年有兩個峰、兩個谷的不對稱曲線。一個近似式為：

$$\text{EoT} \approx 9.87\sin(2B) - 7.53\cos B - 1.5\sin B \quad(\text{分鐘})$$

其中 $B = \dfrac{360°}{365}(N - 81)$，$N$ 是一年中的第幾天。第一項來自黃赤交角，後兩項來自偏心率。

看一個例子：日出最早不在夏至

現在我們能回答開場的謎題了。某地的日出時刻，可以粗略寫成：

$$t_{\text{日出}} \approx 12{:}00 - \frac{H_0}{15°/\text{h}} - \text{EoT} + (\text{時區與經度修正})$$

其中 $H_0$ 是日出時太陽的時角(hour angle)，它只跟太陽赤緯有關，在夏至達到最大（白晝最長）。如果均時差是常數，日出最早就會發生在夏至。

但 EoT 在六月初到六月底之間正快速變化（從約 $+2$ 分鐘減到約 $-3$ 分鐘）。這個逐日漂移的 EoT，疊加到緩慢變化的 $H_0$ 上，使日出時刻的最小值提前到六月初、日落時刻的最大值延後到六月底。白晝長度（日落減日出）的極值仍乖乖落在夏至，但日出、日落各自的極值卻被均時差「拖」開了。

下次有人理直氣壯說「夏至太陽升得最早」，你可以微笑著請他去查日出表。

日射量的積分：一天到底曬進多少能量

入門篇用 $I = I_0\cos\theta$ 描述瞬時的能量密度。進階一層，我們要問：一整天累積的太陽能量（每日日射量，daily insolation）到底是多少？這才是決定季節氣候的真正物理量。

關鍵變數是太陽的赤緯(declination) $\delta$——太陽直射點的緯度，隨季節在 $\pm 23.5°$ 之間擺盪，近似為：

$$\delta \approx 23.5° \times \sin\!\left(\frac{360°}{365}(N - 81)\right)$$

對緯度 $\phi$ 的某地，正午太陽的天頂角 $\theta_z$ 滿足：

$$\cos\theta_z = \sin\phi\sin\delta + \cos\phi\cos\delta\cos H$$

其中 $H$ 是時角（正午時 $H = 0$，每小時增加 $15°$）。日出與日落時太陽位於地平線，$\theta_z = 90°$，由此解出半日長時角 $H_0$：

$$\cos H_0 = -\tan\phi\tan\delta$$

白晝長度就是 $2H_0/15°$ 小時。把瞬時日射率沿白晝對時角積分，得到大氣層頂的每日日射量：

$$Q_{\text{day}} = \frac{S_0}{\pi}\left(\frac{\bar d}{d}\right)^2 \left(H_0\sin\phi\sin\delta + \cos\phi\cos\delta\sin H_0\right)$$

其中 $S_0 \approx 1361\ \text{W/m}^2$ 是太陽常數，$(\bar d/d)^2$ 是日地距離的修正項（近日點時略大於 $1$）。$H_0$ 以弧度代入。

動手算一下：極區夏至的日射弔詭

這條公式藏著一個違反直覺的結果。我們比較夏至日赤道與北極的每日日射量。

夏至時 $\delta = +23.5°$。先看北極（$\phi = 90°$）：此時太陽整天不落（永晝），$H_0 = 180° = \pi$。代入：

$$Q_{\text{極}} = \frac{S_0}{\pi}\left(\frac{\bar d}{d}\right)^2 \left(\pi \cdot \sin 90° \sin 23.5° + \cos 90°\cdots\right) = S_0\left(\frac{\bar d}{d}\right)^2 \sin 23.5°$$

$$Q_{\text{極}} \approx 1361 \times 0.987 \times 0.399 \approx 536\ \text{W/m}^2 \ (\text{日平均})$$

再看赤道（$\phi = 0°$）：$H_0 = 90° = \pi/2$，白晝 $12$ 小時。

$$Q_{\text{赤}} = \frac{S_0}{\pi}\left(\frac{\bar d}{d}\right)^2\left(0 + \cos 23.5°\sin 90°\right) = \frac{S_0}{\pi}\left(\frac{\bar d}{d}\right)^2 \cos 23.5°$$

$$Q_{\text{赤}} \approx \frac{1361 \times 0.987 \times 0.917}{\pi} \approx 392\ \text{W/m}^2$$

結果令人吃驚：夏至這天，北極大氣層頂接收的每日日射量，竟然比赤道還多（約 $536$ vs $392\ \text{W/m}^2$）！原因是北極雖然太陽角度低（$\cos\theta$ 小），卻享有 $24$ 小時不間斷的日照，靠「曬得久」補回了「曬得斜」。

那為什麼北極不熱？因為這是大氣層頂的值。真實地表還要扣除：低角度入射經過更長的大氣路徑而被散射吸收、冰雪極高的反照率(albedo)把多數陽光反射回太空、以及把能量大量耗在融冰與蒸發上。日射量只是氣候方程式的「輸入」，不是「輸出」。

季節遲滯：為什麼最熱的不是夏至

既然夏至日射量最大，為什麼台灣最熱的月份是七、八月，而非六月底的夏至？最冷的也不是冬至，而是一、二月？這是季節遲滯(seasonal lag)，它讓我們把地表氣候看成一個有「熱慣性」的系統。

把地表能量收支寫成一個簡化的熱容方程：

$$C\frac{dT}{dt} = Q(t) - \varepsilon\sigma T^4$$

左邊是熱容 $C$ 乘以溫度變化率，右邊是吸收的日射 $Q(t)$ 減去向外輻射的能量。$Q(t)$ 在夏至達到峰值，但溫度 $T$ 並不會立刻跟上——因為要先把熱量「存」進海洋、土壤、大氣這些有熱容的儲庫裡。只要 $Q(t)$ 仍大於散熱，溫度就持續上升；要等到 $Q(t)$ 降到與散熱平衡，溫度才見頂。

數學上，若把 $Q(t)$ 近似為週期為一年的正弦驅動，這類一階系統的響應會落後驅動一個相位 $\varphi$，滿足：

$$\tan\varphi = \omega\tau$$

其中 $\omega = 2\pi/\text{年}$，$\tau = C/(4\varepsilon\sigma T^3)$ 是系統的熱鬆弛時間(thermal relaxation time)。熱容愈大，$\tau$ 愈長，遲滯愈久。這正是為什麼海洋性氣候（大水體、$C$ 大）的季節遲滯比內陸（$C$ 小）更明顯，也是為什麼一天當中最熱的時刻是午後兩三點而非正午——同一個物理機制，只是把週期從一年換成一天。

重點回顧

• 太陽日比恆星日長約 $4$ 分鐘，源於地球每天公轉前進約 $0.99°$，自轉必須多補這個角度才能讓太陽回到子午線。
• 均時差(EoT) 是視太陽與平均太陽的時間差，由軌道偏心率（週期一年）與黃赤交角（週期半年）兩效應疊加，全年在 $-14$ 至 $+16$ 分鐘間擺動。
• 因 EoT 逐日漂移，日出最早與日落最晚都不發生在夏至，但白晝長度極值仍在二至點。
• 每日日射量 $Q_{\text{day}}$ 由赤緯 $\delta$、緯度 $\phi$ 與半日長時角 $H_0$ 決定；夏至時北極大氣層頂日射量竟超過赤道，靠永晝補足低入射角。
• 季節遲滯來自地表熱容造成的相位落後（$\tan\varphi = \omega\tau$），所以最熱在夏至之後，最冷在冬至之後；同一機制也解釋一日中的午後高溫。

深入探討（研究所視角）

入門篇談到米蘭科維奇循環時，把黃赤交角的擺盪當成一個有固定週期（約 $41{,}000$ 年）、振幅在 $22.1°$–$24.5°$ 之間的規律振盪。但這個「規律」其實建立在一個我們很容易視為理所當然的前提上——月球的存在。一旦把月球拿走，地球的季節故事會變得遠比想像中狂野。

第一，自轉軸的混沌(chaotic obliquity)。 一顆行星的自轉軸進動，會被其他行星攝動所激發的軌道進動共振拉扯。法國天文學家 Laskar 等人在 1990 年代以長期數值積分顯示：若沒有月球，地球的黃赤交角將落入一片混沌帶(chaotic zone)，可在約 $0°$ 到超過 $60°$ 之間大幅、不可預測地漂移。試想黃赤交角擺到 $60°$——那時極區會在夏季被太陽近乎直射、冬季陷入漫長黑暗，回歸線與極圈的位置（入門篇的 $23.5°$ 與 $66.5°$）會劇烈遷移，整個氣候帶的劃分將被徹底改寫。

第二，月球的陀螺穩定作用。 為什麼有月球就不混沌？因為月球（加上太陽）對地球赤道隆起施加的重力矩，使地軸的進動速率被「鎖」在約 $26{,}000$ 年（每年約 $50''$）。這個進動頻率遠離了會引發災難性共振的軌道頻率，於是黃赤交角只在 $\pm 1.5°$ 的窄帶內小幅呼吸。換言之，月球是地球氣候穩定的隱形支柱：它讓四季的劇本年復一年大致照演，沒有突然改寫。這也成為討論系外行星宜居性(habitability)時的一個重要參數——一顆缺乏大衛星的類地行星，其轉軸可能劇烈翻滾，難以維持長期穩定的氣候。

第三，從幾何回到動力學：真近點角與克卜勒方程。 前文的均時差，其偏心率項本質上來自地球在軌道上「真實角位置」與「均勻時間」的差異。要精確計算這個差異，必須求解克卜勒方程(Kepler's equation)：

$$M = E - e\sin E$$

其中 $M$ 是平近點角(mean anomaly，隨時間均勻增加)，$E$ 是偏近點角(eccentric anomaly)，$e \approx 0.0167$ 是地球軌道偏心率。這是個超越方程，無法解析反解，需用牛頓法等數值迭代求出 $E$，再轉成真近點角 $\nu$：

$$\tan\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}}\,\tan\frac{E}{2}$$

$\nu$ 與 $M$ 之差（即中心差，equation of the center）正是均時差偏心率項的根源。當我們把偏心率的長期變化（米蘭科維奇的 $\sim 100{,}000$ 年週期）代入，均時差曲線本身也會在數萬年尺度上緩慢變形——今天日出最早不在夏至這件「天文常識」，在地質時間裡其實也不是常數。

第四，把這些串成一張因果網。 從一條 $4$ 分鐘的恆星日—太陽日差，到均時差的不對稱曲線，再到每日日射量的緯度積分、季節遲滯的相位落後，最後上溯到月球穩定的轉軸與克卜勒方程的數值求解——這些看似零散的進階主題，其實都是同一件事的不同切面：把「地球繞日」從一個圓形示意圖，升級成一個有橢圓、有傾角、有熱慣性、有長期混沌的真實動力系統。 對研究地球氣候、行星宜居性或古氣候(paleoclimate)代理紀錄的人而言，正是這些「不守時」「不對稱」「不穩定」的細節，承載了最豐富的科學訊息。下次抬頭看那顆斜掛南天、似乎準時的太陽時，別忘了它其實天天都在和你的手錶，悄悄地討價還價。