detect_inflation_preeruptive2コマンドマニュアル

(The documentation of detect_inflation_preeruptive2 command)

Last Update: 2023/06/12

◆機能・用途(Purpose)

detect_inflationコマンドの出力を用いて Maeda (2023)に基づき噴火直前の傾斜変動を検知する。
Detect a pre-eruptive tilt change using the output of detect_inflation command on the basis of Maeda (2023).

◆概要(Overview)

detect_inflationコマンドでは波形$s(t)$を用いて全てのタイムウインドウ$[t^{st},t^{en}]$における一方向への変動の指標$M(t^{st},t^{en})$を計算した。 $M(t^{st},t^{en})$の値はそのタイムウインドウが火山の山体膨張・収縮に伴う一方向への傾斜変動を表している可能性の相対的な大小を与えるものである。 $M(t^{st},t^{en})$の値が大きいほど傾斜変動を表す可能性が高いと考えられる。しかしdetect_inflationコマンドの出力はそこまでで止まっており、波形の中に傾斜変動が含まれるのか否かや含まれる場合のタイムウインドウを直接的・確定的に与えるような出力にはなっていない。
The detect_inflation command computes a measure $M(t^{st},t^{en})$ of a unidirectional change in every time window $[t^{st},t^{en}]$ of a waveform $s(t)$. The value of $M(t^{st},t^{en})$ indicates the relative possibility that each time window corresponds to a unidirectional tilt change caused by inflation or deflation of a volcano; the larger the $M(t^{st},t^{en})$, more likely the tilt change. However, the output of the detect_inflation command stops there, and does not consist of direct and definite information on whether the tilt change is present in the waveform and the time of the tilt change window if present.

そこで、detect_inflationコマンドに引き続いて実行する事後処理用のプログラムとしてこのコマンドを用意した。このコマンドでは噴火直前の傾斜変動にターゲットを絞り、 detect_inflationコマンドの出力を用いて傾斜変動の有無やタイムウインドウを具体的に与える。
This command was developed as a postprocessing program for the detect_inflation command. The target signals are limited to tilt changes that immediately precede eruptions. This command definitely determines the existence or absence of the pre-eruptive tilt change and the time window of the tilt change if present.

◆ソースコード(Source code)

$YMAEDA_OPENTOOL_DIR/event_detection/src/detect_inflation_preeruptive2.c

◆使用方法(Usage)

コマンドライン引数でパラメータを指定する。パラメータの一覧を下表に示す。
Specify parameters by command-line arguments. The table below shows a list of parameters.

●「-」から始まらない引数 (Arguments not beginning with “-”)

引数 Argument	与える値 Quantity to be given
第1引数 1st argument	入力ファイル名。全ての$[t^{st},t^{en}]$に対する$M(t^{st},t^{en})$の値を列挙したテキストファイルであり、 detect_inflationコマンドの出力ファイルまたはそれと同様の書式を持つファイルとする。 The input file name. This file must be a text file that consists of $M(t^{st},t^{en})$ for all $[t^{st},t^{en}]$, in the same format as that of the output file of detect_inflation command.
第2引数 2nd argument	出力ファイル名。テキストファイルであり、拡張子は何でも良い。 The name of the output file. This is a text file with an arbitrary extension.

●1つの「-」から始まる引数 (Arguments beginning with a single “-”)

このコマンドでは1つの「-」から始まる引数は存在しない。
This command does not have arguments beginning with a single “-”.

●「--パラメータ名=パラメータ値」の形式の引数 (Arguments of a form “--Parameter name=Parameter Value”)

「--パラメータ名=パラメータ値」の形式の引数は自由な順番で指定できる。「-」から始まらない引数の間に挿入しても良い。相反する指定がなされた場合には後の指定が優先される。デフォルト値を持つパラメータは省略できる。
Arguments of a form “--Parameter name=Parameter Value” can be placed in an arbitrary order. They can even be inserted between arguments not beginning with “-”. In case of conflicting options being specified, the latter option has a higher priority. Parameters that have default values can be omitted.

パラメータ名 Parameter name	意味 Meaning	可能なパラメータ値 Allowed parameter values	デフォルト値 Default value
t_eruption	噴火時刻$t^{eruption}$(s)。時刻の原点は入力ファイル第3列、第4列と一致させること。 Time (s) of the eruption $t^{eruption}$. The time origin must be consistent with that of the 3rd and 4th columns of the input file.	実数。入力ファイルに登場する時刻範囲の先頭を$t_{file}^{min}$、末尾を$t_{file}^{max}$として $t_{file}^{min}<t^{eruption}\le t_{file}^{max}$ でなければならない。 A real number that satisfies $t_{file}^{min}<t^{eruption}\le t_{file}^{max}$, where $t_{file}^{min}$ and $t_{file}^{max}$ are the beginning and end times of the input file, respectively.	$t_{file}^{max}-60$と$(t_{file}^{min}+t_{file}^{max})/2$の大きい方。 The larger one of $t_{file}^{max}-60$ and $(t_{file}^{min}+t_{file}^{max})/2$.
tmin_candidate	傾斜変動の候補となる時刻範囲の下限$t_{candidate}^{min}$(s)。 $t^{st}\le t_{candidate}^{min}$となるタイムウインドウはたとえ$M$が大きくとも傾斜変動ではないと判断する。時刻の原点は入力ファイル第3列、第4列と一致させること。 The lower limit, $t_{candidate}^{min}$, of the time range (s) to regard a time window being potentially a tilt change. A time window that meets $t^{st}\le t_{candidate}^{min}$ is regarded not being a tilt change even if $M$ is large. The time origin must be consistent with that of the 3rd and 4th columns of the input file.	実数。 $t_{file}^{min}\leq t_{candidate}^{min}<t^{eruption}$ でなければならない。 A real number that satisfies $t_{file}^{min}\leq t_{candidate}^{min}<t^{eruption}$.	$t_{file}^{min}+(t^{eruption}-t_{file}^{min})/2$
tmax_candidate	傾斜変動の候補となる時刻範囲の上限$t_{candidate}^{max}$(s)。 $t^{en}\ge t_{candidate}^{max}$となるタイムウインドウはたとえ$M$が大きくとも傾斜変動ではないと判断する。時刻の原点は入力ファイル第3列、第4列と一致させること。 The upper limit, $t_{candidate}^{max}$, of the time range (s) to regard a time window being potentially a tilt change. A time window that satisfies $t^{en}\ge t_{candidate}^{max}$ is regarded not being a tilt change even if $M$ is large. The time origin must be consistent with that of the 3rd and 4th columns of the input file.	実数。 $t_{candidate}^{min}<t_{candidate}^{max}\leq t_{file}^{max}$ でなければならない。 A real number that satisfies $t_{file}^{min}\leq t_{candidate}^{min}<t^{eruption}$.	$t^{eruption}+60$と$t_{file}^{max}$の小さい方。 The smaller one of $t^{eruption}+60$ and $t_{file}^{max}$.
Tsum_max	下記(\ref{eq.M34condition})式における閾値$T_{1-4}^{max}$の値。 A threshold value $T_{1-4}^{max}$ in Eq. (\ref{eq.M34condition}) below.	$t_{file}^{max}-t_{file}^{min}$以下の正の実数。 A positive real number less than or equal to $t_{file}^{max}-t_{file}^{min}$.	$(t^{eruption}-t_{candidate}^{min})/2$
M_thre	下記のアルゴリズムにおける$M_1$-$M_4$の閾値$M^{thre}$。 $M_1$-$M_4$がこの値よりも小さい場合は傾斜変動の候補としない。 The threshold value, $M^{thre}$, for $M_1$-$M_4$ in the algorithm below. Each candidate is regarded as not being a tilt change if $M_1$-$M_4$ is less than this threshold.	1よりも小さな正の実数。 A positive real number less than 1.	0.4
M_ratio_thre	下記のアルゴリズムにおける比$M_1/M_{ave}$-$M_4/M_{ave}$の閾値$M_{ratio}^{thre}$。比がこの値よりも小さい場合は傾斜変動の候補としない。 The threshold value, $M_{ratio}^{thre}$, for the ratios $M_1/M_{ave}$ to $M_4/M_{ave}$ in the algorithm below. Each candidate is regarded as not being a tilt change if the ratio is less than this threshold.	正の実数。 A positive real number.	2.0

◆動作(Behaviour)

第1引数で指定したテキストファイルを読み込み、噴火直前の傾斜変動を探索して結果を第2引数で指定したテキストファイルに出力する。出力ファイルの書式は以下の通り。
Read the text file given by the 1st argument, search the pre-eruptive tilt change, and output the result into the file given by the 2nd argument with the format below.

第1行には噴火直前の傾斜変動が検知された場合は
#Tilt change present
検知されなかった場合は
#Tilt change absent
と出力される。
The 1st line consists of
#Tilt change present
if a pre-eruptive tilt change was detected, and
#Tilt change absent
if not detected.
傾斜変動が検知された場合、次の行には
#Window: 値
と出力される。値は1-4のいずれかであり、 $[t_1^{st},t_1^{en}]$-$[t_4^{st},t_4^{en}]$のうち、どの番号のものが傾斜変動に対応するかを示している。
If the tilt change was detected, the next line consists of
#Window: value
where value is 1-4 that represents which of $[t_1^{st},t_1^{en}]$-$[t_4^{st},t_4^{en}]$ corresponds to the tilt change.
傾斜変動が検知された場合、次の行には
#Start: 値
と出力される。値は傾斜変動の開始時刻(s)が小数点以下3桁まで出力される。
If the tilt change was detected, the next line consists of
#Start: value
where value is the start time (s) of the tilt change with 3 digits in the decimal part.
傾斜変動が検知された場合、次の行には
#End: 値
と出力される。値は傾斜変動の終了時刻(s)が小数点以下3桁まで出力される。
If the tilt change was detected, the next line consists of
#End: value
where value is the end time (s) of the tilt change with 3 digits in the decimal part.
傾斜変動が検知されたか否かに関わらず、次の行には
#Ncandidates: 値
と出力される。値は傾斜変動の候補として検討したタイムウインドウの個数であり、 $[t_1^{st},t_1^{en}]$と$[t_2^{st},t_2^{en}]$ についてのみ検討した場合は2、 $[t_1^{st},t_1^{en}]$-$[t_4^{st},t_4^{en}]$ について検討した場合は4となる。また、どの候補も検討されなかった場合は0となる。これは$M_1=0$である場合(元の波形の全サンプルが0のときに起きる)や、 $[t_1^{st},t_1^{en}]$と$[t_2^{st},t_2^{en}]$だけで全タイムウインドウを占めてしまい $M_3$が定義できない場合に起きる。
Regardless of whether the tilt change was detected or not, the next line consists of
#Ncandidates: value
where value is the number of time windows examined as candidates of the tilt change. The value is 2 if only $[t_1^{st},t_1^{en}]$ and $[t_2^{st},t_2^{en}]$ were examined, and 4 if $[t_1^{st},t_1^{en}]$-$[t_4^{st},t_4^{en}]$ were examined. The value is 0 if no candidate was examined; this occurs when $M_1=0$ (which occurs when all samples of the original waveform are zero), or $[t_1^{st},t_1^{en}]$ and $[t_2^{st},t_2^{en}]$ occupies the entire time window so that $M_3$ was not defined.
次の行には
#Definable non-zero M_i: 値
と出力される。値は定義可能なノンゼロの$M_i$ ($i=1,2,\cdots$)の個数であり、0-8とする (9以上の場合であっても出力される値は8となる)。定義可能とは、$[t_1^{st},t_1^{en}]$から $[t_{i-1}^{st},t_{i-1}^{en}]$までのタイムウインドウで $[t_{file}^{min},t_{file}^{max}]$の全体を占めてしまわないということである。
The next line consists of
#Definable non-zero M_i: value
where value is the number of definable non-zero $M_i$ ($i=1,2,\cdots$) and takes a value from 0 to 8; if more than 8 $M_i$ are definable, the output value is 8. The definable means that the entire $[t_{file}^{min},t_{file}^{max}]$ is not occupied by $[t_1^{st},t_1^{en}]$ to $[t_{i-1}^{st},t_{i-1}^{en}]$.
上記のNcandidatesの値が0でない場合、次の行には
#Mave: 値
と出力される。値は計算に用いた$M_{ave}$である。
If the value of Ncandidates above is not 0, the next line consists of
#Mave: value
where value is $M_{ave}$ used in the computation.
次の行には
#i[TAB]t_i^st(s)[TAB]t_i^en(s)[TAB]M_i
と出力される。ここで[TAB]はタブを表す。
The next line consists of
#i[TAB]t_i^st(s)[TAB]t_i^en(s)[TAB]M_i
where [TAB] indicates a tab.
その次の行以降には $i$, $t_i^{st}$ (s), $t_i^{en}$ (s), $M_i$ の値がタブ区切りで出力される。 1行につき1つの$i$とし、$i$について昇順($i=1,2,\cdots$)に、定義可能なノンゼロの$M_i$の個数分だけ(最大:8)出力される。
The next and following lines consist of the values of $i$, $t_i^{st}$ (s), $t_i^{en}$ (s), and $M_i$, separated by tabs. Each line consists of the values for single $i$, and the data are sorted in an ascending order of $i$ (i.e., $i=1,2,\cdots$) until $i$ reaches the number of definable non-zero $M_i$ (maximum: 8).

以下に出力ファイルのサンプルを示す。
Sample output files are shown below.

#Tilt change present
#Window: 1
#Start: -422.000
#End: 30.000
#Ncandidates: 4
#Definable non-zero M_i: 8
#Mave: 1.530740e-03
#i[TAB]t_i^st(s)[TAB]t_i^en(s)[TAB]M_i
1[TAB]-422.000[TAB]30.000[TAB]9.892279e-01
2[TAB]30.000[TAB]60.000[TAB]8.182390e-02
3[TAB]-451.000[TAB]-422.000[TAB]4.144805e-03
4[TAB]-1739.000[TAB]-1723.000[TAB]2.797923e-03
5[TAB]-704.000[TAB]-655.000[TAB]1.830976e-03
6[TAB]-729.000[TAB]-704.000[TAB]1.542258e-03
7[TAB]-1609.000[TAB]-1585.000[TAB]1.425062e-03
8[TAB]-1723.000[TAB]-1707.000[TAB]1.324663e-03

#Tilt change absent
#Ncandidates: 4
#Definable non-zero M_i: 8
#Mave: 3.746137e-01
#i[TAB]t_i^st(s)[TAB]t_i^en(s)[TAB]M_i
1[TAB]-583.000[TAB]-473.000[TAB]6.106610e-01
2[TAB]-962.000[TAB]-906.000[TAB]4.576497e-01
3[TAB]12.000[TAB]51.000[TAB]4.333082e-01
4[TAB]-1500.000[TAB]-1478.000[TAB]4.119325e-01
5[TAB]-753.000[TAB]-711.000[TAB]4.018677e-01
6[TAB]-853.000[TAB]-778.000[TAB]3.788439e-01
7[TAB]-287.000[TAB]-258.000[TAB]3.596151e-01
8[TAB]-1082.000[TAB]-1045.000[TAB]3.581280e-01

◆使用例(Example)

detect_inflation waveform.sac M.txt --tmin=60.0 --tmax=1860.0
detect_inflation_preeruptive M.txt result.txt --M_thre=0.2

◆アルゴリズム(Algorithm)

以下で説明するアルゴリズムはMaeda (2023)の2.2.4節と同様である。
The algorithm explained below is same as that in Section 2.2.4 of Maeda (2023).

$[t_{file}^{min},t_{file}^{max}]$から$M(t^{st},t^{en})$の最大値を探索し、最大値を$M_1$、それに対応するタイムウインドウを$[t_1^{st},t_1^{en}]$とする。次に$[t_{file}^{min},t_{file}^{max}]$から $[t_1^{st},t_1^{en}]$を除いたタイムウインドウ (すなわち$[t_{file}^{min},t_1^{st}]$および$[t_1^{en},t_{file}^{max}]$)での $M(t^{st},t^{en})$の最大値を探索し、最大値を$M_2$、それに対応するタイムウインドウを$[t_2^{st},t_2^{en}]$とする。以下同様にして$M_3$, $M_4$, …および $[t_3^{st},t_3^{en}]$, $[t_4^{st},t_4^{en}]$, …を求める。
The maximum value of $M(t^{st},t^{en})$ is searched from $[t_{file}^{min},t_{file}^{max}]$. This maximum value is denoted as $M_1$, and the corresponding time window is denoted as $[t_1^{st},t_1^{en}]$. Next, the maximum value of $M(t^{st},t^{en})$ is searched from $[t_{file}^{min},t_{file}^{max}]$ outside of $[t_1^{st},t_1^{en}]$ (i.e., $[t_{file}^{min},t_1^{st}]$ and $[t_1^{en},t_{file}^{max}]$). This maximum is denoted as $M_2$, and the corresponding time window is denoted as $[t_2^{st},t_2^{en}]$. In the same way, $M_3$, $M_4$, … and $[t_3^{st},t_3^{en}]$, $[t_4^{st},t_4^{en}]$, … are defined.

$[t_1^{st},t_1^{en}]$は傾斜変動である可能性が最も高いタイムウインドウである。しかし噴火直前の膨張ではなく噴火時の収縮を表している場合がある。そこで$[t_2^{st},t_2^{en}]$が2つ目の候補となる。典型的な波形においては $[t_1^{st},t_1^{en}]$と$[t_2^{st},t_2^{en}]$が膨張・収縮ペアをなし、 (どちらが膨張に対応するかはケースバイケースであるが) このペアが噴火直前の膨張と噴火時の収縮を表す。
The time window $[t_1^{st},t_1^{en}]$ most likely consists of a tilt change; however, it may be a syn-eruptive deflation. Therefore, $[t_2^{st},t_2^{en}]$ is the secondary candidate. A typical waveform consists of a pair of pre-eruptive inflation and syn-eruptive deflation in $[t_1^{st},t_1^{en}]$ and $[t_2^{st},t_2^{en}]$; which window corresponds to the inflation depends on cases.

イベント的な膨張・収縮が卓越する波形では、 $M_1$, $M_2$は$M_3$以降と比べて有意に大きくなると期待される。一方、同程度の振幅での振動が続く波形では $M_1$, $M_2$と$M_3$以降は近い値になると考えられる。そこで$M_3$と$M_4$の平均値を$M_{ave}$とし、 $M_1/M_{ave}\geq M_{ratio}^{thre}$, $M_2/M_{ave}\geq M_{ratio}^{thre}$の場合に $[t_1^{st},t_1^{en}]$, $[t_2^{st},t_2^{en}]$を傾斜変動イベントの候補とみなす。 $M_1/M_{ave}<M_{ratio}^{thre}$, $M_2/M_{ave}<M_{ratio}^{thre}$の場合にはこの波形に傾斜変動イベントは含まれないと判断する。
In a waveform that consists of significant inflation and deflation, $M_1$ and $M_2$ are expected to be significantly larger than $M_3$ and later values. Contrary, these values are expected to be similar in a waveform that consists of a long-lasting oscillation of nearly constant amplitudes. Based on this consideration, $[t_1^{st},t_1^{en}]$ and $[t_2^{st},t_2^{en}]$ are regarded as candidates of a pre-eruptive tilt change if $M_1/M_{ave}\geq M_{ratio}^{thre}$ and $M_2/M_{ave}\geq M_{ratio}^{thre}$, respectively, where $M_{ave}$ is an average of $M_3$ and $M_4$; the waveform is regarded to have no tilt change event if $M_1/M_{ave}<M_{ratio}^{thre}$ and $M_2/M_{ave}<M_{ratio}^{thre}$.

波形によっては噴火とは無関係の膨張・収縮が $[t_{candidate}^{min},t_{candidate}^{max}]$の中に偶然に入ってしまう場合がある。そこで個々の候補の時間が短い場合には $[t_3^{st},t_3^{en}]$, $[t_4^{st},t_4^{en}]$ も噴火直前の傾斜変動の候補に加える。具体的には \[\begin{equation} \sum_{i=1}^4 (t_i^{en}-t_i^{st}) \le T_{1-4}^{max} \label{eq.M34condition} \end{equation}\] の場合に $[t_3^{st},t_3^{en}]$, $[t_4^{st},t_4^{en}]$ を候補に加える。この場合には$M_{ave}$は$M_5$-$M_8$の平均とし、比$M_1/M_{ave}$, $M_2/M_{ave}$, $M_3/M_{ave}$, $M_4/M_{ave}$ について上記と同様の判断を行う。なお(\ref{eq.M34condition})式が満たされない場合 (個々の膨張・収縮の周期が長くこの4候補だけで波形の大部分を占めてしまう場合) には4つ目の候補まで考えると振動との区別がつかないので上位2つ($[t_1^{st},t_1^{en}]$, $[t_2^{st},t_2^{en}]$)のみを候補とする。
The waveform may consist of occational inflation and deflation in $[t_{candidate}^{min},t_{candidate}^{max}]$ not related to eruptions. To account for it, $[t_3^{st},t_3^{en}]$ and $[t_4^{st},t_4^{en}]$ are included as candidates of the tilt changes if individual time windows are short enough to satisfy Eq. (\ref{eq.M34condition}). In this case, $M_{ave}$ is defined as an average of $M_5$-$M_8$, and the ratios $M_1/M_{ave}$, $M_2/M_{ave}$, $M_3/M_{ave}$, and $M_4/M_{ave}$ are used to determine if each time window is a candidate of the pre-eruptive tilt change, in the same manner as described above. If Eq. (\ref{eq.M34condition}) is not satisfied, which means that the durations of individual inflation and deflation are long enough that the four candidates occupy most parts of the entire waveform, than only the top two candidates ($[t_1^{st},t_1^{en}]$ and $[t_2^{st},t_2^{en}]$) are used because otherwise it would be difficult to distinguish oscillation and transient signals.

噴火直前の傾斜変動と判断するには更にいくつかの条件を要する。まず$M_1<M^{thre}$の場合には比$M_1/M_{ave}$に関わらず噴火直前の傾斜変動ではないと判断する。また$t_1^{st}\geq t^{eruption}$または $t_1^{en}\geq t_{candidate}^{max}$の場合には $[t_1^{st},t_1^{en}]$は噴火時の変動を表している可能性が高いため噴火直前の傾斜変動ではないと判断する。 $t_1^{st}\leq t_{candidate}^{min}$の場合には $[t_1^{st},t_1^{en}]$は波形全体にわたる長周期ノイズを表している可能性が高いため噴火直前の傾斜変動ではないと判断する。 $t_1^{en}+(t_1^{en}-t_1^{st})<t^{eruption}$の場合には噴火直前と見なすには噴火との時間差が大きすぎるため噴火直前の傾斜変動ではないと判断する。 $[t_2^{st},t_2^{en}]$, $[t_3^{st},t_3^{en}]$, $[t_4^{st},t_4^{en}]$ についても同様の判断を行う。
Several additional conditions are needed to determine a pre-eruptive tilt change. If $M_1<M^{thre}$, it is regarded as not being a pre-eruptive tilt change regardless of the ratio $M_1/M_{ave}$. If $t_1^{st}\\geq t^{eruption}$ or $t_1^{en}\geq t_{candidate}^{max}$, $[t_1^{st},t_1^{en}]$ likely represents a syn-eruptive tilt change and thus is regarded as not being a pre-eruptive tilt change. If $t_1^{st}\leq t_{candidate}^{min}$, $[t_1^{st},t_1^{en}]$ likely represents a long-period noise over the waveform, so that it is regarded as not being a pre-eruptive tilt change. If $t_1^{en}+(t_1^{en}-t_1^{st})<t^{eruption}$, the time window is too apart from the eruption time to regard as a pre-eruptive tilt change. Same checks are conducted for $[t_2^{st},t_2^{en}]$, $[t_3^{st},t_3^{en}]$, and $[t_4^{st},t_4^{en}]$.

全ての条件をクリアしたタイムウインドウが複数ある場合には、それらの中で$(t_i^{st}+t_i^{en})/2\le t^{eruption}$(中心時刻が噴火開始以前) となる最初の$i$に対する$[t_i^{st},t_i^{en}]$を噴火直前の傾斜変動と判断する。但し全ての$i$について$(t_i^{st}+t_i^{en})/2> t^{eruption}$の場合は最初の$i$に対する$[t_i^{st},t_i^{en}]$を噴火直前の傾斜変動と判断する。
If there remained multiple time windows that passed all these examinations, $[t_i^{st},t_i^{en}]$ for the first $i$ of them that satisfies $(t_i^{st}+t_i^{en})/2\le t^{eruption}$ (i.e., the central time is before the eruption time) is regarded as the pre-eruptive tilt change; if $(t_i^{st}+t_i^{en})/2> t^{eruption}$ for all $i$, $[t_i^{st},t_i^{en}]$ for the first $i$ of them is regarded as the pre-eruptive tilt change.

◆補足 (Additional notes)

detect_inflationコマンドのマニュアルにMaeda (2023)の計算を再現するスクリプトの例を示してある。
A sample script for reproducing Maeda (2023) is available in the documentation of detect_inflation command.

◆引用文献 (References)

★ は本プログラムを用いた研究成果の発表時に必ず引用すべき論文を表す。
The star (★) indicates the references that must be cited whenever a research that used this program is published or presented.

★ Maeda Y (2023) A systematic survey for precursory tilt changes at all monitored eruptions in Japan, J Volcanol Geotherm Res, 439, 107831. https://doi.org/10.1016/j.jvolgeores.2023.107831

detect_inflation_preeruptive2コマンド マニュアル