(57)【特許請求の範囲】
【請求項１】
今期の初期の緩効性肥料の施肥後に行われる追肥時の施肥量である追肥時施肥量を導出する施肥量導出装置であって、
今期の初期の施肥後から追肥前に至るまでの積算温度を取得し、取得した積算温度に基づいて、追肥時における緩効性肥料の残存量である今期肥料残存量または平年と今期との緩効性肥料の残存量の差である残存量差を導出する残存量導出部と、
前記残存量導出部により導出された前記今期肥料残存量または前記残存量差を反映して前記追肥時施肥量を導出する施肥量導出部と、
を備えることを特徴とする施肥量導出装置。
【請求項２】
追肥前の農作物の植生に関する植生関連値を取得し、植生関連値を入力の１つとし施肥量を出力とする施肥量導出計算式を利用して、取得した植生関連値に基づいて仮の追肥時施肥量を導出する仮施肥量導出部を更に備え、
前記残存量導出部は、今期の初期の施肥後から追肥前に至るまでの積算温度に基づいて養分溶出率または養分残存率を求め、養分溶出率または養分残存率から前記今期肥料残存量を導出し、
前記施肥量導出部は、前記仮施肥量導出部により導出された前記仮の追肥時施肥量から前記残存量導出部により導出された前記今期肥料残存量を減算して最終的な追肥時施肥量を導出する
ことを特徴とする請求項１に記載の施肥量導出装置。
【請求項３】
前記施肥量導出部は、前記仮施肥量導出部により導出された前記仮の追肥時施肥量から前記残存量導出部により導出された前記今期肥料残存量を減算した値が、予め定められた上限値を上回る場合には、上限値を前記最終的な追肥時施肥量として導出することを特徴とする請求項２に記載の施肥量導出装置。
【請求項４】
追肥後から収穫前に至るまでの地力窒素の発現量の予測値である地力窒素予測発現量を取得する地力窒素予測発現量取得部を更に備え、
前記施肥量導出部は、前記仮施肥量導出部により導出された前記仮の追肥時施肥量から前記残存量導出部により導出された前記今期肥料残存量を減算した値に、更に前記地力窒素予測発現量取得部により取得された前記地力窒素予測発現量を反映して前記最終的な追肥時施肥量を導出することを特徴とする請求項２に記載の施肥量導出装置。
【請求項５】
前記施肥量導出部は、前記仮施肥量導出部により導出された前記仮の追肥時施肥量から前記残存量導出部により導出された緩効性肥料の残存量を減算した値に、更に前記地力窒素予測発現量取得部により取得された前記地力窒素予測発現量を反映して導出した値が、予め定められた上限値を上回る場合には、上限値を前記最終的な追肥時施肥量として導出することを特徴とする請求項４に記載の施肥量導出装置。
【請求項６】
前記仮施肥量導出部は、蓄積された過去の実際の追肥時の施肥量および追肥前の植生関連値の組み合わせと農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および追肥前の植生関連値を説明変数とする重回帰分析を行って求められた重回帰式に基づく計算式を前記施肥量導出計算式として利用して前記仮の追肥時施肥量を導出することを特徴とする請求項２から５の何れか１項に記載の施肥量導出装置。
【請求項７】
前記仮施肥量導出部は、農作物の群落に対する遠隔からの観測により得られる第１指標値を前記植生関連値として取得することを特徴とする請求項２から６の何れか１項に記載の施肥量導出装置。
【請求項８】
前記仮施肥量導出部は、ドローンを用いた上空からのセンシングにより得られるＮＤＶＩ、または、衛星写真に基づくＮＤＶＩを前記第１指標値として取得することを特徴とする請求項７に記載の施肥量導出装置。
【請求項９】
前記仮施肥量導出部は、ＮＤＶＩセンサを用いた地上での測定により得られるＮＤＶＩ、葉緑素計を用いた地上での測定により得られるＳＰＡＤ値、または、葉色板を用いた地上での測定により得られる葉色板値を前記植生関連値として取得することを特徴とする請求項２から６の何れか１項に記載の施肥量導出装置。
【請求項１０】
前記残存量導出部は、積算温度に基づいて前記今期肥料残存量を導出し、平年の追肥時における緩効性肥料の残存量である平年肥料残存量と、導出した前記今期肥料残存量との差を求めることによって前記残存量差を導出し、
前記施肥量導出部は、前記残存量導出部により導出された前記残存量差を慣行の追肥時の施肥量に反映させることによって前記追肥時施肥量を導出する
ことを特徴とする請求項１に記載の施肥量導出装置。
【請求項１１】
前記残存量導出部は、平年の積算温度と今期の積算温度との差に基づいて前記残存量差を導出し、
前記施肥量導出部は、前記残存量導出部により導出された前記残存量差を慣行の追肥時の施肥量に反映させることによって前記追肥時施肥量を導出する
ことを特徴とする請求項１に記載の施肥量導出装置。
【請求項１２】
前記残存量導出部に代えて、追肥後から収穫前に至るまでの地力窒素の発現量の予測値である地力窒素予測発現量を取得する地力窒素予測発現量取得部を有し、
前記施肥量導出部は、今期の初期に速効性肥料が施されている場合、前記地力窒素予測発現量取得部により取得された前記地力窒素予測発現量を反映して前記追肥時施肥量を導出する
ことを特徴とする請求項１に記載の施肥量導出装置。
【請求項１３】
今期の初期の緩効性肥料の施肥後に行われる追肥時の施肥量である追肥時施肥量を導出する施肥量導出方法であって、
今期の初期の施肥後から追肥前に至るまでの積算温度を取得し、取得した積算温度に基づいて、追肥時における緩効性肥料の残存量である今期肥料残存量または平年と今期との緩効性肥料の残存量の差である残存量差を導出するステップと、
導出した前記今期肥料残存量または前記残存量差を反映して前記追肥時施肥量を導出するステップと、
を含むことを特徴とする施肥量導出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、農作物についての追肥時の施肥量を導出する施肥量導出装置および施肥量導出方法に関する。
【背景技術】
【０００２】
従来、農作物の栽培に際して行われる肥料の供給に関して、肥料の量に関する有益な情報を提供する技術が存在する。例えば、特許文献１には、今期生育中の稲（水稲）の葉色および茎数を測定し、葉色および茎数から吸収窒素量を求め、吸収窒素量から既に投入した窒素施用量を減算して地力窒素量を求め、地力窒素量と事前に設定された適正窒素量から次期の稲作に必要な窒素施用量（施肥量）を求める技術が記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１８−８２６４８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで稲作、その他の農作物の栽培では、１つの期の中で緩効性肥料を施肥した後に追肥が行われる場合がある。例えば稲作では、緩効性肥料の施肥後、幼穂の形成時期から出穂時期までの所定のタイミングで穂肥と呼ばれる追肥が行われる場合がある。この追肥時に供給される肥料の量は、農作物の品質および収量に大きな影響を与えるため、１つの期の中で緩効性肥料の施肥後に行われる追肥時の施肥量について適切な量を決定することが従来から求められていた。特に施肥時に緩効性肥料を供給するため、この緩効性肥料の養分の溶出に関する要素を踏まえた適切な施肥量を決定することが求められていた。なお、特許文献１に記載の技術は、今期の農作物の栽培に関する要素に基づいて次期の施肥量を求めるものであり、また緩効性肥料の養分の溶出に関する要素を考慮した技術に関して知見を提供するものではなく、このような要求に応えることはできない。
【０００５】
本発明は、このような問題を解決するために成されたものであり、１つの期の中で緩効性肥料を供給した後に行われる追肥の施肥量に関して、施肥時に緩効性肥料が供給されることを踏まえた適切な施肥量を決定できるようにすることを目的とする。
【課題を解決するための手段】
【０００６】
上記した課題を解決するために、本発明では、今期の初期の施肥後から追肥前に至るまでの積算温度を取得し、取得した積算温度に基づいて、緩効性肥料の残存量である今期肥料残存量または平年と今期との緩効性肥料の残存量の差である残存量差を導出し、導出した今期肥料残存量または残存量差を反映して追肥時施肥量を導出するようにしている。
【発明の効果】
【０００７】
緩効性肥料は、施肥後の積算温度によって養分溶出量が変化し、初期の施肥時に供給した緩効性肥料の量を一定とすると、施肥後の積算温度によって土壌に残存する肥料の量が変化する。更に積算温度と養分溶出量とには一定の関係性があり、積算温度に基づいて養分溶出量を導出することが可能である。そして、上記のように構成した本発明によれば、今期の初期の施肥後から追肥前に至るまでの積算温度に基づいて、追肥時における緩効性肥料の残存量、または、平年と今期との緩効性肥料の残存量の差が求められた上で、この残存量或いはこの差が反映されて追肥時の施肥量が導出される。このため、初期の施肥時に緩効性肥料が供給されるという事情、および、追肥時における肥料の残存の状況を踏まえた適切な施肥量を導出できる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１実施形態に係る施肥量導出装置の機能構成例を示すブロック図である。
【図２】測定タイミングにおけるＮＤＶＩと追肥時の施肥量との関係を示す図である。
【図３】窒素系の緩効性肥料について、積算温度と単位面積当たりの窒素溶出量との関係を示す図である。
【図４】本発明の第１実施形態に係る施肥量導出装置の動作例を示すフローチャートである。
【図５】本発明の第２実施形態に係る施肥量導出装置の機能構成例を示すブロック図である。
【図６】本発明の第３実施形態に係る施肥量導出装置の機能構成例を示すブロック図である。
【発明を実施するための形態】
【０００９】
＜第１実施形態＞
以下、本発明の第１実施形態を図面に基づいて説明する。図１は、本実施形態に係る施肥量導出装置１の機能構成例を示すブロック図である。施肥量導出装置１は、稲作に関して、追肥（追肥についての詳細は後述）時に施す肥料の量（以下「追肥時施肥量」という）を導出し、追肥時施肥量を示す情報をユーザに提供する装置である。ユーザは、追肥にあたって施肥量導出装置１により提供された追肥時施肥量を示す情報を参考にして、実際に施す肥料の量を適切に決定することができる。
【００１０】
なお、以下の説明では、対象とする稲、稲作について、稲の品種や、稲作が行われる環境（地域（寒冷地域や、温暖地域、乾燥地域等）、標高、圃場の規模）、稲作の方法、肥料の種類、その他の稲の生育に影響を与える外部的な要素は共通しているものとする。従って、施肥量導出装置１により導出された追肥時施肥量は、特定の環境で特定の方法により栽培される特定の品種の稲について、特定の種類の肥料を施す場合に適切な施肥量ということになる。また、以下の説明では、説明の便宜のため、施肥量導出装置１は、１つの圃場を対象として追肥時施肥量を導出するものとする。ただし、まとまった複数の圃場を対象としてもよいことおよび圃場の一部の領域を対象としてもよいことは勿論である。また、本実施形態で「施肥量」という場合、特に説明がない限り単位面積あたりの肥料の量を意味する。
【００１１】
本実施形態において、「初期施肥」および「追肥」とは以下を意味する。すなわち、一般にある１つの期において、稲は、田植えによって苗代から圃場に移植された後、分げつ期→幼穂形成期→出穂期→開花・受粉期→穂揃期→登熟期→成熟期という段階を経て成長する。通常、稲代から圃場への移植時には、窒素系の緩効性肥料の施肥が行われる。そして本実施形態において「初期施肥」とは稲代から圃場への移植時に行われる窒素系の緩効性肥料の施肥を意味する。初期施肥は、いわゆる基肥に相当する。
【００１２】
また、幼穂形成期前後から出穂期前後に至るまでの期間（以下「追肥対象期間」という。）では、穂に着生する籾の個数の増大や、籾に詰まるデンプンの量の増大等を目的として、所定のタイミングで穂肥と呼ばれる肥料の追加が行われる場合がある。そして本実施形態では、「追肥」とは追肥対象期間における所定のタイミングで行われる肥料の追加を意味する。
【００１３】
図１で示すように、施肥量導出装置１には、液晶ディスプレイや、有機ＥＬディスプレイ等の表示装置２と、マウスやキーボード等の入力装置３とが接続されている。また図１で示すように、施肥量導出装置１は、機能構成として、施肥量導出部１０と、仮施肥量導出部１１と、残存量導出部１２とを備えている。上記各機能ブロック１０〜１２は、ハードウェア、ＤＳＰ（Digital Signal Processor）、ソフトウェアの何れによっても構成することが可能である。例えばソフトウェアによって構成する場合、上記各機能ブロック１０〜１２は、実際にはコンピュータのＣＰＵ、ＲＡＭ、ＲＯＭなどを備えて構成され、ＲＡＭやＲＯＭ、ハードディスクまたは半導体メモリ等の記録媒体に記憶されたプログラムが動作することによって実現される。また、施肥量導出装置１は、記憶手段として、記憶部１３を備えている。記憶部１３に記憶されるデータについては後述する。
【００１４】
本実施形態では、施肥量導出装置１は、ユーザから入力装置３に対して追肥時施肥量を示す情報の提供の指示（以下「情報提供指示」という）があったことをトリガとして、追肥時施肥量を導出し、当該情報を提供する。以下、ユーザから入力装置３に対して情報提供指示があったときの施肥量導出装置１の動作について詳述する。なお、追肥時施肥量を導出するトリガは本実施形態で例示するものに限られず、例えば、予め定められた日時が到来したことをトリガとしてもよい。
【００１５】
施肥量導出部１０は、入力装置３に対して情報提供指示があったことを検出すると、仮施肥量導出部１１に対して仮追肥時施肥量（後述）の導出を要求する。
【００１６】
仮施肥量導出部１１は、施肥量導出部１０から仮追肥時施肥量の導出の要求があった場合、以下の処理を実行する。すなわち、まず、仮施肥量導出部１１は、記憶部１３に記憶された分析用画像データを取得する。分析用画像データとは、追肥を行う前の所定のタイミング（以下「測定タイミング」という）で、ドローンに搭載されたマルチスペクトルセンサにより、追肥を行う対象の圃場（以下「対象圃場」という）の稲の群落を上空から撮影することによって生成された画像データである。なお測定タイミングは、追肥時施肥量の導出に利用する分析用画像データを得るための撮影を行うタイミングとして予め定められたタイミングであり、例えば、追肥を行うタイミングの直前や、出穂の３週間前、出穂の４週間前等とされる。仮施肥量導出部１１は、取得した分析用画像データを分析し、周知の方法でＮＤＶＩを算出し、取得する。なお本実施形態では説明の便宜のため、指標としてのＮＤＶＩの指標値の実際の値を「ＮＤＶＩ」と表現する。ここで取得されたＮＤＶＩは、特許請求の範囲の「植生関連値」および「第１指標値」に相当する。
【００１７】
なお、本実施形態では、ＮＤＶＩの算出、取得にあたり、ドローンに搭載されたマルチスペクトルセンサの撮影結果に基づく分析用画像データを使用する。ここで追肥時施肥量を導出するときに使用する指標値として、ＮＤＶＩに代えて、ＮＤＶＩと同じく葉色に由来する稲の活性度を測る指標として利用可能なＳＰＡＤ値を用いることも可能である。しかしながら、ＳＰＡＤ値を用いる場合、葉緑素計を用いて葉を一枚ずつ測定する必要があり、作業の困難性が高い。特に、圃場が大規模の場合、高い精度でＳＰＡＤ値を算出するためには、適度に分散された地点で葉緑素計を用いた測定を行う必要があり、作業の困難性が高い。ＳＰＡＤ値と同様に葉色板（カラースケイル）を用いた地上でのセンシングにより得られる葉色板値を用いることも可能であるが、この場合もＳＰＡＤ値と同様に困難性が高い。また、携帯型のＮＤＶＩセンサを利用して作業者が地上でＮＤＶＩの測定を行い、これを利用することも可能である。しかしながら、この場合も作業者が圃場を歩き回って測地を行う必要があり、作業の困難性が高い。特に、圃場が大規模、多数の場合は作業の困難性が非常に高い。一方で、本実施形態では、ドローンを利用して、ＮＤＶＩの算出に必要なデータが収集されることになるため、ＮＤＶＩの算出に際して行われる作業が簡易であり、かつ、圃場が大規模であってもスピーディに作業を行うことができる。ただし、ＳＰＡＤ値を用いることを排除するのではなく、ＳＰＡＤ値や、ＳＰＡＤ値に準じた指標を用いるようにしてもよい。また、「ＮＤＶＩセンサを用いた地上での測定により得られるＮＤＶＩ」や、「葉色板（カラースケイル）を用いた地上での測定により得られる葉色板値」を用いることを排除するものではない。
【００１８】
なお、上述した効果は、葉緑素計を用いて葉を一枚ずつ測定するのではなく、稲の群落に対する遠隔からの観測により得られる指標値を、追肥時施肥量を決定するときに使用する指標値として用いることにより得ることが可能である。従って、観測（センシング）に用いる装置は、マルチスペクトルセンサに限られず、例えば、赤外線サーモグラフィであってもよい。また、指標値は、ＮＤＶＩに限られない。また、観測は、必ずしもドローンを用いて行われる必要はなく、塔の上部に設置された装置や、脚立に載った人間が操作する装置を用いて行われたりしてもよい。また、ＮＤＶＩ（またはこれに相当する指標値）の算出に、衛星写真を用いるようにしてもよい。また仮施肥量導出部１１がＮＤＶＩを算出し取得するのではなく、仮施肥量導出部１１が外部の装置により算出されたＮＤＶＩを取得する構成でもよい。
【００１９】
以下、仮施肥量導出部１１により取得されたＮＤＶＩを特に「取得ＮＤＶＩ」という。
【００２０】
取得ＮＤＶＩを取得した後、仮施肥量導出部１１は、施肥量導出計算式を利用して、取得した取得ＮＤＶＩに基づいて仮追肥時施肥量（特許請求の範囲の「仮の追肥時施肥量」に相当）を導出する。施肥量導出計算式とは、ＮＤＶＩを入力とし施肥量を出力とする式であり、単純な数式だけでなく、プログラム上で定義された関数（モデル）や、プログラムから呼び出し可能な関数（モデル）を概念的に表すものである。このことは後述する溶出率計算式についても同様である。
【００２１】
施肥量導出計算式に取得ＮＤＶＩを入力したときに出力される施肥量は、圃場の測定タイミングにおけるＮＤＶＩ（の値）が取得ＮＤＶＩ（の値）である場合に、所定のカテゴリについての農作物（農作物の収穫物も含む）の状態を目標状態とするために必要かつ適切な施肥量である。ここで求められた施肥量は、後述する理由により追肥時に圃場に残存する肥料の量がゼロであることを前提とする。そして、施肥量導出計算式に取得ＮＤＶＩを入力することによって出力として得られる施肥量が「仮追肥時施肥量」に相当する。
【００２２】
本実施形態において、所定のカテゴリについての農作物の状態とは、例えば、収穫物の外観的あるいは食味的な品質や、収穫物の収量、籾数、追肥後の生育過程の指標等である。追肥後の生育過程の指標は、例えば、追肥した後の穂揃期の葉緑素計値（ＳＰＡＤ値）である。なお、穂揃期の葉緑素計値（ＳＰＡＤ値）が一定値以上の場合は、品質が良いということができる。
【００２３】
例えば、施肥量導出計算式の出力は、収穫物の収量（単位当たりの収量）を○○ｋｇとするために必要かつ適切な追肥時施肥量である。この場合、収穫物の収量が「所定のカテゴリについての農作物の状態」に相当し、収量が○○ｋｇであるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって（ただし、追肥時に肥料が残存していない状態であるものとする）、収穫物の収量が○○ｋｇとなることが期待される。また例えば、施肥量導出計算式の出力は、収穫物全体に対する低外観品質の収穫物の割合（以下「低品質割合」という）を○○％とするために必要かつ適切な追肥時施肥量である。なお、低外観品質の収穫物は、例えば、白未熟粒（腹白粒や、背白粒、基部未熟粒、心白粒、乳白粒等）や、その他の未熟粒、被害粒、死米、着色粒等である。この場合、低品質割合が「所定のカテゴリについての農作物の状態」に相当し、低品質割合が○○％であるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって（ただし、追肥時に肥料が残存していない状態であるものとする）、低品質割合が○○％となることが期待される。なお、低品質割合に代えて、収穫物の「整粒」の割合としてもよいことは勿論である。
【００２４】
図２は、測定タイミングにおけるＮＤＶＩを縦軸とし追肥時の施肥量を横軸とする２次元空間に描画された１次式によって、上述した低品質割合（所定のカテゴリについての農作物の状態）を５％（目標状態）とするときの、測定タイミングにおけるＮＤＶＩと追肥時の施肥量との関係を単純化して示している。図２で示すように、低品質割合を５％とするときの追肥時の施肥量と、測定タイミングにおけるＮＤＶＩとの間には、ＮＤＶＩが大きくなるほど施肥量が小さくなる関係がある。このように、所定のカテゴリについての農作物の状態を目標状態とするときの、測定タイミングにおけるＮＤＶＩと追肥時の施肥量とには強い相関関係がある。
【００２５】
本実施形態では、施肥量導出計算式は、蓄積された過去の実際の追肥時の施肥量および測定タイミングにおけるＮＤＶＩ（植生関連値）の組み合わせと農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングにおける植生関連値を説明変数とする重回帰分析を行って求められた重回帰式に基づいて求められる。
【００２６】
詳述すると、事前に、実験や既存の圃場に対する観察により、追肥時の施肥量および測定タイミングにおけるＮＤＶＩの組み合わせと、農作物の状態との関係を示すデータ、つまり、追肥時の施肥量が○○であり、測定タイミングにおけるＮＤＶＩが○○であった場合に、所定のカテゴリについての農作物の状態は○○であったというデータが蓄積される。このとき、全てのデータについて、追肥時に残存する肥料の量はゼロであることが前提とされる。これは、全てのデータについて追肥時の施肥量についての条件を統一するため、および、条件を統一するにあたって残存する肥料の量をゼロ以外の一定量に統一することは非常に困難であるためである。
【００２７】
そして、蓄積されたデータについて、農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングにおけるＮＤＶＩを説明変数とする重回帰分析が行われ、農作物の状態を目標状態とするための追肥時の施肥量と測定タイミングにおけるＮＤＶＩとの関係が定義された重回帰式が算出される。例えば、収穫物の質が上述した低品質割合である場合には、低品質割合を目標値（一例として５％）とするための、追肥時の施肥量と測定タイミングにおけるＮＤＶＩとの関係が定義された重回帰式が算出される。そして重回帰式に基づいて、ＮＤＶＩを入力とし、重回帰式により表されるＮＤＶＩと施肥量との関係から決定される施肥量を出力とする施肥量導出計算式が生成される。この施肥量導出計算式の出力は、追肥時に緩効性肥料が残存していないという条件下で、測定タイミングにおけるＮＤＶＩが入力された値である場合に農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量に相当する。なお重回帰式を、稲の生育時期毎、例えば出穂前３０日、２５日、２０日、１５日、１０日、５日といったように作成するようにしてもよい。この場合において、生育診断する時期が出穂前何日なのかを判定するには、作付け地点、品種、移植日、苗の葉齢及び苗姿、気温や、日長時間等から推定したり、幼穂の長さから推定したりすることができる。なお、農作物の状態は、追肥後の地力窒素発現量（詳しくは後述）に影響を受けるため、上述した方法により生成される施肥量導出計算式の出力値は当然、地力窒素発現量の影響も加味された値となる。
【００２８】
なお、本実施形態の施肥量導出計算式の内容は、あくまで単純化した一例である。施肥量導出計算式は、測定タイミングにおけるＮＤＶＩを入力の１つとし、農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量を出力する計算式であれば、どのようなものであってもよい。一例として、測定タイミングにおけるＮＤＶＩ以外に、農作物の状態に影響を与える要素の要素値を入力とし、当該要素を加味して追肥時の施肥量を出力する式であってもよい。当該要素は、一例として、特定の期間の気温（到来していない期間（例えば登熟期）については予想気温）や、特定の期間における雨量等である。また、施肥量導出計算式は重回帰分析により生成されるものである必要はなく、施肥量導出計算機の生成には様々な機械学習（一例としてニューラルネットワークを用いた機械学習）を応用可能である。
【００２９】
さて、仮施肥量導出部１１は、施肥量導出計算式を利用して、取得した取得ＮＤＶＩに基づいて仮追肥時施肥量を導出する。すなわち、仮施肥量導出部１１は、取得ＮＤＶＩを施肥量導出計算式に入力し、その出力を仮追肥時施肥量とする。仮施肥量導出部１１は、導出した仮追肥時施肥量を施肥量導出部１０に応答する。
【００３０】
また、施肥量導出部１０は、入力装置３に対して情報提供指示があったことを検出すると、残存量導出部１２に対して今期肥料残存量（後述）の導出を要求する。
【００３１】
残存量導出部１２は、施肥量導出部１０から今期肥料残存量料の導出の要求があった場合、以下の処理を実行する。すなわち、まず、残存量導出部１２は、記憶部１３に記憶された積算温度データを取得する。積算温度データとは、初期施肥から現時点に至るまでの積算温度が記録されたデータである。積算温度とは、該当期間における毎日の平均気温を合計したものである。また、本実施形態では、現時点から極めて短い期間内に追肥が行われるものとし、現時点と追肥が行われるタイミングとは同視できるものとする。
【００３２】
積算温度データは、施肥量導出装置１をインターネットに接続可能に構成し、残存量導出部１２が外部サーバから取得するようにしてもよく、また、人為的な観測の下で人為的に積算温度データが作成され、記憶部１３に記憶されるようにしてもよい。
【００３３】
積算温度データを取得した後、残存量導出部１２は、当該データに記録された積算温度を取得し、溶出率計算式により窒素溶出率（養分溶出率）を導出する。以下、残存量導出部１２により取得された積算温度を「取得積算温度」という。図３は、窒素系の緩効性肥料について、積算温度と窒素溶出率との関係を、初期施肥後の積算温度を横軸に取り窒素溶出率を縦軸に取った２次元空間におけるグラフによって示す図である。図３で示すように、窒素系の緩効性肥料について、施肥後の積算温度と窒素溶出率とは比例関係にあり、施肥後の積算温度が大きくなるほど、それに比例して窒素溶出率が大きくなる。そして、溶出率計算式は、積算温度を入力とし、窒素溶出率を出力する計算式である。溶出率計算式に積算温度として取得積算温度を入力したときに出力される窒素溶出率は、初期施肥後、追肥時に至るまでの積算温度が取得積算温度であった場合に予想される窒素溶出率である。
【００３４】
なお、本実施形態の溶出率計算式の内容は、あくまで単純化した一例である。溶出率計算式は、積算温度を入力の１つとし、窒素溶出率を出力する計算式であればどのようなものであってもよい。一例として、積算温度以外に、窒素の溶出に影響を与える要素の要素値を入力とし、当該要素を加味して窒素溶出率を出力する計算式であってもよい。当該要素は一例として、特定の期間における日射量である。また、溶出率計算式の生成には、重回帰分析のほか、既存の機械学習を広く適用できる。また、本実施形態では、積算温度を、該当期間における毎日の平均気温を合計したものとしたが、積算温度はこれに限られない。すなわち、積算温度は、該当期間において積算された温度（気温に限られない）に基づくものであって、窒素溶出率を出力する計算式の入力となるものであれば何でもよい。例えば、水田土壌温度としてもよい。この場合に、気温や日射量等の諸条件や、葉面積指数等の植物体の生育状況から水田土壌温度を推定するようにしてもよい。
【００３５】
さて、残存量導出部１２は、溶出率計算式を利用して、取得した取得積算温度に基づいて窒素溶出率を導出する。すなわち、残存量導出部１２は、取得積算温度を溶出率計算式に入力し、その出力を窒素溶出率とする。
【００３６】
次いで、残存量導出部１２は、記憶部１３に記憶された初期施肥量データを取得する。初期施肥量データは、初期施肥時の施肥量（単位面積あたりの量として表されたもの）が記録されたデータである。次いで、残存量導出部１２は、当該データに記録された施肥量を取得する。次いで、残存量導出部１２は、取得した施肥量（＝初期施肥時の施肥量）に対して導出した窒素溶出率を乗算して窒素溶出量を求め、更に取得した施肥量から、算出した窒素溶出量を減算し、今期肥料残存量（単位面積あたりの量として表されたもの）を導出する。次いで、残存量導出部１２は、導出した今期肥料残存量を施肥量導出部１０に応答する。つまり今期肥料残量は、追肥時における、初期施肥時に施した緩効性肥料の残存量である。
【００３７】
なお、本実施形態では、積算温度から窒素溶出率を求め、初期施肥時の施肥量に窒素溶出率を乗算して窒素溶出量を求め、初期施肥時の施肥量から窒素溶出量を減算することにより今期肥料残存量を求める構成である。しかしながら、今期肥料残存量を求める流れは、この流れと完全に一致している必要はない。例えば、積算温度から窒素残存率を求める計算式を用いて、この計算式に積算温度を入力することによって窒素残存率を求め、初期施肥時の施肥量に窒素残存率を乗算して今期肥料残存量を求めてもよい。また、まず窒素溶出率を求め、「１」から窒素溶出率を減算して窒素残存率を求め、初期施肥時の施肥量に窒素残存率を乗算して今期肥料残存量を求めるようにしてもよい。
【００３８】
施肥量導出部１０は、仮施肥量導出部１１から仮追肥時施肥量の応答を受けると共に、残存量導出部１２から今期肥料残存量の応答を受けると、以下の処理を実行する。すなわち、施肥量導出部１０は、仮施肥量導出部１１により導出された仮追肥時施肥量から残存量導出部１２により導出された今期肥料残存量を引いて求めた値を最終的な追肥時施肥量として導出する。
【００３９】
例えば、仮施肥量導出部１１により導出された仮追肥時施肥量が、１０アール当たり窒素成分量で４ｋｇであったとする。また、残存量導出部１２により導出された今期肥料残存量が、１０アール当たり窒素成分量で２ｋｇであったとする。この場合、施肥量導出部１０は、１０アール当たり窒素成分量で４ｋｇから２ｋｇを引いた２ｋｇを最終的な追肥時施肥量として導出決定する。
【００４０】
施肥量導出部１０がこのようにして最終的な追肥時施肥量を導出する理由は以下である。すなわち、仮追肥時施肥量は、追肥時に緩効性肥料が残存していないという条件の下での適切な追肥時施肥量である。しかしながら現実には、追肥時に緩効性肥料が残存している場合があり、残存量導出部１２により導出された今期肥料残存量が、追肥時に緩効性肥料が残存している場合の肥料の残存量である。従って、仮施肥量導出部１１により導出された仮追肥時施肥量から残存量導出部１２により導出された今期肥料残存量を引いて求めた値を最終的な追肥時施肥量とすることによって、最終的な追肥時施肥量について、緩効性肥料の残存量を反映した適切な値とすることができる。
【００４１】
最終的な追肥時施肥量を導出した後、施肥量導出部１０は、表示装置２を制御して、導出した追肥時施肥量を示す情報を表示装置２の表示領域に表示する。ユーザは、表示装置２を参照することにより、施肥量導出装置１が導出した追肥時施肥量を認識でき、追肥時施肥量を参考とすることができる。
【００４２】
次に本実施形態に係る施肥量導出装置１の動作について図４のフローチャートを用いて説明する。図４（Ａ）は施肥量導出部１０の動作を示し、図４（Ｂ）は仮施肥量導出部１１の動作を示し、図４（Ｃ）は残存量導出部１２の動作を示している。
【００４３】
図４（Ａ）で示すように、施肥量導出部１０は、情報提供指示があったか否かを監視する（ステップＳＡ１）。情報提供指示があった場合（ステップＳＡ１：ＹＥＳ）、仮施肥量導出部１１に対して仮追肥時施肥量の導出を要求する（ステップＳＡ２）。更に施肥量導出部１０は、残存量導出部１２に対して肥料残存量の導出を要求する（ステップＳＡ３）。なおステップＳＡ２とステップＳＡ３との処理の順番は順不同であり、同時に実行されてもよい。
【００４４】
図４（Ｂ）で示すように、仮施肥量導出部１１は、施肥量導出部１０から仮追肥時施肥量の導出の要求があったか否かを監視する（ステップＳＢ１）。当該要求があった場合（ステップＳＢ１：ＹＥＳ）、仮施肥量導出部１１は、記憶部１３に記憶された分析用画像データを取得する（ステップＳＢ２）。次いで、仮施肥量導出部１１は、ステップＳＢ２で取得した分析用画像データを分析し、取得ＮＤＶＩを取得する（ステップＳＢ３）。次いで、仮施肥量導出部１１は、施肥量導出計算式を利用して、取得した取得ＮＤＶＩに基づいて仮追肥時施肥量を導出する（ステップＳＢ４）。次いで、仮施肥量導出部１１は、ステップＳＢ４で導出した仮追肥時施肥量を施肥量導出部１０に応答する（ステップＳＢ５）。
【００４５】
図４（Ｃ）で示すように、残存量導出部１２は、施肥量導出部１０から肥料残存料の導出の要求があったか否かを監視する（ステップＳＣ１）。当該要求があった場合（ステップＳＣ１：ＹＥＳ）、残存量導出部１２は、記憶部１３に記憶された積算温度データを取得する（ステップＳＣ２）。次いで、残存量導出部１２は、当該データに記録された取得積算温度を取得する（ステップＳＣ３）。次いで、残存量導出部１２は、ステップＳＣ３で取得した取得積算温度に基づいて、溶出率計算式を用いて窒素溶出率を導出する（ステップＳＣ４）。
【００４６】
次いで、残存量導出部１２は、初期施肥量データを取得する（ステップＳＣ５）。次いで、残存量導出部１２は、ステップＳＣ５で取得した初期施肥量データに記録された施肥量を取得する（ステップＳＣ６）。次いで、残存量導出部１２は、ステップＳＣ６で取得した施肥量に対して、ステップＳＣ４で導出した窒素溶出率を乗算して窒素溶出量を求め、当該施肥量から窒素溶出量を減算することによって今期肥料残存量を導出する（ステップＳＣ７）。次いで、残存量導出部１２は、ステップＳＣ７で導出した今期肥料残存量を施肥量導出部１０に応答する（ステップＳＣ８）。
【００４７】
図４（Ａ）で示すように、ステップＳＡ３の処理後、施肥量導出部１０は、仮施肥量導出部１１からの仮追肥時施肥量の応答、および、残存量導出部１２からの今期肥料残存量の応答の双方があったか否かを監視する（ステップＳＡ４）。双方の応答があった場合（ステップＳＡ４：ＹＥＳ）、施肥量導出部１０は、仮施肥量導出部１１により導出された仮追肥時施肥量から残存量導出部１２により導出された今期肥料残存量を引いて求めた値を最終的な追肥時施肥量として導出する（ステップＳＡ５）。次いで、施肥量導出部１０は、ステップＳＡ５で導出した追肥時施肥量を示す情報を表示装置２の表示領域に表示する（ステップＳＡ６）。
【００４８】
以上詳しく説明したように、本実施形態に係る施肥量導出装置１は、今期の初期の施肥後から追肥前に至るまでの積算温度を取得し、取得した積算温度に基づいて、緩効性肥料の残存量である今期肥料残存量を導出し、導出した今期肥料残存量を反映して追肥時施肥量を導出する。より詳細には本実施形態に係る施肥量導出装置１は、測定タイミングにおけるＮＤＶＩを取得し、ＮＤＶＩを入力の１つとし施肥量を出力とする施肥量導出計算式を利用して、取得したＮＤＶＩに基づいて仮追肥時施肥量を導出する。更に本実施形態に係る施肥量導出装置１は、今期の初期施肥の後から追肥前に至るまでの積算温度を取得し、取得した積算温度に基づいて養分溶出率を求め、養分溶出率から緩効性肥料の残存量を導出する。そして本実施形態に係る施肥量導出装置１は、仮の追肥時施肥量から残存量を減算して最終的な追肥時施肥量を決定するようにしている。
【００４９】
この構成によれば、今期の初期の施肥後から追肥前に至るまでの積算温度に基づいて、追肥時における緩効性肥料の残存量が求められた上で、この残存量が反映されて追肥時の施肥量が導出される。このため、初期の施肥時に緩効性肥料が供給されるという事情、および、追肥時における肥料の残存の状況を踏まえた適切な施肥量を導出できる。より詳細には、追肥時のＮＤＶＩに基づいて施肥量導出計算式により導出される施肥量を単純に追肥時の最終的な施肥量とするのではなく、積算温度に基づいて養分溶出率および緩効性肥料の残存量が導出された上で、施肥量導出用計算式により導出された仮の追肥時施肥量から、この残存量を引いた値が最終的な追肥時の施肥量として決定される。このため、初期の施肥時に緩効性肥料が供給されることを踏まえた適切な施肥量を決定できる。
【００５０】
＜第１実施形態の変形例＞
次に上記第１実施形態の変形例について説明する。上記実施形態では、施肥量導出部１０は、仮施肥量導出部１１により導出された仮追肥時施肥量から残存量導出部１２により導出された肥料残存量を減算した値（以下「第１値」という）を、最終的な追肥時施肥量として決定していた。一方で、施肥量導出部１０が、第１値を導出した後、この第１値が予め定められた上限値（単位面積あたりの量として表されたもの）を上回る場合には、この上限値を最終的な追肥時施肥量として決定するようにしてもよい。ここで、上限値は、追肥時のＮＤＶＩの多寡にかかわらず、この上限値を超えて肥料を供給した場合に、収穫物の質に悪影響がある可能性があるような値とされる。上限値は、事前の実験や、シミュレーションに基づいて定められる。
【００５１】
＜第２実施形態＞
次に第２実施形態について説明する。図５は、本実施形態に係る施肥量導出装置１Ａの機能構成例を示すブロック図である。図１と図５との比較で明らかな通り施肥量導出装置１Ａは、第１実施形態に係る施肥量導出部１０に代えて施肥量導出部１０Ａを備えている。また施肥量導出装置１Ａは、地力窒素予測発現量取得部２０を備えている。以下、ユーザから入力装置３に対して情報提供指示があったときの施肥量導出装置１Ａの動作について詳述する。
【００５２】
施肥量導出部１０Ａは、情報提供指示に応じて、仮施肥量導出部１１に対して仮追肥時施肥量の応答を要求すると共に、残存量導出部１２に今期肥料残存量の応答を要求し、これに加えて地力窒素予測発現量取得部２０に地力窒素予測発現量（後述）の応答を要求する。
【００５３】
地力窒素予測発現量取得部２０は、施肥量導出部１０Ａからの要求に応じて、地力窒素予測発現量を取得する。地力窒素予測発現量とは、追肥後から収穫前に至るまでの期間における地力窒素の発現量の予測値である。周知の通り、地力窒素とは、人工的に施肥された肥料ではなく、土壌に残っている窒素のことをいう。地力窒素は、微生物に徐々に分解され、農作物に吸収されることにより、農作物の生育に影響を与えることが知られている。
【００５４】
地力窒素予測発現量取得部２０の動作について詳述すると、施肥量導出部１０Ａから地力窒素予測発現量の応答の要求があると、地力窒素予測発現量取得部２０は、記憶部１３に記憶された地力窒素予測発現量データを取得する。地力窒素予測発現量データとは、事前に導出された地力窒素予測発現量が記録されたデータである。地力窒素予測発現量は、追肥後から収穫前に至るまでの期間における地力窒素の発現量に影響を与える種々の要素（例えば、予測される積算温度（気温や水田土壌温度を含む地温等）や、農作物の品種、土壌に関する諸条件、栽培方法に関する諸条件、農作物が所定の状態となるタイミング（例えば、幼穂が形成されたタイミングや、出穂のタイミング、登熟すると予測されるタイミング等）等）に基づいて事前に導出される。地力窒素予測発現量取得部２０は、取得した地力窒素予測発現量データに基づいて地力窒素予測発現量を認識し、施肥量導出部１０に応答する。
【００５５】
なお地力窒素予測発現量取得部２０が、追肥後から収穫前に至るまでの期間における地力窒素の発現量に影響を与える種々の要素の要素を入力し、地力窒素予測発現量を出力するモデルを利用して、地力窒素予測発現量を導出する構成でもよい。この場合、当該種々の要素は事前に記憶部１３に記憶されるか、都度ユーザに入力される。
【００５６】
施肥量導出部１０Ａは、仮施肥量導出部１１から仮追肥時施肥量を取得し、残存量導出部１２から今期肥料残存量を取得し、地力窒素予測発現量取得部２０から地力窒素予測発現量を取得する。次いで施肥量導出部１０Ａは、取得した地力窒素予測発現量のレベルを判定する。本実施形態では、地力窒素予測発現量の大きさに応じて、大レベル、中レベル、小レベルの３つのレベルが事前に設定されており、各レベルを区分けする閾値が事前に設定されている。次いで施肥量導出部１０Ａは、レベルに応じて事前に設定された対応値（対応値の意義については後述）を設定する。各レベルの対応値は、レベルが大きいほど、大きさが大きい。すなわち、「大レベルに対応する対応値＞中レベルに対応する対応値＞小レベルに対応する対応値」である。
【００５７】
次いで施肥量導出部１０Ａは、仮追肥時施肥量および今期肥料残存量に基づいて、仮追肥時施肥量から今期肥料残存量を引いた値（以下「暫定値」という。この暫定値は、第１実施形態で導出される最終的な追肥時施肥量に相当）を導出する。次いで施肥量導出部１０Ａは、暫定値から、地力窒素予測発現量のレベルに対応する対応値を引いた値を求め、この値を最終的な追肥時施肥量として導出する。施肥量導出部１０Ａは、第１実施形態と同様、導出した追肥時施肥量を示す情報を表示装置２に表示する。
【００５８】
ここで、上述したように仮追肥時施肥量は、過去の実際の地力窒素発現量が反映して生成された施肥量導出計算式を利用して導出される値であり、当然、暫定値にも、過去の実際の地力窒素発現量が反映されている。しかしながら、暫定値は、今期において発現することが予測される地力窒素発現量（＝地力窒素予測発現量）が的確に反映されているとは言えない。そして、「暫定値から、地力窒素予測発現量のレベルに対応する対応値を引く」処理は、最終的に導出される追肥時施肥量に、今期の地力窒素予測発現量を的確に反映させるための処理である。より具体的には、当該処理は、地力窒素予測発現量の大きさに応じて、地力窒素予測発現量の大きさが大きいほど、最終的に導出される追肥時施肥量が適切に小さくなるよう調整し、これにより最終的に導出される追肥時施肥量の精度を向上させる処理である。なお、各レベルの対応値は、地力窒素予測発現量に応じて最終的に導出される追肥時施肥量を調整するという観点の下、テストやシミュレーションの結果に基づいて適切に設定される。
【００５９】
本実施形態によれば、以下の効果を奏する。すなわち、追肥後、農作物に吸収される窒素には、追肥時に人工的に提供された肥料に由来する窒素だけではなく、地力窒素も含まれる。そして本実施形態によれば、地力窒素予測発現量も考慮して追肥時施肥量が導出されるため、より精度の高い追肥時施肥量を導出できる。
【００６０】
なお第２実施形態について、第１実施形態の変形例を応用して以下の構成としてもよい。すなわち施肥量導出部１０Ａが、暫定値から対応値を引いた値が予め定められた上限値を上回っているときに、この上限値を最終的な追肥時施肥量として決定する構成としてもよい。
【００６１】
また第２実施形態について、暫定値から地力窒素予測発現量のレベルに応じた対応値を引いた値を最終的な追肥時施肥量としたのは単純化した一例であり、他の方法で最終的な追肥時施肥量が導出される構成でもよい。すなわち施肥量導出部１０Ａが、暫定値（仮追肥時施肥量から今期肥料残存量を引いた値）に対して、地力窒素予測発現量を反映して最終的な追肥時施肥量を導出する構成であればよい。一例として、施肥量導出部１０Ａが、地力窒素予測発現量に対して所定の係数を乗じた値を導出し、暫定値からこの値を引いて求めた値を最終的な追肥時施肥量とする構成でもよい。この場合、所定の係数は、最終的な追肥時施肥量を適切な値とするために、地力窒素予測発現量の値を調整するような係数とされ、事前のテストやシミュレーションにより適切に決定される。また、施肥量導出計算式が、どの程度の地力窒素発現量を前提として出力値を導出しているのかを把握し、施肥量導出計算式が前提とする地力窒素発現量と、今期の予測値である地力窒素発現量との差分を、暫定値に反映させる構成でもよい。
【００６２】
また第２実施形態について以下の構成としてもよい。すなわち、第２実施形態体は、初期施肥時（基肥時）に緩効性肥料を施用することを前提としていた。一方、初期施肥時に速効性肥料を施用するケースがある。このケースでは、追肥時の追肥量を決定するにあたって、肥料の残存量を考慮する必要はない。一方、追肥後の地力窒素発現量を考慮することによって、精度を高めることが可能である点は上述した通りである。そして上記ケースを想定して、第２実施形態に係る施肥量導出装置１Ａを以下の構成とすることができる。
【００６３】
第２実施形態に係る施肥量導出装置１Ａについて、残存量導出部１２を有さない構成とする。その上で、ユーザから入力装置３に対して情報提供指示があった場合、施肥量導出装置２Ａは以下の処理を実行する。すなわち、施肥量導出部１０Ａは、情報提供指示に応じて、仮施肥量導出部１１に対して仮追肥時施肥量の応答を要求すると共に、地力窒素予測発現量取得部２０に地力窒素予測発現量の応答を要求する。
【００６４】
仮施肥量導出部１１は、施肥量導出計算式を利用して、仮追肥時施肥量を導出する。この施肥量導出計算式は、測定タイミングにおけるＮＤＶＩを入力し、初期施肥時に速効性肥料が施用されるという環境において、農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量を出力する計算式とされる。仮施肥量導出部１１は、導出した仮追肥時施肥量を施肥量導出部１０Ａに応答する。更に地力窒素予測発現量取得部２０は、地力窒素予測発現量を取得し、施肥量導出部１０Ａに応答する。
【００６５】
施肥量導出部１０Ａは、仮施肥量導出部１１から仮追肥時施肥量を取得すると共に、地力窒素予測発現量取得部２０から地力窒素予測発現量を取得する。施肥量導出部１０Ａは、仮追肥時施肥量に地力窒素予測発現量を反映させて最終的な追肥時施肥量を導出する。仮追肥時施肥量に地力窒素予測発現量を反映させる方法の例は、第２実施形態で説明した通りであり、例えば、仮追肥時施肥量から、地力窒素予測発現量のレベルに対応する対応値を減算し、また例えば、仮追肥時施肥量から、地力窒素予測発現量に所定の係数を乗じた値を減算する。以上のような構成によれば、初期施肥時に初期施肥時に速効性肥料を施用するケースにおいても、予想される地力窒素の発現量を反映して、導出される追肥時の施肥量の精度を向上できる。
【００６６】
＜第３実施形態＞
次に第３実施形態について説明する。図６は、本実施形態に係る施肥量導出装置１Ｂの機能構成例を示すブロック図である。図１と図６との比較で明らかな通り、施肥量導出装置１Ｂは、第１実施形態に係る施肥量導出部１０に代えて施肥量導出部１０Ｂを備えている。また施肥量導出装置１Ｂは、第１実施形態に係る残存量導出部１２に代えて残存量導出部１２Ｂを備えている。また施肥量導出装置１Ｂは、第１実施形態に係る仮施肥量導出部１１を備えていない。以下、ユーザから入力装置３に対して情報提供指示があったときの施肥量導出装置１Ｂの動作について詳述する。
【００６７】
施肥量導出部１０Ｂは、情報提供指示に応じて、残存量導出部１２Ｂに対して残存量差（後述）の応答を要求する。残存量導出部１２Ｂは、施肥量導出部１０Ｂからの要求に応じて、以下の処理を実行する。すなわち、まず残存量導出部１２Ｂは、第１実施形態と同様の方法で積算温度に基づいて今期肥料残存量を導出する。次いで残存量導出部１２Ｂは、記憶部１３に記憶された平年肥料残存量データを参照し、当該データに記録された平年肥料残存量を認識する。平年肥料残存量とは、平年の追肥時における緩効性肥料の残存量（初期施肥で施された肥料の残存量）を意味する。平年肥料残存量は、過去の実際の圃場における残存量の実際値や、緩効性肥料の残存量に関するテスト／シミュレーションに基づいて事前に導出される。
【００６８】
次いで残存量導出部１２Ｂは、平年肥料残存量と、導出した今期肥料残存量との差（平年肥料残存量−今期肥料残存量）を残存量差として導出する。このようにして導出される残存量差は、平年の肥料の残存量に対する、今期の肥料の残存量の過不足量を示している。残存量差は、「平年肥料残存量＞今期肥料残存量」であればプラスになり、「平年肥料残存量＜今期肥料残存量」であればマイナスになる。残存量導出部１２Ｂは、導出した残存量差を施肥量導出部１０Ｂに応答する。
【００６９】
残存量導出部１２Ｂから残存量差を取得すると、施肥量導出部１０Ｂは、以下の処理を実行する。すなわち施肥量導出部１０Ｂは、記憶部１３に記憶された慣行追肥時施肥量データを参照し、当該データに記録された慣行追肥時施肥量を認識する。慣行追肥時施肥量とは、慣行の追肥時施肥量であり、追肥時の肥料の残存量が上述した平年肥料残存量であることを想定して設定された追肥時施肥量の基準値を意味している。次いで施肥量導出部１０Ｂは、慣行追肥時施肥量に残存量差を加算して、最終的な追肥時施肥量を導出する。
【００７０】
残存量差がプラスの場合は、今期は平年に比して、初期施肥において施した緩効性肥料の残存量が少ないということである。そしてこの場合、慣行追肥時施肥量に残存量差（＞０）を加えた値が最終的な追肥時施肥量とされるため、最終的な追肥時施肥量は、平年に比して不足する分の肥料が適切に補われた値となる。一方、残存量値がマイナスの場合は、今期は平年に比して、初期施肥において施した緩効性肥料の残存量が多いということである。そしてこの場合、慣行追肥時施肥量に残存量差（＜０）を加えた値（つまり残存量差の絶対値を引いた値）が最終的な追肥時施肥量とされるため、最終的な追肥時施肥量は、平年に比して過剰な分の肥料が適切に削減された値となる。
【００７１】
以上の通り、本実施形態によれば、今期の初期の施肥後から追肥前に至るまでの積算温度に基づいて平年と今期との緩効性肥料の残存量の差である残存量差が求められた上で、この残存量差が反映されて追肥時の施肥量が導出される。このため、初期の施肥時に緩効性肥料が供給されるという事情、および、追肥時における肥料の残存の状況を踏まえた適切な施肥量を導出できる。
【００７２】
＜第３実施形態の変形例＞
次に第３実施形態の変形例について説明する。上記第３実施形態では、残存量導出部１２Ｂは、今期肥料残存量を導出し、平年肥料残存量から導出した今期肥料残存量を引くことによって残存量差を導出した。この点に関し、以下の構成でもよい。
【００７３】
すなわち、平年の積算温度（第１実施形態で説明したように、初期施肥から追肥時に至る期間の積算温度）と、今期の積算温度とを入力とし、残存量差（平年肥料残存量と今期肥料残存量との差）を出力とするモデルが事前に生成される。なお、「平年の積算温度と今期の積算温度との差」と、「平年肥料残存量と今期肥料残存量との差（残存量差）」とには強い相関関係があるため、相関関係を分析することにより、上述したモデルを生成することができる。なお、モデルに、平年の積算温度および今期の積算温度が入力されるのではなく、これら積算温度の差が入力される構成としてもよい。また、モデルの入力に、平年の積算温度および今期の積算温度（或いはこれらの積算温度の差）以外の、肥料の残存量に影響を与える１つ以上の要素の要素値が含まれ、モデルは、１つ以上の要素値を反映して、残存量差を導出するものでもよい。
【００７４】
残存量導出部１２Ｂは、残存量差の導出に際し、以下の処理を実行する。すなわち残存量導出部１２Ｂは、平年の積算温度および今期の積算温度を認識する。これらは、記憶部１３に記憶されたデータに事前に記録される。ただし、ユーザが都度入力する構成でもよい。そして残存量導出部１２Ｂは、モデルに平年の積算温度および今期の積算温度を入力し、その出力を残存量差とする。以上の構成でもよい。
【００７５】
以上、本発明の一実施形態（変形例を含む）を説明したが、上記各実施形態は、本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその要旨、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００７６】
例えば上記各実施形態では、施肥量導出装置１は、稲を対象として、稲作の追肥時の施肥量を導出した。しかしながら、対象とする農作物は稲に限られない。すなわち、本発明は、栽培に際して緩効性肥料の施肥が行われた後に追肥が行われる農作物に広く適用可能である。一例として小麦、大麦、大豆その他の穀物や、野菜、果物等に適用可能である。
【００７７】
また上記第１実施形態（変形例を含む）では、緩効性肥料は窒素系の肥料であった。しかしながら、本発明が対象とする緩効性肥料は窒素系の肥料に限らず、積算温度と養分溶出率との間に何らかの関係があり、積算温度に基づいて養分溶出率が予測できるものであれば何でもよい。一例として、水田土壌の積算温度で推定可能な地力窒素であってもよい。
【００７８】
また、記憶部１３は、施肥量導出装置１とは異なる外部記憶装置が有する構成でもよい。
【００７９】
また、施肥量導出装置１を、インターネット上のサーバとして構成し、インターネットを介してクライアントから追肥時施肥量を示す情報の提供の要請があったときに、上述した手法で追肥時施肥量を導出し、クライアントに提供する構成としてもよい。
【００８０】
また、施肥量導出装置１は、単体のコンピュータである必要はなく、複数のコンピュータにより構成されていてもよい。例えば、インターネットを介して接続された端末とクラウドサーバとが協働して施肥量導出装置１として機能し、端末が適宜クラウドサーバと協働して処理を実行する構成でもよい。
【００８１】
また本実施形態では、初期施肥を「圃場への移植時に行われる窒素系の緩効性肥料の施肥」と定義したが、これはあくまで初期施肥の一例である。初期施肥は、ある一つの期において、追肥の前の所定のタイミングで行われる施肥であればよい。
【符号の説明】
【００８２】
１、１Ａ、１Ｂ  施肥量導出装置
１０、１０Ａ、１０Ｂ  施肥量導出部
１１  仮施肥量導出部
１２、１２Ｂ  残存量導出部
２０  地力窒素予測発現量取得部
【図１】
0006916558.tif 000002
【図２】
0006916558.tif 000003
【図３】
0006916558.tif 000004
【図４】
0006916558.tif 000005
【図５】
0006916558.tif 000006
【図６】
0006916558.tif 000007