Otsuji's homepage

2006/12/16

NOAA10798の研究(9):Velocityの時間変化(I)

ここでは、arch filamentの上昇速度の時間変化を調べる。
その前に、5波長（Hα center、±0.5Å、±0.8Å）の観測点で、何km/sまでのDoppler shiftをcloud modelから求めることが可能かを調べた。

Cloud modelの適用限界

Source function、optical depth、Doppler widthを固定してDoppler shiftのみを0.0Åから2.0Å(約90km/s相当)まで変化させた場合のcloud modelのあてはめの良さを調べた。
以下の図において、横軸はcontrastを計算するために用いたDoppler shift量、縦軸はそのcontrastを用いてcloud modelから求めたDoppler shift量。
またfittingの様子を表した動画も載せる。点線が実際のContrast、実線はfitting結果。

初期値を固定した場合

fittingの初期値を[0.15,1.0,0.5,0.0]で固定した場合。

Source function : 0.15 optical depth : 1.0 Doppler width : 0.5 Fitting movie (GIF)	Source function : 0.25 optical depth : 1.0 Doppler width : 0.5 Fitting movie (GIF)
Source function : 0.15 optical depth : 0.5 Doppler width : 0.5 Fitting movie (GIF)	Source function : 0.15 optical depth : 1.0 Doppler width : 0.8 Fitting movie (GIF)

見る限り、大体±0.8Å以上ではfittingがうまくいかない。
ただし、Doppler widthが大きければ多少あてはめは良くなる。

初期値を変化させた場合

実際のパラメータの小数第2桁目に乱数を加えたものをFittingの初期値とした場合。
（周囲の状況からパラメータの初期値を推測することに相当）

Source function : 0.15 optical depth : 1.0 Doppler width : 0.5 Fitting movie (GIF)	Source function : 0.25 optical depth : 1.0 Doppler width : 0.5 Fitting movie (GIF)
Source function : 0.15 optical depth : 0.5 Doppler width : 0.5 Fitting movie (GIF)	Source function : 0.15 optical depth : 1.0 Doppler width : 0.8 Fitting movie (GIF)

先ほどよりは適用限界は広がったが、やはりshift量が大きくなるとあてはめが悪くなることは否めない。
ただ、正確なfitting初期値を推測することさえできれば、100km/s程度までの速度をcloud modelから導くことは可能であると考える。