Page 1

hmiya@icrr.u-tokyo.ac.jp

~7000 year old cedar tree

東京大学　宇宙線研究所

中世の温暖期と近世の小氷期における

太陽活動と気候変動

宮原ひろ子

－樹木年輪中の炭素同位体の分析から－

topics

◇ Introduction

太陽活動は気候に作用するのか？

太陽活動の長期変動と小氷期

太陽活動、太陽黒点数、宇宙線の11年変動

◇ Measurements and Results

樹木年輪中の炭素14濃度測定による太陽活動の復元

太陽活動11年周期の変遷と気候変動への影響　－小氷期と中世温暖期を例に－

◇ Appendix

太陽活動度の長期変動を正確に理解するために

日射量の可変性を正確に理解するために

太陽活動の予測に向けて

IPCC第4次報告書（2007年）より

※ただし、1750AD以降の日射量変動分しか含まれていない

太陽活動は気候に作用するのか？

気候変動に見られる周期　（赤枠：太陽周期と近似しているもの）

From Hoyt and Schatten, 1998

The role of the Sun in Climate Change

過去1200年間における気温変化と太陽活動

小氷河期

中世の温暖期

黒点

年輪中炭素14による

南極氷床中ベリリウム同位体による

全球氷床中ベリリウム同位体による

西暦年

太陽活動活発期

太陽活動極小期（静穏期）

（Usoskin, 2004）

（IPCC2007より）

炭素14は太陽活動の

指標であるが

気温との相関が良いため、

歴史的に気温の指標として

用いられていたこともある

日本における中世温暖期と小氷期

炭素13による

屋久杉年輪中炭素13による気温復元値

（Kitagawa and Matsumoto, 1995）

琵琶湖沿岸の熱伝導計測による気温復値

(Goto et al., 2005)

Year AD

←北半球の気温

（IPCC 2007より）

日本のローカルな気温の復元値⇒

熱伝導計測による

日射量, 紫外線

銀河宇宙線

あるいは

気候システムの内的要因

気候の十年変動を駆動しているのは何か？

(Svensmark, 2007)

Svensmark, 2007

宇宙線

雲量

0.1%

3%

15%

銀河宇宙線（宇宙から飛来する荷電粒子）

【大気をイオン化し、雲凝結核の生成を

促進する】という間接的な影響を持つと

考えられているが、

詳細なメカニズムはまだ解明されていない

1975

1980

1985

1990

1995

2000

Year AD

0

50

100

150

200

250

太陽活動の11年周期

ようこう衛星による太陽のX線画像

太陽活動極大期

太陽活動極小期

太陽黒点数

衰退

活発化

ひので衛星がとらえた巨大黒点

日射量変動

Total Solar Irra

dia

n

ce

[W/m

2

]

Year

Mon

th

ly Sunspot Nu

mbe

r

ACRIM TSI composite

太陽圏

太陽から噴き出した

磁場とプラズマの風の

ひろがり

太陽圏の磁場は宇宙を飛びかう放射線（荷電粒子）から

地球を防護する役割を果たしている。

しかし、そのバリアの強度は、太陽活動度に依存して変化する

⇒ 地球に飛来する宇宙線は時間とともに変化する

（たとえば太陽の11年周期変動にともなう増減）

↓コロナグラフでみた太陽風

太陽

太陽圏

日射量変動と根本的に異なる宇宙線量の振る舞い … 太陽の磁場の影響

2000

2500

3000

3500

4000

4500

1950 1960 1970 1980 1990 2000

Year AD

0

100

200

300

400

黒点数と宇宙線量

Ｎ

Ｓ

Ｎ

Ｓ

太陽圏の磁場は複雑な構造をしているため、

太陽の双極子磁場の向きが上向きの時と下向きの時では

宇宙から飛来する荷電粒子（宇宙線）に対する遮蔽効果が変わってくる

そのため、地球に飛来する宇宙線の量には

11年周期に加えて22年周期（太陽磁場の反転の周期）も見られる

太陽黒点数

地球に飛来する

宇宙線量

太陽は双極子磁場を持ち、

11年に一度、黒点数の極大で

正⇒負、あるいは負⇒正に反転する

太陽活動度、日射量、黒点数には

11年周期しか見られないが、

太陽磁場の向きと宇宙線変動には

22年周期が見られる

太陽圏磁場の構造

大きさ：100AU （天文単位）

1AUは、太陽から地球までの距離

中心に太陽

Ｎ

Ｓ

太陽圏の磁場の形＆太陽の磁場の向き　が宇宙線の来やすさに与える影響

太陽活動の活発さ

宇宙線量

(Kota, 1983; 2003)

宇宙線

100AU

100AU

100AU

100AU

Ｎ

Ｓ

Ｎ

Ｓ

Ｎ

Ｓ

宇宙線

宇宙線

宇宙線

太陽の磁石が上向き

太陽の磁石が下向き

太陽活動が活発なとき

太陽活動が静穏なとき

太陽圏磁場

太陽圏磁場

太陽圏磁場

太陽圏磁場

Ｎ

Ｓ

Ｎ

Ｓ

太陽活動が静穏なほど

宇宙線は来やすい。

太陽の磁石が上向きの方が

宇宙線は来やすい。

0

50

100

150

200

250

1600

1700

1800

1900

2000

マウンダー極小期

(西暦1645-1715 年)

太陽黒点数の長期変動とマウンダー極小期

太陽黒点数

太陽活動の

長期変動周期

～88年周期

～208年周期

望遠鏡の発明

テムズ川が凍ったとされる

小氷期の絵画

西暦年

マウンダー極小期のような黒点の消失は過去に度々

発生していたことが年輪中炭素14により明らかになっている

たとえば

シュペーラー極小期　（西暦1416-1534年）

ウォルフ極小期

（西暦1282-1342年）

など

太陽活動の復元方法　（炭素14を指標とする手法）

太陽活動が活発化

⇒ 太陽（圏）の磁場が宇宙線をバリア

⇒ 地球に飛んでくる宇宙線（陽子）量が減少

⇒ 宇宙線（陽子）が作る炭素14の量が減少

⇒ 樹木年輪に取り込まれる炭素14の濃度が減少

二酸化炭素として循環（ 14CO

2

）

● 陽子

○ 中性子

（陽子）

大気分子

地球

光合成により年輪に取り込まれる

太陽圏（磁場）

宇宙線

窒素原子核

炭素14

中性子

高エネルギー宇宙線

(a) 太陽圏（太陽の磁場が広がる空間）が

宇宙を飛びかう高エネルギー粒子をさえぎる

(b) 宇宙線によって炭素14が生成される

Carbon cycle

surface ocean

deep ocean

troposphere

marine sediment

soil

stratosphere

biosphere

炭素14の循環

宇宙線の変動に起因する炭素14の変動は

炭素循環のなかで減衰する

-12

-9

-6

-3

0

3

6

1880

1900

1920

1940

Year AD

Δ

14C

anom

aly (per m

ill)

0

100

200S

unspot N

um

ber

2000

2500

3000

3500

4000

4500

1950 1960

1970 1980 1990 2000

Year AD

N

eutron Flux at C

lim

ax

0

100

200

300

400

Sunspot N

u

m

b

er .

黒点数と炭素14濃度

（大気循環による3年のタイムラグ補正済み）

黒点数⇒宇宙線量⇒年輪中炭素14

Blue : 14C data

Red : 3-yr running ave

黒点数と宇宙線量

太陽活動活発（黒点数増）

= 14C濃度低下

太陽活動静穏（黒点数減）

= 14C濃度上昇

Ｎ

Ｓ

Ｎ

Ｓ

Ｎ

Ｓ

Ｎ

Ｓ

Ｎ

Ｓ

年輪中炭素14濃度を測定すると

太陽活動の変動の周期性を

復元することができる。

しかし、濃度自体は活動度レベルの

正確な指標とはならない。

ただし、周期長が太陽活動度の指標として

使える可能性がある。

Motivation

◇ 太陽活動の11年周期は本当に気候に影響するのか？

◇ 太陽活動と気候変動を媒介しているものは何なのか？

日射量？　宇宙線？　あるいは・・・？

樹木年輪中の炭素14濃度を1年ごとに測定し、

黒点数のデータの無い中世の太陽活動活発期やマウンダー極小期

における太陽活動周期（～11年周期）、磁場反転周期（～22年周期）

の特徴を明らかにすれば、

太陽活動と気候変動の関係性とそのメカニズムを過去にさかのぼって

調べることができる

炭素14濃度測定法

年輪の剥離から炭素14濃度測定まで

１．年輪の絶対年代の決定

年輪数のカウント

核実験による炭素14濃度のピーク（西暦1964年）の検出

年輪幅パターンによる年代決定

２．年輪の剥離

３．セルロース抽出　（年輪間を移動しない成分）

酸、アルカリによる洗浄

亜塩素酸ナトリウムによる漂白

４．セルロース試料からグラファイトを合成

セルロース燃焼 ⇒ 二酸化炭素 ⇒ 水素還元 ⇒ グラファイト

５．グラファイト中の炭素14濃度を測定

加速器質量分析計による高精度分析　（東大、名古屋大ほか所有）

樹齢~1950年の

屋久杉の切り株

30N

40N

136E

144E

128E

東京

屋久島

室生寺杉

1998年の台風

により倒木

奈良県室生寺

０．試料採取

0

200

400

600

800

1952 1958 1964 1970 1976 1982 1988

0

200

400

600

800

1952 1958 1964 1970 1976 1982 1988

樹齢712年の屋久杉

樹齢382年の室生寺杉

１．年輪の絶対年代の決定 – 核実験による炭素14濃度ピークの検出 –

1963年の部分的核実験禁止条約の施行を前に、各国により相次いで大気圏内での核実験が行われ、

その結果、大気中の炭素14濃度は、1964年に従来の約2倍に達した。そのピークの痕跡は年輪中に残されている。

西暦年

西暦年

年輪中炭素

14

濃度

（千分率）

年輪中炭素

14

濃度

（千分率）

木材からの収率　約40％

２．年輪の剥離　＆　３．セルロース抽出

カッターを用いて

年輪を１枚ずつ剥離

薬品による洗浄・漂白

セルロース

年輪

木材ブロック

（2cm ×2cm）

年輪幅　～1mm

酸化剤

セルロース

４．セルロース試料からグラファイトを合成

（セルロース燃焼 ⇒ 二酸化炭素）

計量し

酸化剤とともに

ガラス管へ

真空ラインを用いて

ガラス管内を真空に引き

バーナーでシーリング

950℃のオーブンで燃焼

４．セルロース試料からグラファイトを合成

（二酸化炭素 ⇒ 水素還元 ⇒ グラファイト）

二酸化炭素を

真空ラインへ導入

二酸化炭素に水素ガスを足して

630℃に加熱し、水素還元

⇒グラファイト（炭素）の完成

1.5mgのグラファイトを

ターゲットホルダーに詰めて

加速器のイオン源へ

５．　加速器質量分析計によるグラファイト中の炭素14濃度の測定

グラファイトをイオン化 ⇒ 炭素イオンの電流を加速器へ

電磁石で軽い炭素12と重い炭素14を分離 ⇒ 14C量、12C量を測定

加速管

電磁石

イオン源

測定結果

太陽活動11年周期の変遷と

気候変動への影響

0

500

1000

1500

2000

2500

3000

3500

4000

4500

5000

0.02

0.03

0.04

0.05

0.06

0.07

0.08

0.09

周波数

時間(t)

時系列データ y(t)

ウェーブレットφ(x)

フーリエ変換とは

W(f,t) = y(t) ・φ(x) dx

振幅大

振幅小

ウェーブレット変換とは

周波数

振幅

平面波 Φ*

(x)

を用いた周波数解析

フーリエ変換は平面波を

用いた周波数解析であるため

時間方向の情報を失うが、

ウェーブレット変換は

波束を用いているため、

周波数が時間とともに変化する

様子をとらえることができる

F(f) = y(t) ・φ*(x) dx

測定結果と周期解析結果：　過去1100年間における太陽”11年”周期の変遷

(10年値)

マウンダー極小期

シュペーラー極小期

ダルトン極小期

オールト極小期

~9年周期

11~13 年周期

10~11年周期

9~11年周期

11~14年周期

~14年周期

~13年周期

Δ14

C (pe

rmil)

~11 年周期

（青点: Miyahara et al., 2004, 2006, 2007

黒点: Stuiver, 1998; 灰点：Damon, 2003

）

中世の活発期

中世の活発期

西暦年

-30

-10

10

30

800

1100

1400

1700

2000

-4

0

4

8

1600

1650

1700

1750

1800

Year AD

Δ14

C

(per m

ill)

0

100

200

p

Miyahara et al., Sol. Phys. 2004,

JGR, 2006,

EPSL, 2008

14C

黒点

周期解析結果⇒

Temperature

(Vinther, 2003)

(Miyahara, et al., EPSL, 2008)

活発

静穏

マウンダー極小期における太陽活動11年周期の気候への影響

通常時

マウンダー極小期

太陽極性が正のときのほうが

寒くなりやすい

太陽極性が負のときのほうが

寒くなりやすい

Solar activity

太陽活動の活発さ

宇宙線量

(Kota, 1983; 2003)

Ｎ

Ｓ

Ｎ

Ｓ

+

-

Correlation between d18O and

relative humidity in June in Japan

65

70

75

80

85

90

95

1970

1980

1990

2000

2010

25.5

26

26.5

27

27.5

28

28.5

29

29.5

Correlation coeff.

年輪中δ18O による

日本の相対湿度の復元

δ18O

相対湿度

Year AD

・梅雨前線の活発度

（＝相対湿度）に依存した

蒸散時の同位体分別

グリーンランドの寒冷化と日本の湿度変化の比較

グリーンランドの氷床δ18O

による北半球の気温

（Vinther et al., 2003）

室生寺スギ中δ18O

による奈良の湿度

（≒降水量？）

(Our results, 2008)

グリーンランド、日本とも、気候変動に

太陽磁場極性依存がみられる

青：太陽11年周期の極小 + 太陽極性負

赤：太陽11年周期の極小 + 太陽極性正

マウンダー極小期

測定誤差0.2%

Our results (Yamaguchi et al.)

Vinther, 2003

太陽の磁石の向きが変わるだけでも、地球の気候は変わる

＝　太陽の磁場が地球に飛来する宇宙線の量を左右することで、

地球の気候に影響している可能性が大きい

－

＋

Solar max min

max min max

湿潤

乾燥

マウンダー極小期における湿潤イベントの太陽磁場極性依存

極性負

極性正

活発

衰退

太陽活動度

相対湿度

Our results (Yamaguchi, in prep)

Ｎ

Ｓ

Ｎ

Ｓ

・Δt は、活動度が低いほど長くなる

気温

黒点

気温

気温

マウンダー極小期

通常期（活発期）

気候の”11年/22年”周期のモード

移行期

Miyahara et al., Earth & Planetary

Science Letters, 2008

太陽活動の活発さ

宇宙線量

(Kota, 1983; 2003)

Ｎ

Ｓ

Ｎ

Ｓ

Maunder Minimum

Present

気温（年輪幅）変動にみられる”22年”周期　－太陽の双極子磁場の向きの影響ー

ダルトン極小期

気候の周期 = ~28年

0.06

0.04

0.02

0

[℃]

マウンダー極小期

気候の周期 = ~19 年

中世の温暖期

(太陽11年周期 = ~14 年)

(太陽11年周期 = ~9 年)

気候の周期 = ~26年

太陽活動の活発さ

宇宙線量

Ｎ

Ｓ

Ｎ

Ｓ

Maunder Minimum

Present

(Kota, 1983; 2003)

(太陽11年周期 = ~13 年)

太陽活動レベルの長期変動を

正確に理解するために

(Vonmoos, 2006)

(Solanki, 2004;

2005)

400

0

Φ (MeV

)

用いる核種によって異なる太陽活動の長期変動に関する見解

14Cによる

10Beによる

年輪中炭素14濃度（データ） → [炭素循環] → 宇宙線による炭素14生成量 → 太陽活動度の推定

気候変動による炭素循環の変化は考慮されていない

人為起源CO2放出による炭素14の希釈効果が完全には補正できない

20世紀の太陽活動は異常？！

20世紀の太陽活動は異常なほど活発だったのか？

“11年周期変動の伸び縮み”に基づいた太陽活動度レベルの復元の可能性

(Decadal)

Maunder

Spoerer

Wolf

Dalton

Oort

~9 yrs

11~13 yrs

10~11 yrs

9~11 yrs

11~14 yrs

~14 yrs

~13 yrs

Δ14

C (permil)

~11 yrs

中世の太陽活動活発期

太陽活動11年周期の実際の長さ（年）

次のサイクルの最大黒点数

Watari, 2008

太陽活動衰退期

周期解析結果⇒

-30

-10

10

30

800

1100

1400

1700

2000

Year AD

中世の太陽活動活発期では太陽周期は約9年であった。

20世紀の太陽活動よりも活発であった可能性が高い。

人為起源の温暖化ガスやヒートアイランド現象等に加えて

どの程度太陽が温暖化に起因したのかを探るには

中世と現代との気候の比較がカギ。

日射量の可変性を

正確に理解するために

Total Solar Irradiance Database

日射量

[W/m

2

]

Year

1975 1980 1985

1990 1995 2000 2005

1375

1370

1365

1360

250

0

太陽黒点数

日射量のデータベース

日射量に長期トレンドはあるのか？

1975

1980

1985

1990

1995

2000

2005

1375

1370

1365

1360

?

0

50

100

150

200

250

1600

1700

1800

1900

2000

Year AD

太陽黒点数

1363

1364

1365

1366

1367

1600

1700

1800

1900

2000

Year

T

S

I (W

/m

2)

Wang (2005)による

Lean (2000)による

Maunder Minimum

Maunder Minimum

?

400年前のマウンダー極小期では

日射量はどれだけ減っていたのか？

日射量に長期トレンドはあるのか？

月面10mの温度測定による過去500年間の日射量復元

日射により月面が暖められる

0 m

10 m

約 500年前の月面温度（=日射量）の情報

時間とともに熱は深部に伝播

大気がほとんどない月では、月面の温度は日射量のみによって決定される。

月面の熱は拡散により深部に伝播する。そのため、深さ方向に順次、温度を測定して

いくことで、過去の日射量の変化を復元することができる。

岩盤

微隕石の衝突

によって

積もった砂

（レゴリス）

↑月面の

熱収支

中緯度

（赤道，極域）

について

(W/m^2)

月面の熱伝導シミュレーション結果

1363

1364

1365

1366

1367

1600

1700

1800

1900

2000

Year

T

S

I (W

/m

2)

Wang, 2005

Lean, 2000

マウンダー極小期における

日射量のモデルの差異 ~2 W/m2

→ 月面深度～8mにて ~0.03 K の温度差

温度偏差

Wangのシナリオ – Leanのシナリオ (K)

月の赤道

月の南極

(Miyahara et al.,

Geophysical Research Letters, 2008)

太陽活動の予測に向けて

太陽はマウンダー極小期に向かっているのか？

1996年の最後の極小期から

すでに12年が経過している

（NASAによる予測

in 2004）

spaceweather.com

0

100

200

300

400

1975

1980

1985

1990

1995

2000

2005

2010

Year AD

太陽黒点数

science.nasa

太陽はマウンダー極小期に向かっているのか？

←

過去50年間で最低レベルの太陽風

↓宇宙線量の増加

次の太陽活動周期における最大黒点数の予測

Pesnell, 2007 より

太陽黒点数の予測

0

100

200

300

400

1975

1980

1985

1990

1995

2000

2005

2010

太陽黒点数

太陽活動を予測する物理モデルは

まだ確立されていない

Hathaway, Solar Physics, 2004

Hathaway, ApJ, 2003

黒点データにより導出された黒点数と周期長の関係性

太陽活動の11年周期が伸びると

次のサイクルの最大黒点数（振幅）

が少なくなる（小さくなる）

Watari, 2008

太陽活動11年周期の実際の長さ（年）

次のサイクルの最大黒点数

0

100

200

300

400

1975

1985

1995

2005

太陽活動度（黒点数）の増減の経験則

1996年に始まった最後の太陽サイクルは

終わるまでに12年以上かかった

⇒ 次のサイクルは活動が弱くなる

可能性が大きい？！

黒点数

-8

-5

-2

1

4

1750

1770

1790

1810

1830

Year AD

Δ

14

C

(p

e

r m

ill)

0

100

200

300

S

unspot N

um

ber

いずれの型の太陽活動極小期の場合も、

発生の直前にサイクル周期が～１３年程度に伸びている

Year AD

Year AD

シュペーラー極小期

マウンダー極小期

11yr

9yr

~13yr

13yr

13yr

ダルトン極小期

Frequency (1/yr)

Frequency (1/yr)

太陽活動極小期の前兆現象？

気候モデルへのインプットとして

- 日射量の可変性に関する正確な理解

- 太陽活動度レベルの長期変動　（周期長の復元にもとづく）

宇宙線が気候を左右するメカニズムは未だ解明されていない

- イオン化率の変化が雲形成に与える影響の素過程を調べる実験に期待

・スイスCERNでの基礎実験（SKY、CLOUD実験）

・国内でも名古屋大、甲南大などによりエアロゾル生成実験がスタート

まとめと展望

気候が持つ数十年スケールの複雑な変動は、

太陽と宇宙線が持つ複雑な周期的変動の影響で説明できる可能性がある。

～200年スケールで変化する太陽活動の”11年”周期（9年～15年）

太陽磁場反転の影響による”22年周期” （18～30年）