Summary
この作品は、水質改善目標の指定された一連の農業保全慣行の最適な（最も低コストの）配置を解くために、進化的アルゴリズムを用いた最適化コンポーネントを使用した水質モデルの統合を示します。ソリューションは、トレードオフの明確な定量化を可能にする、マルチ客観的なアプローチを使用して生成されます。
Abstract
風景渡る特定の水質目標を達成するためのターゲティング保全実践投資のコスト効率の高い（ すなわち 、最も低コストの）方法を見つけることは流域管理における最も重要である。流域コンテキスト（ 例えば 、5,12,20）で最も低コストの解決策を見つけるの伝統的な経済学の方法は、外部への影響を正確に生成された敷地内の汚染の割合として記述することができることを想定しています。このようなアプローチは、汚染源の影響はしばしば複雑な生物物理学的プロセスによって決定される流域における実際の汚染プロセスを代表することはほとんどありません。現代の物理ベース、空間的に分散水文学的シミュレーションモデルの使用は、プロセス表現の面でリアリズムの大きい程度を可能にしますが、モデルは最適化の不可欠な一部となるシミュレーション最適化フレームワークの開発が必要である。
進化進アルゴリズムは流域シミュレーション·最適化問題のコンビナトリアル性質に対処するとフル水質モデルの使用を可能にすることが特に有用な最適化ツールであるように見えます。進化的アルゴリズムは、候補解として、流域の保全慣行の特定の空間的配分を扱い、最適化の目標に対して改善を見つけるために反復的に選択、交叉および突然変異の確率的演算子を適用する解候補の集合（集団）を利用する。この場合の最適化の目的は同時に保全慣行のコストを最小限に抑え、流域における非点、源汚染を最小にするためである。研究の最近の拡大とセットは同様の方法を使用しようとすると広義の進化的最適化手法3,4,9,10,13-15,17-19,22,23,25と水質モデルを統合しています。このアプリケーションでは、Rabotyagovらの適切に従うプログラムを実証ACHと保全の実践と保全の実践と、ユーザーが指定した水質のコストとの間の完全なトレードオフのフロンティアを検索するため、その費用の現代および一般SPEA2 26多目的進化アルゴリズムとSWATの水質モデル7を使用し、ユーザが指定したセットを統合目的。フロンティアは、様々な水質改善の目標に関連するコストの完全な範囲を提示することによって、流域管理者が直面してトレードオフを定量化する。プログラムは、指定された水質改善の目標と保全慣行の最適化された配置のマップの生産を達成する流域の構成を選択することができます。
Protocol
1。流域モデルを準備して、最適化のための入力データを提供
- i_SWATデータベースを作成
- "回転子"と呼ばれるプログラムを使って、土壌、気象、管理や肥料などの複数の入力データベースからデータベースを構築します。
- あるいは、既存のSWATの実行は（おそらくArcSWATまたはAVSWATで作成された）i_SWAT.exeでインポートすることができます。この場合、プログラムは "swat_rewrite"が管理またはフィールド·レベルのデータに基づいて、他のHRU情報を交換するために使用することができます。
- SWATモデルのキャリブレーションと検証は、この時点で実行する必要があります。 Jhaはらによって説明されているように、このアライグマ川流域のEAモデリングフレームワーク内に組み込まれ、SWAT（バージョン2005）のモデルはもともと最大総デイリー負荷（TMDL）研究のためにキャリブレーションされ、検証されました。 （2010）。 SWATモデルの更なる校正および検証はdescriとして、ラクーン川流域マスタープランの開発を支援するために実施されたこの研究のために使用されたSWATのモデルですAgren社（2011）、のベッド。
- SWAT2005GA.exe呼ばSWAT2005.exeの修正版を使用してください。
- "流域presets.csv"ファイルの準備 - これはアライグマ流域に固有の設定の多くを保存するテキストファイルです。それはワンクリックで分水嶺のために複数のコントロールとデータフィールドを設定するには、次のGeneticiSWATとMapSWATによって読み取られます。
- 保全の実践要素の費用を準備します。この例では、これらは、テーブルに格納されています[実践費用 - カウンティーアライグマ]データベース内の "実践コストsubbasin Josh.mdbによって"。
候補解の合計コストは、流域単位（ "水文応答ユニット"、またはHRUs）に適用保全慣行のコストの合計です。最適化プログラムは、流域内のすべての農地HRUの保全慣行の特定のセットから単一の保全の実践の最適割り当てを考慮します。のHRUに割り当てられている可能性保全の実践のETSは、対立遺伝子セットと呼ばれます。
- SWATのフォルダを作成します。この実行では、16個のCPUがSWAT2005GA.exeの16コピーは16個のフォルダで実行された、つまり、使用された（同じ数が少ないCPUを、 例えばを持つシステムに適用され、4つのフォルダは、 "クアッドコア"プロセッサ用に作成されるべきです。）
2。最適化パラメータを選択する
- 最適化は、 "GeneticISWAT"と呼ばれるプログラムによって制御されます。最適化、オープンGeneticISWAT.exeを実行します。
- 次に "開く"、 "ファイル"に移動し、i_SWATデータベース "アライグマGA.mdb"を選択します。
- SWATモデルの実行可能ファイル（SWAT2005GA.exe）にパスを割り当てるには、 "ファイル"を選び、 "設定"に進みます。
- "実行"を選び、 "対立遺伝子セットを選択"に進みます。このステップでは、最適化で使用される保全の実践の組み合わせを決定します。この実行では、対立遺伝子のセット＃14は、保全実務の23の組み合わせを持っている、使用されていました。 AVailable対立遺伝子セットは、コンフィギュレーションファイル "Alleles.csv"に格納されています。
- "実行"に進み、その後、SPEA2進化アルゴリズムを用いた多目的最適化を実行するために "SPEA2アーカイブのベースラインを意識サブセット"を選択します。
図1は、最適化の目的とパラメータの設定。
最適化パラメータを選択する：
プリセット：最適化する流域を選択します。 "Apply"をクリックすると、この画面上のコントロールの値を埋めるためにプリセットファイル "流域presets.csv"からエントリを選択します。
出力変数：最適化のための環境目標を選択します。 （Nアウトレット、アウトレットP）を選択したように、これは3次元の目的関数を定義しています：窒素（有機N + NO3 + NH4 + NO2）はリン（オルガ、出口での5年間の平均NICのP +ミネラルP）は出口で5年間の平均化され、保全実務の総費用。これは3次元のトレードオフのフロンティアを作成することに注意してください。多目的プログラム（{出力変数}、トータルコスト）を最小化するためにある別の出力変数は、選択することができます。
人口規模：設定初期人口規模。これは、解候補数の初期値を決定します。オプションの "各アレルによるシード"が選択された場合、流域内のすべての耕作地HRUsに設定アレルで指定された各保全実践の統一的適用を表す解候補が最初に作成されています。残りの候補解は耕作地HRUsに設定アレルから保全慣行の無作為割り当てによって作成されます。オプションの "各アレルによるシード"を選択した場合は、1は、初期集団のサイズが対立遺伝子セット内の対立遺伝子の数（このデモでは23）と少なくとも同じ大きさであることを確認する必要があります。
世代の > 数：最適化のラン（ランが再開できることに注意してください）の世代の希望の番号を（反復）を設定します。
交叉確率：2つの候補解が新たな解候補を作成するために選択された場合は、交叉確率は異なる新たな解が（このデモのために1に設定）作成されている確率を指定します。
一時的な人口の大きさ：これは、作成した新しい解候補の数を決定します。この値はプロセッサスレッドの数の整数倍（このデモで16）である場合にプロセッサのリソースを最も効率的に使用されています。
突然変異確率：対立遺伝子セットから別の省エネ実践にHRU割り当てのランダムな変化の確率を指定します。 （それはこのデモのために0.03に設定されています）。
数スレッド：使用プロセッサまたはスレッドの数を選択します。 16は、このデモで使用されています。
カーブ番号校正係数：これは、SWATのモデルのキャリブレーションから提供される。
テキストファイル内の人口を保存：これは、指定した数だけ繰り返された後に実行する最適化を再起動するように一つ願いが完了したかどうか選択することが重要です。このオプションをチェックすると、すべての存続候補解（個体）のすべてのHRUのアレルの値を指定したテキストファイルを生成します。これは実行を再開し、続行します後ろに読み取ることができます。
二次最適化パラメータ
最初の年は：過去の気象情報の開始年後に設定されており、遅くとも7年以上このデータの終わりの前にしなければなりません。
トウモロコシの価格：肥料削減のコストを見積もるために歩留まり損失の式と組み合わせて使用します。
パージ方法：支配した。すべて3次元で悪化している個人が支配し、パージされます。
HUCソース： "指定した場所"に設定し、次のフィールドから"7100006"の値を意味する"流域HUC"はHUCゾーンテーブルの行を検索するために使用されます。値は、 "07100006"はアライグマ流域のための8桁のHUCコードです。
コストソース：CRP以外の費用は上記の練習のコストテーブルに郡のFIPSコードによって決定されるであろうことを示すために、 "郡（HRUロケーション·コード）"に設定されています。
コストソースCRP：CRPのコストは上記の練習のコストテーブルに郡のFIPSコードによって決定されるであろうことを示すために、 "1場所"に設定してください。
3。代表的な結果 GeneticiSWAT.exeはすべての解候補（個人）だけでなく、最適化の実行を再開するために使用することができます最終的なアルゴリズムの反復との結果をコードする "保存"ファイルの設定と結果を示すログファイルを生成します。 この時点で、1次のステップに従うことで、パレート効率的なソリューション（トレードオフのフロンティア）のセット全体を視覚化することができます： 出力は、POV-Rayのプログラムを使用して、次に "レンダリング"、 "ファイルキュー"を選択して画像ファイルに一度にレンダリングされたすべてのことができるファイルのシリーズです。画像は自分で使用したり、アルゴリズムの進行を示すムービーにまとめることができます。 必要であれば、アルゴリズムの進行を示すムービー "はFramescanner.exe"を実行し、次の手順を実行することによって作成できます。 フロンティア内の各ポイントは、流域構成（景観上の保全慣行の特定の割り当て）を表します。これらの構成のマップは、次の手順を実行して、全体のフロンティアとして見ることができる： 関心のある特定の流域構成（個人）をエクスポートします。 多くの場合、目的の質問は特定の水質目標を達成する具体的な流域の構成（個人）を選択することです。例えば、人は個人ラージヒルを見つけたがっているかもしれませんベースラインの負荷に対して20％の30％、リンによって窒素を低減フロンティアでidual。 MapSWATは1つが指定された目標に対する最小ユークリッド距離を使用して個々のためにフロンティアを検索することができます。これは、次の操作を実行して行うことができます。 files/ftp_upload/4009/4009fig7.jpgの "alt ="図7 "/> さらに分析するために、地図データをエクスポートすると、以下の手順に従うことにより可能である： 表1に必要なプログラムの表。 表2に必要なサンプルファイルの表。 
図2：Scが3次元可視化のためのフロンティア "スナップショット"を作成するためのreenshot。 
図3はトレードオフのフロンティアの静的可視化。 
図4：最後のフロンティア内の個々のマップを作成するのスクリーンショット。 
図5：水質目標に基づくフロンティアで特定の個人を検索のスクリーンショット。 
図6：検索出力のスクリーンショット
図7。フロンティアで選択された個々のサンプルを記述したマップのスクリーンショット。 拡大図を表示するには、ここをクリックしてください 。 プログラムの名称 ソース 説明 回転子 カード土壌、気象、流域の管理データとi_SWATデータベースを作成し、塗りつぶします。 Swat2005GA.exe 米国農務省草地、土壌·水研究所 流域シミュレーションモデル i_SWAT.exe カード SWATのデータベースインターフェイス GeneticISWAT.exe カード 進化アルゴリズムのSWATのコントローラ。からGALib組み込んhttp://lancet.mit.edu/ga/ 。 MapSWAT.exe カード i_SWATデータベースとシェープファイルを読み込んで、世代や個人の画像を生成します。 POV-Rayは Povray.org ビジョンレイトレーサーの永続性。 Framescanner.exe トッドキャンベル AVIファイルコンバータへのPNG画像 Windows Liveムービーメーカー マイクロソフト AVを圧縮するために使用私はWMVへ ファイルの名前 タイプ 説明 アライグマGA.mdb Accessデータベース アライグマ流域の構造と管理の説明。 GeneticiSWATとMapSWATによる読み込み。 流域presets.csv テキスト GeneticiSWAT.exeとMapSWATの設定プリセット Alleles.csv テキスト 進化的アルゴリズムのための対立遺伝子セットの一覧。 アライグマアレルHRU.txt テキスト GeneticISWATが農地であると判断、それらの対立を記載することによって作成されたファイル。 MapSWATによる読み込み。 subbasin Josh.mdbの訓練費 Accessデータベース 管理の実践や郡によってコストを削減できます。 テラスZones.mdb Accessデータベース 表[HUCデータ]流域用のテラスと歩留まりゾーン番号を保持しています。 NRI Budgets.mdb Accessデータベース この実行で使用されていない作物＆マシンのテーブルのGeneticISWAT.exeによる読み込み。 phucrp 2008-12-15.dat テキスト この実行で使用されていない植物の熱ユニットルックアップテーブル、。 Management.mdb Accessデータベース この実行で使用されていない回転ルックアップテーブル、。 アライグマジョージア2011年9月28日1313.log、アライグマジョージア2011年9月29日0732.log、アライグマジョージア2011年10月7日0644.log テキスト GeneticISWAT実行のログファイル。 アライグマGA.wmv アニメーション 世代別の個人の3Dディスプレイ Subbasin.shp ESRIシェイプファイル 流域におけるsubbasinsの概要。 アライグマMap.wmv アニメーション フロンティア上の個々の各subbasinための支配的な対立の表示。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
私たちは、流域レベルの栄養素削減目標の範囲を達成するために農業保全慣行の最低コストミックスと場所を含む流域構成のパレート効率的なセットを検索するための統合されたシミュレーション·最適化フレームワークを構築します。シミュレーション·最適化システムの概念図を図8に示されている。農業保全慣行の水質への影響をシミュレートすることを含む、流域シミュレーションは、Windowsベースのデータベース管理システム、i_SWAT 6,8と相まって水文モデル、SWAT2005によって処理されます。最適化コンポーネントは、SWATの水文応答単位（HRUs）上で動作し、同時に栄養負荷（N、P、または両方）と保全慣行のコストを最小限に抑える保全慣行の配分を見つけるために進化アルゴリズム26のロジックを採用しています。アルゴリズムの反復が、個々を存続のセットを終えられた後にualsは、おおよそのトレードオフのフロンティアを表しています。 2つの栄養素を同時に（硝酸塩-N及び全リン）が検討されているので、我々は3次元のトレードオフのフロンティアを入手してください。トレードオフのフロンティア上の個々の点は、流域内の各意思決定ユニット（HRU農地）の保全慣行の特定の構成を規定しています。保全の実践が選択されているかを確認するために、我々は栄養素のターゲットを指定してから、栄養削減基準を満たす個々の構成にはトレードオフのフロンティアを検索する必要があります。選択された保全の実践の場所とミックスは、（このようなデータがHRUs作成時に利用可能であるなら）流域のフィールドレベルの空間的意思決定単位にマッピングすることができます。保全慣行の特定の混合と分配を指定する当社のアプローチは、水質の改善に向けた省エネルギー政策のより良いターゲティングのためのツールと政策立案者を提供することができます。実装の観点から、アルゴリズムの処方で武装して、政策立案者は（11で提案された方法）をターゲットの支払いを提供する、または入札を引き出すとガイダンスとしてモデリングの結果を使用してそれらを受け入れるか、または拒否することができます。もちろん、ターゲットとプラクティスの特定のセットは、特定の水質目標と学び特定流域に依存します。しかし、水文モデルの今後の改善と経済的コストの見積もりが容易にシミュレーション·最適化システムに組み込むことができます。開発したフレームワークは、容易に一般化可能であり、非点源汚染削減の複雑な問題に有用と政策関連の洞察を提供することが可能である。
図8実験の全体の流れ。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
特別な利害関係は宣言されません。
Acknowledgments
この研究は、米国環境保護庁のターゲット流域グラント·プログラム（プロジェクト＃WS97704801）、結合自然と人間系の国立科学財団のダイナミクス（プロジェクト＃DEB1010259-CARD-クリン）、及び米国司法省から受け取ったサポートから一部で賄われていたFoodand農業協定の農業プロジェクトの農研究所（プロジェクト＃20116800230190-CARD-）の。
References
- American Society of Agricultural and Biological Engineers. Design, layout, construction, and maintenance of terrace systems. ASAE Standard. , S268 (2003).
- Raccoon River Watershed water quality master plan. , Agren, Inc. Carroll, IA. (2011).
- Arabi, M., Govindaraju, R. S., Hantush, M. M. Cost-effective allocation of watershed management practices using a genetic algorithm. Water Resour. Res. 42, W10429 (2006).
- Bekele, E. G., Nicklow, J. W. Multiobjective management of ecosystem services by integrative watershed modeling and evolutionary algorithms. Water Resour. Res. 41, W10406 (2005).
- Carpentier, C. L., Bosch, D. J., Batie, S. S. Using spatial information to reduce costs of controlling agricultural nonpoint source pollution. Am. J. Agr. Econ. 61, 404-413 (1998).
- CARD: Resource and Environmental Policy Interactive Software Programs [Internet]. , Center for Agricultural and Rural Development, Iowa State University. Ames, IA. Available from: http://www.card.iastate.edu/environment/interactive_programs.aspx (2007).
- Gassman, P. W., Reyes, M., Green, C. H., Arnold, J. G. The Soil and Water Assessment Tool: Historical development, applications, and future directions. Trans. ASABE. 50 (4), 1211-1250 (2007).
- The i_SWAT software package: a tool for supporting SWAT watershed applications. Gassman, P. W., Campbell, T., Secchi, S., Jha, M., Arnold, J. G. SWAT2003: The 2nd International SWAT Conference, 1-4 July, Bari, Italy, , Instituto di Ricerca sulle Acque, IRSACNR. Bari, Italy. 66-69 Forthcoming.
- Gitau, M. W., Veith, T. L., Gburek, W. J., Jarrett, A. R. Watershed level best management practice selection and placement in the Town Brook watershed, New York. J. Am. Water Resour. Assoc. 42, 1565-1581 (2006).
- Jha, M., Rabotyagov, S., Gassman, P. W. Optimal placement of conservation practices using genetic algorithms with SWAT. Int. Agr. Eng. J. 18 (1-2), 41-50 (2009).
- Jha, M. K., Wolter, C. F., Schilling, K. E., Gassman, P. W. Assessment of total maximum daily load implementation strategies for nitrate impairment of the Raccoon River, Iowa. J. Envir. Qual. 39, 1317-1327 (2010).
- Kanna, M., Yang, W., Farnsworth, R., Onal, H. Cost effective targeting of CREP to improve water quality with endogenous sediment deposition coefficients. Am. J. Agr. Econ. 85, 538-553 (2003).
- Kramer, R. A., McSweeny, W. T., Kerns, W. R., Stravros, R. W. An evaluation of alternative policies for controlling agricultural nonpoint source pollution. Water Resour. Bull. 20, 841-8446 (1984).
- Lant, C. L., Kraft, S. E., Beaulieu, J., Bennett, D., Loftus, T., Nicklow, J. Using GIS-based ecological-economic modeling to evaluate policies affecting agricultural watersheds. Ecol. Econ. 55, 467-484 (2005).
- Muleta, M. K., Nicklow, J. W. Evolutionary algorithms for multiobjective evaluation of watershed management decisions. J. Hydroinform. 4 (2), 83-97 (2002).
- Muleta, M. K., Nicklow, J. W. Decision support for watershed management using evolutionary algorithms. J. Water Resour. Plan. Manag. 131 (1), 35-44 (2005).
- Rabotyagov, S. S., Campbell, T., Jha, M., Gassman, P. W., Arnold, J., Kurkalova, L., Secchi, S., Feng, H., Kling, C. L. Least Cost Control of Agricultural Nutrient Contributions to the Gulf of Mexico Hypoxic Zone. Ecol. Appl. 20, 1542-1555 (2010).
- Rabotyagov, S., Jha, M., Campbell, T. Impact of crop rotations on optimal selection of conservation practices for water quality protection. J. Soil. Water Conserv. 65 (6), 369-380 (2010).
- Rabotyagov, S. S., Jha, M., Campbell, T. D. Nonpoint-source pollution for an Iowa watershed: An application of evolutionary algorithms. Canadian Journal of Agricultural Economics/Revue canadienne d'agroeconomie. , (2010).
- Rabotyagov, S. S., Jha, M., Campbell, T. Searching for Efficiency: Least Cost Nonpoint Source Pollution Control with Multiple Pollutants, Practices, and Targets. J. Nat. and Environ. Sci. 1, 75-90 (2010).
- Ribaudo, M. O. Consideration of off-site impacts in targeting soil conservation programs. Land Econ. 62, 402-411 (1986).
- Ribaudo, M. O. Targeting the conservation reserve program to maximize water quality benefits. Land Econ. 65, 320-332 (1989).
- Srivastava, P., Hamlett, J. M., Robillard, P. D., Day, R. L. Watershed optimization of best management practices using AnnAGNPS and a genetic algorithm. Water Resour. Res. 38 (3), 1-14 (2002).
- Veith, T. L., Wolfe, M. L., Heatwole, C. D. Development of optimization procedure for cost-effective BMP placement. J. Am. Water Resour. Assoc. 39 (6), 1331-1343 (2003).
- GAlib: Matthew's Genetic Algorithms Library [Internet]. , Available from: http://lancet.mit.edu/ga/ (1996).
- Whittaker, G., Confesor, R., Griffith, S. M., Färe, R., Grosskopf, S., Steiner, J. J., Mueller-Warrant, G. W., Banowetz, G. M. A Hybrid Genetic Algorithm for Multiobjective Problems with Activity Analysis-based Local Search. Eur. J. Operat. Res. , (2007).
- Zitzler, E., Laumanns, M., Thiele, L. SPEA2: Improving the Strength Pareto Evolutionary Algorithm for Multiobjective Optimization. CIMNE. , 95-100 (2002).
