豊葦原中津谷のニニギ

食糧自給率の向上を目指して! (2100年の日本へワープ)

カテゴリ:気象予測 > サクラ開花予想

桜の開花時期を予想するためには、桜の開花を決める法則を探す必要がありました。
それについては、三つの候補を見つけました。
(他にもいろいろありそうですが、今は三つで充分です)
 
次には、未来の気温の変化を予測する必要があります。
これについても、桜開花予想の記事の中でも、一度触れています。
今回は、もう少し踏み込んでみましょう。
 
以前の記事では、二次関数で気温の変化を回帰したい旨を書いています。
ですが、気温の変化は、本来なら正弦曲線を描くはずですから、
二次関数で描ける範囲は、変曲点から次の変曲点までの期間に限られます。
そこで、
7月1日から3月31日までを三次式で回帰し、変曲点を探してみました。

東京1953年回帰

1952年7月1日から翌年3月31日までを三次関数で回帰したグラフです。
係数を基に、変曲点を調べたところ、1952年11月1日でした。

東京2013年回帰

2012年7月1日から翌年3月31日までを三次関数で回帰したグラフです。
係数を基に、変曲点を調べたところ、こちらも2012年11月1日でした。
東京では、11月1日頃に変曲点がありそうです。
 
次回は、札幌の気象データで調べてみることにします。

前回にも書いたように、
本来なら正弦曲線になる気温の変化を二次関数で回帰することを考えています。
しかし、二次関数で描ける範囲は、変曲点から次の変曲点までの期間に限られます。
そこで、
前回は、東京の日平均気温について
7月1日から3月31日までを三次式で回帰し、変曲点を探してみました。
 今回は、札幌の日平均気温で変曲点を探してみました。

札幌1953年回帰


1952年7月1日から翌年3月31日までを三次関数で回帰したグラフです。
係数を基に、変曲点を調べたところ、1952年10月28日でした。


札幌2013年回帰



2012年7月1日から翌年3月31日までを三次関数で回帰したグラフです。
係数を基に、変曲点を調べたところ、こちらも2012年11月5日でした。
札幌では、11月1日頃を中心に変曲点がありそうです。
 
 
東京と札幌の例を見ると、
変曲点を避けるには、11月上旬は回帰式のデータに含めない方が良さそうです。
桜の開花予想は、3月1日に発表したいので、
使用できる気象庁のデータは、2月26日(来年は27日)までとなります。
つまり、
前年の11月11日から当年の2月26日までの気象庁のデータを用い、
桜の開花時期の予測を行うことになります。
 
次回は、実際に過去のデータで検証してみることにします。

桜の開花予想、忘れていません!
昨日は、桜の開花予想からちょっと外れていましたが、元に戻ってきました。
 
前回は、気温の変化を二次関数で回帰するために、変曲点を探しました。
その結果を踏まえ、
11月11日から2月26日までのデータを基に回帰式を求めることを決めました。
今回は、実際に過去のデータを基に回帰式を求めてみました。
使用したデータは、11月11日から2月26日までの東京の日平均気温です。
この期間のデータを1953年から2013年までの61年分について
回帰式を計算してみました。
その計算結果から、4730Kの法則に基づいて7.8℃になる日付を推定し、
その16日後を桜の開花日として算出してみました。
 
その結果は、全く芳しくないものでした。
一例を上げると、
1988年の桜の開花日の推定値は、なんと9月5日となりました。
実績は、4月2日でしたから、5ヶ月もずれています。
(オーストラリアの桜開花日を予測してしまった?
 
その原因を探るため、
気温の実測値と回帰式で推定した気温をグラフにしてみました。
 
まずは、上手くいかなかった例です。

1988年桜開花予想(東京)

次は、上手くいった方です。

1963年桜開花予想(東京)


こちらは、実際の開花日を推定することができました。
 
 
上手くいく場合もありますが、
気温の推定方法を検討しないと、桜の開花予想はできないと思っています。

昨日は、東京の過去の気温を用いて、
開花時期を推定するための7.8℃となる日付を推定してみました。
しかし、
その結果は、惨憺たるものでした。
気温変化の推定が、誤差という言葉では済まされない値が出てしまいました。
 
では、札幌ではどうなるのか?
実は、札幌で桜が開花する5月中旬頃は、変曲点を超えた頃になります。
なので、札幌は、東京以上にズレを生じやすいのです。
でも、
念のため、試してみました。

1978年桜開花予想(札幌)


予想した通り、札幌は東京よりもズレが大きく現れました。
上のグラフは、1978年春の計算結果です。
これから予想された1978年の札幌の桜開花は、12月14日でした。
笑うしかありません。
実際の開花は、5月10日でした。
 
 
さて、対策ですが、いくつかのアイデアがあります。
ですが、
その前に解決しなければならない問題があります。
私は、EXCELシートとマクロの組み合わせで計算しています。
61年分のデータを一気に処理するために、
上記の組み合わせを繰り返して計算しています。
この時、EXCELシートの計算がマクロの速さに追いつけないのです。
この問題も解決しないと、検証に時間が掛かりすぎて話になりません。
 
色々、障害だけはたくさんあります。
 

ウェザーマップから、早くも桜開花予想が発表されました。
私の「サクラ開花予想」は、今月末の発表予定なのに!

ところで、桜には、「休眠打破」と呼ばれる性質があると言われています。つまり、ある期間以上に寒さに晒されないと開花しない性質です。
私は、7.8℃より下がらないと開花しないと考えています。
鹿児島県では、一月下旬の平均気温でも7.9℃もあり、真っ当に開花するか、ギリギリのところです。

ウェザーマップでも、充分に気温が下がっていないため、鹿児島県や熊本県の開花時期は遅れると判断しています。


さて、私の「サクラ開花予想」ですが、現時点までに発表した内容では、2月下旬以降でなければ開花時期を予想することはできません。
この対策は、次々回で発表する予定です。
ただ、「お知らせ」にも掲載しましたように
、ちょっと面倒な式を解く必要があり、時間がかかりそうです。
なんとか、2月末日に間に合うように頑張ってみます。

乞うご期待!

↑このページのトップヘ