Presentation is loading. Please wait.

Presentation is loading. Please wait.

天体物理学 I : 授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。

Similar presentations


Presentation on theme: "天体物理学 I : 授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。"— Presentation transcript:

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

2 K 輻射の方程式 (II) 今回の内容 (K.1) 線形大気 ソースファンクション S(τ)が S(τ) =aτ+b の時の大気を調べます。
K 輻射の方程式 (II) 今回の内容 (K.1) 線形大気   ソースファンクション S(τ)が S(τ) =aτ+b の時の大気を調べます。 (K.2) モーメント方程式   各点での輻射方程式に cosnθ をかけて全立体角で積分した式とは?    (K.3) エディントン近似 (Eddington approximation)   無限に開いたモーメント方程式に K=J/3 という関係を仮定して強引に   閉じさせてしまう近似です。 (K.4) Rossland mean opacity κR   ただ全波長で平均した吸収係数はあまり使えないのです。役に立つのは調   和平均です。それはなぜ。 (K.5) 恒星大気のエディントンモデル   線形大気の方程式にエディントン近似を組み合わせると、色々な星のスペク   トルを作り出す事が出来ます。 (K.6) 温室効果    地球の大気を星の大気として扱うと、温暖化に何が効いているか判ります。 授業の内容は下のHPに掲載されます。 A: 水素原子

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

4 Iλ (μ=1,τλ=0) τλ=0 Y Z τλ X Iλ (μ=1,τλ)
下の図を見て下さい。Y,Z軸は星の大気の表面に乗っています。そして、X軸は星の内側(下側)方向が正の方向としています。X=0が表面です。  水平大気では密度、温度などの物理量 A(X,Y,Z) が高さ ( というより深さ X ) のみの関数 A(X)であると仮定します。  (1) マイナスX軸方向(上向き)の輻射強度 I に対する方程式は、               dI/dx= k I-ε Iλ (μ=1,τλ=0) τλ=0 τλ Iλ (μ=1,τλ) K: 輻射の方程式

5 μdI(μ,x)/dx=I(μ,x)κ(x)-ε(x)
 (2) 次に、X軸に角度 θ を成す直線 (μ=cosθ) に沿って、輻射方程式を考え     ましょう。直線に沿っての長さを t とすると、(1)と同じく             dI/dt= κI-ε      となります。図を見ると分かるようにdt と X軸に沿っての深さdXとの関係     は、dt=dx/μ なので、書き直して μdI(μ,x)/dx=I(μ,x)κ(x)-ε(x) Iλ (μ,τλ=0)  垂直に進む (1) が上の式の μ=1 の場合に相当することは言うまでもありません。 τλ=0 以降、Iλ 、 κλ等を I、κと省略します。 τλ θ Iλ (μ,τλ) K: 輻射の方程式

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

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

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

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

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

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

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

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

14 2次モーメント Iλ (μ,τλ=0) τλ=0 (表面) 斜め方向の輻射方程式 t X Iλ (μ,τλ) τ 1次モーメント
 M1(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次モーメント  M2(x,λ)=(1/4π)∫ (cosθ)2I(cosθ, x,λ) dΩ   = (1/2)∫μ2 I(μ, x,λ)dμ   =K (x,λ) Iλ (μ,τλ=0) τλ=0  (表面) 光圧力 P(ν) = (4π/c) K(ν) 斜め方向の輻射方程式 θ X軸に沿って光学的深さτを定めます。μ方向の光線に沿っては、以前にやったように Iλ (μ,τλ) τ L: エディントン近似

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

16 K.3. エディントン近似 (Eddington approximation)
   μdI/dτ=I-S  をΩで積分してモーメント方程式に変えます。 × ∫dΩ/4π   : × ∫μdΩ/4π  :  I(X,μ)に対する方程式では変数が2種類(波長まで入れると3種類)あって大変でしたが、J(X),H(X),K(X)なら解けそうな気がしてきます。 この系列はμ2  μ3  と上げても閉じません。 しかし、上の式を見るとJ(X),H(X),K(X)で3つの変数を求めなければいけないのに、式の数は2つです。このように、      式の数<変数の数 のような場合、もう一つ式がないと解が求まらないのです。 モーメント方程式をどこかでむりやり閉じないと解けません。     エディントン近似    L: エディントン近似

17 J=Io, K=(1/2)∫1-1 Io μ2 dμ= Io / 3 = J / 3
エディントン近似が正確に成り立つ例 (i) 完全等方輻射 I(Ω)=Ioの場合 J=Io, K=(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‐ 4H L: エディントン近似

18 エディントン近似モデル(iii) τ=0 Io 4Ho 4Ho 4Ho L: エディントン近似

19 K.4.Rossland mean opacity κR
今やっていること モーメントやらエディントン近似やらいろいろ出てきて、何のためにこんな事を、 と思い始めてきた頃でしょう。今いる場所を整理しましょう。   出発点 元々、出発点の式は輻射の方程式です。    光学的深さを、 dτλ=κλdX で定義すると、下式が解きたい式なのです。  角度θを消す しかし、これでは難し過ぎるので、まずモメントを取って角度θを消しました。 こうして J(X,λ)、 H(X,λ),  K(X,λ) のモーメント方程式になりました。 L: エディントン近似

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

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

22 で定義される平均吸収係数、ロスランド平均吸収κR を導入します。 κRを使うと、
( ここに、B(T)=∫Bλ)(T)dλ  ちょっと変ですが、方程式系の形を整えるため、B(T)から K(x) に戻し、 改めて整理すると、モーメント方程式の一つ を積分し、 L: エディントン近似

23 この式にロスランド平均吸収 (Rossland mean opacity) を導入して、
( ここに、B(T)=∫Bλ(T)dλ ) ちょっと変ですが、方程式系の形を整えるため、B(T)から K(x) に戻し、 というわけです。 以上で平らな大気に対するモーメント方程式は、エディントン近似と局所熱平衡の仮定の下で、平均輻射強度 J(x), 総フラックスH(x), 2次モーメント、K(x)に関する次の3つの方程式にまとめられました。 L: エディントン近似

24 κRに効くのは、κiが小さい所とΔBiが大きい所でκiが大きい所は効かない。
κRの意味 κRに効くのは、κiが小さい所とΔBiが大きい所でκiが大きい所は効かない。 F∝∑ΔBi /κi=ΔB /κR Fi ∝ΔBi(T) /κi Bi Bi+ΔBi L: エディントン近似

25 K.5. 恒星大気のエディントンモデル エディントン大気の解
K.4.で出した平らな層状構造に対する輻射の式を恒星大気の場合について解いてみましょう。 恒星大気では核反応は起きていないので E=0 です。従って の左式から、    H (x) = Ho ( = 一定) となります。 次に、 ロスランド平均光学的深さ τR を dτR = κR dx で定義すると、 第2式は、 となるので、K(τR ) = Ho・ τR +C と簡単に解けます。Cは積分定数です。 すると、第3式から、J (τR ) = 3Ho・ τR +3C です。 さらに、局所熱平衡の仮定、J(x)=S(x)=B[T(x)]から、 S( τR ) = 3Ho・ τR +3C、 B[T(x)] = 3Ho・ τR +3C です。 L: エディントン近似

26 F=π[3Ho・ (2/3) +3C] = π(2Ho+3C) 一方、フラックスFは一次モーメントHとの関係は、 F=4πHo
積分定数 C の決定 ここで、「K.1.線形モデル」最後の問題を思い出して下さい。S(τ)=aτ+ b の時、表面からの放射フラックス F は、F=πS(τ=2/3) でしたから、            F=π[3Ho・ (2/3) +3C] = π(2Ho+3C) 一方、フラックスFは一次モーメントHとの関係は、            F=4πHo  です。この二つから C=(2/3)Ho と決まります。 温度分布 「G.3.黒体輻射のその他の性質」で全輻射強度 B(T)=(σ/π)T4 と学びました。前ページの B(T) の式と合わせると、温度 T と光学的深さ τR の関係が判ります。             (σ/π)T4 == 3Ho・ τR +3C  です。  L: エディントン近似

27 一方、モーメントの定義のところでやったように、フラックスFは一次モーメントHとF=4πH という関係があります。
温度分布 「G.3.黒体輻射のその他の性質」で全輻射強度 B(T)=(σ/π)T4 と学びました。前ページの B(T) の式と合わせると、温度 T と光学的深さ τR の関係が判ります。             (σ/π)T4 == 3Ho・ τR +3C  積分定数 C の決定 ここで、「K.1.線形モデル」最後の問題を思い出して下さい。S(τ)=aτ+ b の時、表面からの放射フラックス F は、F=πS(τ=2/3)= でした。 一方、モーメントの定義のところでやったように、フラックスFは一次モーメントHとF=4πH という関係があります。 K(τR ) = Ho・ τR +C と簡単に解けます。Cは積分定数です。 すると、第3式から、J (τR ) = 3Ho・ τR +3C です。 さらに、局所熱平衡の仮定、J(x)=S(x)=B[T(x)]から、 S( τR ) = 3Ho・ τR +3C、 B[T(x)] = 3Ho・ τR +3C です。 L: エディントン近似

28 K(τ) = τHo+ (2/3)Ho =Ho (τ+ 2/3)
ここまでの結果は、エディントン近似モデルの      (iii)  I(τ)= I+ (τ)  μ>0 = I‐ (τ)  μ<0 でも考えられるます。  H(τ)=Ho=一定=(I+ – I‐)/4   K(τ)=τHo+ C=(I+ + I‐)/6 を解いて、    I+ (τ) =2H(τ)+3K(τ)=2 Ho +3(τHo+ C)    I‐(τ)=3(τHo+ C)- 2 Ho 仮定 : 表面τ=0で、I=Io (μ>0) =0 (μ<0) とすると、C=(2/3)Ho , Io=4Ho    H(τ)=Ho=一定    K(τ) = τHo+ (2/3)Ho =Ho (τ+ 2/3) で、前ページと同じになります。 L: エディントン近似

29 有効温度 Te エディントンモデルに入るパラメターはHoだけです。 パラメターHoを温度で表現する為、F= 4πHo =σTe 4 で
有効温度 Te エディントンモデルに入るパラメターはHoだけです。 パラメターHoを温度で表現する為、F= 4πHo =σTe 4 で   有効温度 Te  を導入します。すると、     Ho=σTe 4/4π    J(τ)=S(τ)=B(τ)=(3σTe 4/4π) (τ+2/3)=(σ/π)T4 (τ) 表面(τ=0)温度 To はTeよりやや低く、      To4 = (1/2)Te4、 (To=0.84Te) また、 T(τ=2/3)=Te  ここにも、τ=2/3 が現れています。 L: エディントン近似

30 J,H,Kのτによる変化 温度Tのτによる変化 4H 0 2/3 1 2 3 τ 0 2/3 1 2 3 τ J 3H K 2H H H
/ τ 1.5 / τ J 3H K Te To 2H H H 表面 L: エディントン近似

31 K.6.地球大気の温室効果 F(λ) 太陽 地球 λ 可視 赤外 F σTg4 地表 地球表面の温度は基本的には、
   太陽輻射による熱流入(主に可視域)=地表からの熱放射(主に赤外域) で決まります。 F(λ) 太陽         地球 λ 可視 赤外 F σTg4  地表 この時の熱平衡の式は、地表温度=Tgとおくと、F(太陽)=σTg4 です。 L: エディントン近似

32 F(太陽)=Fo 地球の表面が太陽から浴びる平均のフラックスはいくつになるでしょう?
To=太陽有効温度=5780K、 Ro=太陽半径、  D=1AU=215Ro とします。すると、   F=σTo4(Ro/D)2 : 地球軌道での太陽フラックス 地球の吸収断面積はΣ=π・Re2 なので、地球の吸収輻射 P =Σ・F これを地球の全表面積 S=4・π・Re2 で平均すると、 地表での平均太陽フラックス Fo=P/S=(1/4)F= (1/4)σTo4(Ro/D)2   To 断面積=Σ 表面積= S L: エディントン近似

33 地表の反射率を可視光Bとします。すると地表への熱流入は(1-B)・Fo です。 前ページに与えられた数値を入れると、
地表の反射率を可視光Bとします。すると地表への熱流入は(1-B)・Fo です。 前ページに与えられた数値を入れると、   σTg4 =(1-B)Fo=(1-B) (1/4)σTo4(Ro/D)2     Tg =(1-B) 1/4 (Ro/2D)1/2 To となります。これが、真空中にむき出しにさらされた場合の地球表面温度です。 上の第1式左辺に(1-B)が掛かっていないことに注意して下さい。これは地表からの放射は赤外域で起きていて、そこでの反射率をゼロと仮定していることを意味します。 実際に数値を入れて計算すると下の表のようになります。      B    0     0.3     0.5     0.7    0.9  Tg   279    255    234     206    156 氷河期 この表を見ると、反射率Bが上がると地表温度Tgが低下することがはっきりと判ります。これは、氷河期のメカニズムの基本です。 氷が増えると、反射率が増加し、すると温度が下がってさらに氷が増えというサイクルが回り出すと地球が氷河に覆われるようになることが上の表には含まれているのです。  L: エディントン近似

34 温室効果 (1)単層モデル ここまで地球大気の効果を無視してきました。大気にはどんな効果があるのでしょう?そこで、次のような仮定をおきます。  1)大気は一様な温度Ta を持つ。    2)太陽光は可視、地上からは赤外のみ放射  3)大気は可視で透明、赤外は不透明で黒体     4)可視太陽光の地表反射率=B 定常状態で、雲と太陽の間、雲と地球の間の輻射のやり取りを考えると、下の図のようになります。図中のσ=ステファン・ボルツマン定数です。 太陽 Fo  Fo=σTa4 +BFo   BFo  σTa4 大気 Ta σTa4 σTg4     Fo+ σTa4 = σTg4 +BFo   大地 Tg L: エディントン近似

35 1)前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 地表の間を流れる輻射のつり合いの二つの式を立てて下さい。
 1)前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と    地表の間を流れる輻射のつり合いの二つの式を立てて下さい。  2) 太陽のフラックスFoと地表反射率Aとが与えられたとして、    Ta と Tg をFoとAを使って求めて下さい。    L: エディントン近似

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

37 大気なしでは反射ゼロでやっと0°Cを上まわる温度だったのが、大気ありになると反射率50%以下なら0°Cを上まわることが判ります。
これまでの結果をまとめると、 大気なし      Tg= (1-B) 1/4 (Ro/2D)1/2 To 大気あり      Ta= (1-B) 1/4 (Ro/2D)1/2 To Tg=21/4 (1-B) 1/4 (Ro/2D)1/2 To 大気なし   B    0     0.3     0.5     0.7    0.9   Tg   279    255    234     206    156 大気あり   Ta   279    255    234     206    156   Tg   332    303    278     245    186 大気なしでは反射ゼロでやっと0°Cを上まわる温度だったのが、大気ありになると反射率50%以下なら0°Cを上まわることが判ります。 L: エディントン近似

38

39

40

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

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

43 大気上端 To TG 地表 L: エディントン近似

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

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

46 CO2 H2O CO2 H2O L: エディントン近似

47

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

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

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

51

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

53

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

55

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

57 リム・ダークニング ( limb darkening )と表面輻射強度
θ α0 I(θ) α 天体表面で輻射強度が鉛直方向からの角度θにより、 I(cosθ)で表されるとする。I(cosθ)は星の表面輝度分布F(α)に反映される。 上の図で なので、 逆に、F(α)が求まったら、 ところで、表面輝度分布I(cosθ)は源泉関数S(τ)と関係している。 仮に、S(τ)=a+b・τ+c・τ2+...と展開されたとすると、 L: エディントン近似

58 ここでも逆に I(cosθ) からS(τ)を以下のように求められるます。 I(cosθ)=A+B・cosθ+C・cos2θ+... なら、
S(τ)=A+B・τ+(C/2)・τ2+(D/6)・τ3+... τ S(τ) 恒星大気内でLTEが成立していると、源泉関数S(τν)=B(T, ν) から、τνの深さでの温度が決まります。 結局、星の表面の輝度分布がある波長(周波数)で決まると、大気内の温度変化がその波長での光学的深さの関数T(τν)として求まることが判ります。 L: エディントン近似

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

60 Fλ=∫μIλ(μ,τ=0)dΩ= 2π∫10μ( aλ+ bλμ)dμ=2π(aλ/2 +bλ/3)
では、そのように方向により大きさが異なる輻射強度全体のフラックス F はどうなるでしょうか?Fの計算式は、  F=∫ I(θ) cos θ dΩ=∫ I(θ) cosθ sinθdφdθ = 2π∫01 I(μ)μdμ   Fλ=∫μIλ(μ,τ=0)dΩ= 2π∫10μ( aλ+ bλμ)dμ=2π(aλ/2 +bλ/3)       Source Function  Sλ (τ)=aλ+bλτ だったから、 Fλ=π[aλ+(2/3)bλ]=πSλ(τ=2/3) です。  温度Tの黒体表面からのフラックスがπBλ(T),ここにBλ(T)は輻射強度、  だったことを考えると、線形大気では、τλ=2/3の深さの所を見て  いると言えます。  I(τ=0)  a  0  τλ=0  1/3  τλ=μ=cosθ S(τ=2/3)  2/3  1  τλ=1 S  a+b  a+bμ K: 輻射の方程式

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

62 1)前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と 地表の間を流れる輻射のつり合いの二つの式を立てて下さい。
 1)前ページの図を参考にして、太陽と雲の間を流れる輻射のつり合い、雲と    地表の間を流れる輻射のつり合いの二つの式を立てて下さい。  2) 太陽のフラックスFoと地表反射率Bとが与えられたとして、    Ta と Tg をFoとBを使って求めて下さい。    1)  Fo=σTa4 +B・Fo    Fo+ σTa4 = σTg4 +BFo 2) 2式を辺々足して、  2Fo+ σTa4 = σTa4 +2・B・Fo+ σTg4        2式を辺々引いて、  σTa4 = ーσTa4 +σTg4    したがって、       σTg4 = 2・(1-B)Fo                   σTa4 =(1/2)σTg4 = (1-B)Fo L: エディントン近似


Download ppt "天体物理学 I : 授業の内容 天文学は天体からの光を研究する学問です。 そこでこの授業では、「光」をどう扱うかの基礎を学びます。"

Similar presentations


Ads by Google