温暖化がサクラの開花期間に及ぼす影響
Inﬂuence of the global warming on blooming periods of cherry
塚原 あずみ1＊

1＊Azumi TSUKAHARA

・林陽生

22and Yousay HAYASHI

NPO 法人環境カウンセラー東京城北 ECA
2筑波大学生命環境系
1Tokyo Johoku Environment Counselor Association, NPO
2
Graduate School of Life and Environmental Sciences, University. of Tsukuba
1摘要

従来の研究では，地球温暖化とサクラの開花期の関係が多く取り扱われている。本
研究では，桜の開花日と満開日について，気温との関係に基づき，将来の気候変動が桜の開花期間にどのような変化をもたらすか評価した。開花から満開の期間を 1 つの生育ステージと考え，高緯度地域ほど開花期間が短いという地域特有の傾向の補正を行ったのち，発育速度 DVR の概念を適用し，開花期間を推定するモデルを構築した。
このモデルをもとに，IPCC-AR4 で提供された将来の気候を用いた予測を行った。その結果，全国的な傾向として，温室効果ガス排出シナリオや予測対象とする年代および地点によらず，今世紀中ごろまでに開花期が遅くなり，現在より 2.1 日減の 6.2 日で開花から満開に至ることが明らかになった。この結果は，サクラを観光資源としてみた場合，お花見の日数が減ることになり人びとが楽しむ期間が減り，観光資源が減少することを示唆している。

キーワード：開花－満開期間，観光産業，気候変動，桜，発育速度（DVR）
Key words：blooming period, tourism, climate change, cherry ﬂower,
DVR（developmental rate）

1．はじめに：研究の背景と目的
古くから人びとは，生物季節を観察して季節の移り変わりを知るとともに，生活の中に生物相が変化する情報を取り込み，農作業の暦を決める指標として用いてきた。特にサクラの開花は生活と密着しており，人びとの関心が高いため，こうした観察の対象とされることが多い。地球温暖化によるサクラの季節の変化について，これまで多数の研究が行われ1）ているが，近年では丸岡・伊藤 が，将来の地球温暖化がソメイヨシノ開花に与える影響を予測している。
しかし，これらの研究の多くは開花日に注目したものであり，満開日に関しての予測の研究はほとんど行われていない。本稿では，開花日と満開日の双方に焦点をあて，両者と気温との関係性を解析し，その結果を利用して将来の気候変動が桜の開花期間にどのような変化をもたらすかを明らかにすることを目的とする。
サクラの開花期には，お花見や観光行事が企画される。満開日および開花期間の予測は，観光産業にとって重要な情報となっている。本稿では，満開日ばかりでなくお花見ができる期間が将来どのように変化するかを明らかにする。
2．方法
2.1 開花－満開期間推定モデル2），3）
杉浦らの研究結果を参考にして，サクラの開花から満開を 1 つの生育ステージとし，開花から満開までの期間に DVI（Developmental Index）＝1 になると考える。次に，生育は日々進むことから，1 日ごとの発育速度 DVR（Developmental Rate）を定義する。換言すると，DVR は開花から満開までの期間の気温に応じて発育する割合の指標である。この関係を式（1）に示す。
（B・Ti ）
DVR＝A・exp（1）－1

ここで，DVR：発育速度［day ］，T i：i 日平均気－1
温［℃］，A：定数［day ］，B：定数［℃］である。
推定モデルを適用する地点は，本州・四国・九州の全 52 地点（表 1）で，解析期間は 1985～2009 年の－1

受付；2011 年 10 月 3 日，受理：2012 年 2 月 14 日＊
〒 305-0044 茨城県つくば市並木 3-660，e-mail：waxyleaf@gmail.com
2012 AIRIES31塚原・林：温暖化がサクラの開花期間に及ぼす影響表1 析対象地点一覧解（右図は対象地点の分布を示す）．

1

青森14前橋27名古屋40

徳島

2秋田

15水戸

28

京都41

下関

3盛岡

16熊谷

29静岡

42

松山

4酒田

17福井

30

大島43

福岡

5仙台

18

銚子

31津44

高知

6山形

19東京

32

神戸

45

潮岬

7新潟

20甲府

33奈良

46

佐賀

8福島

21

鳥取

34大阪

47

大分

9輪島

22松江

35

岡山48

熊本

10富山

23舞鶴

36

広島49

長崎

11長野

24横浜

37高松

50

福江

12金沢

25岐阜

38

和歌山

51

宮崎

13宇都宮

26

彦根39

厳原

52

鹿児島

25 年間とした。また，気象庁が作成している日平
均気温と生物季節データのうち，サクラ開花日と満［1］［2］開日を使用した。ここで，対象とするサクラの品種は，各地点ともソメイヨシノである。なお，モデルを構築するにあたり，気温と開花期間の日数に一定の関係が認められた。すなわち，高緯度の地点ほど開花期間が短いという地域に依存した傾向が見られたため，発育速度の概念を適用する前に補正を4）行った。詳細は塚原 に記述した。
2.2 開花－満開期間の将来予測2.1 で構築した開花－満開期間推定モデルと，青5）
野・守屋 が求めた開花日推定モデル B の式（2）および式（3）を組み合わせて将来予測を行った。
青野・守屋のモデルは，温度変換日数（the number of days
［3］
transformed to standard temperature：DTS） を用いており，開花と密接に関係する標準的な温度を基準にして開花日を推定するもので，式（2）および式（3）を用いる。この関係を用いて，将来の気候変動下での開花日の変化を推定した。開花日の指定にモデル B を選定した理由は，次の通りである。まず，モデルそのものの形が複雑ではない。さらに，開花［4］をもたらす休眠効果 の考慮要素が，気候変動とは独立な緯度や海岸からの距離であり，将来予測に適していると判断したためである。
E（Taij－T s）（ts）（2）
ij＝exp
RTijTs｛23.8＝Σj ＝Dme（ts）ijD｝（3）

ここで，（ts）ij；温度変換日数，Tij；i 年 j 日の日平均気温［K］，T s；標準温度［K］（＝288.2 K＝15℃），－1－1（＝70 kJ mol ），R；モル
Ea；温度特性値［J mol ］－1－1
気体定数（＝8.314）［J K mol ］，D；推定開花日である。また，D me は起算日であり，式（4）のように地点ごとの緯度φや海岸からの距離 L などを用い32て求める。
Dme＝136.75－7.689φ＋0.133φ －1.307lnL＋2
0.144TF＋0.285TF（4）
2

ここで，Dme：使用モデルの起算日，φ；緯度［°N］，
L：海岸からの距離［km］，T F：1，2，3 月の平均＊気温の平年値［℃］，また上式の相関係数は r ＝5）0.927 である。（4）式は，青野・守屋 が説明変数として，さまざまな期間について求めた平均気温の平年値，緯度・経度，暖かさの指数，加えて国土数値［5］
情報 による標高や海岸距離や地形開放度などの地形因子などを用いた重回帰分析を行い，説明変数＊ごとに自由度調整済み重相関係数（R ）を考慮して求めたものである。
将来予測に使用する気候変化シナリオは，IPCC［6］が 2007 年に発表した AR4 における温暖化予測デ6）ータを用いた。具体的には，Okada ら が作成した日本付近の将来の気候変化を表すデータセット（気候変化メッシュデータ Ver.2）を使用した。このデー
タセットに採用されている気候変化シナリオは，表2 に示す 6 つの大気海洋結合モデル（CGCM：Atmosphere-Ocean Coupled General Circulation Model）
による合計 8 つの気候変化シナリオである。
本研究では，温室効果ガス排出シナリオとして［7］SRES -A1B と A2 シナリオを想定して考察を行った。これらのデータセットの日平均気温を利用し，累積密度分布（CDF：Cumulative Distribution Function）補正法によって系統的な誤差の補正を行った
のちに，開花日などの将来予測の解析に用いた。影響評価の対象とする時期は，2046～2065 年（今世紀中ごろ），2081～2100 年（今世紀末）の 2 つの期間に分けて行った。

地球環境 Vol.17 No.1

表231－36（2012）

将来予測シナリオで使用された CGCM モデル ．
9）

モデルの名称

数分解度（経度×緯度）

国

大気海洋

3.8 ﾟ× 3.8 ﾟ

1.9 ﾟ× 1.9 ﾟ

カナダ

1.9 ﾟ× 1.9 ﾟ

1.9 ﾟ× 0.8 ﾟ

オーストラリア

CGCM3_1

3

CSIRO_Mk3_01

CSIRO_Mk3_51

GFDL_CM2_01

2.5 ﾟ× 2.0 ﾟ

1.0 ﾟ× 1.0 ﾟ

アメリカ合衆国

MIROC3_2_MEDRES

12.8 ﾟ× 2.8 ﾟ

1.4 ﾟ× 1.4 ﾟ

日本

MRI CGCM2 3 212.8 ﾟ× 2.8 ﾟ

2.5 ﾟ× 2.0 ﾟ

日本

3．結果
3.1 将来開花日
構築した発育速度モデルは，次式のようになった。
（－11.6・Tre ）DVRr＝0.345exp－1（5）

ここで，T は r 日目の平均気温［℃］（地点ごとに補正を行っているので区別するため T re と表記）である。式（5）の標準誤差は 1.85 日である。
式（2），（3）のモデルを利用して，将来の開花日を予測し，現在との偏差の全国空間分布図を図 1 に示す。これを見ると，将来の開花日は北日本などでは早まる傾向にあるが，西南日本では逆に遅くなる傾向が見られた。この傾向は，先行研究の丸岡・伊1）
藤 の結果と一致する。この傾向が現れる理由は，西南日本では，将来の地球温暖化時に冬期の必要な低温量が満たされず，休眠が十分にされなくなるためと考えられる。
3.2 将来開花期間
気候シナリオおよび排出シナリオを違えた全ての組み合わせに対して，式（5）を使い将来の予測を行った結果，今世紀中ごろおよび今世紀末には，地球温暖化に伴い開花から満開までに必要な日数は短くなる結果となった。
推定した開花から満開までの期間（開花期間）の日数を，地点別に現在および 2 つの将来の期間（2046～2065 年と 2081～2100 年）別に平均した結果を図2 に示す。2 つの期間の将来の値は，6 つの CGCM（計 8 つのシナリオ）でほとんど差がなかった。
また，
SRES-A1B 排出シナリオ下における，現在（1985～2009 年）の平均開花期間の日数との偏差割合を図 3に示す。
これらの結果は，将来的にすべてのケースで開花から満開までに要する日数は短くなることを示している。地域的な特徴を見ると，九州地方（地点番号46～52）などは，比較的偏差が小さかった。最も満開までの期間が短くなったのは，SRES-A2 シナリオの 2081～2100 年における松江（地点番号 22）の－4.7日（現在平均期間の 55.3％短くなる）である。
同じ期間での排出シナリオの違いによる日数の差はあまり大きくない。また，同じ排出シナリオの場

図1

現 在の気候条件における平均開花日からの偏差（A1B シナリオの場合）．
（上図：2046～2065 年，下図：2081～2100 年）

合に，今世紀中ごろと今世紀末の予測期間による日数の差も，排出シナリオの違いによる差と同程度である。開花期間が短縮する割合は今世紀末に近づくほど大きくなる傾向にあり，同じ予測時期では，SRES-A2 シナリオでの方が SRES-A1B より大きくなる。

33

塚原・林：温暖化がサクラの開花期間に及ぼす影響

図2

現在と将来の気候条件における解析地点ごとの開花から満開までの期間の日数．

黒実線：SRESa1b（2046 ～ 2065 年）の平均期間日数，黒破線：SRESa1b（2081 ～ 2100 年）の平均期間日数
灰色実線：SRESa2（2046 ～ 2065 年）の平均期間日数，灰色破線：SRESa2（2081 ～ 2100 年）の平均期間日数
赤点線：1985 ～ 2009 年の実測平均期間日数

4．考察
4.1 開花から満開までの期間の短縮
将来の地球温暖化で日平均気温が高くなることによって，開花から満開までの期間は現在よりも短くなることが考えられる。予測結果の現在平均からの偏差割合を示した図 3 を見ると，現在の気候条件における開花期間の 6.4％～55.3％の幅で開花から満開までに要する日数が，将来，短くなることが予測された。また，地域を比較すると，東北ほかの地域で開花期間が短くなる割合が大きく，九州地方などは小さい。この理由として，九州地方の気温の昇温幅が他地域に比べて小さいことが考えられる。開花から満開までの期間の変化は，現在の気候から将来の気温の変化とよく対応しているため，気温変化の地域差が開花から満開までの期間の変化の地域差に影響を及ぼしていると言える。
現在の気候条件では，全国平均して 8.3 日で開花から満開に至るが，2046～2065 年には 2.1 日減の6.2 日で満開に至るようになる。すなわち，それだけ花見ができる日数が減ることになり，市民が楽しむ機会が減ることが予想される。さらに観光産業にとっても，観光資源が減少するという負の影響が現れるかもしれない。
4.2 地域による特徴
例えば宮崎は，現在約 10 日で開花から満開に至るが，SRES-A2 シナリオでの 2081～2100 年の間に8 日強となり，現在の松江と同程度の期間になると予測される。また同じ条件で，現在も 6 日弱と期間の短い北日本の盛岡は約 3 日になり，現在の半分近くに短縮してしまうと予測される。盛岡は現在，5月のゴールデンウィークごろにサクラが咲くことが多く，この時期に「さくら祭り」が開催され，欠か34

図3現在と A1B 排出シナリオの気候条件における開花から満開までの日数の偏差割合（％）．

（上図：2046 ～ 2065 年，下図：2081 ～ 2100 年）
マイナスは期間が短縮することを示す．

地球環境 Vol.17 No.1

せない観光地要素となっている。しかし近年，次第に開花日が早まり，ゴールデンウィークには咲き終わってしまう年が増えていることが指摘されている。
今世紀末ごろの盛岡を例に取ると，将来の地球温暖化に伴う開花日の早まりと期間の短期化により，現在の満開日である 4 月下旬が 4 月上中旬に移動することが予測される。すなわち，ゴールデンウィークまでには桜は咲き終わってしまう可能性が高まると考えられる。
4.3 将来の変化予測
将来の満開日の変化についてもう一つの興味深い点は，日本の中で南北の地域差が不明瞭になる点である。北日本では将来開花日は早まり，西南日本では遅くなる傾向はあるが，西南日本の方が北日本よりも早期に咲くことには変化がない。しかし，このような開花日の変動により，将来の満開日の北日本と西南日本との差は開花日の差ほど明瞭でなくなると予想される。例えば長崎と仙台を比較すると，開花日は長崎が先だが満開日は仙台とほぼ同じになる。このような予測結果が将来の状況を適切に再現しているとすれば，北と南でほぼ同日に満開の桜が見られる状況になるだろう。
次に図 2 に示したように，今世紀の中ごろと末ころにおける開花から満開までの期間の変化の傾向について，今世紀中ごろ（2046～2065 年）には現在との差は大きいが，今世紀中ごろから今世紀末ごろ（2081～2100 年）の差は比較的小さくなった。この要因は気温変動自体やそれに応答する開花日の変動を反映しているためである。すなわち，現在から今世紀中ごろにかけての，開花が予想される期間の気温上昇について，SRES-A1B と A2 のそれぞれについて比較すると，今世紀中ごろから今世紀末ごろまでの気温上昇の割合は，約 50％と小さな規模である。こうした傾向が，全国の開花期間の短縮に大きく影響していると考えられる。

31－36（2012）の期間が他の地域と比較して高い気温の時期にあたるため，補正の効果が大きく現れる可能性が考えられる。今後はこの補正法も課題になると思われる。
2010 年から気象庁は開花日予測の業務を行わないことにした。これに対し，黒岩は，サクラの開花予［8］
想には生物季節の観察が不可欠であると述べた 。
開花予測精度の向上には生物観察が基本であること
7）
は言うまでもない 。このような意見を参考に，サクラの開花・満開の研究がより進歩することを望みたい。
謝辞

本研究を進めるにあたり，独立行政法人果樹研究所の杉浦俊彦博士には，研究初期からさまざまな助言をいただきました。また，大阪府立大学の青野靖之准教授に多くの文献を提供していただき，参考にさせていただきました。厚く御礼申し上げます。さらに筑波大学生命環境科学研究科の岡田将誌氏には，本研究で使用した将来予測の気候データの利用を許可していただくとともに，解析方法についてご教示いただきました。このほか，解析結果について議論していただいた多数の方々に対して，深く感謝いたします。
資料および注［1］気象庁 全国日平均気温データ（CD-RM）.
［2］気象庁 全国桜開花日・満開日データ情報（CD-RM）.
［3］DTS：Konno and Sugihara がアレニウスの法則を7）用いて，生物活性への温度影響の指標として提案
し，「ある温度で，ある日数置かれた条件が標準温度に変換すると，何日に相当するかを表したもの」と定義されている。生物の成長度合いを示すもの
で，この積算した値がある値を超えると開花に至るとされている。
［4］休眠効果：秋季に葉を落とし，翌春まで外見的な生長をしない状態のことであり，動物でいう冬眠

5．今後の課題に類似する。休眠には，自発休眠期と他発休眠期があり，先に迎える自発休眠期には一定量の低温

本研究に用いた予測モデルは，現在のサクラの樹木に固有の気温に対する特性は将来も変わらないことが条件であった。もしサクラの樹木が，その地域の環境，特に気温の上昇に適応して特性を変化させる性質を持つ場合には，環境適応性を含め，生物的・遺伝子的考察を加えたモデルの構築が必要になるだ
ろう。
また，モデル構築の際に行った補正法によって，その地点の気温が高いほど補正が負の方向に働く特徴を持っている。このため，東北地方と比べて将来の気温上昇幅が小さい九州地域などでは，開花期間の偏差も小さくなることが考えられる。特に九州では，現在より開花日が遅くなり，開花から満開までが必要であり，その量に達すると休眠打破となり，その後の他発休眠期では，一定量の高温が必要になってくる。他発休眠を終えると発芽となる。
［5］気象庁 全国観測所地理データ（気象庁 HP）〈http://www.jma.go.jp/jma/index.html〉［6］AR4：気候変動に関する政府間パネル（The Intergovernmental Panel on Climate Change：IPCC）が2007 年に発表した第 4 次評価報告書（Four th Assessment Report：AR4）を言う。
［7］SRES（Special Report on Emission Scenarios）シナリ
オ：IPCC が 2000 年にまとめたもので，温室効果ガスの排出シナリオを大きく，経済成長速度の程度とグローバル化の進展程度の組み合わせによる 435

塚原・林：温暖化がサクラの開花期間に及ぼす影響つのストーリーラインに分類している。本稿で使（2009）Mesh climate change data of Japan, Ver.2 for

用した A1B シナリオは，経済発展が重視されてグ

Climate Change Impact Assessments under IPCC

ローバル化が進む社会であり，化石エネルギー，SRES A1B and A2. Journal of Agricultural Meteorolo-

原子力エネルギー，自然エネルギーをバランス良gy, 65, 97-109.

く利用する社会を想定している。
［8］毎日新聞，2007 年 4 月 15 日（日曜日）朝刊，発言席「桜開花予測に，生物観察も」

7） Konno, T. and S. Sugihara（1986）
Temperature index
for characterizing biological activity in soil and its application to decomposition of soil organic matter.
Bulletin of National Institute of Agroenvironmental

引用文献1） 丸岡知浩・伊藤久徳（2009）わが国のサクラ（ソメイヨシノ）の開花に対する地球温暖化の影響．農業気象，65，283-296.
2） 杉浦俊彦・小野祐幸・鴨田福也・朝倉利員・奥野隆・浅野聖子（1991）ニホンナシの自発休眠覚醒期から開花期までの発育速度モデルについて．農業気象，46，197-203.
3） 杉浦俊彦・本條 均（1997）ニホンナシの自発休眠覚
醒と温度の関係解明およびそのモデル化．農業気象，53，285-290.

Sciences, 1, 51-68.

塚原あずみ

Azumi TSUKAHARA

東京都出身。筑波大学第一学群自然学類在籍時に，気候学・気象学を専攻。環境問題に興味を持ち，生物気象について研究し，主にサクラを用いた気候変動の影響評価を行った。現
在は，観光業界を経て，環境問題や自然について専門知識のない人びとに伝える橋渡し役になるべく，環境カウンセラーやネイチャーゲームリーダーなど環境教育分野の NPO に所属し活動している。

林陽生4） 塚原あずみ（2010）桜の開花から満開までの期間とYousay HAYASHI

気温と関係性及び将来予測，筑波大学第 1 学群自神奈川県出身。独立行政法人農業環境技術研究所地球環境部長を経て，国立大学法人筑波大学生命環境科学研究科教授。
専門は気候学・気象学および農業気象学。地球温暖化の農業生産への影響の評価および植物群落内外のエネルギー収支などを研究テーマとしている。著書に『地球温暖化で日本の農業はどう変わる』（家の光協会），『日本の気候Ⅱ』（二宮書店，共著）などがある。

然学類平成 21 年卒業論文，1-64.
5） 青野靖之・守屋千晶（2003）休眠を考慮したソメイ
ヨシノの開花日推定モデルの一般化．農業気象，
59，167-177.
6） Okada, M., T. Iizumi, M. Nishimori and M. Yokozawa36