知財求人 - 知財ポータルサイト「IP Force」
7681303養分供給量推定方法、養分供給体の利用方法、プログラム、および推定システム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2025-05-14
(45)【発行日】2025-05-22
(54)【発明の名称】養分供給量推定方法、養分供給体の利用方法、プログラム、および推定システム
(51)【国際特許分類】
A01C 21/00 20060101AFI20250515BHJP
【ＦＩ】
A01C21/00 Z
【請求項の数】
7
(21)【出願番号】P 2021144059
(22)【出願日】2021-09-03
(65)【公開番号】P2023037363
(43)【公開日】2023-03-15
【審査請求日】2024-05-20
【新規性喪失の例外の表示】特許法第３０条第２項適用 ・刊行物名「日本土壌肥料科学雑誌、第９１巻、第５号、第３６６～３７３頁」 発行日 ２０２０年１０月５日 ・ウェブサイトのアドレス ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔａｇｅ．ｊｓｔ．ｇｏ．ｊｐ／ａｒｔｉｃｌｅ／ｄｏｊｏ／９１／５／９１＿９１０５０３／＿ａｒｔｉｃｌｅ／－ｃｈａｒ／ｊａ／ 掲載日 ２０２０年１０月１０日
(73)【特許権者】
【識別番号】501203344
【氏名又は名称】国立研究開発法人農業・食品産業技術総合研究機構
(74)【代理人】
【識別番号】110000338
【氏名又は名称】弁理士法人 ＨＡＲＡＫＥＮＺＯ ＷＯＲＬＤ ＰＡＴＥＮＴ ＆ ＴＲＡＤＥＭＡＲＫ
(72)【発明者】
【氏名】原 嘉隆
【審査官】大澤 元成
(56)【参考文献】
【文献】特開２０２０－０５４３９５（ＪＰ，Ａ）
【文献】特開２０１５－０２７２９６（ＪＰ，Ａ）
【文献】米国特許出願公開第２０２０／００６８７９７（ＵＳ，Ａ１）
【文献】特開２００４－３５２７２７（ＪＰ，Ａ）
【文献】国際公開第９９／０４９８６３（ＷＯ，Ａ１）
【文献】Yoshitaka Hara，Application of the Richards function to nitrogen release from coated urea at a constant temperature and relationships among the calculated parameters，Soil Science and Plant Nutrition，第４６巻第３号，2000年，pp.683-691
【文献】石橋英二ら，堆肥等有機質資材からの窒素無機化率の推定における反応速度論的手法の新たな解析方法の提案，日本土壌肥料科学雑誌，第８５巻，第４号，2014年08月，pp.362-368
【文献】杉原進ら，土壌中における有機態窒素無機化の反応速度論的解析法，農業環境技術研究所報告，第１号，1986年03月，pp.127-166
【文献】小林新ら，ガウス補正法による溶出モデル式の改良，日本土壌肥料学雑誌，第６８巻第５号，日本，1997年，pp.487-492
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ０１Ｃ ２１／００
Ａ０１Ｇ ７／００
Ｇ０６Ｑ ５０／０２
(57)【特許請求の範囲】
【請求項1】
土壌、有機物および有機質肥料から選択される養分供給体から供給される養分であって、無機態窒素である前記養分の量を推定する養分供給量推定方法であって、養分供給量を推定するためのパラメータ群を含む下記式（ｉ）、又は当該式（ｉ）を数学的に変換した式を参照して、任意の時点における前記養分供給量を推定する養分供給量推定工程を包含し、
ｎ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］ －１／ｄ ＋ｎ ｉ 式（ｉ）
前記パラメータ群には、
基準となる時点における前記養分供給量を示す第１のパラメータｎ ｉ と、
前記養分供給量の最大変化量を示す第２のパラメータＮと、
前記養分供給量の推移を表す式の形状を示す第３のパラメータｄと、
前記養分供給量の変化における速度定数を示す第４のパラメータｋと、
が含まれ、
式（ｉ）中、
ｎは前記養分供給量を示し、
ｔは前記基準となる時点からの経過時間を示す、
養分供給量推定方法。
【請求項2】
前記第４のパラメータｋは、基準となる温度における、前記養分供給量の変化における速度定数を示す第５のパラメータｋ ｓ と、
当該速度定数の温度依存性を示す第６のパラメータＥと、
を含む下記式（ｉｉ）を参照して、算出され、
ｋ＝ｋ ｓ ・ｅｘｐ｛Ｅ／Ｒ・（１／Ｔ ｓ －１／Ｔ）｝ 式（ｉｉ）
式（ｉｉ）中、
Ｒは気体定数を示し、
Ｔ ｓ は前記基準となる温度を示し、
Ｔは前記養分供給体の温度を示す、
請求項１に記載の養分供給量推定方法。
【請求項3】
前記養分供給量推定工程において、前記式（ｉ）における前記養分供給量ｎに所望の養分供給量を入力することにより、任意の時点における前記養分供給量ｎの代わりに、前記所望の養分供給量に達するための時間ｔを推定する、請求項１または２に記載の養分供給量推定方法。
【請求項4】
（ｉ）請求項１または２に記載の養分供給量推定方法によって推定した養分供給量ｎ、ならびに、（ｉｉ）当該養分供給量ｎの推定に用いた前記パラメータ群に含まれる少なくとも１つのパラメータ、の少なくとも一方を参照して養分供給体の状態を推定する工程を包含し、前記養分供給体の状態は、前記養分供給量ｎの過不足、および前記養分の供給時期の遅速から選択される、
養分供給体の利用方法。
【請求項5】
前記推定する工程において、前記パラメータ群に含まれる少なくとも１つのパラメータから算出される、前記養分供給体から植物に供給される養分の任意の養分供給量が得られるまでの時間を参照して、前記養分供給体の状態を推定する、請求項４に記載の養分供給体の利用方法。
【請求項6】
請求項１～３のいずれか１項に記載の養分供給量推定方法、または、請求項４または５に記載の養分供給体の利用方法をコンピュータに実行させるためのプログラム。
【請求項7】
（ｉ）請求項１または２に記載の養分供給量推定方法によって推定した養分供給量ｎ、（ｉｉ）請求項４または５に記載の養分供給体の利用方法によって推定した養分供給体の状態、ならびに、（ｉｉｉ）当該養分供給量ｎの推定に用いた前記パラメータ群に含まれるパラメータの値のいずれか、のうちの少なくとも１つを、養分供給体と関連付けて記憶するデータベースと、前記データベースに記憶されたデータを参照して、被験養分供給体の養分供給量ｎを推定する推定装置と、
を備えた、推定システム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、養分供給量推定方法、養分供給体の利用方法、プログラム、および推定システムに関する。
【背景技術】
【0002】
近年、アプリケーションなどの営農支援サービスが数多く開発されている。営農支援サービスのうち、栽培支援サービスにおいて、土壌、有機物および肥料などの養分供給体から作物への養分供給量を把握することは、精密に肥料施用を行い、作物生育の予測および制御を行うスマート農業の推進のために重要であることが多い。これまで、養分供給体からの養分供給量を推定する技術が研究されてきた。
【0003】
このような技術として、例えば非特許文献１には、土壌に含まれる窒素の無機化および有機化が並行して進行するというモデルに基づき、２つの一次反応式の差を用いる解析方法が記載されている。
【先行技術文献】
【非特許文献】
【0004】
【文献】石橋英二、藤原（芝）宏子、鷲尾建紀、大家理哉，「堆肥等有機質資材からの窒素無機化率の推定における反応速度論的手法の新たな解析方法の提案」，日本土壌肥料科学雑誌，第８５巻，第４号，２０１４年，ｐ．３６２－３６８
【発明の概要】
【発明が解決しようとする課題】
【0005】
養分供給体からの養分供給量を推定するためのモデルを複雑化させパラメータの種類を増加させると、養分供給量の推定結果は実験的測定結果に良好に合致する一方で、パラメータの値を一定に定めることが難しい傾向がある。非特許文献１に記載の技術は、最大で７つのパラメータが設定されたモデルに基づくため、養分供給体それぞれのパラメータ設定に必要とするデータが多い。また、そのデータの取得に手間がかかるため、生産現場で多用するには至っていない。したがって、養分供給体からの養分供給量を簡易に推定する技術が求められている。
【0006】
本発明の一態様は、上記の問題に鑑みてなされたものであり、その目的は、養分供給体からの養分供給量を簡易に推定する養分供給量推定方法、およびその関連技術を提供することを目的とする。
【課題を解決するための手段】
【0007】
本発明の一態様に係る養分供給量推定方法は、養分供給体から供給される養分であって、植物が利用できる形態である養分の量を推定する養分供給量推定方法であって、養分供給量を推定するためのパラメータ群を参照して、任意の時点における前記養分供給量を推定する養分供給量推定工程を包含し、前記パラメータ群には、基準となる時点における前記養分供給量を示す第１のパラメータと、前記養分供給量の最大変化量を示す第２のパラメータと、前記養分供給量の推移を表す式の形状を表す第３のパラメータと、前記養分供給量の変化における速度定数を示す第４のパラメータと、が含まれる。
【0008】
本発明の一態様に係る養分供給体の利用方法は、（ｉ）本発明の一態様に係る養分供給量推定方法によって推定した養分供給量、ならびに、（ｉｉ）当該養分供給量の推定に用いた前記パラメータ群に含まれる少なくとも１つのパラメータ、の少なくとも一方を参照して養分供給体の状態を推定する工程を包含する。
【0009】
本発明の一態様に係るプログラムは、本発明の一態様に係る養分供給量推定方法、または、本発明の一態様に係る養分供給体の利用方法をコンピュータに実行させるためのプログラムである。
【0010】
本発明の一態様に係る推定システムは、（ｉ）本発明の一態様に係る養分供給量推定方法によって推定した養分供給量、（ｉｉ）本発明の一態様に係る養分供給体の利用方法によって推定した養分供給体の状態、ならびに、（ｉｉｉ）当該養分供給量の推定に用いた前記パラメータ群に含まれるパラメータの値のいずれか、のうちの少なくとも１つを、養分供給体と関連付けて記憶するデータベースと、前記データベースに記憶されたデータを参照して、被験養分供給体の養分供給量を推定する推定装置と、を備える。
【発明の効果】
【0011】
本発明の一態様によれば、養分供給体からの養分供給量を簡易に推定する養分供給量推定方法、およびその関連技術を提供することができる。
【図面の簡単な説明】
【0012】
【図1】実施例の無機態窒素量の推定結果を示す図である。
【発明を実施するための形態】
【0013】
本発明者らは、養分供給体の状態を参照して設定された特定のパラメータを用いることで、任意の時点における養分供給量を容易に推定することができることを見出した。また、本発明者らは鋭意研究を重ね、特定のパラメータそれぞれは、当該養分供給体の状態に強く関連しており、一の養分供給体の状態と他の養分供給体の状態との間の類似性を参照することによって、他の養分供給体におけるパラメータを推定することができることをさらに見出し、本発明を完成させた。
【0014】
〔養分供給量推定方法〕
本発明の一実施形態に係る養分供給量推定方法は、養分供給体から供給される養分であって、植物が利用できる形態である養分の量を推定する養分供給量推定方法であって、養分供給量を推定するためのパラメータ群を参照して、任意の時点における前記養分供給量を推定する養分供給量推定工程を包含し、前記パラメータ群には、基準となる時点における前記養分供給量を示す第１のパラメータと、前記養分供給量の最大変化量を示す第２のパラメータと、前記養分供給量の推移を表す式の形状を表す第３のパラメータと、前記養分供給量の変化における速度定数を示す第４のパラメータと、が含まれる。
本発明の一実施形態に係る養分供給量推定方法は、養分供給体から供給される養分であって、植物が利用できる形態である養分の量を推定する養分供給量推定方法であって、養分供給量を推定するためのパラメータ群を参照して、任意の時点における前記養分供給量を推定する養分供給量推定工程を包含し、前記パラメータ群には、基準となる時点における前記養分供給量を示す第１のパラメータと、前記養分供給量の最大変化量を示す第２のパラメータと、前記養分供給量の推移を表す式の形状を表す第３のパラメータと、前記養分供給量の変化における速度定数を示す第４のパラメータと、が含まれる。
【0015】
〔養分〕
本明細書中で使用される場合、用語「養分」とは、根を介して植物に吸収され、植物の生育に対して有益な作用を与え得る成分、もしくは当該成分の前駆体、またはそれらの組み合わせを意味する。養分の例には、窒素、カリウム、カルシウム、マグネシウム、リン、硫黄、塩素、鉄、ホウ素、マンガン、亜鉛、銅、ニッケルおよびモリブデン、ならびにこれらの任意の組み合わせが含まれるが、これらに限定されない。また、養分は、１種類単独であってもよく、２種類以上の任意の組み合わせであってもよい。推定の対象である養分は、所望の種類が適宜に選択され得る。
本明細書中で使用される場合、用語「養分」とは、根を介して植物に吸収され、植物の生育に対して有益な作用を与え得る成分、もしくは当該成分の前駆体、またはそれらの組み合わせを意味する。養分の例には、窒素、カリウム、カルシウム、マグネシウム、リン、硫黄、塩素、鉄、ホウ素、マンガン、亜鉛、銅、ニッケルおよびモリブデン、ならびにこれらの任意の組み合わせが含まれるが、これらに限定されない。また、養分は、１種類単独であってもよく、２種類以上の任意の組み合わせであってもよい。推定の対象である養分は、所望の種類が適宜に選択され得る。
【0016】
本発明の一実施形態において、その量を推定する対象となる養分は、植物が利用できる形態である。植物が利用できる形態とは、当該植物が根を介して吸収できる形態を意味する。なお、植物の根域外の、植物が即座に吸収できない領域に存在する養分であっても、当該養分に植物の根が近づいた場合に植物が吸収できる場合には、当該養分も、推定の対象となる養分である。植物が利用できる形態の例には：アンモニウム態窒素および硝酸態窒素などの無機態の交換性窒素、ならびにアミノ酸態などの小分子量の一部の有機態窒素などの窒素；カリウムイオンを含む交換性カリウムなどのカリウム；リン酸イオンを含む可溶性リンなどの可給態リンなどのリンが含まれる。
【0017】
〔養分供給体〕
本発明の一実施形態において、その量を推定する対象となる養分は、養分供給体から供給される。本明細書中で使用される場合、用語「養分供給体」とは、主に植物へ、養分を供給する主体となり得る物体を意味する。養分供給体が植物へ養分を供給することは、当該養分供給体が含有する成分が、溶出などの物理的変化、ならびに、他の物質、微生物または熱によって媒介される分解または合成などの化学的変化に適宜に供された後、植物が利用できる形態になることによって、達成され得る。
本発明の一実施形態において、その量を推定する対象となる養分は、養分供給体から供給される。本明細書中で使用される場合、用語「養分供給体」とは、主に植物へ、養分を供給する主体となり得る物体を意味する。養分供給体が植物へ養分を供給することは、当該養分供給体が含有する成分が、溶出などの物理的変化、ならびに、他の物質、微生物または熱によって媒介される分解または合成などの化学的変化に適宜に供された後、植物が利用できる形態になることによって、達成され得る。
【0018】
なお、養分供給体が供給する「養分」には、植物が利用できる形態である「養分」、および植物が利用できる形態になり得る「養分」の両方が含まれる。すなわち、養分供給体が供給する養分の形態は、（ｉ）植物が利用できる形態そのものであってもよいし、（ｉｉ）そのままでは植物が利用できないが、分解などの状態変化によって植物が利用できる形態へと変換される形態であってもよい。
【0019】
養分供給体の例には、土壌、有機物、肥料および根域担体、ならびにこれらの任意の組み合わせが含まれる。養分供給体は、固体であってもよいし、液体であってもよい。
【0020】
土壌の具体例には：岩、礫、砂、シルトおよび粘土などの鉱物；土および泥；土に含まれる有機物、すなわち土壌有機物；土に含まれる有機物と無機物との複合体、すなわち有機無機複合体；が含まれるが、これらに限定されない。
【0021】
有機物の具体例には、糞、動物残渣、植物残渣およびこれらの堆肥が含まれるが、これらに限定されない。
【0022】
肥料の具体例には：成分を基準に分類すると、無機質肥料および有機質肥料；供給する養分を基準に分類すると、窒素肥料、カリ肥料、リン酸肥料、マグネシウム肥料、カルシウム肥料；養分を供給する速さを基準に分類すると、速効性肥料および緩効性肥料；が含まれるが、これらに限定されない（ただし、被覆尿素肥料を除く）。
【0023】
根域担体とは、植物の根域を担持する物体を意味する。根域担体の具体的な例には、スポンジ、ロックウール、ヤシガラ、および水が含まれるが、これらに限定されない。
【0024】
本明細書中で使用される場合、用語「養分供給量」とは、養分供給体から供給される養分の量を意味する。当該量は、経時的積算量、経時的差分量または微分量であり得るが、これらに限定されない。当該量は、適宜に所望の種類が選択され得る。また、選択された種類に応じて、後述する数学的処理は、変形され得る。このような変形は、当業者には公知である。以下では、非限定的な例として、経時的積算量が養分供給量として選択された場合について説明する。
【0025】
本発明の一実施形態に係る養分供給量推定方法は、後述するパラメータ群を定数として含む式に対して、時間の値を変数として入力し、任意に温度の値を変数としてさらに入力することによって、養分供給量を数学的に推定する方法であり得る。パラメータ群に含まれるパラメータの種類は、少なく、最小で４つである。さらに、複数のパラメータはそれぞれ独立して、養分供給体が有する特性のうち１つ以上に関連する。したがって、複数のパラメータのうち１つが変化したときに、他のパラメータが変化することによって最終的に算出される推定値が維持されるという現象、すなわち異なるパラメータ間の干渉、が本発明の一実施形態においては防がれる、または減じられる。ゆえに、本発明の一実施形態において、各パラメータを独立して取り扱うことが可能である。具体的には、複数の養分供給体の状態の間の類似性を参照することによって、各パラメータを独立して推定することができる。そのため、すべての養分供給体において実験的方法により各パラメータを設定することが難しい場合であっても、推定により各パラメータを設定でき、次いでは養分供給量を推定できる。
【0026】
なお、本明細書中で単に「パラメータ」と記載される場合、特に明示しない限り、当該記載は「パラメータ」、「パラメータの値」および「パラメータの値を表す情報」のうち少なくとも１つを意味する。
【0027】
以下に、本発明の一実施形態に係る養分供給量推定方法において参照されるパラメータ群に含まれる各パラメータの非限定的な特性を説明する。なお、説明を簡潔にするため、「養分供給量」を「供給量」と略記することがある。同様に、「養分供給体」を「供給体」と略記することがある。
【0028】
〔第１のパラメータ〕
第１のパラメータｎｉは、基準となる時点における養分供給量を示す。本発明の一実施形態において、基準となる時点は、適宜に選択され得る。基準となる時点の例には、肥料を施用する時点、および当該時点から特定期間の前または後、が含まれるが、これらに限定されない。
第１のパラメータｎｉは、基準となる時点における養分供給量を示す。本発明の一実施形態において、基準となる時点は、適宜に選択され得る。基準となる時点の例には、肥料を施用する時点、および当該時点から特定期間の前または後、が含まれるが、これらに限定されない。
【0029】
一例として、第１のパラメータｎｉは、供給体および植物の少なくとも一方の収集およびこれらが含有する成分の化学的分析を伴う実験的方法を用いて、設定され得る。供給体の複数が種類である場合には、種類ごとに成分の化学的分析を行ってもよい。実験的方法は、当技術分野で公知の方法であってよい。実験的方法の例には、無機態窒素量、またはアンモニウム態窒素量を測定することができる方法が含まれるが、これらに限定されない。他の例として、第１のパラメータｎｉを、過去の知見に基づく養分供給体の種類毎の代表値（平均値など）を用いて設定することもできる。また、肥料を施用する場合、第１のパラメータｎｉは、投与される肥料中の養分の量とみなせる場合がある。
【0030】
〔第２のパラメータ〕
第２のパラメータＮは、養分供給量の最大変化量を示す。最大変化量とは、基準となる時点から十分に時間が経過した時点における供給量と、基準となる時点における供給量との間の差を意味する。一般的に、化学反応ならびに物理的状態の変化によって、供給量は変化し得る。しかしながら、外的要因による、成分の添加および減少ならびに場の状態の変化が無い場合には、十分に時間が経過した時点で供給量は定常状態に達し得る。したがって、第２のパラメータＮは、基準となる時点から定常状態に達する時点までの供給量の変化量に相当し得る。
第２のパラメータＮは、養分供給量の最大変化量を示す。最大変化量とは、基準となる時点から十分に時間が経過した時点における供給量と、基準となる時点における供給量との間の差を意味する。一般的に、化学反応ならびに物理的状態の変化によって、供給量は変化し得る。しかしながら、外的要因による、成分の添加および減少ならびに場の状態の変化が無い場合には、十分に時間が経過した時点で供給量は定常状態に達し得る。したがって、第２のパラメータＮは、基準となる時点から定常状態に達する時点までの供給量の変化量に相当し得る。
【0031】
一例として、第２のパラメータＮは、当技術分野で公知の実験的方法を用いて、供給体が含有する成分の種類および量を決定することによって、設定され得る。別の一例として、第２のパラメータＮは、施用される供給体の量を参照することによって、設定され得る。
【0032】
〔第３のパラメータ〕
第３のパラメータｄは、養分供給量の推移を表す式の形状を示す。本発明の一実施形態において、供給量の推定に用いられる式は、時間の値を変数として供給量の推移を表す。本発明の一実施形態において、この式の形状は特定の１種類に限定されず、複数の種類の数学的式の特性を合わせて有し得る。第３のパラメータｄは、式の形状の決定に関連するものである。典型的な一例として、式の形状は、積分形一次反応速度式、シグモイド式、またはこれらの変形によって得られる式である。なお、変形の様式の例には、積分形一次反応速度式およびシグモイド式の間を内挿する式となること、またはこれらを外挿する式となることの両方が含まれる。なお、これら２つの式の間を内挿する式となることとは、数式に含まれる少なくとも１つのパラメータが所定の値をとることによって、数式の形状がこれら２つの式の形状それぞれに近似する場合に、これら２つの式の形状それぞれに対応する２つの所定の値の間に含まれる値に当該パラメータを設定できることを指す。また、これら２つの式を外挿する式となることとは、数式に含まれるパラメータの少なくとも１つが所定の値をとることによって、数式の形状がこれら２つの式の形状それぞれに近似する場合に、これら２つの式の形状それぞれに対応する２つの所定の値の間に含まれない値に当該パラメータを設定できることを指す。第３のパラメータｄは、変形の様式および程度を決定し得る。
第３のパラメータｄは、養分供給量の推移を表す式の形状を示す。本発明の一実施形態において、供給量の推定に用いられる式は、時間の値を変数として供給量の推移を表す。本発明の一実施形態において、この式の形状は特定の１種類に限定されず、複数の種類の数学的式の特性を合わせて有し得る。第３のパラメータｄは、式の形状の決定に関連するものである。典型的な一例として、式の形状は、積分形一次反応速度式、シグモイド式、またはこれらの変形によって得られる式である。なお、変形の様式の例には、積分形一次反応速度式およびシグモイド式の間を内挿する式となること、またはこれらを外挿する式となることの両方が含まれる。なお、これら２つの式の間を内挿する式となることとは、数式に含まれる少なくとも１つのパラメータが所定の値をとることによって、数式の形状がこれら２つの式の形状それぞれに近似する場合に、これら２つの式の形状それぞれに対応する２つの所定の値の間に含まれる値に当該パラメータを設定できることを指す。また、これら２つの式を外挿する式となることとは、数式に含まれるパラメータの少なくとも１つが所定の値をとることによって、数式の形状がこれら２つの式の形状それぞれに近似する場合に、これら２つの式の形状それぞれに対応する２つの所定の値の間に含まれない値に当該パラメータを設定できることを指す。第３のパラメータｄは、変形の様式および程度を決定し得る。
【0033】
本発明の一実施形態において、供給量を数学的に推定するために用いられる式は、積分形一次反応速度式、シグモイド式、またはこれらの変形によって得られる式である。用いられる式の違いは、供給体が養分を供給する反応（溶出および分解など）の反応系の複雑さに起因すると考えられる。例えば、供給体の種類は、１種類に限られず、２種類以上であることがある。一般的に、供給体の種類が少ないほど反応系がより単純になり、反応の反応速度式は、一次反応速度式により近似していく。このため、供給量の推移を表す式は、時間の値を変数とする積分形一次反応速度式に対して、より高い類似性を示すようになる。逆に、供給体の種類が多いほど反応系がより複雑になり、反応の反応速度式は、高次反応速度式により近似していく。このため、供給量の推移を表す式は、時間の値を変数とするシグモイド式に対して、より高い類似性を示すようになる。
【0034】
本発明の一実施形態において、供給量を数学的に推定するために用いられる式が第３のパラメータｄを含むことにより、当該式に対して反応系の複雑さが反映される。このため、当該式が用いられることにより、供給体の種類および組み合わせによらず、本発明の一実施形態に係る養分供給量推定方法は、包括的で多種多様な種類の養分および供給体を対象として、供給量を推定することができる。
【0035】
〔第４のパラメータ〕
第４のパラメータｋは、養分供給量の変化における速度定数を示す。供給量の変化における速度定数とは、供給体が養分を供給する反応の速度定数を意味する。典型的には、供給体が養分を供給する反応が速いほど、第４のパラメータｋは大きくなり得る。
第４のパラメータｋは、養分供給量の変化における速度定数を示す。供給量の変化における速度定数とは、供給体が養分を供給する反応の速度定数を意味する。典型的には、供給体が養分を供給する反応が速いほど、第４のパラメータｋは大きくなり得る。
【0036】
第４のパラメータｋは、供給体の温度に依存して変化するパラメータであり得る。供給体温度が変化し得る場合には、第４のパラメータｋは、第５のパラメータｋｓと、第６のパラメータＥとを参照して、算出されてもよい。このようにして、第４のパラメータｋが算出されることにより、本発明の一実施形態に係る養分供給量推定方法は、温度変化に応じて、養分供給量を推定することができる。
【0037】
すなわち、本発明の一実施形態に係る養分供給量推定方法において、第４のパラメータｋは、基準となる温度における、養分供給量の変化における速度定数を示す第５のパラメータｋｓと、当該速度定数の温度依存性を示す第６のパラメータＥと、を参照して、算出されてもよい。
【0038】
本発明の一実施形態において、供給体の温度は、適宜に設定することができ、推定が行われる時点を含む期間において可変であってもよい。供給体の温度は、好ましくは５℃～４０℃である。このように、温度域により低温の領域まで含めることで、後述する第６のパラメータＥが供給体の状態をより強く示すようになるため、第６のパラメータＥを推定することが容易になることに加え、供給量の推定精度が向上する。
【0039】
〔第５のパラメータ〕
第５のパラメータｋｓは、基準となる温度における、養分供給量の変化における速度定数を示す。本発明の一実施形態において、基準となる温度は、適宜に選択され得る。基準となる温度は、これらに限定されないが、５℃、１０℃、１５℃、２０℃、２５℃、３０℃、または４０℃であり得る。
第５のパラメータｋｓは、基準となる温度における、養分供給量の変化における速度定数を示す。本発明の一実施形態において、基準となる温度は、適宜に選択され得る。基準となる温度は、これらに限定されないが、５℃、１０℃、１５℃、２０℃、２５℃、３０℃、または４０℃であり得る。
【0040】
〔第６のパラメータ〕
第６のパラメータＥは、速度定数の温度依存性を示す。本発明の一実施形態において、第６のパラメータＥは、供給体が養分を供給する反応の反応系のうち律速反応の活性化エネルギーに対応する値と考えられる。供給体の種類が多い場合、または供給体が養分を供給する反応が多段階反応である場合、反応系は複雑さを増す。しかしながら、当該反応系のうち、全体としての反応速度は律速反応の反応速度が支配的であるため、全体としての反応速度の温度依存性は律速反応の温度依存性、すなわち律速反応の活性化エネルギーが支配的であると考えられる。
第６のパラメータＥは、速度定数の温度依存性を示す。本発明の一実施形態において、第６のパラメータＥは、供給体が養分を供給する反応の反応系のうち律速反応の活性化エネルギーに対応する値と考えられる。供給体の種類が多い場合、または供給体が養分を供給する反応が多段階反応である場合、反応系は複雑さを増す。しかしながら、当該反応系のうち、全体としての反応速度は律速反応の反応速度が支配的であるため、全体としての反応速度の温度依存性は律速反応の温度依存性、すなわち律速反応の活性化エネルギーが支配的であると考えられる。
【0041】
〔代替パラメータ〕
本発明の一実施形態において、パラメータ群に含まれるパラメータは、上述の第１、第２、第３、第４のパラメータ、追加の第５、第６のパラメータのセットに限定されない。パラメータ群には、上述のパラメータのうち少なくとも１つのパラメータを数学的変換式で変換して得られるパラメータが、前記少なくとも１つのパラメータの代替として、含まれてもよい。
本発明の一実施形態において、パラメータ群に含まれるパラメータは、上述の第１、第２、第３、第４のパラメータ、追加の第５、第６のパラメータのセットに限定されない。パラメータ群には、上述のパラメータのうち少なくとも１つのパラメータを数学的変換式で変換して得られるパラメータが、前記少なくとも１つのパラメータの代替として、含まれてもよい。
【0042】
すなわち、パラメータ群には、パラメータ群に含まれる少なくとも１つのパラメータを数学的変換式で変換したパラメータが、パラメータ群に含まれる少なくとも１つのパラメータの代替として含まれてもよい。本発明の一実施形態に係る養分供給量推定方法を実施するユーザは、数学的変換式を適宜に選択することにより、得られる代替パラメータが示す特性を適宜に変換または追加することができる。これにより、ユーザは、供給体の状態および特性の概要を容易に理解できる、または直感的に知覚できるようになる。これにより、例えば供給体に関する情報の発信側と利用者側との間において、供給体の特徴を表す情報の伝達が容易であり、このような情報を供給体の流通上の指標とすることができる。このような代替パラメータを得るための数学的変換式の例は、後述する。
【0043】
以下に、本発明の一実施形態に係る養分供給量推定方法が包含し得る各工程の非限定的な構成を説明する。
【0044】
〔パラメータ設定工程〕
本発明の一実施形態に係る養分供給量推定方法は、後述する養分供給量推定工程の前に、パラメータ設定工程をさらに包含してもよい。パラメータ設定工程は、パラメータ群に含まれる少なくとも１つのパラメータを設定する工程であり得る。
本発明の一実施形態に係る養分供給量推定方法は、後述する養分供給量推定工程の前に、パラメータ設定工程をさらに包含してもよい。パラメータ設定工程は、パラメータ群に含まれる少なくとも１つのパラメータを設定する工程であり得る。
【0045】
少なくとも１つのパラメータを設定する方法の例には：（ｉ）ユーザが測定した供給量を参照して、パラメータを設定する方法；および（ｉｉ）任意のデータベースに記憶されたパラメータをユーザが参照して、パラメータを設定する方法；が含まれるが、これらに限定されない。パラメータのそれぞれは、同一の方法によって設定されてもよいし、別々の方法によって設定されてもよい。
【0046】
パラメータを設定する方法のうち、（ｉ）ユーザが、測定した供給量を参照して、パラメータを設定する方法について、説明する。以下、この方法を「方法（ｉ）」と略記することがある。
【0047】
方法（ｉ）は、供給体の収集および当該供給体が含有する成分の分析によって供給量を測定する、実験的方法であり得る。収集および分析は、当技術分野で公知の方法から適宜に選択され得る。
【0048】
方法（ｉ）において、測定された供給量を、パラメータ群を定数として含む式に対して入力して、パラメータをフィッティングすることによって、パラメータを設定することができる。フィッティングは、当技術分野で公知の方法で、達成され得る。フィッティングにおいて、パラメータ群に含まれるパラメータがいずれも未知である場合には、これらのパラメータすべてをフィッティングしてもよい。また、フィッティングにおいて、パラメータ群に含まれるパラメータのうちいずれか１つ以上が既知である場合には、既知であるパラメータの値を固定して、未知である残りのパラメータのみをフィッティングしてもよい。
【0049】
なお、フィッティングの条件は、当業者に公知の条件を設定することができる。例えば、フィッティング結果が発散する、または複数存在するもしくは存在しない場合には、パラメータのいずれかの初期値を所定の定数、例えば０または１として、フィッティングを行ってもよい。
【0050】
パラメータを設定する方法のうち、（ｉｉ）任意のデータベースに記憶されたパラメータをユーザが参照して、パラメータを設定する方法について、説明する。以下、この方法を「方法（ｉｉ）」と略記することがある。
【0051】
方法（ｉｉ）は、任意のデータベースに記憶されたパラメータをユーザが選択し、そのまま、または適宜に数学的変換を行い、パラメータを設定する方法であり得る。データベースの例には、コンピュータ読み取り可能な記憶媒体、およびインターネットなどの一般に閲覧可能な媒体が含まれる。
【0052】
代替的には、方法（ｉｉ）においては、記憶されたパラメータを選択するための情報がユーザにより入力され、次いで、入力された情報に基づいて自動的に選択されたパラメータが、パラメータとして設定されてもよい。一例として、ユーザは記憶されたパラメータを選択するための情報として、場の状態に関する情報を入力し、次いで、入力された情報を参照して、対応するパラメータが自動的に呼び出されてもよい。
【0053】
〔パラメータ推定工程〕
パラメータは、上述したように設定する方法の他にも、下記に示すように、推定してもよい。すなわち、本発明の一実施形態に係る養分供給量推定方法は、養分供給量推定工程の前に、一の養分供給体の状態と他の養分供給体の少なくとも１つの状態との間の類似性、および、前記一の養分供給体における前記パラメータ群に含まれる少なくとも１つのパラメータを参照して、他の養分供給体の少なくとも１つにおける前記パラメータ群に含まれる少なくとも１つのパラメータを推定するパラメータ推定工程をさらに包含してもよい。
パラメータは、上述したように設定する方法の他にも、下記に示すように、推定してもよい。すなわち、本発明の一実施形態に係る養分供給量推定方法は、養分供給量推定工程の前に、一の養分供給体の状態と他の養分供給体の少なくとも１つの状態との間の類似性、および、前記一の養分供給体における前記パラメータ群に含まれる少なくとも１つのパラメータを参照して、他の養分供給体の少なくとも１つにおける前記パラメータ群に含まれる少なくとも１つのパラメータを推定するパラメータ推定工程をさらに包含してもよい。
【0054】
パラメータ群に含まれるパラメータそれぞれは、養分供給体の状態に強く関連している。したがって、他の養分供給体が養分供給量推定方法の対象である場合、一の養分供給体の状態および対象である他の養分供給体の状態との間の類似性、および、一の養分供給体において設定された既知であるパラメータを参照することにより、実験的方法を伴わなくとも、推定により他の養分供給体におけるパラメータを設定することができる。
【0055】
パラメータ推定工程において、参照に用いられる養分供給体は、合計３つ以上であってもよい。例えば、パラメータが既知である２つ以上の養分供給体と、パラメータが未知である１つの養分供給体と、を参照に用いて、未知であるパラメータを推定してもよい。
【0056】
パラメータ推定工程において、ユーザは、任意の方法を用いて、一の養分供給体において設定された既知であるパラメータを取得することができる。一例として、ユーザは、コンピュータ読み取り可能な記憶媒体に記憶されるか、またはインターネットなどの一般に閲覧可能な媒体に記載された、既知であるパラメータを取得することができる。
【0057】
パラメータ推定工程において、参照される養分供給体の状態の例には：養分供給体の種類；養分供給体が含有する成分、および量または含有率；養分供給体が含まれる場が含む、養分供給体以外の成分；当該場の土壌の種類；当該場の水分条件；当該場に存在する微生物の活性が含まれる。ユーザは、当技術分野で公知の方法を用いて、一の養分供給体の状態、および他の養分供給体の状態に関する情報を取得してよい。一例として、過去の知見に基づいて、類似する養分供給体の種類毎の代表値（平均値など）を、第１のパラメータｎｉとして設定することができる。
【0058】
パラメータ推定工程において、パラメータの推定は、一の養分供給体の状態と他の養分供給体の少なくとも１つの状態との間の類似性が高いほどパラメータの類似性も高い、という仮定に基づき、行われ得る。一例として、養分供給体の種類に関して一の養分供給体の状態と他の養分供給体とを比較し、類似性が高い場合には、他の養分供給体のパラメータのうち第６のパラメータＥおよび第３のパラメータｄのうち少なくとも１つを一の養分供給体において設定された既知であるパラメータと類似する値に設定することができる。
【0059】
パラメータ推定工程において、推定され得るパラメータは、第１～第６のパラメータに限定されない。例えば、上述の代替パラメータが推定されてもよい。すなわち、パラメータ推定工程においては、一の養分供給体におけるパラメータ群に含まれる代替パラメータを参照して、他の養分供給体の少なくとも１つにおけるパラメータ群に含まれる代替パラメータを推定してもよい。
【0060】
〔取得工程〕
本発明の一実施形態に係る養分供給量推定方法は、養分供給量推定工程の前に、取得工程をさらに包含してもよい。取得工程は、パラメータを定数として含む式に対して変数として入力される値を取得する工程であり得る。
本発明の一実施形態に係る養分供給量推定方法は、養分供給量推定工程の前に、取得工程をさらに包含してもよい。取得工程は、パラメータを定数として含む式に対して変数として入力される値を取得する工程であり得る。
【0061】
変数として入力される値は、ユーザが所望する任意の時点を示す時間の値を含み、任意に供給体温度の値をさらに含む。
【0062】
〔養分供給量推定工程〕
本発明の一実施形態に係る養分供給量推定方法は、養分供給量を推定するためのパラメータ群を参照して、任意の時点における養分供給量を推定する養分供給量推定工程を包含する。
本発明の一実施形態に係る養分供給量推定方法は、養分供給量を推定するためのパラメータ群を参照して、任意の時点における養分供給量を推定する養分供給量推定工程を包含する。
【0063】
本発明の一実施形態において、パラメータ群を参照することは、パラメータ群を定数として含む式を用いることにより、達成され得る。当該式は、ユーザが所望する任意の時点を示す時間の値を変数として含み、当該変数が入力されることにより、所望される任意の時点における、推定される供給量を出力し得る。
【0064】
養分供給量推定工程において用いられる式は、上述したパラメータ群に含まれるパラメータが示す供給体の状態を考慮して、ユーザが適宜に設定することができる。一例として、用いられる式は、第１のパラメータｎｉを切片とする式であり、第２のパラメータＮに対して、第３のパラメータｄ、第４のパラメータｋ、および変数ｔを含む式によって算出される値を乗算することにより得られる値が前記切片に加算される式であり得る。本発明の一実施形態において用いられる式の非限定的な例である式（ｉ）を以下に示す。
【0065】
ｎ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ＋ｎｉ 式（ｉ）
式（ｉ）中、ｎ（単位：ｍｇ／１００ｍｇ 養分供給体）は、時間ｔ（単位：ｄａｙ）において養分供給体１００ｍｇから供給される養分供給量（単位：ｍｇ）を示す。ｎｉ（単位：ｍｇ／１００ｍｇ 養分供給体）は第１のパラメータであり、ｔ＝０において養分供給体１００ｍｇから供給される養分供給量を示す。Ｎ（単位：ｍｇ／１００ｇ）は第２のパラメータである。ｄ（単位：－）は第３のパラメータである。ｋ（単位：ｄａｙ－１）は第４のパラメータである。
式（ｉ）中、ｎ（単位：ｍｇ／１００ｍｇ 養分供給体）は、時間ｔ（単位：ｄａｙ）において養分供給体１００ｍｇから供給される養分供給量（単位：ｍｇ）を示す。ｎｉ（単位：ｍｇ／１００ｍｇ 養分供給体）は第１のパラメータであり、ｔ＝０において養分供給体１００ｍｇから供給される養分供給量を示す。Ｎ（単位：ｍｇ／１００ｇ）は第２のパラメータである。ｄ（単位：－）は第３のパラメータである。ｋ（単位：ｄａｙ－１）は第４のパラメータである。
【0066】
式（ｉ）は、第１のパラメータｎｉを切片とする式であり、第２のパラメータＮに対して、第３のパラメータｄ、第４のパラメータｋ、および変数ｔを含む項（ｉＡ）によって算出される値を乗算することにより得られる値が前記切片に加算され、ｎが算出される式である。
【0067】
［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ 項（ｉＡ）
【0068】
式（ｉ）は、積分形一次反応速度式、シグモイド式、またはこれらの変形によって得られる式であり、選択される式および変形の程度は、ｄの値に依存する。
【0069】
一例として、式（ｉ）において、第３のパラメータｄがｄ＝－１を満たす場合、ｎによって表される供給量の推移は変数ｔを含む積分形一次反応速度式である。このとき、供給体の種類が少なく、供給体が養分を供給する反応の反応系が単純であることが推察される。別の一例として、式（ｉ）において、第３のパラメータｄがｄ＝１を満たす場合、ｎによって表される供給量の推移は変数ｔを含むロジスティック式（シグモイド式の特殊形）である。このとき、供給体の種類が多く、供給体が養分を供給する反応の反応系が複雑であることが推察される。
【0070】
別の一例として、式（ｉ）において、第３のパラメータｄが－１＜ｄ＜１を満たす場合、第３のパラメータｄが小さいほど、積分形一次反応速度式に対する供給量の推移の類似性が高くなり、反応系が単純であることが推察される。また、第３のパラメータｄが大きいほど、ロジスティック式に対する供給量の推移の類似性が高くなり、反応系が複雑であることが推察される。
【0071】
なお、式（ｉ）において、第３のパラメータｄがとり得る値の範囲は、－１≦ｄ≦１の範囲に限定されず、ｄ＜－１、１＜ｄの範囲であってもよい。反応系の複雑さの観点からは、反応系を構成する反応それぞれに対して第１、２、４、５および６のパラメータが設定されることが従来一般的である。しかしながら、本発明の一態様によれば、第３のパラメータｄがとり得る値の範囲が広い、具体的にはｄ＜－１、１＜ｄの範囲も取り得ることにより、複雑な反応系を少ない種類のパラメータを用いて簡易に取り扱うことができる。
【0072】
また、ｄ＝０の場合は、式（ｉ）の項（ｉＡ）の値は数学的には計算できない。しかしながら、ｄを正または負の値から０に近づけた極限において項（ｉＡ）の値が一致するため、ｄ＝０の場合の項（ｉＡ）の値は、その極限における値と等しいものとして扱ってもよい。また、項（ｉＡ）の下記底（ｉＢ）が負の値である場合、項（ｉＡ）の値は実数でないことがある。この場合には、底（ｉＢ）の値を０として扱い、項（ｉＡ）を含む式（ｉ）を計算するか、または、底（ｉＢ）の値に、まず－１を掛けてから正の値としてから項（ｉＡ）の値を算出し、次いで、項（ｉＡ）値に－１を掛けて、式（ｉ）を計算してもよい。
【0073】
［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］ 底（ｉＢ）
【0074】
本発明の一実施形態において、第４のパラメータｋは、第５のパラメータｋｓと、第６のパラメータＥとを参照して、算出されてもよい。このようにして、第４のパラメータｋが算出されることにより、本発明の一実施形態に係る養分供給量推定方法は、温度変化に応じて、供給量を推定することができる。第５のパラメータｋｓと、第６のパラメータＥとを参照することは、これらのパラメータを定数として含む式を用いることにより、達成され得る。当該式は、供給体温度を示す値を変数として含み、当該変数が入力されることにより、入力された供給体温度における、第４のパラメータｋを出力し得る。本発明の一実施形態において第４のパラメータｋを算出するために用いられる式の非限定的なアレニウス式を利用した例を下記式（ｉｉ）に示す。
【0075】
ｋ＝ｋｓ・ｅｘｐ｛Ｅ／Ｒ・（１／Ｔｓ－１／Ｔ）｝ 式（ｉｉ）
式（ｉｉ）中、ｋｓ（単位：ｄａｙ－１）は第５のパラメータであり、基準となる温度２５℃における、供給量の変化における速度定数を示す。Ｅ（単位：ｋＪ・ｍｏｌ－１）は第６のパラメータである。Ｒ（８．３１ｋＪ・ｍｏｌ－１・Ｋ（絶対温度））は気体定数を示す。Ｔｓ（単位：Ｋ（絶対温度））は、基準となる温度２５℃に対応する絶対温度を示す。Ｔ（単位：Ｋ（絶対温度））は、供給体温度を示す。
式（ｉｉ）中、ｋｓ（単位：ｄａｙ－１）は第５のパラメータであり、基準となる温度２５℃における、供給量の変化における速度定数を示す。Ｅ（単位：ｋＪ・ｍｏｌ－１）は第６のパラメータである。Ｒ（８．３１ｋＪ・ｍｏｌ－１・Ｋ（絶対温度））は気体定数を示す。Ｔｓ（単位：Ｋ（絶対温度））は、基準となる温度２５℃に対応する絶対温度を示す。Ｔ（単位：Ｋ（絶対温度））は、供給体温度を示す。
【0076】
なお、養分供給量推定工程においては、パラメータ群に含まれるパラメータのうち少なくとも１つの代わりに、代替パラメータを参照して、供給量を推定してもよい。ここで、代替パラメータは、パラメータ群に含まれる少なくとも１つのパラメータを数学的変換式で変換したパラメータである。
【0077】
代替パラメータを得るための数学的変換の例として、温度が基準となる温度２５℃の条件、すなわちｋ＝ｋｓの条件の例を挙げる。このとき、式（ｉ）に対して、ｎ＝０．２５Ｎ＋ｎｉを入力すると、下記式（ｉｉｉ）が得られる。
【0078】
０．２５Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ２５％｝］－１／ｄ＋ｎｉ 式（ｉｉｉ）
【0079】
式（ｉｉｉ）は、第３のパラメータｄおよび第５のパラメータｋｓから、代替パラメータｔ２５％を得るための変換式である。代替パラメータｔ２５％は、供給量の最大変化量に対する供給量の変化量の割合が２５％に達する時間を示す。したがって、ユーザは、代替パラメータｔ２５％を参照することにより、供給体から供給される養分が最大変化量の２５％に達する時間を容易に理解することができる。
【0080】
同様にして、式（ｉ）に対して、ｎ＝０．５Ｎ＋ｎｉ、ｎ＝０．７５Ｎ＋ｎｉをそれぞれ入力すると、供給体から供給される養分が最大変化量の５０％、７５％に達する時間を示す代替パラメータｔ５０％、ｔ７５％が得られる。ユーザはｔ２５％、ｔ５０％、およびｔ７５％をそれぞれ得ることにより、供給体から供給される養分が最大変化量の２５％、５０％、および７５％に達する時間を理解でき、植物への養分供給の特徴を容易に理解することができる。
【0081】
さらに、ｔ２５％、ｔ５０％、ｔ７５％のうち少なくとも２つを用いた数学的変換を行うことにより、第３のパラメータｄおよび第５のパラメータｋｓが算出され得る。すなわち、ユーザの理解を促進する代替パラメータは、養分供給量推定工程において用いられる式に含まれるパラメータとしても、機能し得る。
【0082】
代替パラメータは、上記の例に限定されず、任意の数学的変換によって設定される。別の例として、ｔ２５％、ｔ５０％、およびｔ７５％から、中央値（ｔ５０％）および四分位偏差｛（ｔ７５％－ｔ２５％）／２｝を算出できる。これらを用いれば、ユーザは、中央値±四分位偏差の期間に半分の養分が供給されることを理解でき、大まかな養分の供給時期を理解できる。これら２つのパラメータからも、可逆的に、第３のパラメータｄおよび第５のパラメータｋｓを算出することができる。
【0083】
〔養分供給体の利用方法〕
本発明の一態様は、養分供給体の利用方法に関する。本発明の一実施形態に係る養分供給体の利用方法は、（ｉ）本発明の一実施形態に係る養分供給量推定方法によって推定した養分供給量、ならびに、（ｉｉ）当該養分供給量の推定に用いたパラメータ群に含まれる少なくとも１つのパラメータ、の少なくとも一方を参照して養分供給体の状態を推定する工程を包含する。
本発明の一態様は、養分供給体の利用方法に関する。本発明の一実施形態に係る養分供給体の利用方法は、（ｉ）本発明の一実施形態に係る養分供給量推定方法によって推定した養分供給量、ならびに、（ｉｉ）当該養分供給量の推定に用いたパラメータ群に含まれる少なくとも１つのパラメータ、の少なくとも一方を参照して養分供給体の状態を推定する工程を包含する。
【0084】
本発明の一実施形態に係る養分供給体の利用方法においては、上述した供給量およびパラメータの少なくとも一方を参照することにより、供給体の状態を推定することができる。ここで、供給体の状態の例には、限定するものではないが、供給量の過不足、および供給時期の遅速が含まれる。したがって、本発明の一実施形態に係る養分供給体の利用方法によれば、ユーザは、推定された供給体の状態に基づいて、最適な供給体の選択、ならびに、供給体の施用量、施用時期および施用期間等の設定が可能である。一例として、ユーザは、推定された供給体の状態に基づいて、場に含まれる供給体をそのまま利用するか、または、当該場に対して供給体の追加を行うか、等の計画を設定することができる。
【0085】
本発明の一実施形態に係る養分供給体の利用方法において、推定された養分供給量および推定に用いられたパラメータ群に含まれる少なくとも１つのパラメータは、任意の媒体に記載または記憶された養分供給量およびパラメータが参照されてもよい。媒体の例には、養分供給体の包装体、および説明書が含まれるが、これらに限定されない。また、任意の媒体には、養分供給量およびパラメータそのものではなく、これらを閲覧可能な媒体へのアクセス方法、例えばホームページアドレスまたはＱＲコード（登録商標）等、が記載または記憶され、参照されてもよい。
【0086】
本発明の一実施形態に係る養分供給体の利用方法は、推定する工程において、パラメータ群に含まれる少なくとも１つのパラメータから算出される、養分供給体から植物に供給される養分の任意の養分供給量が得られるまでの時間を参照して、養分供給体の状態を推定してもよい。本発明の一実施形態によれば、任意の養分供給量が得られるまでの時間を参照することにより、所望の量の養分が植物に供給されるまでに必要とされる時間を考慮して、養分供給体の状態を推定することができる。したがって、ユーザは、より適切な計画を設定することができる。
【0087】
〔ソフトウェアによる実現例〕
本発明の一実施形態に係る養分供給量推定方法、または、本発明の一実施形態の養分供給体の利用方法、の各工程を実行する装置についても本発明の範疇に含まれる。そして、装置の機能は、当該装置としてコンピュータを機能させるためのプログラムであって、当該装置の各制御ブロックとしてコンピュータを機能させるためのプログラムにより実現することができる。すなわち、本発明の一実施形態に係る養分供給量推定方法、または、本発明の一実施形態の養分供給体の利用方法、をコンピュータに実行させるためのプログラムについても本発明の範疇に含まれる。
本発明の一実施形態に係る養分供給量推定方法、または、本発明の一実施形態の養分供給体の利用方法、の各工程を実行する装置についても本発明の範疇に含まれる。そして、装置の機能は、当該装置としてコンピュータを機能させるためのプログラムであって、当該装置の各制御ブロックとしてコンピュータを機能させるためのプログラムにより実現することができる。すなわち、本発明の一実施形態に係る養分供給量推定方法、または、本発明の一実施形態の養分供給体の利用方法、をコンピュータに実行させるためのプログラムについても本発明の範疇に含まれる。
【0088】
この場合、上記装置は、上記プログラムを実行するためのハードウェアとして、少なくとも１つの制御装置（例えばプロセッサ）と少なくとも１つの記憶装置（例えばメモリ）を有するコンピュータを備えている。この制御装置と記憶装置により上記プログラムを実行することにより、上述した養分供給量推定方法、または養分供給体の利用方法の各工程が実現される。
【0089】
上記プログラムは、一時的ではなく、コンピュータ読み取り可能な、１または複数の記録媒体に記録されていてもよい。この記録媒体は、上記装置が備えていてもよいし、備えていなくてもよい。後者の場合、上記プログラムは、有線または無線の任意の伝送媒体を介して上記装置に供給されてもよい。
【0090】
また、上記各制御ブロックの機能の一部または全部は、論理回路により実現することも可能である。例えば、上記各制御ブロックとして機能する論理回路が形成された集積回路も本発明の範疇に含まれる。この他にも、例えば量子コンピュータにより上記各制御ブロックの機能を実現することも可能である。
【0091】
さらに、上述した養分供給量推定方法、または養分供給体の利用方法の各工程は、ＡＩ（Artificial Intelligence：人工知能）に実行させてもよい。この場合、ＡＩは上記制御装置で動作するものであってもよいし、他の装置（例えばエッジコンピュータまたはクラウドサーバ等）で動作するものであってもよい。
【0092】
〔推定システム〕
さらに、上述のデータを、養分供給体と関連付けて記憶するデータベースと、養分供給量を推定する推定装置と、を備えた、推定システムも本発明の範疇に含まれる。すなわち、本発明の一実施形態に係る推定システムは、（ｉ）本発明の一実施形態に係る養分供給量推定方法によって推定した養分供給量、（ｉｉ）本発明の一実施形態に係る養分供給体の利用方法によって推定した養分供給体の状態、ならびに、（ｉｉｉ）当該養分供給量の推定に用いたパラメータ群に含まれるパラメータの値のいずれか、のうちの少なくとも１つを、養分供給体と関連付けて記憶するデータベースと、データベースに記憶されたデータを参照して、被験養分供給体の養分供給量を推定する推定装置と、を備える。
さらに、上述のデータを、養分供給体と関連付けて記憶するデータベースと、養分供給量を推定する推定装置と、を備えた、推定システムも本発明の範疇に含まれる。すなわち、本発明の一実施形態に係る推定システムは、（ｉ）本発明の一実施形態に係る養分供給量推定方法によって推定した養分供給量、（ｉｉ）本発明の一実施形態に係る養分供給体の利用方法によって推定した養分供給体の状態、ならびに、（ｉｉｉ）当該養分供給量の推定に用いたパラメータ群に含まれるパラメータの値のいずれか、のうちの少なくとも１つを、養分供給体と関連付けて記憶するデータベースと、データベースに記憶されたデータを参照して、被験養分供給体の養分供給量を推定する推定装置と、を備える。
【0093】
本発明の一実施形態に係る推定システムにおいて、データベースは、ユーザによって、アクセスされ得、（ｉ）養分供給量、（ｉｉ）養分供給体の状態、ならびに（ｉｉｉ）パラメータの値のいずれか、のうちの少なくとも１つが、書き込みおよび読み込みの少なくとも１つが行われ得る。したがって、ユーザは、自身または他のユーザが実施した推定によるデータやパラメータの値を参照して、所望の被験養分供給体の養分供給量を推定することができる。
【0094】
本発明の一実施形態に係る推定システムは、ユーザの操作を参照して、またはユーザが実施した推定の結果を参照して、供給量を数学的に推定するために用いられる式を変形する手段を備えてもよい。一例として、データベースは、ユーザによって実施された供給量の推定結果と、対応する供給量の実際の結果とを記憶し、推定結果と実際の結果の差異が小さくなるように、用いられる式を適宜に変形する手段を備えてもよい。このような構成によれば、本発明の一実施形態に係る推定システムは、推定の精度をより改善することができる。
【0095】
〔変形例〕
本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
【0096】
例えば、供給量を数学的に推定するために用いられる式は、養分供給量に関する情報の所望の様式に応じて、適宜に設定することができる。一例として、以下に表す式（ｉｖ）が用いられてもよい。
【0097】
ｎ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ 式（ｉｖ）
式（ｉｖ）は、式（ｉ）においてｎｉ＝０と設定した場合の式であり、ｎは養分供給量の経時的差分量を示す。すなわち、ｎｉの値を適宜に設定することにより、経時的差分量を推定することも、本発明の一態様の範疇である。
式（ｉｖ）は、式（ｉ）においてｎｉ＝０と設定した場合の式であり、ｎは養分供給量の経時的差分量を示す。すなわち、ｎｉの値を適宜に設定することにより、経時的差分量を推定することも、本発明の一態様の範疇である。
【0098】
別の一例として、以下に表す式（ｖ）または式（ｖｉ）が用いられてもよい。
【0099】
ｎ／Ｎ＝［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ＋ｎｉ／Ｎ 式（ｖ）
ｎ／Ｎ＝［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ 式（ｖｉ）
式（ｖ）は、供給された養分の割合ｎ／Ｎを算出する式である。また、式（ｖｉ）は、式（ｖ）においてｎｉ＝０と設定した場合の式であり、ｎ／Ｎは養分供給量の経時的差分量の割合を示す。式（ｖ）および式（ｖｉ）において、Ｎ＝１と設定し、ｎが割合に対応するように、各式を設定してもよい。すなわち、式を適宜に変形することにより、養分供給体あたりの供給量の比率を推定することも、本発明の一態様の範疇である。
ｎ／Ｎ＝［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ 式（ｖｉ）
式（ｖ）は、供給された養分の割合ｎ／Ｎを算出する式である。また、式（ｖｉ）は、式（ｖ）においてｎｉ＝０と設定した場合の式であり、ｎ／Ｎは養分供給量の経時的差分量の割合を示す。式（ｖ）および式（ｖｉ）において、Ｎ＝１と設定し、ｎが割合に対応するように、各式を設定してもよい。すなわち、式を適宜に変形することにより、養分供給体あたりの供給量の比率を推定することも、本発明の一態様の範疇である。
【0100】
別の一例として、用いられる式において、変数として入力される値は、時間の値に加えて、場の水分ポテンシャルおよびｐＨなどの、養分供給に影響を与える可能性がある条件に対応する補正を含んでもよい。
【実施例】
【0101】
本発明の一実施例について以下に説明する。
【0102】
本実施例において、養分供給量として、土壌および堆肥から供給される無機態窒素の量を推定した。上記の非特許文献１に記載され、実験的方法によって測定された、土壌および堆肥に含まれる無機態窒素量を参照して、パラメータ群を設定した。次いで、設定されたパラメータ群を用いて無機態窒素量を推定し、推定された無機態窒素量（推定値）と参照された無機態窒素量（実験値）とを図１に示すグラフ中にプロットし、比較した。図１は、実施例の無機態窒素量の推定結果を示す図である。
【0103】
〔実施例１〕
非特許文献１の「表１」に記載された温度ごとの土壌および堆肥を施用した土壌に含まれる無機態窒素量の値と、式（ｉ）および式（ｉｉ）による推定値との差の残差平方和が最小になるように、第１のパラメータ、第２のパラメータ、第３のパラメータ、第５のパラメータ、および第６のパラメータを設定した。なお、本実施例において、無機態窒素量とは、アンモニウム態窒素の量と硝酸態窒素の量との合計量を指す。
非特許文献１の「表１」に記載された温度ごとの土壌および堆肥を施用した土壌に含まれる無機態窒素量の値と、式（ｉ）および式（ｉｉ）による推定値との差の残差平方和が最小になるように、第１のパラメータ、第２のパラメータ、第３のパラメータ、第５のパラメータ、および第６のパラメータを設定した。なお、本実施例において、無機態窒素量とは、アンモニウム態窒素の量と硝酸態窒素の量との合計量を指す。
【0104】
ｎ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋ・ｔ｝］－１／ｄ＋ｎｉ 式（ｉ）
式（ｉ）中、ｎは時間ｔにおいて土壌に含まれる無機態窒素量の推定値を示す。ｎｉは第１のパラメータであり、ｔ＝０において土壌に含まれる無機態窒素量を示す。Ｎは第２のパラメータであり、無機態窒素量の最大変化量、すなわち十分に時間が経過した時点における無機態窒素量の変化量を示す。ｄは第３のパラメータであり、無機態窒素量の推移を表す式（ｉ）の形状を示す。ｋは第４のパラメータであり、土壌温度Ｔ（単位：Ｋ（絶対温度））における無機態窒素量の変化における速度定数を示す。
式（ｉ）中、ｎは時間ｔにおいて土壌に含まれる無機態窒素量の推定値を示す。ｎｉは第１のパラメータであり、ｔ＝０において土壌に含まれる無機態窒素量を示す。Ｎは第２のパラメータであり、無機態窒素量の最大変化量、すなわち十分に時間が経過した時点における無機態窒素量の変化量を示す。ｄは第３のパラメータであり、無機態窒素量の推移を表す式（ｉ）の形状を示す。ｋは第４のパラメータであり、土壌温度Ｔ（単位：Ｋ（絶対温度））における無機態窒素量の変化における速度定数を示す。
【0105】
ｋ＝ｋｓ・ｅｘｐ｛Ｅ／Ｒ・（１／Ｔｓ－１／Ｔ）｝ 式（ｉｉ）
式（ｉｉ）中、ｋｓは第５のパラメータであり、２５℃における無機態窒素量の変化における速度定数を示す。Ｅは第６のパラメータであり、速度定数の温度依存性を示す。Ｔｓは基準となる温度、すなわち２５℃に対応する絶対温度を示す。Ｒは気体定数（８．３１ｋＪ・ｍｏｌ－１・Ｋ（絶対温度））である。
式（ｉｉ）中、ｋｓは第５のパラメータであり、２５℃における無機態窒素量の変化における速度定数を示す。Ｅは第６のパラメータであり、速度定数の温度依存性を示す。Ｔｓは基準となる温度、すなわち２５℃に対応する絶対温度を示す。Ｒは気体定数（８．３１ｋＪ・ｍｏｌ－１・Ｋ（絶対温度））である。
【0106】
各パラメータの設定は、条件（１）肥料が施用されていない土壌のみ、（２）牛豚鶏糞堆肥が施用された土壌、（３）牛糞堆肥が施用された土壌、（４）鶏糞が施用された土壌、について行われた。また、各パラメータの設定は、堆肥の効果を理解するため、堆肥の効果の差が大きかった条件（１）および条件（４）を用いて、（１）の無機態窒素量の値と、（４）の無機態窒素量の値との差、すなわち、条件（５）土壌の効果を除いた鶏糞堆肥のみ、についても行われた。なお、各条件（１）～（５）について、土壌温度は、１０℃、２０℃、および３０℃に設定された。
【0107】
下記表１に、非特許文献１の「表１」に記載され、本実施例の各パラメータの設定に用いられた温度別無機態窒素量の値（２連）を示す。
【0108】
【表1】
【0109】
設定されたパラメータを参照した推定の結果として得られた推定値は、実験値に近い値であった。また、第６のパラメータは、土壌および堆肥などの種類に依らず、条件（１）～（５）において、近い値となった。
【0110】
〔実施例２〕
そこで、さらに、すべての条件（１）～（５）において第６のパラメータの値が同一であると条件付けして、式（ｉ）および式（ｉｉ）による推定値と実験値との差の残差平方和が最小になるように、第１のパラメータ、第２のパラメータ、第３のパラメータ、第５のパラメータ、および第６のパラメータを設定した。
そこで、さらに、すべての条件（１）～（５）において第６のパラメータの値が同一であると条件付けして、式（ｉ）および式（ｉｉ）による推定値と実験値との差の残差平方和が最小になるように、第１のパラメータ、第２のパラメータ、第３のパラメータ、第５のパラメータ、および第６のパラメータを設定した。
【0111】
表２に、実施例２において設定された各パラメータの結果を示す。
【0112】
【表2】
【0113】
図１に、無機態窒素量の実験値と、実施例２において設定された各パラメータを含む式（ｉ）による推定値とを示す。図１において、点は実験値のデータを示し、曲線は式（ｉ）による推定値のデータを示す。ただし、条件（５）では、条件（４）の推定値から条件（１）の推定値を引いた値を実験値として示している。
【0114】
〔実施例３〕
さらに、無機態窒素量の変化量の、最大変化量に対する割合が、２５％、５０％または７５％に達する時間ｔ２５％、ｔ５０％またはｔ７５％を示す代替パラメータを、第３のパラメータｄおよび第４のパラメータを用いて算出した。具体的には、代替パラメータｔ２５％、ｔ５０％またはｔ７５％のそれぞれは、下記式（ｉｉｉ）、式（ｖｉｉ）、または式（ｖｉｉｉ）から、算出された。
さらに、無機態窒素量の変化量の、最大変化量に対する割合が、２５％、５０％または７５％に達する時間ｔ２５％、ｔ５０％またはｔ７５％を示す代替パラメータを、第３のパラメータｄおよび第４のパラメータを用いて算出した。具体的には、代替パラメータｔ２５％、ｔ５０％またはｔ７５％のそれぞれは、下記式（ｉｉｉ）、式（ｖｉｉ）、または式（ｖｉｉｉ）から、算出された。
【0115】
０．２５Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ２５％｝］－１／ｄ＋ｎｉ 式（ｉｉｉ）
０．５０Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ５０％｝］－１／ｄ＋ｎｉ 式（ｖｉｉ）
０．７５Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ７５％｝］－１／ｄ＋ｎｉ 式（ｖｉｉｉ）
０．５０Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ５０％｝］－１／ｄ＋ｎｉ 式（ｖｉｉ）
０．７５Ｎ＋ｎｉ＝Ｎ・［１＋ｄ・ｅｘｐ｛ｅｘｐ（ｄ＋１）－１－ｋｓ・ｔ７５％｝］－１／ｄ＋ｎｉ 式（ｖｉｉｉ）
【0116】
表３に、実施例３において算出された代替パラメータを示す。
【0117】
【表3】
【0118】
〔比較例１〕
非特許文献１の「表１」に記載された温度ごとの土壌に含まれる無機態窒素量の値と、非特許文献１の「モデル式２」による推定値との差の残差平方和が最小になるように、モデル式２に含まれるパラメータを設定した。設定された各パラメータは、非特許文献１の「表２」に記載された結果と同様であった。非特許文献１の「表２」に記載された各パラメータを、下記表４に示す。各パラメータの詳細な意味は、非特許文献１を参照されたい。
Ｎ＝Ｎ０［１－ｅｘｐ（－ｋ・ｔ）］－Ｎ０ｉｍ［１－ｅｘｐ（－ｋｉｍ・ｔ）］＋ｂ （非特許文献１の「モデル式２」）
非特許文献１の「表１」に記載された温度ごとの土壌に含まれる無機態窒素量の値と、非特許文献１の「モデル式２」による推定値との差の残差平方和が最小になるように、モデル式２に含まれるパラメータを設定した。設定された各パラメータは、非特許文献１の「表２」に記載された結果と同様であった。非特許文献１の「表２」に記載された各パラメータを、下記表４に示す。各パラメータの詳細な意味は、非特許文献１を参照されたい。
Ｎ＝Ｎ０［１－ｅｘｐ（－ｋ・ｔ）］－Ｎ０ｉｍ［１－ｅｘｐ（－ｋｉｍ・ｔ）］＋ｂ （非特許文献１の「モデル式２」）
【0119】
ただし、非特許文献１の「モデル式２」中、ｋｉｍ、ｋ、ｔは、以下のように算出される。
ｋｉｍ＝Ｂｅｘｐ（－Ｅａｉｍ／ＲＴ）
ｋ＝Ｃｅｘｐ（－Ｅａ／ＲＴ）
ｋｉｍ＝Ｂｅｘｐ（－Ｅａｉｍ／ＲＴ）
ｋ＝Ｃｅｘｐ（－Ｅａ／ＲＴ）
【0120】
【表4】
【0121】
なお、非特許文献１の「表２」または「表４」に合わせて、表４中のＮ０、およびＮ０ｉｍは、条件（１）～（４）については、５０ｇ中の土壌に対する量（単位：ｍｇ）として表され、条件（５）については、施用した鶏糞堆肥の重量百分率として表される。
【0122】
〔議論〕
表２に示すように、速度定数の温度依存性を示す第６のパラメータＥについて、実施例１では、各条件（１）～（５）において近い値となったため、同じ値を設定することができた。このことは、土壌に含まれる微生物によって媒介されると考えられる反応系において、無機態窒素の生成のための律速反応が堆肥の種類に依らず同一であるため、当該律速反応の活性化エネルギーに対応する第６のパラメータもまた同一の値になることを示唆する。対して、表４に示すように、比較例１では、各条件（１）～（５）において、温度依存性を示すＥａおよびＥａｉｍのばらつきが大きかった。このことは、類似する意味を有するパラメータが複数含まれるモデル式においては、当該パラメータ同士の干渉が生じるため、最終的に算出される推定値Ｎの値が正確であっても、個々のパラメータの値自体に対して意味付けすることが難しいことがあることを示す。
表２に示すように、速度定数の温度依存性を示す第６のパラメータＥについて、実施例１では、各条件（１）～（５）において近い値となったため、同じ値を設定することができた。このことは、土壌に含まれる微生物によって媒介されると考えられる反応系において、無機態窒素の生成のための律速反応が堆肥の種類に依らず同一であるため、当該律速反応の活性化エネルギーに対応する第６のパラメータもまた同一の値になることを示唆する。対して、表４に示すように、比較例１では、各条件（１）～（５）において、温度依存性を示すＥａおよびＥａｉｍのばらつきが大きかった。このことは、類似する意味を有するパラメータが複数含まれるモデル式においては、当該パラメータ同士の干渉が生じるため、最終的に算出される推定値Ｎの値が正確であっても、個々のパラメータの値自体に対して意味付けすることが難しいことがあることを示す。
【0123】
本実施例で示したように、第６のパラメータＥは、複数の反応系の間で一定値とすることができることがある。第６のパラメータＥは、律速反応の活性化エネルギーと対応するパラメータであるから、類似の反応系を取り扱う場合には、第６のパラメータＥを一定値として扱えると推察される。したがって、このように温度依存性を第６のパラメータ１つを用いて示し、温度の影響を第６のパラメータによって除くことにより、温度の影響がない定温条件下のように、養分供給量の推定を容易に取り扱うことができる。
【0124】
次いで、表２に示した第１のパラメータｎｉ、または「第２のパラメータＮ＋第１のパラメータｎｉ」の値は、図１に示したように、それぞれｔ＝０のときの窒素供給量、または十分に時間が経ったときの窒素供給量に近い値となっている傾向があった。このことは、第１のパラメータｎｉ、および第２のパラメータＮによって、供給量に関する情報を大まかに把握できることを示す。
【0125】
さらに、条件（１）、（４）および（５）の比較について、議論する。表２に示すように、実施例１では、条件（４）において設定された第１のパラメータｎｉ、第２のパラメータＮの値から、条件（１）において設定された当該パラメータの値を減算した値が、条件（５）において設定された当該パラメータの値と近似するという傾向があった。このことは、第１のパラメータｎｉがｔ＝０において土壌または堆肥に含まれる無機態窒素量に好適に対応し、第２のパラメータＮが無機態窒素量の最大変化量に好適に対応していることを示す。このことはさらに、一の養分供給体に対応する既知のパラメータから、他の養分供給体に対応する未知のパラメータを推定可能であることを示している。
【0126】
なお、比較例１の結果を示した表４において、ｂ、Ｎ０、およびＮ０ｉｍは、土壌５０ｇあたりの量として表されており（ただし、条件（５）においては鶏糞堆肥の重量百分率として表されているため、条件（５）はここでの議論の対象から除かれる）、土壌１００ｇあたりの量で表している表２および図１と比較するには、表４の値を２倍する必要があることに留意されたい。表４のパラメータｂについては、対応する第１のパラメータｎｉと同様に、ｔ＝０のときの窒素供給量を示す傾向が見られた。しかしながら、第２のパラメータＮと対応する表４のパラメータＮ０を用いて算出される「パラメータｂ＋パラメータＮ０」の値は、十分に時間が経ったときの窒素供給量との関係性が不明瞭であった。そのため、パラメータｂおよびパラメータＮ０を用いても、供給量に関する情報を大まか把握しにくかった。
【0127】
また、無機態窒素量の変化量の、最大変化量に対する割合が、所定の値に達する時間を代替パラメータが示すように、設定することができた。一例として、表３に示すように、条件（５）土壌を除いた鶏糞堆肥のみ、において、鶏糞堆肥から無機態窒素量へと変化する割合は、２５℃の場合、施用から６０日後に５０％に達することが、直感的に理解できる。
【産業上の利用可能性】
【0128】
本発明は、営農支援に利用することができる。
【図1】