天体物理学 Ｉ ： 授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。 天体物理学 Ｉ　：　授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。 授業計画は、 Ａ．水素原子　　Ｂ．エネルギー準位　　Ｃ．熱平衡　　Ｄ．線吸収　　Ｅ．連続吸収　　 Ｆ．光のインテンシティ　Ｇ．黒体輻射　Ｈ．等級　　Ｉ．色等級図　　 Ｊ．光の伝達式　Ｉ　　Ｋ．光の伝達式　ＩＩ　Ｌ．星のスペクトル という順で進めます。 最後まで行くと、星のスペクトルがどんな仕組みで決まっているかが判る、 というのが目標です。 ＡからＥまでは光の吸収に関係する物理の話です。Ｆでは光の強さをきちん と定義します。ＧからＩは光の強さを天文学でどう使うかを示します。ＪからＬは 光がガス中を伝わる様子を式に表わし、その式を解いて星のスペクトルを導き ます。それでは、始めましょう。 Ａ：　水素原子

Ｋ 輻射の方程式 （ＩＩ） 今回の内容 （Ｋ．１） 線形大気 ソースファンクション Ｓ（τ）が Ｓ（τ） ＝aτ+b の時の大気を調べます。 Ｋ　輻射の方程式　（ＩＩ） 今回の内容 （Ｋ．１）　線形大気 　　ソースファンクション Ｓ（τ）が　Ｓ（τ） ＝aτ+b　の時の大気を調べます。 （Ｋ．２）　モーメント方程式 　　各点での輻射方程式に cosｎθ をかけて全立体角で積分した式とは？　　　 （Ｋ．３）　エディントン近似 （Eddington approximation） 　　無限に開いたモーメント方程式に　Ｋ＝Ｊ／３　という関係を仮定して強引に 　　閉じさせてしまう近似です。 （Ｋ．４）　Rossland mean opacity κR 　　ただ全波長で平均した吸収係数はあまり使えないのです。役に立つのは調 　　和平均です。それはなぜ。 （Ｋ．５）　恒星大気のエディントンモデル 　　線形大気の方程式にエディントン近似を組み合わせると、色々な星のスペク 　　トルを作り出す事が出来ます。 （Ｋ．６）　温室効果　 　　地球の大気を星の大気として扱うと、温暖化に何が効いているか判ります。 授業の内容は下のHPに掲載されます。 http://www.ioa.s.u-tokyo.ac.jp/kisohp/STAFF/nakada/intro-j.html Ａ：　水素原子

Ｋ．１．線形大気 半径Ｒの星を考えます。その星の厚みＨの大気を考えます。左の図のように星が球形なので、大気の構造も全体としては球殻状です。しかし、Ｈ/Ｒ＜＜１なので、右の図のように大気構造の曲率は無視できます。 大気構造の曲率を無視できる時、大気は水平な層構造と考えて構いません。　 Ｈ Ｒ そのように、平らな層が重なっている大気を水平大気と呼ぶことにしましょう。 Ｋ：　輻射の方程式

Iλ (μ=1,τλ=0) τλ＝０ Ｙ Ｚ τλ Ｘ Iλ (μ=1,τλ) 下の図を見て下さい。Ｙ，Ｚ軸は星の大気の表面に乗っています。そして、Ｘ軸は星の内側（下側）方向が正の方向としています。Ｘ＝０が表面です。　 水平大気では密度、温度などの物理量 Ａ（Ｘ，Ｙ，Ｚ） が高さ ( というより深さ Ｘ ) のみの関数 Ａ（Ｘ）であると仮定します。 　（１）　マイナスＸ軸方向（上向き）の輻射強度 Ｉ に対する方程式は、 　　　　　　　　　　　　　　ｄＩ/ｄｘ＝ k Ｉ－ε Iλ (μ=1,τλ=0) τλ＝０ Ｙ Ｚ τλ Ｘ Iλ (μ=1,τλ) Ｋ：　輻射の方程式

μdI(μ,x)/dx＝I(μ,x)κ(x)－ε(x) 　（２）　次に、Ｘ軸に角度 θ を成す直線 (μ＝cosθ) に沿って、輻射方程式を考え 　　　　ましょう。直線に沿っての長さを ｔ とすると、（１）と同じく 　　　　　　　　　　　　ｄＩ/ｄｔ＝ κＩ－ε 　　　 　となります。図を見ると分かるようにｄｔ と Ｘ軸に沿っての深さｄＸとの関係 　　　　は、ｄt＝ｄｘ／μ なので、書き直して μdI(μ,x)/dx＝I(μ,x)κ(x)－ε(x) Iλ (μ,τλ=0) 　垂直に進む　（１）　が上の式の μ＝１ の場合に相当することは言うまでもありません。 τλ＝０ 以降、Ｉλ 、 κλ等を Ｉ、κと省略します。 ｘ ｔ τλ Ｘ θ Iλ (μ,τλ) Ｋ：　輻射の方程式

形式解 μ>0 μ<0 μdI (μ,Ｘ)/dＸ＝I (μ,x) k(x)－ε(x) を書き直し、 　　　　　　　　　μdI (μ,Ｘ )/[ k(Ｘ) dＸ]＝I　(μ,Ｘ) －[ε(Ｘ)/k(Ｘ)　] 　ｄτ＝ k ｄＸ　とおきます。　Ｓ ＝ε(Ｘ)/ k(Ｘ)　は源泉関数 (Source Function)、 τは大気表面から真下に向かって測った光学的深さです。ですから、 光線に沿っての光学的深さではないことに注意して下さい。すると、 μdI （ μ, τ）/ dτ＝I （ μ, τ） －S （τ） 　 I （ μ<0 , τ=0） = 0　 　τ＝０ （表面） この方程式は、もしS （τ） が与えられていれば解くことができます。解は光線が 上向きの場合（ μ>0 ） と 下向きの場合（ μ<0 ）とで ちょっと違います。 それは、大気では表面、つまりτ＝０で下向きの輻射がゼロという境界条件があるからです。　 μ>0 τ μ<0 Ｋ：　輻射の方程式

= eτ/μ ∫∞τ S ( t ) e－t/μ dt/μ μ>0 μ>0　の時：　Ａ点での輻射強度は ｔ＝０ 　 I( τ, μ ) = ∫∞τ S(t) exp[－(t－τ)/μ] dt/μ = eτ/μ ∫∞τ S ( t ) e－t/μ dt/μ μ>0 Ａ この式で光学的深さτ, t は表面に対し垂直に測っています。右の図では点Ａの光学的深さはτ、点Ｂの光学的深さはｔです。ＢからＡへの経路は傾いているのでその間の光学的深さは（ｔ－τ）/μになります。 τ S(ｔ) τＢＡ＝ （ｔ－τ）/μ Ｂ ｔ Ｂ点ではそこでの源泉関数 (Source Function) Ｓ(ｔ) で光が生まれています。その光はＡ点に達するまでに exp[－(ｔ－τ)/μ ] に減光されて I (τ,μ ) に寄与するというのが積分式　S(t) exp[－(t－τ)/μ] dt/μ　の意味です。 Ｋ：　輻射の方程式

=∫τ0 S(t,λ) e－[ (τ－t) / (－μ) ] dt / (－μ) μ<0　の時： A点での輻射強度は I(τ＝０, μ＜０) =０ ｔ＝０ 　　I(τ,μ)= －∫τ0 S ( t, λ) exp[ (τ－t)/μ ] dt/μ = －eτ/μ ∫τ0 S( t, λ) e－t/μ dt /μ =∫τ0 S(t,λ) e－[ (τ－t) / (－μ) ] dt / (－μ) S(ｔ) ｔ B （ｔ－τ）/μ τ μ<0 の場合、 I (τ,μ ) に寄与するのは、τより浅い点 t 、つまりτ＞ｔ　です。ソースファンクションＳ(ｔ) で生まれた光が 　　　　 exp[－(ｔ－τ )/( - μ )] 　　　　= exp[ (ｔ－τ)/μ ] に減光されて 寄与することはμ>0　の場合と同じです。 ｔ＝０　の表面で外から星の表面に降り注ぐ輻射はゼロと仮定しているので、積分区間は　ｔ＝０　から　τ　までの間です。 Ａ Ｋ：　輻射の方程式

表面からの輻射強度 表面から角度θで出る輻射Ｉの解は下のように与えられます。 I(τ=0 、μ) = (1/μ)∫∞0 S(τ) exp(－τ/μ) dτ 上式を見ると、星の表面から放出される光は表面で生まれたものではなく、 星の内部、様々な深さで生まれた S(τ) の光が表面までの光学的深さτ/μ に応じて exp(－τ/μ) の減光を受けながら表面に辿りつき、宇宙空間に 放射されていくことが判ります。 I(τ=0 、μ) = S１ exp(－τ１ /μ) dτ/ μ + S2 exp(－τ2 /μ) dτ/ μ S3 exp(－τ3 /μ) dτ/ μ ．．． S１ exp(－τ１ /μ) S2 exp(－τ2 /μ) S3 exp(－τ3 /μ) τ１ S１ τ2 S2 τ3 /μ τ3 S3 Ｋ：　輻射の方程式

S(τλ)＝一定のスラブ表面からの輻射強度 簡単な例として、S(τλ)＝Ｓ＝一定（０ <τ<τo） の板を考えます。 板を垂直に見た時の光学的深さが τo です。 板の表面での Ｉ (τ=0 , μ) を計算すると、 I (τ=0 , μ) =　(1/μ)∫τo0 S exp ( -τ/μ) dτ　＝　S [ 1－exp (－τo /μ) ] 右図から放射方向θによる変化を見てとって下さい。 θ I(τ=0 , μ) τo S(τ） Ｋ：　輻射の方程式

線形解のフラックス 星の大気では内部に向かって高温になっていきます。ですから、ソースファンクション一定より、τの増加関数にする方が現実にマッチします。そこで、 　　　　　　　　　　S(τ)= a + bτ という形を仮定して、表面からの輻射強度を調べましょう。すると、 I ( τ=0 , μ>0 ) = (1/μ) ∫∞0 S (t) exp(‐t / μ) dt 　　　　　　＝(1/μ) ∫∞0 (ａ＋ｂt ) exp (－t / μ ) dt 　 ＝ (1/μ) [ a ∫∞0 exp (‐t /μ) dt + b ∫∞0 t exp(‐t /μ) dt ] ＝ a+ b μ = S ( τ = μ )　　 I ( τ=0 , μ<0 ) = 0 θ 　τ＝０ つまり、 I ( τ=0 , μ>0 ) は光線が来る方向にτ= １分伸ばした点のS ( τ = μ )です。　　 　τ＝μ 　τ＝１ Ｋ：　輻射の方程式

Ｆ＝∫ Ｉ(θ) cos θ ｄΩ＝∫ Ｉ(θ) cosθ sinθdφdθ = 2π∫01 Ｉ(μ)μｄμ です。 では、そのように方向により大きさが異なる輻射強度全体のフラックス Ｆ はどうなるでしょうか？　フラックスＦの計算式は　「Ｆ：　輻射強度」　でやったように、 　Ｆ＝∫ Ｉ(θ) cos θ ｄΩ＝∫ Ｉ(θ) cosθ sinθdφdθ = 2π∫01 Ｉ(μ)μｄμ です。 ｄΩ θ 星表面 Ｋ：　輻射の方程式

Ｋ．２． モーメント方程式 ソースファンクションＳ（τ）が判れば星の表面での輻射強度もフラックスも計算できることは、ここまでに示した通りです。ですから、問題はＳ（τ）を決めることです。これは、簡単には解けません。 そのため、様々な工夫が凝らされてきました。ここではモーメントを使う方法を紹介します。 普通、輻射は軸対称で、 I (x,θ,φ）＝ I ( x,θ） です。 μ＝cosθとし、Ｎ次モメント　MＮ　というものを以下のように定義します。 ｎ Ω θ ＭＮ（x, λ）＝(1/4π)∫(cosθ)Ｎ　I (θ, x) dΩ 　　　　　　　＝(1/4π) ∫∫ (cosθ)Ｎ　I (θ, x) (sinθ) dφdθ 　　　　　　　＝(1/2)∫μＮ I (μ, x)dμ ０次モーメント　　 Ｍ0（x,λ）＝ (1/4π)∫I (μ, x, λ) dΩ 　 ＝ (1/2)∫I (μ, x, λ) dμ = J （x,λ）＝ 平均輻射強度 (mean intensity) φ Ｌ：　エディントン近似

２次モーメント Iλ (μ,τλ=0) τλ＝０ (表面) 斜め方向の輻射方程式 ｔ Ｘ Iλ (μ,τλ) τ １次モーメント 　Ｍ1（x,λ）＝ (1/4π)∫cosθI (θ,x,λ) dΩ ＝ (1/2)∫μI (μ, x,λ) dμ = H（x,λ） エネルギーフラックス F（n, x ,λ) ＝∫ cosθ I (θ,x,λ) dΩ＝ 4πH ( x, λ) ２次モーメント 　Ｍ2（x,λ）=(1/4π)∫ (cosθ)2I(cosθ, x,λ) dΩ 　 = (1/2)∫μ2 I(μ, x,λ)dμ 　 =K （x,λ） Iλ (μ,τλ=0) τλ＝０ 　(表面) 光圧力 P(ν) = (4π/c) K(ν) 斜め方向の輻射方程式 ｔ θ Ｘ軸に沿って光学的深さτを定めます。μ方向の光線に沿っては、以前にやったように Ｘ Iλ (μ,τλ) τ Ｌ：　エディントン近似

dＨλ/dτλ= Jλ – Sλ dK λ/dτλ= Hλ μdI/dτ＝I－S　 ( i ) 両辺をdΩ/4πで積分する。(0次モーメント) ∫[μdI/dτ]dΩ/4π＝∫IdΩ/4π－ ∫ＳdΩ/4π = d[∫μIdΩ/4π]/dτ dＨλ/dτλ= Jλ – Sλ (ii) 両辺にμをかけてdΩ/4πで積分(１次モーメント) 　　d[∫μ2IdΩ/4π]/dτ ＝∫μIdΩ/4π－∫μSdΩ/4π 　　　∫1-1μdμ=0 に注意すると、　　　 dK λ/dτλ= Hλ Ｌ：　エディントン近似

Ｋ．３. エディントン近似 （Eddington approximation） 　　　μdI/dτ＝I－S　　をΩで積分してモーメント方程式に変えます。 × ∫dΩ/4π　　　： × ∫μdΩ/4π 　：　 Ｉ（Ｘ，μ）に対する方程式では変数が２種類（波長まで入れると３種類）あって大変でしたが、Ｊ（Ｘ），Ｈ（Ｘ），Ｋ（Ｘ）なら解けそうな気がしてきます。 この系列はμ2 　μ3 　と上げても閉じません。 しかし、上の式を見るとＪ（Ｘ），Ｈ（Ｘ），Ｋ（Ｘ）で３つの変数を求めなければいけないのに、式の数は２つです。このように、 　　　　　式の数＜変数の数 のような場合、もう一つ式がないと解が求まらないのです。 モーメント方程式をどこかでむりやり閉じないと解けません。  　　　エディントン近似　　　 Ｌ：　エディントン近似

J=Io, Ｋ＝(1/2)∫1-1 Io μ2 dμ= Io / 3 = J / 3 エディントン近似が正確に成り立つ例 (i) 完全等方輻射 I(Ω)=Ioの場合 J=Io, Ｋ＝(1/2)∫1-1 Io μ2　dμ=　Io / 3 = J / 3 Jλ=(1/2)∫1-1I dμ=(1/2)∫1-1(Io+I1μ)dμ=Io(λ) Hλ=(1/2)∫1-1Iμdμ=(1/2)∫1-1(Io+I1μ)μdμ＝(1/3)I1(λ) Kλ=(1/2)∫1-1Iμ2dμ=(1/2)∫1-1(Io+I1μ)μ2dμ=(1/3)Io(λ) θ (ii) 　I(τ,λ,μ)=Io(λ)+I1(λ)μ (iii) 　I(τ,λ,μ)= I+ (λ) μ>0 = I‐(λ) μ<0 J=(I+ + I‐)/2 H=[I+ /2 – I‐/2]/2=(I+ – I‐)/4 K=[I+ /3+ I‐ /3]/2=J/3 I+ I‐ ４Ｈ Ｌ：　エディントン近似

エディントン近似モデル(iii) τ＝０ Io 4Ho 4Ho 4Ho Ｌ：　エディントン近似

Ｋ．４．Rossland mean opacity κR 今やっていること モーメントやらエディントン近似やらいろいろ出てきて、何のためにこんな事を、 と思い始めてきた頃でしょう。今いる場所を整理しましょう。　　 出発点 元々、出発点の式は輻射の方程式です。　　　 光学的深さを、　ｄτλ＝κλｄＸ　で定義すると、下式が解きたい式なのです。　 角度θを消す しかし、これでは難し過ぎるので、まずモメントを取って角度θを消しました。 こうして　Ｊ（Ｘ,λ）、　Ｈ（Ｘ,λ）, 　Ｋ（Ｘ,λ）　のモーメント方程式になりました。 Ｌ：　エディントン近似

この系は変数Ｊ，Ｈ，Ｋの数３が方程式の数２より多いので困ります。そこで、 エディントン近似 　　　　エディントン近似 という仮定を立てて、変数の数と方程式の数を３つずつにしました。 次に減らしたいのが波長λです。前と同じ考え方だと、モーメント方程式を波長で積分することになります。今いるのはこの地点です。 角度θを消す 波長積分　Ｉ. モーメント方程式の一つ、 はκλを両辺にかけて としてからλで積分を掛け、 この式は、E(x) = 深さ ｘ での光吸収　∫κλ Ｊλｄλと放出∫ελｄλとの差、がフラックスの増加になることを示しています。 Ｌ：　エディントン近似

とすると、　 の中の　Ｋ(X, λ)　は、 エディントン近似により、　（1/3) J(X, λ）　に置き換えられます。局所熱平衡の近似を追加して、J(X, λ）＝ Ｂλ[Ｔ（Ｘ）]　（プランク関数） モーメント方程式の一つ、 波長積分　ＩＩ． λで積分すると、　 の辺は直接積分でき、 となります。 左辺はもう少し変形が必要です。 に着目すると、 となりますから、 Ｌ：　エディントン近似

で定義される平均吸収係数、ロスランド平均吸収κＲ を導入します。 κＲを使うと、 （　ここに、Ｂ（Ｔ）＝∫Ｂλ）（Ｔ）ｄλ　 ちょっと変ですが、方程式系の形を整えるため、Ｂ（Ｔ）から　Ｋ（ｘ）　に戻し、 改めて整理すると、モーメント方程式の一つ を積分し、 Ｌ：　エディントン近似

この式にロスランド平均吸収 （Rossland mean opacity） を導入して、 （　ここに、Ｂ（Ｔ）＝∫Ｂλ（Ｔ）ｄλ　） ちょっと変ですが、方程式系の形を整えるため、Ｂ（Ｔ）から　Ｋ（ｘ）　に戻し、 というわけです。 以上で平らな大気に対するモーメント方程式は、エディントン近似と局所熱平衡の仮定の下で、平均輻射強度　Ｊ（ｘ),　総フラックスＨ（ｘ），　２次モーメント、Ｋ（ｘ）に関する次の３つの方程式にまとめられました。 Ｌ：　エディントン近似

κＲに効くのは、κiが小さい所とΔBiが大きい所でκiが大きい所は効かない。 κＲの意味 κＲに効くのは、κiが小さい所とΔBiが大きい所でκiが大きい所は効かない。 Ｆ∝∑ΔBi /κi＝ΔB /κＲ Fi ∝ΔBi(T) /κi Bi Bi＋ΔBi Ｌ：　エディントン近似

Ｋ．５. 恒星大気のエディントンモデル エディントン大気の解 Ｋ．４．で出した平らな層状構造に対する輻射の式を恒星大気の場合について解いてみましょう。 恒星大気では核反応は起きていないので　Ｅ＝０　です。従って の左式から、　　　　H (x) = Ho ( = 一定) となります。 次に、　ロスランド平均光学的深さ τR を dτR = κR dx で定義すると、 第２式は、 となるので、Ｋ(τR ) = Ho・ τR ＋Ｃ　と簡単に解けます。Ｃは積分定数です。 すると、第３式から、Ｊ (τR ) = ３Ho・ τR ＋３Ｃ　です。 さらに、局所熱平衡の仮定、Ｊ（ｘ）＝Ｓ（ｘ）＝Ｂ［Ｔ（ｘ）］から、 Ｓ（ τR ) = ３Ho・ τR ＋３Ｃ、　Ｂ［Ｔ（ｘ）］ = ３Ho・ τR ＋３Ｃ　です。 Ｌ：　エディントン近似

Ｆ＝π[３Ho・ (2/3) ＋３Ｃ] = π（２Ho＋３Ｃ） 一方、フラックスＦは一次モーメントＨとの関係は、 Ｆ＝４πＨo 積分定数　Ｃ　の決定 ここで、「Ｋ．１．線形モデル」最後の問題を思い出して下さい。Ｓ（τ）＝aτ＋ ｂ の時、表面からの放射フラックス　Ｆ　は、Ｆ＝πＳ（τ＝２/３） でしたから、 　　　　　　　　　　　Ｆ＝π[３Ho・ (2/3) ＋３Ｃ]　= π（２Ho＋３Ｃ） 一方、フラックスＦは一次モーメントＨとの関係は、 　　　　　　　　　　　Ｆ＝４πＨo　 です。この二つから　Ｃ＝（２/３）Ｈo と決まります。 温度分布 「Ｇ．３．黒体輻射のその他の性質」で全輻射強度　Ｂ（Ｔ)＝（σ/π）Ｔ４　と学びました。前ページの　Ｂ（Ｔ)　の式と合わせると、温度　Ｔ　と光学的深さ　τR　の関係が判ります。 　　　　　　　　　　　　（σ/π）Ｔ４　＝= ３Ho・ τR ＋３Ｃ　 です。　 Ｌ：　エディントン近似

一方、モーメントの定義のところでやったように、フラックスＦは一次モーメントＨとＦ＝４πＨ という関係があります。 温度分布 「Ｇ．３．黒体輻射のその他の性質」で全輻射強度　Ｂ（Ｔ)＝（σ/π）Ｔ４　と学びました。前ページの　Ｂ（Ｔ)　の式と合わせると、温度　Ｔ　と光学的深さ　τR　の関係が判ります。 　　　　　　　　　　　　（σ/π）Ｔ４　＝= ３Ho・ τR ＋３Ｃ　 積分定数　Ｃ　の決定 ここで、「Ｋ．１．線形モデル」最後の問題を思い出して下さい。Ｓ（τ）＝aτ＋ ｂ の時、表面からの放射フラックス　Ｆ　は、Ｆ＝πＳ（τ＝２/３）＝　でした。 一方、モーメントの定義のところでやったように、フラックスＦは一次モーメントＨとＦ＝４πＨ　という関係があります。 Ｋ(τR ) = Ho・ τR ＋Ｃ　と簡単に解けます。Ｃは積分定数です。 すると、第３式から、Ｊ (τR ) = ３Ho・ τR ＋３Ｃ　です。 さらに、局所熱平衡の仮定、Ｊ（ｘ）＝Ｓ（ｘ）＝Ｂ［Ｔ（ｘ）］から、 Ｓ（ τR ) = ３Ho・ τR ＋３Ｃ、　Ｂ［Ｔ（ｘ）］ = ３Ho・ τR ＋３Ｃ　です。 Ｌ：　エディントン近似

K(τ) = τＨo+ (２/３)Ｈo ＝Ｈo (τ+ ２/３) ここまでの結果は、エディントン近似モデルの　 　　　　(iii) 　I(τ)= I+ (τ)　 μ>0 = I‐ (τ) 　μ<0 でも考えられるます。 　H(τ)=Ho＝一定＝(I+ – I‐)/4 　　K(τ)=τＨo+ C＝(I+ + I‐)/６　を解いて、 　　 I+ (τ) ＝２Ｈ(τ)＋３K(τ)=２ Ho ＋３(τＨo+ C) 　　 I‐(τ)＝３(τＨo+ C)－ ２ Ho 仮定 ： 表面τ＝０で、Ｉ＝Io (μ＞０) ＝0 (μ＜０) とすると、Ｃ＝(２/３)Ｈo , Io=4Ho 　　 H(τ)=Ho＝一定 　　 K(τ) = τＨo+ (２/３)Ｈo ＝Ｈo (τ+ ２/３) で、前ページと同じになります。 Ｌ：　エディントン近似

有効温度 Ｔｅ エディントンモデルに入るパラメターはＨｏだけです。 パラメターＨｏを温度で表現する為、Ｆ＝ 4πＨｏ ＝σＴe ４ で 有効温度　Ｔｅ エディントンモデルに入るパラメターはＨｏだけです。 パラメターＨｏを温度で表現する為、Ｆ＝ 4πＨｏ ＝σＴe ４　で 　　有効温度　Ｔｅ　 を導入します。すると、 　　　 Ｈｏ＝σＴe ４／４π 　　 J(τ)＝S(τ)=B(τ)=(３σＴe ４／４π) (τ+2/3)＝(σ/π)T4 (τ) 表面(τ＝０)温度 To はＴｅよりやや低く、　 　　　　To4 = (1/2)Te4、 (To=0.84Te) また、　T(τ=2/3)＝Ｔｅ　 ここにも、τ=2/3　が現れています。 Ｌ：　エディントン近似

Ｊ，Ｈ，Ｋのτによる変化 温度Ｔのτによる変化 ４H 0 2/3 1 2 3 τ 0 2/3 1 2 3 τ J ３H K ２H H H 0 2/3 1 2 3 τ １.５ 0 2/3 1 2 3 τ J ３H K Te To ２H H H 表面 Ｌ：　エディントン近似

Ｋ．６．地球大気の温室効果 Ｆ（λ） 太陽 地球 λ 可視 赤外 F σＴｇ４ 地表 地球表面の温度は基本的には、 　　　太陽輻射による熱流入（主に可視域）＝地表からの熱放射（主に赤外域） で決まります。 Ｆ（λ） 太陽　　　　　　　　　地球 λ 可視 赤外 F σＴｇ４　 地表 この時の熱平衡の式は、地表温度＝Ｔｇとおくと、Ｆ（太陽）＝σＴｇ４　です。 Ｌ：　エディントン近似

Ｆ（太陽）＝Ｆｏ 地球の表面が太陽から浴びる平均のフラックスはいくつになるでしょう？ To=太陽有効温度＝５７８０Ｋ、　Ｒｏ＝太陽半径、　　Ｄ＝１ＡＵ＝２１５Ｒｏ とします。すると、 　　F=σTo4(Ro/D)2 ：　地球軌道での太陽フラックス 地球の吸収断面積はΣ＝π・Ｒｅ２　なので、地球の吸収輻射　Ｐ　＝Σ・Ｆ これを地球の全表面積 Ｓ＝４・π・Ｒｅ２　で平均すると、 地表での平均太陽フラックス　Ｆｏ＝Ｐ/Ｓ＝（１/４）Ｆ＝ （１/４）σTo4(Ro/D)2 　 Ｆ Ｔｏ 断面積＝Σ 表面積＝ Ｓ Ｌ：　エディントン近似

地表の反射率を可視光Ｂとします。すると地表への熱流入は（１－Ｂ）・Ｆｏ です。 前ページに与えられた数値を入れると、 地表の反射率を可視光Ｂとします。すると地表への熱流入は（１－Ｂ）・Ｆｏ　です。 前ページに与えられた数値を入れると、 　　σＴｇ４　＝（１－Ｂ）Ｆｏ＝（１－Ｂ） （１/４）σTo4(Ro/D)2 　 　　Ｔｇ　＝（１－Ｂ） １/４ (Ro/２D)１/２ To となります。これが、真空中にむき出しにさらされた場合の地球表面温度です。 上の第１式左辺に（１－Ｂ)が掛かっていないことに注意して下さい。これは地表からの放射は赤外域で起きていて、そこでの反射率をゼロと仮定していることを意味します。 実際に数値を入れて計算すると下の表のようになります。　　　　 　Ｂ　　　　０　　　　 ０．３　　　 　０．５　　　　　０．７　　　　０．９ 　Ｔｇ　　　２７９　　　 ２５５　 　　２３４　　　 　２０６　　　　１５６ 氷河期 この表を見ると、反射率Ｂが上がると地表温度Ｔｇが低下することがはっきりと判ります。これは、氷河期のメカニズムの基本です。 氷が増えると、反射率が増加し、すると温度が下がってさらに氷が増えというサイクルが回り出すと地球が氷河に覆われるようになることが上の表には含まれているのです。　 Ｌ：　エディントン近似

温室効果　（１）単層モデル ここまで地球大気の効果を無視してきました。大気にはどんな効果があるのでしょう？そこで、次のような仮定をおきます。 　１）大気は一様な温度Ｔａ を持つ。　　　 ２）太陽光は可視、地上からは赤外のみ放射 　３）大気は可視で透明、赤外は不透明で黒体　　　 　４）可視太陽光の地表反射率＝Ｂ 定常状態で、雲と太陽の間、雲と地球の間の輻射のやり取りを考えると、下の図のようになります。図中のσ＝ステファン・ボルツマン定数です。 太陽 Fo　 Fo=σTa4 ＋ＢＦo　　 ＢFo　 σTa4 大気 Ta σTa4 σTg4 　　 　Ｆo＋ σTa4 ＝ σTg4 ＋ＢＦo　　 大地 Tg Ｌ：　エディントン近似

１）前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 地表の間を流れる輻射のつり合いの二つの式を立てて下さい。 　１）前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 　　　地表の間を流れる輻射のつり合いの二つの式を立てて下さい。 　２）　太陽のフラックスＦｏと地表反射率Ａとが与えられたとして、 　　　Ｔａ と Ｔｇ　をＦｏとＡを使って求めて下さい。　　　 Ｌ：　エディントン近似

下に示す二つは「大気なし」と「あり」とでの輻射の勘定です。 赤外 可視 地球外 Fo Ｂ・Fo (1- Ｂ)・Fo 地表　Tg 大気あり 赤外 可視 地球外 Ｂ・Fo (1- Ｂ)・Fo Fo 大気　Ta (1- Ｂ)・Fo 2(1- Ｂ)・Fo 地表　Tg Ｌ：　エディントン近似

大気なしでは反射ゼロでやっと０°Ｃを上まわる温度だったのが、大気ありになると反射率５０％以下なら０°Ｃを上まわることが判ります。 これまでの結果をまとめると、 大気なし 　　　　 Tg＝ （１－Ｂ） １/４ (Ro/２D)１/２ To 大気あり 　　　　 Ta＝ （１－Ｂ） １/４ (Ro/２D)１/２ To Tg＝21/4 （１－Ｂ） １/４ (Ro/２D)１/２ To 大気なし 　 Ｂ　　　　０　　　　 ０．３　　　 　０．５　　　　　０．７　　　　０．９ 　 Ｔｇ　　　２７９　　　 ２５５　 　　２３４　　　 　２０６　　　　１５６ 大気あり 　 Ｔa　　　２７９　　　 ２５５　 　　２３４　　　 　２０６　　　　１５６ 　 Ｔｇ　　　３３２　　　 ３０３　 　　２７８　　　 　２４５　　　　１８６ 大気なしでは反射ゼロでやっと０°Ｃを上まわる温度だったのが、大気ありになると反射率５０％以下なら０°Ｃを上まわることが判ります。 Ｌ：　エディントン近似

Ｔ２/3 ＝ Ta ＝To(Ro/D)1/2 (1-Ｂ) 1/4 温室効果　（２）エディントン大気モデル 大気は赤外では不透明で、温度勾配を持ちます。そこで、平面大気中を地表から大気表面へ向け F = (1-A)Fo の赤外フラックスが流れると考えます。 可視 赤外 大気表面 　　τ＝０ Fo=太陽の可視フラックス ＢFo=反射。 (1-Ｂ)Fo=地表で吸収。 － エディントン近似では、大気温度は大気表面から地表にかけて上昇し、ロスランド平均光学深さτＲ＝2/3 での温度Ｔ２/3が単層モデルの大気温度Taと等しい。 Ｔ２/3 ＝ Ta ＝To(Ro/D)1/2 (1-Ｂ) 1/4 τ＝2/3 地表面 τ＝τＧ Ｌ：　エディントン近似

τλG＜２/３ だと 地表（T=TG）が直接見えてしまう。 地球大気の特異性は、後に示すように大気の上端から地表までのロスランド平均光学的深さτRGが2/3より小さいことです。このため大気温度がTeまで達することは起きないのです。 そこで、通常のエディントン大気モデルを次のように変更します。 １：　地表までの大気温度分布はエディントン大気モデルを採用 　　　　　大気上端温度＝To　　　　　　　地表温度＝TG　　とすると、 ２：　地表までの光学的深さτλGは波長による。 　　　τλG＞２/３　ならば　Fλはτλ＝２/３、すなわちτR=（２/３）（τRG/τλG） 　　　の深さを見る。その深さでの温度 T は、 τλG＜２/３　だと　地表（T=TG）が直接見えてしまう。 Ｌ：　エディントン近似

大気上端 To TG 地表 Ｌ：　エディントン近似

大気から地表までの光学的厚みτλG Ｌ：　エディントン近似

大気吸収の近似式 大気吸収を取り扱い易くするために前々ページのグラフから、 吸収＝水蒸気連続吸収＋水蒸気バンド吸収+炭酸ガスバンド吸収 と考え、以下のように近似します。 次ページの上図は前と同じ大気の光学的深さの測定値で、下図は上の近似式をグラフにしたものです。バンドの形等細かいところでの違いはありますが、今後はこの近似式で話を進めます。 Ｌ：　エディントン近似

CO2 H2O CO2 H2O Ｌ：　エディントン近似

ロスランド平均光学的深さ τR ロスランド平均の定義は、 大気のκλは温度T の関数ですが、温度依存を無視すると、大気表面から地表までの平均光学深さ τR＝∫ κR ρｄｘ は上と同じ形、 で計算できます。 なので とおき、 となります。 Ｌ：　エディントン近似

簡単に前頁の積分を、下式のような和で置き換え、T=２８０Kて計算すると、 大気上端の温度として（ちょっと高いが）To=260Kを採用します。　すると、 Ｌ：　エディントン近似

吸収の強い波長帯では T=To （大気上端）の黒体輻射 吸収の弱い波長帯では T=TG(地表)の黒体輻射が卓越します。 Ｌ：　エディントン近似

炭酸ガス増加の影響 大気中CO２が地表温度を上げる効果を持つとよく言われます。そのメカニズムをエディントンモデルを使って調べてみましょう。 炭酸ガス増加により、１０、１５μｍ付近の吸収バンドが強くなります。他は水蒸気が卓越するので変化ありません。

新しいτでロスランド平均光学的深さを計算すると、 To=260K、τR＝０．３９３なので、 CO2が現在値ではTG=２８９．６Kだったから、２．４Kの温度上昇に相当します。 注意しておくと、炭酸ガスバンドのピーク波長付近ではもともとτが深いのでτRの増加に寄与しません。寄与が大きいのは炭酸ガス増加前はτが小さかったバンドの縁のあたりなのです。 つまり、炭酸ガスの増加により以前は透明だった吸収バンドの縁（ウイングと呼ばれます）が不透明になり、大気の平均光学的深さが大きくなることが炭酸ガスによる温暖効果の強化メカニズムです。 グラフを見ると分かるように、ロスランド平均光学的深さには水蒸気吸収も大きな役割を果たしています。したがって、現在はそれ程問題にされていませんが、地球の乾燥化も地球温暖化に影響があるはずです。

赤線スペクトルが炭酸ガスを倍に増やした時のものです。それを元と比べると、 F(λ)が　　Ａ：変わらない、Ｂ：増加、Ｃ：減少の３パターンあります。 その簡単な解釈を下にしめします。 A：　CO増加前に既にτ（λ）＞１だった波長はどちらにせよToを見るのでF（λ）不変 B:　CO増加前に既にτ（λ）＜＜１だったCOバンドがかぶらない波長は上昇したTGを 　　　見るのでF（λ）増加 C:　COバンドのウイング部分で弱い吸収を受け、2/3<τ（λ）<1、大気の下部を見て 　　　いた波長ではCO吸収が強くなると、大気の上部を見るようになり、温度が低下 　　　してF(λ)が減少する。

リム・ダークニング ( limb darkening )と表面輻射強度 θ α０ Ｉ（θ） α １ 天体表面で輻射強度が鉛直方向からの角度θにより、 Ｉ（cosθ）で表されるとする。Ｉ（cosθ）は星の表面輝度分布Ｆ（α）に反映される。 上の図で なので、 逆に、Ｆ（α）が求まったら、 ２ ところで、表面輝度分布Ｉ（cosθ）は源泉関数Ｓ（τ）と関係している。 仮に、Ｓ（τ）＝a+b・τ＋ｃ・τ２＋．．．と展開されたとすると、 Ｌ：　エディントン近似

ここでも逆に Ｉ(cosθ) からＳ（τ）を以下のように求められるます。 Ｉ(cosθ)＝Ａ＋Ｂ・cosθ＋Ｃ・cos２θ＋．．． なら、 Ｓ（τ）＝Ａ＋Ｂ・τ＋（Ｃ/２）・τ２＋(D/6)・τ３＋．．． τ S(τ) ３ 恒星大気内でＬＴＥが成立していると、源泉関数Ｓ（τν）＝Ｂ（Ｔ, ν） から、τνの深さでの温度が決まります。 結局、星の表面の輝度分布がある波長（周波数）で決まると、大気内の温度変化がその波長での光学的深さの関数Ｔ（τν）として求まることが判ります。 Ｌ：　エディントン近似

Ｔ（τν）を逆に表現すると、ある周波数（波長）での光学的深さτνが温度Ｔの関数τν（Ｔ）として ４ 前頁の変換をグラフで示すと下図の ようになります。 Ｔ（τν）を逆に表現すると、ある周波数（波長）での光学的深さτνが温度Ｔの関数τν（Ｔ）として 表されます。 このプロセスを二つの波長λ１と λ２で繰り返して、τ１（Ｔ）と τ２（Ｔ）を得ます。 I(α） I(θ） ５ これは、同じ温度Ｔの地点までの 光学的深さが波長によって異な るためです。したがって、各波 長での吸収係数をｋ１、ｋ２とする と、 １ ０ １ cosθ １ sinθ Ｌ：　エディントン近似

Ｆλ＝∫μIλ(μ,τ=0)ｄΩ＝ 2π∫１0μ( aλ+ bλμ)dμ=２π(aλ/2 +bλ/3) では、そのように方向により大きさが異なる輻射強度全体のフラックス Ｆ はどうなるでしょうか？Ｆの計算式は、 　Ｆ＝∫ Ｉ(θ) cos θ ｄΩ＝∫ Ｉ(θ) cosθ sinθdφdθ = 2π∫01 Ｉ(μ)μｄμ 　　Ｆλ＝∫μIλ(μ,τ=0)ｄΩ＝ 2π∫１0μ( aλ+ bλμ)dμ=２π(aλ/2 +bλ/3) 　 　　　　Source Function 　Sλ (τ)=aλ+bλτ　だったから、 Ｆλ＝π[aλ+(2/3)bλ]=πSλ(τ=2/3)　です。 　温度Ｔの黒体表面からのフラックスがπＢλ(T),ここにＢλ(Ｔ)は輻射強度、 　だったことを考えると、線形大気では、τλ＝２/３の深さの所を見て 　いると言えます。 　Ｉ（τ＝０） 　ａ 　０ 　τλ＝０ 　1/３ 　τλ＝μ＝ｃｏｓθ S(τ＝２/３) 　２/３ 　１ 　τλ＝１ S 　ａ＋ｂ 　ａ＋ｂμ Ｋ：　輻射の方程式

Ｆ＝∫ Ｉ(θ) cos θ ｄΩ＝∫ Ｉ(θ) cosθ sinθdφdθ = 2π∫01 Ｉ(μ)μｄμ です。 では、そのように方向により大きさが異なる輻射強度全体のフラックス Ｆ はどうなるでしょうか？　フラックスＦの計算式は　「Ｆ：　輻射強度」　でやったように、 　Ｆ＝∫ Ｉ(θ) cos θ ｄΩ＝∫ Ｉ(θ) cosθ sinθdφdθ = 2π∫01 Ｉ(μ)μｄμ です。 前ページで求めた、 　　Ｉ（μ）＝ a+ b μ 　を上式に代入します。 すると、 　　Ｆ＝ 2π∫01 Ｉ(μ)μｄμ 　　　＝ ２π ∫01 (a + bμ)μｄμ ＝２ π[(1/2)aμ2+(1/3)bμ3] 01 ＝２ π[(1/2)aμ2+(1/3)b] = π[a+b(2/3)] すなわち、 　　Ｆ＝πＳ[τ＝(2/3)] 源泉関数 S が線形、Ｓ＝a + b τ　の場合には、表面からのフラックスは表面が 輻射強度　Ｉ＝Ｓ[τ＝(2/3)]で等方的に放射している時と同じです。　 言い方を変えると、星からのフラックスは星の表面からの光ではなく、 τ＝(2/3) の所からの光を見ているとも考えられます。 ｄΩ θ 星表面 Ｋ：　輻射の方程式

１）前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 地表の間を流れる輻射のつり合いの二つの式を立てて下さい。 　１）前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 　　　地表の間を流れる輻射のつり合いの二つの式を立てて下さい。 　２）　太陽のフラックスＦｏと地表反射率Bとが与えられたとして、 　　　Ｔａ と Ｔｇ　をＦｏとBを使って求めて下さい。　　　 １）　 Fo=σTa4 ＋B・Ｆo 　　　Ｆo＋ σTa4 ＝ σTg4 ＋BＦo ２）　２式を辺々足して、　　２Ｆｏ＋ σTa4 ＝ σTa4 ＋２・B・Ｆo＋ σTg4 　　　 　　　２式を辺々引いて、　　σTa4 ＝ ーσTa4 ＋σTg4 　　　したがって、　　　　　　 σTg4 ＝ ２・（１－B）Ｆｏ 　　　　　　　　　　　　　　　　　 σTa4 ＝（１/２）σTg4 ＝ （１－B）Ｆｏ Ｌ：　エディントン近似