土壌水分データの分析方法

How to analyze soil moisture data
見積依頼

土壌水分データの分析について 土壌水分データの分析について

貢献者

あなたは土壌水分量と水ポテンシャル(土壌吸引力)センサーを地面に埋め、ATMOS 41をフィールドに設置し、ZL6 データロガーをセットアップしました。機器のネットワークは数週間、あるいはシーズン中ずっとデータを収集しています。さてどうする?データから意味のある推論を推定する方法を知り、何が起こっているのかについて全体像の結論を出し、問題を特定してトラブルシューティングすることは、時に困難な場合があります。

この記事では、土壌水分量、土壌温度、土壌水ポテンシャル、および大気の測定値をどのように使用すれば、痕跡の背後にある意味を発見できるかを理解するために、複数のデータセットを順を追って説明します。この記事では、データ中の以下の事象を特定する方法を学びます:

  • 異なる土壌タイプにおける土壌水分センサーの挙動
  • 浸潤
  • 洪水
  • 土壌のひび割れ
  • 凍結
  • 空間的変動性
  • 温度効果
  • 水力再分配による日周パターン
  • 壊れたセンサー
  • インストールの問題

各例はグラフで表される。これらのグラフの情報のすべての側面を理解する必要はない。それぞれのグラフは、一般的な土壌水分データのパターンと、最も有用な情報を推定する方法を説明するためのものである。どのグラフにも、右上隅に土壌の種類と作物の種類を記したボックスがあり、このボックスによって変動要因をよりよく理解することができる。

提供されたデータはすべて、次のようなデータロガーによって収集された。 ZL6シリーズなどで収集し ZENTRA Cloud.すべてのデータセットはMETERの研究サイトにある計測器から得られたものか、データ所有者から提供されたもので、許可を得て掲載しています。

 

A photograph of a ZL6 next to a tablet showing ZENTRA Cloud data
図1ZL6 Basic データロガーとZENTRA Cloud プラットフォーム内に収集・保存されたデータ

 

図2は、芝草を被覆作物として用いた人工ローム砂のデータである。この例で実験を行った目的は、芝草の灌漑を改善することでした。この芝草の根域はかなり浅く、中央部の深さは約6cm、底部は約10cmであった。時間の経過とともに、この例では、まず6月と7月にかけて比較的湿った状態が始まり、7月と8月には一定の乾燥期間状態となり、8月と9月には水の取り込みが停止するまで乾燥が続いた。

 

土壌の種類による影響
A graph showing water content and water potential measurements for a turf grass in loamy sand in wet conditions
図2.湿潤状態のローム質砂地における芝草の含水率と水ポテンシャルの測定結果

 

図2は、2種類の土壌水分データも示している: 体積含水率左のY軸に体積含水率,右のY軸にマトリックポテンシャル,または 水ポテンシャルを右のY軸に示している。X軸は初夏から秋の始まりまでの時間である。これらのデータ群から何がわかるかを理解するには、それぞれのデータセットを個別に見なければならない。

A graph showing the turf grass in loamy sand shown in Figure 2 in wet conditions
図3.図2のローム質砂地における湿潤状態の芝草

 

図3は、前述の芝草の湿潤状態における含水量と水ポテンシャルの両方を示したものである。この芝草はローム質の砂の中にあった。グラフ上部の点線で表された水ポテンシャルセンサーはまったく反応しなかった。しかし、土壌含水量センサーは、日単位での灌漑イベントを含め、驚くほど詳細なデータを示している。

 

A photograph of the TEROS 10, 11, and 12 volumetric soil water content sensors
図4.TEROS 10、11、12体積土壌水分センサー

 

圃場は毎晩灌漑され、水がセンサーに当たったとき、6cmレベルのセンサーに見えるスパイクがあった。また、根域の最下部である15cmのところにも小さなスパイクがあった。30cmでも含水量の増加が見られたが、カーブは15cmのセンサーよりも丸みを帯びていた。水ポテンシャルはまったく変化を示さなかった。粒径が大きすぎて、センサーが粒子に保持されている水分を検出できなかったのだ。代わりに、このローム質砂の最適条件で何が起こったかを見てみると、かなり興味深い詳細が見えてくる。

A graph showing the turf grass in loamy sand shown in Figure 2 in optimal conditions
図5.図2に示したローム質砂地での芝草の最適条件

 

図5を見ると、6cmの土壌水分量データは、夜間は平坦になり、日中は低下している。これは連日見られ、植物が根域の最下部である6センチメートルでどれだけの水分を取っているかを知ることができた。15センチでも日ごとに低下しているが、根が水を吸い上げる最下部であるため、それほど顕著ではない。

 

A photograph of the TEROS 21 matric water potential sensor
図6.TEROS 21マトリック水ポテンシャルセンサー

 

また、図5では、プロファイルを伝って排水される水はそれほど多くなかった。7月14日に30cmレベルセンサーから少しピークが見られたが、その後の灌漑イベントでは変動はなかった。このローム質の砂は、散水された水に非常に反応する。水ポテンシャルのデータは、6cmレベルで小さな反応を示した。これは-200~-400kPaの範囲に低下しただけで、この芝草のストレス範囲をまだ超えているため、ストレスを示すものではなかった。

A graph of the turf grass in loamy sand shown in Figure 2 in drought conditions
図7.図2に示したローム質砂地における干ばつ条件下での芝草

 

干ばつ条件(図7)では、6cmレベルでも徐々に最適な水の吸収が見られた。このデータセットで問題となったのは15cmレベルのセンサーで、水位は6cmレベルのセンサーと同じように高く、水が吸収されずに土壌中に浸出していることを示していた。9月5日頃に突然止むまで、毎日の取り込みが見られた。この時点で、草は土壌から水を取り込むことができなくなり、活発な生育から休眠へと変化した。

この土壌水分データセットでは、水ポテンシャルが実に興味深い曲線を示しており、マイナスから-1500kPa、つまり永久萎凋点まで下がっている。この草は休眠状態にあり、草が取り込める水分がなかったからである。含水量と水ポテンシャルの測定値はいずれも、このような干ばつ条件下での減少を明確に示していた。残念なことにこのケースでは、土壌が非常に乾燥するまで、農家はデータの兆候に反応しなかった。

A graph showing an example of water content data for seed potatoes in a well-irrigated clay loam
図8.よく灌漑された粘土壌土における種イモの含水量データの例

 

土壌がローム状の砂ではなく、粘土質のロームというきめの細かい土壌になったらどうなるだろうか?図8は、アイダホ州南部で種イモを栽培していた直径700m近い粘土壌土に、6ヵ所のセンサーを設置したものである。このグラフは、シーズンを通して含水率がほとんど変化せず、2~3%程度しか変動していないことを表している。生産者はこのデータを見て、この畑の水を止めるべきタイミングをどうやって判断すればいいのだろうと考えた。前回のウェビナー「水資源の捕捉」で説明したように、このデータだけではその判断は非常に難しい:水をバイオマスに変える含水量データは水の存在と量を判断するのには非常に役立つが、植物にストレスがかかっているかどうかや、いつ十分な水が得られたかを理解するのにはほとんど役に立たない。

A graph showing the water potential measurements for the same field represented in Figure 8
図9.図8と同じフィールドの水ポテンシャル測定結果

 

植物のストレスレベルと水分更新能力を理解するためには、マトリックポテンシャル(水ポテンシャルとも呼ばれる)を調べる必要がある(図9)。土壌水分のデータからは、シーズンを通して一貫した、特筆すべきことのない水やりが行われていたことがうかがえるが、6カ所のうち3カ所では、マトリックポテンシャルがストレス域に低下しており、1カ所では永久枯死点近くまで低下していた。IRT赤外線温度計のような機器で測定したこれらの場所の植物の葉温は、気温よりもはるかに高い温度で記録された。そして、これらの場所での収量は、ストレスが測定されなかった場所での収量よりもはるかに低く、マトリックポテンシャルの測定とその問題の兆候の妥当性を示している。

 

浸潤
A graph showing volumetric water content for a native system in a desert sandy loam soil in the fall
図10.秋の砂漠砂質壌土における自生系の体積含水率

 

図10の砂質壌土は、ラッシュ・バレーの砂漠に見られる。この例では、作物を調査しているのではなく、在来系の中にある外来種を調査している。この計測器を設置した理由は、チートグラスのような外来種がこの地域の在来植生を駆逐している理由を理解するためであった。緑色の線は降水イベントを表し、水平線は異なる深さに設置された3つの水分センサーを表している。

 

A photograph of the ATMOS 41 all-in-one weather station
図11.ATMOS 41一体型ウェザーステーション

 

図10を見ると、4~5mmの湿潤イベントでも、5cmのレベルセンサーにはほとんど影響がなく、10cmと20cmのセンサーには目立った影響がないことがわかる。なぜ降水イベントは含水量データに現れなかったのでしょうか?いくつかの要因がある。この測定期間の前には、非常に長く、暑く、乾燥した夏があり、土壌の温度はほぼ毎日40℃を超え、土壌を疎水性にしていた。加えて、土壌が粉状に乾燥していたため、水分がすべて吸収され、表面で保持された後、再び蒸発し、水分がシステムの深部へ移動する機会がなかった。

A graph showing volumetric water content for a native system in a desert sandy loam soil in the summer
図12.夏の砂漠砂質壌土における自生系の体積含水率

 

図12は、同じ砂漠の砂質壌土の土壌水分データである。この図には降水量データは含まれていないが、深さ5cmの水分センサーに見られる各スパイクに相関して、湿潤イベントが発生している。5cm水位センサーでは、5月28日前後に湿潤現象が起きているが、10cm水位センサーにも20cm水位センサーにも反映されていないことに注意。6月2日前後の湿潤は、10cm深度で針を動かした証拠を示したが、20cmのマークには達していない。さらに驚くべきことに、最大の湿潤現象は6月14日頃に発生し、10cmの深さではまったく現れなかったが、20cmの深さでは小さなスパイクが生じた。これらの土壌水分データの分析から何がわかるだろうか?

図10と同様に、土壌の表層は1回目の湿潤から多くの水分を吸収し、土壌を下方に排出する機会を与えることなく、高い蒸発需要を供給した。2回目の湿潤時に水位が上昇し、一部の水が深さ10cmのセンサーに到達したが、20cmの地点には到達しなかった。より大きな難問は、最後の湿潤イベントであった。なぜ5cmと20cmのセンサーは水量の増加を記録したのに、10cmのセンサーは水量の増加を示さなかったのだろうか?

雨が降ると、土壌表面に水が一様に分布し、均等に浸透すると考えがちだが、必ずしもそうではない。水は土の中を1つの巨大な塊として進むのではなく、枝分かれした指のように土の中を進み、必ずしもすべての土の粒子に触れるとは限らない。今回のケースで最も可能性が高いのは、水の指の1本が10cmの深度センサーの周囲を回り、20cmのレベルセンサーまで移動し続けたことである。これがこのデータ異常の最も可能性の高い説明であるが、それでも、そのエリアの浸透に問題がないことを確認するために、そのエリアをモニターすることは重要であろう。

 

洪水の兆候

 

気象モニタリングが土壌の中で起こっていることを理解する上で違いを生み出すことができる使用例をいくつか探ってみよう。数年前、METER本社のあるワシントン州プルマンで洪水が発生しました。町の主要道路に平行して流れるミズーリ・フラット・クリークと呼ばれる非常に小さな川が、通りや複数の企業に浸水し、地域社会の中心部に大きな破壊をもたらしました。なぜこのようなことが起こったのか?このデータは、差し迫った洪水の警告サインを理解するのに役立つのだろうか?

(気象モニタリングの詳細については、気象モニタリング・マスタークラスを参照のこと)

A graph showing precipitation in Pullman, WA for the week proceeding the flood
図13.洪水が発生した週のワシントン州プルマンの降水量

 

図13の降水量データでは、なぜ洪水が起きたのか説明できない。最大の降水量は約3mmにしか達していない。複数のセンサーが同じ降水量を示したが、これはあの小さな小川で発生した規模の洪水の可能性を示唆していないようだ。何が起こったのかをよりよく理解するためには、この地域の土壌水分量の測定に目を向ける必要がある(図14)。

A graph showing water content measurements leading up to the flooding event in Pullman, WA
図14.ワシントン州プルマンで発生した洪水までの含水量測定結果。

 

土壌の水分レベルは、4月6日に湿潤現象が起こるまでは、非常に安定していた。4月7日の終わり頃、水は土壌の上層部から排出され、下層部へとろ過され始めた。その後の降雨により、9日には60cmの水位が満水となり、他の水位と同様に平坦になり始めた。4月9日末、さらに雨が降り、ミズーリ・フラット・クリークは堤防を越えた。

A graph showing water depths on the bank of the Missouri Flat Creek just before and at the beginning of the flooding event
図15.洪水発生直前と発生当初のミズーリ・フラット・クリークの堤防の水深。

図15は、クリークの側面に沿って水深を示したもので、水位が1mに保たれていることがわかる。最終的な降雨は、水位が堤防を越えて2.5m近くまで達し、水量のプラトーと相関した。

 

クラッキング

 

テキサス州南部に位置する高収縮膨潤粘土である船舶用粘土(図16)では、土壌水分データの分析が土壌亀裂の発生時にどのようになるかを説明するために、複数の含水率センサーを挿入した。

A graph showing water content sensors in a high shrink-swell clay in south Texas
図16.テキサス州南部の高収縮膨潤粘土中の含水量センサー

 

センサーの挙動を比較すると、興味深いパターンが見られた。深さ20cmのセンサーのひとつは、各湿潤イベントの後に緩やかな低下を示したが、同じ深さのもうひとつのセンサーは、一般的に砂質土壌で予想される急な低下を示した。このような測定値を出すために、粘土の中で何が起きていたのだろうか?このデータセットであなたが目撃している現象は、急カーブのセンサーが粘土の収縮率の高い部分に埋め込まれていたため、土壌がセンサーの周囲から収縮して空隙が生じ、電磁センサーが期待されたほど高い値を読み取れなかったということです。図16は、この急激な低下が土壌のひび割れを示していることを示す完璧な例である。

 

凍結
A graph showing water content and precipitation during freezing in sandy loam soil with invasive species growing on top during winter
図17.冬季に外来種が生育した砂質壌土の凍結時の含水量と降水量

 

図17は砂質壌土の圃場を示しており、気温のために枯死した外来種が含まれている。雨が降るたびに、含水量の測定値が跳ね上がり、その後階段状にギザギザに下がっているのがわかる。これらのセンサーは埋もれなくなったのだろうか?このデータの変動を説明できるものは何だろうか?温度測定値をグラフに加えると、何が起きていたかが明らかになる。

A graph showing water content, precipitation, and temperature during freezing in sandy loam soil with invasive species growing on top during winter
図18.冬季に外来種が生育した砂質壌土の凍結時の含水量、降水量、気温。

 

図18の太く黒い水平線は0℃、つまり凍結を示す。温度測定と凝固点の追加により、凍結すると含水率の測定値が下がり、変動することが明らかになった。気温が氷点下より上がると、含水率の測定値は再び上昇し、予想される範囲に戻った。この関係は理にかなっている。水が凍れば凍るほど、センサーが水の存在を検出するために使っている電気的磁場に対して水分子が消えていくのだ。すべての水が凍ったわけではないので、測定された水の量はゼロにはならなかったが、大幅に減少した。解凍が始まると、データは土壌の水分レベルを平滑化し、凍結前の状態に戻した。これは、凍結時に予想されるパターンである。図19は、同じデータセットを1年全体で見たものである。凍結現象が、夏のデータの滑らかな干満とは明らかに異なっていることに注目してほしい。冬季の点線は、凍結現象が発生しない場合の含水量データに期待されるパターンです。

A graph showing the same data set from Figure 18 expanded to include data from the entire year
図19.図18と同じデータセットを拡大し、通年のデータを含むようにしたもの。
空間的変動性

 

図20では、粘土壌土に種イモを植え、3~4km以内にある7つのエリアにセンサーを設置した。図20のセンサーは、雪の後に雨が降った雨の多い冬の後に設置された。

A graph showing water content readings for seven locations within a clay loam soil, all planted with seed potatoes
図20.粘土壌土内の7地点の含水量測定値(すべて種イモを植えたもの

 

含水量データを読む経験レベルにもよるが、データセットの最初の方でさえ、なぜこれら7つのセンサーの間にこれほどのばらつきがあるのか不思議に思うかもしれない。このばらつきは土壌の種類によるものである。土壌の種類が同じように分類されたとしても、各スポットには独自のベースラインがあるため、土壌水分量の測定値を同じスポットでの過去の測定値と比較することが重要であり、どんなに近い場所であっても、ある場所から次の場所まで同じ測定値を期待することはできない。しかし、同じ期間の同じ圃場の水ポテンシャルデータを調べると(図21)、これらのセンサーの水ポテンシャルはすべて±10kPa以内でスタートしており、驚くほど近いことがわかる。これが、土壌の種類に関係なく圃場全体の水ポテンシャル測定値を使用することの威力である。

A graph showing water potential measurements in the same seven locations shown in Figure 20
図21.図20と同じ7箇所の水ポテンシャル測定結果
温度感受性

 

夏場は、特に土壌表面に近い場所での測定において、ある程度の水分量の変動を考慮することが重要である。このように、すべての含水量センサーは温度センサーとしても機能する。

A graph showing water content and temperature measurements across the summer months showing the effect of temperature fluctuations on water content readings
図22.夏期における水分量と気温の測定値から、気温の変動が水分量の測定値に及ぼす影響を示す。

 

図22では、1日の温度変化が±14℃であった。このような温度変化に関する知識がなければ、水分量測定値のギザギザは水の水文学的な動きと誤って解釈されるかもしれないが、実際には土壌水分量に対する熱変化の影響を示している。各水分含量の小さな変化は、実際には0.0003m3 m-3/℃に過ぎない。

水力再分配

 

ここまで、当初は植生による水の取り込みの兆候に見えたかもしれないが、データをさらに調べると説明がつかなくなるような例をいくつか紹介してきた。では、あなたのデータでは、水力再分配はどのように見えるだろうか?土壌から植物への水の取り込みの存在を証明できるかどうか、4つのグラフを連続して調べてみよう。各グラフは、互いに500m以内の同じ灌漑小麦畑の同じ期間のデータを強調している。

A graph showing water content measurements at a depth of 15 cm across six sites on the same field combined with watering and precipitation data
図23.同一圃場内の6地点における深さ15cmの含水量測定値と散水および降水量データの組み合わせ。

 

図23は、圃場全体の深さ15cmの含水量を示している。灌漑は7月末までオン・オフを繰り返し、灌漑停止後に降水が数回あった。各水分センサーは典型的な日周パターンを示しており、一見すると気温の変動に見えるかもしれないが、この時期の圃場は小麦の葉キャノピーが充実しており、Leaf Area Index (LAI)は4~5程度で、土壌表面まではほとんど放射されていなかった。このため、日周変動が温度変化に起因するとは考えにくい。このパターンが止まったのは、灌漑が停止され、植物が可能な限り吸収した時だけであった。

A graph showing water content measurements for the same field shown in Figure 23, but at a depth of 45 cm
図24.図23と同じ圃場での含水量測定(ただし深さ45cm)。

 

図24の水深45cmでは、6月上旬には日周パターンは見られず、6月下旬になって顕著になった。7月下旬に水を止めた頃、日周パターンが最も顕著になり、日中は下がり、夜間はプラトーになった。土壌の奥深くにあるため、このデータは、温度変化とは対照的に、植物による水の吸収によく一致していた。

A graph showing water content measurements for the same field shown in Figure 23, but at a depth of 65 cm
図25.図23と同じ圃場での含水量測定(ただし深さ65cm)。

 

土壌の奥深くに進むと(図25)、7月中旬から下旬までは同じような日周の階段状の変化は見られず、8月に入っても見られた。

センサー故障

 

精度と信頼性は、私たちが製造するすべての計測器に関して最も重要な要素の2つですが、設置やメンテナンスなど、センサーの故障を引き起こす可能性のある他の要素もあります。

A graph showing an example of data that indicates the potential of a failed sensor
図26.センサーの故障の可能性を示すデータの例

 

図 26 は、故障したセンサーがデータの中でどのように見えるかを示す好例である。いくつかのセンサーは非常に安定したコースを描き、スムーズにデータを流していたが、突然一 つのセンサーだけがほぼ瞬時に落下し、大量の変動データを提供し始めた。このユーザーにとって幸運だったのは、センサーがZENTRA Cloud に接続されていたことで、センサー故障を自動的に警告し、連続測定の中断を最小限に抑えてタイムリーに対処することができた。この場合、センサーはプラグから外れていたが、データロガーに再び差し込むと、センサーは完璧に作動し続けた。

インストールの問題

 

この記事で取り上げる最後のデータ異常は、インストール・プロセス中の問題に起因する欠陥データである。

A graph showing the water content measurements of a silt loam soil without any cover crop
図27.被覆作物のないシルトローム土壌の含水量測定結果

 

すべての科学的な取り組みと同様、測定値を分析する前に、どのような 情報が得られるかを知っておくことが重要である。図27に示すようなシルトローム土壌の場合、土壌はかなり湿っており、体積含水率30%以上の値を示すと予想される。その代わりに、2つのセンサーが10%以下の値を示していた。これは、センサーを検査し、設置し直す価値があると思われる好例である。

 

土壌水分データ解釈の鍵

 

それぞれの状況やデータセットは異なる。データから正しい推測と推論を導き出すことは、結論の妥当性を高めるために極めて重要です。このことを念頭に置いて、あなたのユニークなデータを解釈する際に留意すべきことをいくつかまとめてみましょう。

  • 予期せぬものを見つけることを期待するセンサーを設置するのは、土壌で起こっていることをすべて把握しているわけではないからだ。知っているつもりでいると、重要な兆候が無視されることになる。
  • 土壌水分の挙動は土壌によって異なる。粗い土壌の土壌水分のデータが、細かい土壌のデー タと同じ挙動を示すことはない。時間をかけて土壌のタイプを理解し、専門家と協力して正しい測定値を得るようにしましょう。
  • 土壌が水を取り込み、貯蔵する能力は、水の浸透を止めることがある。データでテーブルトッピングに気づいたら、浸透が止まっているか、土壌が水浸しになっているかに注意してください。
  • ひび割れしやすい土壌に注意してください。一部の土壌、特に粘土は、センサーの周囲に亀裂を生じさせ、結果に影響を与えることがあります。どのような土壌で作業しているかを把握し、データからこのような挙動が見られないか注意してください。
  • 凍結によって土壌の水分がセンサーから消えてしまうことがある。凍った水はセンサーから見えなくなり、奇妙な挙動を引き起こし、連続性を犠牲にして精度を維持するために、しばしばデータの一部をデータセットから削除する必要がある。
  • 土壌水分量は 空間変動が大きい。同じ圃場全体に設置された複数のセンサーがばらつきを示すはずである。ばらつきがまったくない場合は、心配すべき兆候である。一方、土壌水分ポテンシャルセンサーに一貫性があれば、何か正しいことをしている証拠である。
  • 日周パターンは複数の意味を持ちうる。水分含量の測定値の日ごとの干満は、土壌表層付近の温度 変動や、根の水分や再分布によって引き起こされる可能性があ る。日周パターンの原因を特定するためには、水分量だけでなく、すべての変数を調べることが重要である。
  • 設置不良はセンサーの性能に大きく影響します。電気的に故障したセンサーだけでなく、設置が不十分なセンサーからも問題が生じる。研究への影響を最小限に抑えるためには、ZENTRA Cloud から提供されるアラートに注意し、迅速に対処することが重要です。
質問は?

当社の科学者は、研究者や生産者が土壌-植物-大気の連続体を測定するのを何十年も支援してきた経験がある。

測定に関する洞察

すべての記事を見る

土壌水分センサー-その仕組み。研究グレードでないものがある理由

TDR、FDR、静電容量、抵抗:一般的な土壌水分センシング方法の比較、その長所と短所、独自のアプリケーション。

測定インサイトを読む

土壌水分センサー:どの土壌センサーが最適か?

METER土壌センサーを使用した何千もの査読付き出版物の中で、どのタイプが最も好ましいというものはありません。したがって、センサーの選択はニーズとアプリケーションに基づいて行う必要があります。あなたの研究に最適なセンサを特定するために、これらの考察をお役立てください。

測定インサイトを読む

土壌水分とは何か?測定の背後にある科学

たいていの人は、土壌の水分を1つの変数(含水量)だ けから見ている。しかし、土壌中の水の状態を表すには、2種類の変数が必要である。

測定インサイトを読む

A photo of a METER publication in book form open on a flat surface

ケーススタディ、ウェビナー、記事

最新のコンテンツを定期的にお届けします。

アイコン角 アイコンバー アイコンタイムズ