収量予測
2020年収量結果
８月は晴れの日が続いたのですが、９月に入ってから集中豪雨が頻繁に襲ってきました。そのため、刈取りは天気予報を見ながら雨の降らない合間を狙って急遽実施しました。
今まで使っていたコンバインは、令和元年東日本台風（台風19号）による水没によって廃車となり、今年からは新型のコンバインになります。初運転となるので、ヤンマーの担当者さんから操作のレクチャーを受けました。また、実際に刈り取る圃場でも丁寧に教えていただきました。担当者さんも言っていたのですが、梅雨の長雨による不十分な中干し、刈取り前の豪雨で土壌を乾燥することができないといった・・・コンバインの操縦にはハードな土壌状態での刈取りとなりました。
新しいコンバインでの刈取り
今年の出穂期までの積算日照時間は、下図のようにモニタリングを始めてから最も低い値を示しています。そのため、日照不足であった2019年の収量予測パラメータ（日照不足用）を用いて収量を予測しました。
（アメダス：鳩山地点を用いて作成）
2014～2020年の移植日から収穫日まで積算日照時間
それでは、１ヶ月前に投稿した収量予測（予測値：1464 kg、474kg/10a）の答え合わせです。
試験サイト全体の精玄米収量 ：1491 kg（483kg / 10a）
収量予測との誤差 約27kg （1.8％）
今年は、ほぼほぼピッタリの予測ができました（ニアピン賞ぐらい）。また、過去の結果から幅を持たせて予測した1464～1512kg（474～489kg/10a）の間には収まりました。ドローン水稲モニタリングを始めてから７年目になるので、いろいろな知見から対応できつつあります。
ドローン運用開始からのコシヒカリ収量（10a当たりの精玄米収量）およびタンパク質含有率の結果
今年はドローン水稲モニタリングを始めて歴代２位の収量となりました。日照時間から見ると2020年は2019年と同様だったので、10a当たり465kgの収量でもおかしくなかったのですが、農閑期に行った土づくりの見直しが功を奏したのかも知れません。
出穂期・収量予測(2020年)
今年も昨年と同様に梅雨の長雨で梅雨明けが、平年より11日（昨年より８日）遅い８月１日となりました。梅雨明けが遅かったので、試験サイトでも２日遅い（５年間の平均）８月５日に出穂期を迎えました。
イネの開花（2020年８月５日撮影）
2014年からの水稲モニタリングの記録から、2020年は2016、2019年の積算日照時間・積算温度に近い傾向となっています。特に2020年の積算日照時間は過去のデータと比べると短く、最も長かった2018年の58％となっています。2020年の積算温度は1775.8℃と平均に近い値を示しました。これらの結果から、この地域におけるコシヒカリの出穂期までの積算温度の目安は1700℃ぐらいだと考えられます。
| 年 | 移植日 | 出穂期 | 日数 | 積算日照時間（h） | 積算温度（℃） |
| 2020 | |||||
| 2019 | |||||
| 2017 | ５月21日 | 7月30日 | 71 | 399.8 | 1697.6 |
| 2016 | ５月21日 | ８月４日 | 76 | 381.3 | 1774.5 |
| 2015 | ５月23日 | ８月３日 | 73 | 427.8 | 1742.1 |
| 2014 | ５月24日 | ８月３日 | 72 | 400.4 | 1714.0 |
| Average | ５月21日 | ８月１日 | 72.0 | 387.2 | 1736.1 |
気象データは最寄りのアメダス地点：鳩山を用いて算出
出穂期のモニタリングも終えたので、2020年収量予測マップを作成してみました。今年は2019年に近似しているので、2019年の収量パラメータを使ってみたところ、以下のような結果となりました。
試験サイト全体の玄米収量予測 （ドローン）： 1464 kg
10aあたり
速報値：2019年収量
連休を使って、稲刈り＆籾摺りを行いました。天気予報では雨となっていたので、刈り取りができるか心配でしたが、なんとか無事に収穫できました。
雨が降りそうな天気での稲刈り
さて、１ヶ月前に投稿した収量予測の答え合わせです。収量予測では、前年並みの収量と予測したのですが、結果は…以下の通りになりました。
試験サイト全体の精玄米収量 ：1438.2 kg（465kg / 10a）
収量予測との誤差 約64kg (4.3％)
2018年と比べると減収になりました。やはり、幼穗形成期の日照不足が大きく影響していると考えられます。ちなみに、2018年（猛暑パラメータ）で計算すると、精玄米収量：1610 ㎏と大きく外れてしまいます。今年は2016年の天候と傾向が似ていたことから、2016年パラメータを使用しましたが、さまざまな天候に対応するにはデータの蓄積が必要です。ドローンモニタリングを2014年から始めていますが、まだ５年分のデータしか蓄積ができていないので、予測モデルの精度が満足できる状態になるまではもう少し時間がかかりそうです。
ドローン運用開始からのコシヒカリ収量（10a当たりの精玄米収量）およびタンパク質含有率の結果
これまで増収を続けていましたが、2019年で記録はストップしてしまいました。来年は増収できるように栽培方法をちょっと見直してみたいと思います。
収量予測(2019年)
2019年は梅雨の影響で、分げつ期〜幼穂形成期〜穂ばらみ期の日
（アメダス：鳩山地点を用いて作成）
移植日からの積算日照時間（2014～2019）
2016年の収量パラメータを使って、2019年の
試験サイト全体の玄米収量予測 （ドローン）： 1502 kg
10aあたり
結果はもうすぐです。
日射量
千葉大の濱氏は、ドローンを用いたモニタリングデータと日射量を
2018年の収量は、日射量を組み合わせた濱モデルが精度良く推定することができました。詳しくはこちら。
その研究の一環として、2017年から自宅の屋上に全天日射計を
影の影響を受けない屋上の南側に設置
日平均の日射量（2017～2018年）
オレンジ色:日単位に集計した埼玉県坂戸市の全天日射量（20
紺色:東京の全天日射量（平年値）
2017年の６月は空梅雨だったので、日射量が高
次に、2018年は早すぎる梅雨明けで７～８月は猛暑が続きました。梅雨の期間は日射量が少なくなるのですが、2018年は日射量が多い結果を示してい
時間解像度は2017年の観測は１時間単位で行ってきましたが、電池・記憶媒体の余裕もあることから2018年からは
平成最後の収穫
2018年の栽培期における天候は、一言で言うと「猛暑が続いた年」でした。８月末から９月上旬になると、台風20号（Cimaron）の接近・ゲリラ豪雨によって、この地域は瞬間的に強風・大雨が襲いました。その影響で、圃場内の一部のイネが倒伏程度３（出穂期２週間前の倒伏リスク診断でリスクがあった箇所）になりましたが、この程度であればコンバインによる作業には全く問題はありません。
下図は、2014～2018年の移植日から収穫日までの積算日照時間になります。図からわかるように今年の天候は例年と比べて日照時間が大きく外れています。本当に暑かったです。
（アメダス：鳩山地点を用いて作成）
2014～2018年の移植日から収穫日まで積算日照時間
収穫時における2018年の積算日照時間は例年の約1.3倍となりました。ちなみに、2017年も出穂期までは2018年と同様な空梅雨・猛暑でしたが、８月に入ると途端に冷夏となり日照不足となっています。
さて、４年間の水稲モニタリングで得たNDVIによる収量予測の結果では、今年の収量予測は2017年（471kg/10a）と比べると減収になりました。収量予測は出穂期のNDVIと前年までに得た地上サンプルとの相関式で求めます。４年間分のパラメータがあるのですが、今年はどのパラメータを用いても減収予測でした。その理由は、2018年の出穂期のNDVIが前年比85％と低かったためです。
NDVIのみ
収量予測：462kg/10a
共同研究者でもある濱ほか（2018）の論文では、NDVIと日射量（コシヒカリの場合は出穂期から20日間）を用いて収量を予測するモデルを示しています。そのモデルによる収量予測は増収となります。
濱モデル
収量予測：488kg/10a
今年の収穫結果は・・・ 498kg/10a で増収となりました。農家としては嬉しい結果ですが、研究者の立場からだと外れてしまった・・・という感じです。また、今年の収量は世代交代してから初の反収あたり８俵を超えることができました。この地域のコシヒカリの目標は８俵（JAいるま野広報誌2018年３月号 pp.５）なので、満足な結果です。
ドローン導入からのコシヒカリ収量（10a当たりの精玄米重量）およびタンパク質含有率の結果
ドローン導入から５年目で30%の増収になりました。籾摺り → 袋詰め → 冷蔵管理までの一連の作業は終わりましたので、まもなく2018年度産「どろーん米」をネットショッピングで販売開始いたします。
意外な結果!?
収穫の時期になると集落のあちらこちらから「ピーピーピー」と音が響いてきます。音の正体はコンバインの動作音なので、この音を聞くと稲刈りに向けて気持ちがソワソワしてしまいます。
先日、2017年産「どろーん米」の収穫が何とか終了しました。台風18号（TALIM）が来る前に刈り取れたので、一安心です。
こうべを垂れる「どろーん米」
しかし、今年は使用しているコンバインの調子がすこぶる悪く、JAさんに何回も修理を依頼したりと予定通りにできませんでした。修理をお願いしたら、すぐに対応してくれるJAさんには感謝です。ありがとうございます。
刈取りは、「コンバイン → 乾燥機 → 籾摺り → 袋詰め」 が一連の流れとなっているので、この工程でどれかひとつの農機具が故障してしまうと、作業全体がストップしてしまいます。故障が長引いてしまうと、稲刈りの適期を逃してしまう可能性があります。現代の農業は如何に機械に依存しているか身をもって実感しました。
さて、2017年の収量ですが、事前予測では８月の日照不足のため、３反の試験サイト圃場の収量は
試験サイト全体の精玄米収量 ： 1440.0 kg （465kg / 10a）
※未熟米 ：110.2 kg
日照不足の影響はそんなに大きくなかったのか、2016年の精玄米収量 : 1437.5kg とほぼ同量の結果となりました。事前予測で求めた日照不足の影響を考慮に入れていない推定値 1474 kg（粗玄米）の方が良い結果を示しています。
（2017年度の精玄米収量＋未熟米） / 2016年度の収量予測式（粗玄米） で求めると誤差が5.2％となりました。
日照時間が多く、気温が高い状態が続いた2016年のパラメータが冷夏の場合にも使えるとは思っていなかったので、ますます水稲モニタリングの難しさを感じます。
未熟米に分類された玄米
※2017年は日照不足の影響によって未熟米が例年より多くなりました。そこで、未熟米と分類された玄米を再度ふるいにかけたところ、110.2kgのうち14.1kgが精玄米として再分類されたので、その分は「2ndチャンス米」として、自宅で消費します。
収量予測・・・減収の見込み
８月に入ってから、オホーツク海高気圧からの「やませ」が関東地方まで流れ込んでいるため、埼玉県も日照不足となっています。東北の太平洋側（岩手県、宮城県、福島県）では、平年を大きく下回っていることから、いもち病の心配があるそうです。
試験サイト周辺の日照不足を定量的に見るために、最寄りのアメダス地点：鳩山の日照時間を８月１日～20日までまとめてみました。今年の８月１日～20日間の日照時間は36.1時間と平年の34.5％と大きく下回っています。出穂期からの日照時間は収量・食味に大きく影響します。気象庁によると８月下旬からは平年並みに戻る見込みだそうなので、晴れることを祈ります。
８月１日～20日までの積算日照時間（2014～2017）
今年も収量予測をしてみました。使用するのは７月30日の出穂期のデータになります。ただし、生育が順調に進んだ出穂期のデータなので、それ以降の日照不足を反映していません。そのため、ここで推定する値は日照不足がなかった場合の値になります。ちなみに、平年並みの日照時間があった2016年の収量結果はこちらから閲覧できます。
試験サイト全体の玄米収量予測 （ドローン）： 1474 kg
８月の日照不足を考慮にいれると、この求めた推定値（玄米収量1474kg）の約20%の減収（玄米収量1180kg）になると考えています。
今年の収量から解析して得られる「単位面積あたりの収量とNDVIの相関式」は冷夏用のパラメータとして、今後の栽培に活かせるはずです。農業技術が進んでも、天候次第で収量・品質が大きく左右されるのは昔から変わりません。
収量計算確定（2016年）
先日掲載した地上サンプリングからメッシュごとに単位面積あたりの玄米重量（g/㎡）を求めました。
単位面積当たりの玄米重量は、メッシュ内の株数・１株当たり平均茎数・１穂あたりの玄米重量から計算できます。
求めた玄米重量とNDVI（出穂期）の回帰分析から線形回帰式を求めます。
以下の式は、収量予測に用いた推定式になります。
【2014年】
Y = 968.42 × NDVI – 33.50
【2016年】
Y = 1822.2 × NDVI – 341.43
Y：単位面積あたりの玄米重量（g/㎡）
2016年推定式で計算した収量分布が下図になります。
収量マップ（2016年）
８月中旬に予測したマップと比較すると，分布傾向は似ていますが，多くのメッシュで数量が異なりました。
収穫１カ月前に予測した収量マップ（2016年）
※収穫１カ月前に予測した収量マップ（2016年）は、２０１５年の推定式を用いて作成しています。
収量検証（2016年）
計算で求めた収量：1461 kg
実際の収量：1437.5 kg
この結果から、2016年の推定式は1.6%の誤差で推定玄米収量を求めることができました。
2014年は5.3%の誤差がありましたが、年々精度は向上しています。
収量結果（2016年）
2016年の収量は10a当たり465 kgとなりました。圃場全体では、1440 kg（屑米を除く）です。
ドローン水稲モニタリングを始めて３年目になりますが、１年目と比べると収量は約20%の増加になりました。また、お米の美味しさの指標となる玄米タンパク質含有率も2015年では約6.1%と日本人が好むやわらかいご飯になっています。2016年はこれから分析です。
ちなみに、一般的な玄米のタンパク質含有率は6.8%とされており、タンパク質含有率が低いほどやわらかいご飯となり、数値が高いとしっかりとした硬いご飯となります。
ドローン水稲モニタリングの成果（2014－2016）
ドローン水稲モニタリングの導入１年目は、これまで行っていた「勘と経験」の水稲栽培の問題点を洗い出しました。２年目以降は浮かび上がった問題点を改善するような栽培を行うことによって、収量・食味の向上に結び付けることができました。収量・食味UPにつなげるためには、単年だけではなく、複数年のモニタリングが必要です。
収穫前に予測した収量(８月19日)は、地上観測による予測1520kg、ドローンによる予測1570kgでした。地上観測の予測値が実測と近い結果となりましたが、両者とも約10%の誤差の範囲内に収まっています。
ドローンによる予測は、2015年のパラメータをそのまま利用したのが誤差を大きくした要因ではないかと考えられます。
今年は株間（16cm→18cm）を変更した影響も考えられるので、様々な状況下のパラメータを取得することによって、今後の予測精度の向上を図りたいと思います。
農作業も一段落した農閑期は、2016年のモニタリングで得たデータを詳細に解析する期間です。