ブドウ樹の積算温度指標、温暖化による生育促進を遅延と誤認させる可能性

新たな統計研究によると、ブドウ樹の研究で一般的に使われる気象条件の測定方法は、誤解を招く結果をもたらすことがある。実際には生育を促進するはずの温暖な条件が、植物の生育を遅らせるように見えてしまう場合もある。

学術誌Discover Dataに掲載されたこの研究は、開花や果実の成熟といった段階の時期を指す、ブドウ樹のフェノロジー（生物季節）のモデル化手法を検証した。著者のS. Behnamian氏とF. Fogh氏は、多くの分析では、現象が観察された日までの累積熱量を算出することで、モデルに隠れたバイアスが持ち込まれていると指摘する。季節の進行に伴って熱量は自然に蓄積するため、観察日が遅いほど累積熱量も大きくなるのがほぼ常だ。その結果、説明に用いるはずの気象変数が、説明対象である結果によって部分的に決まってしまうという統計上の問題が生じる。

著者らは、この問題によって気温の影響の見かけ上の方向が逆転することがあると述べた。具体的には、実際の生物学的影響は逆であるにもかかわらず、モデルが気温上昇によって開花や成熟が遅れると示す可能性がある。

著者らが提案する解決策は単純だ。現象が起きた日までの気象条件を合計するのではなく、現象を検知できると考えられる時期より前に終了する、あらかじめ固定した期間内だけで気象条件を算出する。ブドウ樹の事例では、開花については1月1日から4月30日まで、成熟については1月1日から6月29日までを対象期間とした。通常の現象発生日より前に期間が終了するため、気象指標が観察された現象の時期と機械的に結び付くことはない。

研究では、まずシミュレーションを通じて問題を検証した。1,000回の反復計算で、研究者らはシーズン前の熱量が及ぼす真の影響を既知としたフェノロジーの発生時期を生成した。この設定では、熱量が標準偏差1つ分増えると現象の発生が15％早まるはずであり、時間比では0.85となる。正しい外生的な気象指標をモデルに用いた場合、この結果はほぼ正確に再現され、推定時間比の平均は0.851だった。一方、現象が起きる日までの累積熱量を使うと、平均推定値は1.237にずれた。これは45％の上方バイアスに相当し、促進ではなく遅延を示す、誤った方向の結果だった。

次に著者らは、USA National Phenology Networkの栽培ヨーロッパブドウ（Vitis vinifera）の観察記録にこの手法を適用した。記録は2012年から2025年までを対象とし、米国内5地点にある7本の個体について、16人の観察者が提出した41,274件の状態記録を含んでいた。対象を開花と成熟果実の2段階に絞った結果、分析に使用できる記録は8,274件となり、これらを植物・地点・年ごとの50観察単位に整理した。

研究では、野外観察に伴う基本的な問題に対応するため、生存時間分析を用いた。観察者が毎日植物を訪れることは通常ないため、ブドウ樹が段階を移行した正確な日付は分からないことが多い。分かるのは、前回の訪問時には開花や成熟が確認されず、初めて記録された次の訪問までの間に起きたということだけだ。著者らは、各現象を単一の観察日に無理に当てはめず、区間打ち切りデータとして扱った。また、現象が一度も記録されなかった年も除外せず、右打ち切り観察として扱った。

これは、より単純な手法では、遅い現象や未観察の現象について有用な情報が含まれているにもかかわらず、そうした年を除外することが多いため重要だ。研究によると、この生存時間分析の枠組みは、観察記録に基づくフェノロジーに伴う不確実性をより忠実に捉えられる。

ブドウ樹のデータでは、補正後のモデルで開花よりも成熟の方にはっきりとした傾向が示された。シーズン前の気温が高いほど、一貫して成熟が早まった。主要なワイブルモデルでは時間比は0.864で、信頼区間は0.813から0.917だった。著者らによると、このデータセットでは、シーズン前の熱量が標準偏差1つ分（約150生育度日）増えるごとに、成熟日の中央値が約29日早まった。対数ロジスティック分布に基づく第2のモデルでは、時間比は0.962と影響が小さくなった。方向は同じだったものの、推定値の大きさはモデルの設定にも左右されることが示された。

開花については異なる結果が得られた。シーズン前の気温上昇に対する推定反応は弱く、統計的にも不確実で、すべてのモデルで信頼区間が1をまたいでいた。また、このサンプル規模では、降水量の明確な影響も確認されなかった。基準蒸発散量を差し引いた水収支指標に累積降水量を置き換えると、わずかに見られた降雨の影響も消えた。

補正後のモデルと、問題のあるモデルとの対比は、論文の中心的な論点だった。研究者らが意図的に、現象の発生日まで累積した気象変数を使うと、気温の上昇によって開花も成熟も遅れるという結果がモデルから示された。シミュレーション検証によれば、これは生物学的なシグナルではなく、統計上のアーティファクトだった。

著者らは、ブドウ樹の事例を、あらゆるブドウ園の気候応答について一般化する主張ではなく、手法を示すための例と位置付けた。サンプルは小規模で、偏りもあった。当初の記録数は多かったものの、対象は5地点のわずか7本であり、固定期間モデルの分析に寄与したのはカリフォルニア州とニューヨーク州の3地点だけだった。50の植物・地点・年単位のうち43では品種情報が欠けていた。また、ワイン用、食用、観賞用のブドウが混在し、大半は灌漑されていたため、生育に影響しうる栽培管理の違いも加わっていた。

論文では、一連の感度分析も行われた。研究者らは、生育度日を算出する際の基準温度、対象期間の長さ、終点のより厳格な定義を変えたほか、より均質な単一地点のサブセットも検証した。こうした変更によって点推定値は動いたものの、主な結論は変わらなかったと報告している。また、打ち切りが情報を含むかどうかも調べたところ、観察間隔の中央値は4日だった。地点単位のクラスター化は成熟について中程度で、クラスターを考慮したモデルでは標準誤差がおよそ2倍になる可能性があったが、全体的な解釈は変わらなかった。

この手法を他の種や観測ネットワークに適用するには、一定の条件を満たす必要があると研究は指摘する。観察頻度は、年ごとの現象発生時期のばらつきに対して十分でなければならない。右打ち切りが過度に多く、パラメトリック推定が不安定になるような状況も避ける必要がある。選択したモデルに十分なサンプル数が必要であり、固定する気象期間は、現象が検知される可能性が最も早い日より前に終了しなければならない。そうでなければ、同じバイアスが再び生じる可能性がある。

飲料業界にとって、この研究は、将来的にブドウ園の作業や収穫時期の予測精度向上につながれば重要な意味を持つ可能性がある。ブドウの開花や成熟時期をより正確に推定できれば、ワイン生産における人員計画や果実の出荷時期の調整に役立ちうる。ただし、今回のサンプルは限られており、実際の運用上の価値が明らかになるには、より幅広い商業栽培の条件で検証する必要がある。

この論文は、気候応答研究におけるより広い問題も示している。フェノロジー研究では、積算温度などの累積指標が、簡便で広く用いられているため頻繁に使われる。しかし、予測対象となる現象との関係でこうした指標を定義すると、回帰モデルや生存時間モデルの基本的な前提に反する可能性があると、本研究は示唆している。ブドウ樹の場合、温暖化で生育が早まる反応が、遅れとして見えてしまうほどの影響が生じうる。

著者らによると、シーズン前の期間を固定し、区間打ち切りを考慮した生存時間モデルを使うことで、この落とし穴を避け、実際の野外観察の収集方法により即した気候シグナルを取り出せる。