DMSTAr: A reproducible R framework for network-based simulation of hydrology and phosphorus dynamics in treatment wetlands

技術要約：DMSTAr – 水文学およびリン動態のための再現可能なRフレームワークとしてのDMSTAr

問題提起
湿地の水質モデルは、管理された景観、特にフロリダ・エバーグレーズズにおけるリン（P）除去と復元計画を評価するために不可欠である。広く使用されているレガシーモデルである「Dynamic Model for Stormwater Treatment Areas (DMSTA)」は、Microsoft ExcelワークブックとVisual Basic for Applications (VBA) コードとして実装されている。機能的ではあるものの、このスプレッドシートベースのアーキテクチャは、透明性、再現性、および現代的な科学的ワークフローへの統合を制限している。ロジックが不透明な数式や分散したセル依存関係の中に隠れていることが多く、体系的な検証、バージョン管理、およびスケーラブルなシナリオ分析を困難にしている。本研究は、確立された数値的挙動や概念的定式化を変更することなく、このレガシーツールを近代化する必要性に対処するものである。

手法
著者らは、レガシーなDMSTAモデルの数値的挙動を再現しつつ、モジュール化されたテスト可能なフレームワークを提供する、オープンソースのRパッケージであるDMSTArを開発した。開発プロセスでは、Excel/VBAのロジックを単なる行単位のコード変換ではなく、明示的なR関数へと系統的に翻訳する手法が用いられた。

• アーキテクチャ: パッケージは、元のDMSTAの概念構造を反映するように階層的に構成されている：
  1. コア関数: 単一の処理セル内における水文学および全リン（TP）動態をシミュレートし、日次の水収支、P質量収支、およびサブデイ（1日未満）の数値積分をカプセル化する。
  2. ケースレベルのシミュレーション: 共通の強制力（forcing）とパラメータを持つ複数のセルを調整する。
  3. ネットワーク構成: 複数のケース間における上流・下流の接続性、流量ルーティング、および質量移動をサポートする。
• 忠実度の維持: DMSTArは、元の支配方程式、数値積分スキーム（Euler法およびRunge-Kutta 4法をサポート）、およびバージョン固有の挙動（例：DMSTA バージョン2E対2C2Bのリリース・ゲーティングおよびオフライン休止期間）を明示的に保持している。
• 検証戦略: 実装の検証は、4つの主要な計画アプリケーション（Restoration Strategies [RS]、Central Everglades Planning Project [CEPP]、およびWestern Everglades Planning Project [WERP]）から得られた、同一の入力、パラメータ、初期条件、およびネットワーク構造を用いて、レガシーなDMSTAの出力と比較することで行われた。
• 指標: 一致度は、符号付きバイアス、中央絶対差、平方根平均二乗誤差（RMSE）、対称パーセント差、およびBland–Altmanプロットを用いたグラフィカルな診断を用いて評価された。また、10年間にわたるステップサイズ感度分析も実施され、両方の実装間の数値的安定性と収束性を比較した。

主な貢献

• オープンソースの実装: DMSTArは、ソースコード、ドキュメント、およびビネット（vignettes）への完全なアクセスを提供するオープンソースのRパッケージとして、CRANおよびGitHubで公開されている。
• モジュール化されたワークフロー: 本パッケージは、モデルの構成、検証、およびシミュレーション実行を、再利用可能な個別の関数（例：dmsta_flowP_case(), run_network_of_cases()）に分離しており、ユーザーがデータ入力からシナリオ評価へとシームレスに移行することを可能にする。
• 透明性と拡張性: 暗黙的なスプレッドシートの依存関係を明示的な状態空間更新へと再定式化することで、診断やテストのために中間状態やフラックス（通量）を直接取得できるフレームワークとなっている。
• レガシーの継続性: 本フレームワークは、元のモデルとの厳格な数値的等価性を維持しており、過去の計画結果との比較可能性を確保しながら、将来的な拡張を可能にしている。

結果

• 数値的等価性: DMSTArは、テストされたすべてのアプリケーションにおいて、無視できる程度の系統的バイアスと高い一致度をもって、レガシーなDMSTAの挙動を再現した。水文学的状態（深度、容積、流量）およびP指標（濃度、質量、負荷量）の時系列軌跡は、ほぼ同一であった。
• 誤差分析: Bland–Altman分析によれば、差はゼロ付近に密に集まっており、大きさに応じた発散は見られなかった。ほとんどのケース固有のバイアスおよび誤差推定値は1%以内であった。
• 感度と安定性: ステップサイズ感度分析により、DMSTArとDMSTAは、積分解像度の変化（1日から16ステップ/日）に対してほぼ同一の応答を示すことが明らかになった。両方の実装は、近似的に負のサブリニアな収束（傾き約-0.5）を示したが、これはアルゴリズムの数値的不安定性ではなく、モデルの閾値ベースの運用ロジックと不連続な強制力に起因すると著者らは考えている。
• 差異の要因: 両実装間の残留差は、積分ステップに関わらず安定していることが判明しており、これは根本的なアルゴリズムの変化ではなく、スプレッドシート/VBAモデルをRへ翻訳する際に生じる、微細な実装レベルの要因（例：丸め処理、条件評価の順序）に起因することを示唆している。

意義および主張
本論文は、DMSTArが、監査が困難な広く利用されている処理湿地モデルを、透明性が高く、モジュール化された、再現可能な生態学的モデリングフレームワークへと正常に変換したことを主張している。主な意義は、レガシーモデルの挙動を実質的に変えることなく、透明性、拡張性、およびアクセシビリティを向上させた点にある。

著者らは、本研究の範囲は元の実装とRへの翻訳との間の同等性を実証することに限定されており、基礎となるDMSTA手法の厳密さを評価したり、モデルの再較正を行ったりするものではないことを強調している。代わりに、DMSTArは、シナリオ分析、診断評価、および長期的な科学的应用に適した、忠実な実行エンジンとして機能する。本フレームワークは、レガシーな生態学的意思決定支援モデルを、過去の適用例との連続性を維持しながら、いかにしてオープンでテスト可能な研究ソフトウェアへと翻訳できるかを示すものであり、より強固で再現可能な湿地管理計画を支援するものである。