Non-invasive imaging of solute redistribution below evaporating surfaces using 23Na-MRI

本研究は、$^{23}\text{Na-MRI}$を用いることで、高透水性多孔質媒体における密度駆動型のフィンガリング不安定性が蒸発溶質の顕著な下方再分配を引き起こし表面の飽和を防ぐ一方で、低透水性媒体では地表付近での極端な濃度蓄積を許容することを実験的に実証するものである。

要約

研究論文の解説：塩、砂、そして「豪雨」効果

想像してみてください。あなたには2つのカップの砂があります。一つは粗くてゴツゴツした砂利（これを「速い砂」と呼びましょう）が入っており、もう一つは非常に細かく粉末状の砂（「遅い砂」）が入っています。両方のカップに塩水を注ぎ、完全に浸るまで満たします。それから、水を蒸発させるために、温かく乾燥した風の中に置いておきます。

研究者たちが問いかけたのは、水が消えていくとき、塩はどこへ行くのか？ ということです。

常識では、水たまりに膜ができるように、塩は一番上の表面に積み上がるだけだと考えがちです。しかし、この研究では特別な「魔法のカメラ」（MRI）を使用して、砂の内部で実際に何が起きているかを観察しました。すると、驚くべき事実が判明しました。それは、砂の穴の大きさによって全く異なるのです。

ツール：「魔法のカメラ」

砂を掘り起こすことなく湿った砂の内部を見るために、科学者たちは病院のMRIスキャナーに似た装置を使用しました。ただし、膝の写真を撮る代わりに、彼らはナトリウム（食塩の「Na」）が見えるように調整しました。

これは、湿った砂の中で塩が光っている様子を見ることができるカメラのようなものです。彼らは3Dスナップショットを時系列で撮影することで、まるで部屋の中を移動する群衆のタイムラプス動画を見るかのように、塩がどのように動くかを観察することができました。

実験：2つの異なる物語

研究者たちは、透水性（水がどれくらい通りやすいか）が大きく異なる2種類の砂を用いて実験を行いました。

1. 「遅い砂」（細かい粒子）

• 何が起きたか： 上部から水が蒸発するにつれ、塩は上へ行く以外に道がなくなりました。塩は表面近くの、非常に狭くタイトな隙間に挟まり込んでしまいました。
• 結果： 塩は表面のすぐ下に、濃度の高い厚い層として積み上がり、結晶化（固体の岩になること）に近い状態にまで達しました。
• 例え： 混雑した廊下を想像してください。人々（水分子）が外に出ようとしていますが、ドアがとても小さいのです。人々は出口のすぐそばで立ち往生し、「塩」（重いバックパック）がその場所に積み上がります。そのせいで、他の誰も外に出ることが非常に困難になります。塩が上部を詰まらせたため、蒸発は著しく遅くなりました。

2. 「速い砂」（粗い粒子）

• 何が起きたか： 最初は、遅い砂と同じように、塩は上部に積み上がろうとしました。しかし、この砂の穴は大きかったため、異なる現象が起きました。上部の塩水が非常に重くなった（高密度になった）のです。
• 結果： 重力が支配権を握りました。重い塩水は上部に留まることができず、プールの底に落ちる重い石のように、砂の中へと沈んでいきました。これにより、塩を柱の深くまで運ぶ「プルーム（柱）」や「指」のような動きが生まれました。
• 例え： 開放的なスタジアムにいる群衆を想像してください。人々が去っていく際、重いバックパックを背負ったグループがあまりに重いため、列の先頭に留まることができなくなります。代わりに、彼らは群衆の脇をすり抜け、部屋の奥へと沈んでいきます。上部は比較的クリアな状態に保たれ、「塩」は砂の深い部分へと再分配されます。

なぜこれが重要なのか（論文による説明）

この研究は、科学者たちがこれまで推測するしかなかった理論を裏付けるものです。塩は必ずしも表面に留まるわけではないということです。

• タイトな砂の場合： 塩が表面に留まり、極度に濃縮され、蒸発を遅らせます。これは建築材料（風化の原因）や土壌（塩類集積の原因）にとって好ましくない状況です。
• 緩い砂の場合： 塩が沈み込みます。これは、表面の塩分濃度が低く保たれることを意味し、蒸発は一定のペースで進み続け、塩は表面に膜を作るのではなく、地面の深い場所に移動します。

「指」か「地殻」か

研究者たちは、実際の観察結果とコンピュータ・シミュレーションを比較しました。コンピュータ・モデルは、MRIで見たものと全く同じ結果を予測しました。

• 遅い砂は、上部に塩の「地殻（クラスト）」を形成しました。
• 速い砂は、下方へと滴り落ちる塩の「指（フィンガー）」を形成しました。

また、彼らは「レイリー数」（流体が十分に重くなって沈むかどうかを測るための、高度な概念）を用いて数学的な検証も行いました。数学的予測では、速い砂は不安定になって沈み、遅い砂はそのまま留まることが示されました。MRIカメラは、その数学が正しいことを証明したのです。

まとめ

この論文は、科学者が特別なカメラを使って、塩と砂に関する謎を解いた探偵物語のようなものです。彼らは、砂の粒の大きさが「交通整理」の役割を果たしていることを発見しました。

• 小さな粒は、塩を入り口で捕まえ、激しい交通渋滞（地殻）を作り出します。
• 大きな粒は、重い塩を地下室へと沈ませ、より多くの水が蒸発するための道を切り開きます。

これは、湖が干上がるとき、建物が損傷するとき、あるいは庭の土壌に影響を与えるときなど、自然界で塩がどのように移動するかを理解する助けとなります。

技術概要

技術要約：$^{23}\text{Na-MRI}$を用いた蒸発面下における溶質再分配の非侵襲的イメージング

問題提起
多孔質媒体からの塩水蒸発は、土壌の塩類集積や建築材料の風化に関連する極めて重要なプロセスである。先行研究では、蒸発の段階（SS1–SS3）が特定されているが、表面付近における溶質蓄積のダイナミクスは依然として複雑である。近年の理論モデル（例：Bringedal et al., 2022）は、完全飽和状態の高透水性多孔質媒体において、蒸発が密度駆動型の不安定性を誘発する可能性を示唆している。塩分が表面に蓄積するにつれ、結果として生じる流体密度の増加が下方へのフィンガリング流（指状流）を引き起こし、蒸発フロントから溶質を遠ざけるように再分配する可能性がある。しかし、この現象は主に理論的、あるいは二次元のヒーレ・ショーセル（Hele-Shaw cell）において観察されてきた。自然界の多孔質媒体におけるウィッキング条件下において、密度駆動の下方流が発生するかどうか、特に固有浸透率が表面濃縮から下方への再分配への転移にどのように影響するかについての三次元的な実験的証拠は不足している。

手法
この空白を埋めるため、著者らは対照的な固有浸透率を持つ2つの多孔質媒体サンプル、F36（中程度の砂、$k \approx 7.83 \times 10^{-12} \, \text{m}^2$）およびW3（非常に細かい砂、$k \approx 3.2 \times 10^{-14} \, \text{m}^2$）を用いて蒸発実験を行った。両サンプルは当初、1.0 mol L⁻¹のNaCl溶液で飽和させ、蒸発はタービュレントな$\text{N}_2$ガス流によって駆動される一方で、底部のリザーバーから同一の溶液を供給して飽和状態を維持した（ウィッキング条件）。

本研究の核心的な革新は、NaClの時空間的な再分配をモニタリングするための非侵襲的な$^{23}\text{Na-MRI}$（ナトリウム23磁気共鳴イメージング）の適用である。

• イメージング・プロトコル： 研究では4.7 TのBrukerスキャナーを使用した。より細かいW3の細孔で見られる速い$T_2$緩和に対応するため、異なるパルスシーケンス（F36には3D RARE、W3にはMSSE）を採用した。
• **校正：**既知のNaCl濃度（0から5 mol L⁻¹）を持つサンプルから得られた校正曲線を用いて、定量的なNa濃度マップを作成した。F36については線形関係が確立されたが、表面緩和効果のため、W3については二次多項式が必要となった。
• 検証： サンプルが実験を通じて完全に飽和状態に保たれていることを確認するために、$^1\text{H-MRI}$を並行して使用した。MRIから得られたNa質量と蒸発による質量損失データを比較することで、質量収支を算出した。
• 数値モデリング： 実験結果は、密度駆動流を組み込み、不安定性を誘発するための小さなランダム摂動で初期化されたDuMuxシミュレータを用いた2次元数値シミュレーションと比較された。

主な結果

1. 分岐する溶質挙動： 実験により、浸透率に基づいた根本的に異なる溶質再分配パターンが明らかになった。

   • W3（低浸透率）： 表面付近での緩やかで浅いNaの濃縮を示した。上部1 mmの濃度は約6.7 mol L⁻¹に達した（溶解度限界である6.14 mol L⁻¹に接近または超過）。下方へのフィンガリングは観察されなかった。その結果、表面濃度の上昇に伴い、蒸気圧が低下したため、蒸発速度は〜9 mm d⁻¹から〜4 mm d⁻¹へと急激に低下した。
   • F36（高浸透率）： 初期段階では表面濃縮が見られたが、最初の1時間以内に密度駆動の下方プルーム（フィンガリング）が発達した。このプルームはNaClをカラムの深部へと再分配し、極端な表面蓄積を防いだ。表面濃度は中程度に維持され（ピーク時で約2.5 mol L⁻¹）、実験全体を通して蒸発速度は比較的一定（約6 mm d⁻¹）であった。

2. 安定性解析： 観察された挙動は、Woodingら(1997a)によって定義された安定性基準と一致している。F36のレイリー数（$R_\delta$）は33.5と計算された（不安定性を示す）のに対し、W3は0.59（安定）であった。Bringedalら(2022)による線形安定性解析は、F36における不安定性の発生時間を0.61時間と予測しており、これは実験でフィンガリングが出現した約1時間という観察結果と一致した。

3. モデルとの比較： DuMuxを用いた数値シミュレーションは、実験データと定性的な一致を示した。モデルは、W3における表面濃縮と、F36における下方フィンガーの形成をうまく再現した。しかし、後半段階の濃度プロファイルには不一致が見られ、これは境界効果（有限のサンプル高さ）およびモデルにおける簡略化された境界条件に起因すると考えられる。

意義と主張
本研究は、蒸発中の飽和多孔質媒体における密度駆動型の溶質再分配に関する、三次元スケールでの初の実験的確認を提供するものである。著者らは、固有浸透率が、蒸発が表面の地殻形成につながるか、あるいは深い溶質の再分配につながるかを決定する決定的な要因であるという理論的予測を、本研究が検証したと主張している。

具体的には、論文は以下のことを断言している：

• 高浸透率の媒体は、上方への蒸発フラックスに対抗する密度駆動の下方流を許容し、それによって表面の飽和と地殻形成を遅らせる。
• 低浸透率の媒体は溶質を表面付近に閉じ込め、それによって急速な濃度上昇、蒸発速度の低下、および早期の地殻形成をもたらす。
• $^{23}\text{Na-MRI}$は、信号対雑音比（S/N比）の課題はあるものの、密度駆動型の不安定性を解像することが可能な、多孔質媒体における溶質再分配を定量化するための実行可能な非侵襲的ツールである。

著者らは、これらの知見が、蒸発速度、塩類地殻形成のタイミング、およびプラヤやソルトフラットに見られるような塩水地下水系の水文学に対する示唆を持つと結論付けている。彼らは、今後の研究として異なる種類の塩や不均質な媒体を調査すべきであることを控えめに示唆しているが、主要な貢献は、これまで理論のみであった密度駆動流のメカニズムを実験的に検証した点にある。