地中の温度は一次元の熱伝導方程式を解いて計算する.
 |
(A.55) |
ここで は地中温度 (K), は地面物質の密度
(kgm-3), は地面物質の比熱 (Jkg-1K-1),
は熱伝導率 (Wm-1K-1) である. 地表面温度 は
で与えられる.
地表面での境界条件は,
 |
(A.56) |
である. ここで は地表面での太陽放射フラックス
(鉛直下向きを正),
は地表面での反射能 (アルベド), は地表面から射出される正味の赤外放射フラックス,
は地表面から大気へ輸送される顕熱フラックスである.
赤外放射, 顕熱フラックスは鉛直上向きを正にとっている.
計算領域下端における境界条件は断熱壁を置く.
アルベド, 地面物質の密度と比熱, 熱伝導率は
Kieffer et al. (1977) の
標準モデルにおけるパラメータを用いる. すなわち
= 0.25,
= 1650 kgm-3,
= 588 JK-1kg-1,
= 7.63×10-2
JK-1m-1sec-1 とする. このとき熱慣性
は 272
Wm-2sec1/2K-1, 日変化の周期 に対する表皮深さ
はおよそ 8.2 cm となる.
|