日本の二酸化炭素濃度の時系列推移(月次)

R, RStudio, Quarto and STEM(Science, Technology, Engineering and Mathematics)

日本国内に設けられた3つの観測所(大気環境観測所(大船渡市綾里)、南鳥島気象観測所および与那国島特別地域気象観測所 )において測定されています二酸化炭素濃度の時系列推移を確認します。

なお、与那国島での観測は2024年3月末で終了しています。

1 用語の定義

  • ppmは大気中の分子100万個中にある対象物質の個数を表す単位です。
  • 大気中の温室効果ガス濃度は、濃度があらかじめわかっている標準ガスと比較して求めますが、標準ガスの濃度は使用期間中に変化することがあります。このため、観測に使用した標準ガス濃度の変化を確認した段階で、観測値にその値を反映させ確定値とします。速報値と確定値の濃度差は、通常二酸化炭素で0.1ppm以下、メタンで1ppb以下です。

    出典 : https://www.data.jma.go.jp/ghg/kanshi/obs/co2_monthave_ryo.html

  • 気象庁の二酸化炭素やメタンなどの温室効果ガスの観測では、広い水平空間を代表する平均的な濃度（バックグランドデータ）を求めています。観測所で観測される濃度値には、気象状況などにより生じる局地的で時間スケールの短い濃度変動が含まれているため、これらの条件のデータを選別し除去することにより、広い水平空間を代表する平均的な濃度と考えられるバックグランドデータを求めることができます。具体的には以下の手順によってバックグランドデータを求めます。
  1. 測器の点検、故障時などを除いた、全ての観測値を1時間平均して時別値を求めます。
  2. 時別値を求める際、算出に用いる観測値が本来得られるべき数の半数未満の場合は時別値を求めません。
  3. 時別値を求める際、算出に用いる観測値に一定以上のばらつきがある場合は時別値を求めません。
  4. 時別値を求める際、その前後の時別値との差が両方とも一定以上の値を超えた場合は時別値を求めません。
  5. 上記の手順により残った時別値を、バックグランドデータとします。
  6. 日別値は、バックグランドデータとして選択された時別値を平均します。
  7. 月別値は、バックグランドデータとして選択された時別値を平均します。
  8. 年平均濃度は、求められた月別値を平均します。

出典 : https://www.data.jma.go.jp/ghg/kanshi/obs/month_calculation.html

2 観測所の所在地

3 時系列データ

3.1 対象期間

観測所 データ開始 データ終了
1 綾里 1987年01月 2025年10月
2 南鳥島 1993年02月 2025年10月
3 与那国島 1997年01月 2024年03月
Table 1

3.2 時系列チャート

Figure 1
Figure 2
Figure 3

3.3 時系列テーブル

Table 2

4 線形回帰モデル

4.1 綾里

Observations 465 (1 missing obs. deleted)
Dependent variable 月平均値(綾里)[ppm]
Type OLS linear regression
F(1,463) 8497.2239
R² 0.9483
Adj. R² 0.9482
Est. S.E. t val. p
(Intercept) 346.8493 0.4982 696.1788 0.0000
Month 0.1705 0.0018 92.1804 0.0000
Standard errors: OLS
Table 3

4.2 南鳥島

Observations 391 (75 missing obs. deleted)
Dependent variable 月平均値(南鳥島)[ppm]
Type OLS linear regression
F(1,389) 16586.3007
R² 0.9771
Adj. R² 0.9770
Est. S.E. t val. p
(Intercept) 341.7883 0.4093 835.1560 0.0000
Month 0.1798 0.0014 128.7878 0.0000
Standard errors: OLS
Table 4

4.3 与那国島

Observations 326 (140 missing obs. deleted)
Dependent variable 月平均値(与那国島)[ppm]
Type OLS linear regression
F(1,324) 9130.1925
R² 0.9657
Adj. R² 0.9656
Est. S.E. t val. p
(Intercept) 341.3519 0.5838 584.7449 0.0000
Month 0.1865 0.0020 95.5520 0.0000
Standard errors: OLS
Table 5

5 月別中央値

5.1 ボックスプロット

Figure 4
Figure 5
Figure 6

5.2 テーブル

  • 各観測所の二酸化炭素濃度月平均値の月別中央値(ppm)
月 綾里
1 5 390.4
2 3 389.5
3 4 389.4
4 2 388.6
5 1 387.7
6 12 387.4
7 11 386.0
8 6 384.9
9 10 384.2
10 9 379.1
11 7 377.8
12 8 377.8
(a)
月 南鳥島
1 5 391.5
2 3 391.1
3 4 390.5
4 1 390.1
5 2 389.8
6 6 389.0
7 12 387.8
8 7 387.3
9 8 386.1
10 10 386.1
11 11 385.9
12 9 385.5
(b)
月 与那国島
1 3 396.7
2 4 396.6
3 2 396.1
4 5 395.3
5 12 395.3
6 1 395.1
7 6 393.2
8 11 392.4
9 7 390.8
10 10 389.0
11 8 388.8
12 9 387.9
(c)
Table 6

6 データ出所

  • 気象庁

7 最終更新

Sys.time()
[1] "2025-12-28 11:43:56 JST"

8 R & Quarto

R.Version()$version.string
[1] "R version 4.5.1 Patched (2025-08-22 r88678 ucrt)"
quarto::quarto_version()
[1] '1.8.26'

9 著者