Dependence of coefficient of diffusivity of potential vorticity on wind stress and topographic roughness in an eddy saturation regime

技術要約：渦飽和領域における風応力および海底地形の粗さに対する潜在渦拡散係数の依存性

問題提起
南大洋などの海洋領域において、風応力の増加が必ずしも体積輸送の比例的な増加につながらない現象があり、「渦飽和（eddy saturation）」として知られている。高解像度の渦解像モデル（FR）は、この非感受性を再現できる一方で、低解像度モデル（CR）は、解像されていない渦の統計的効果を表現するためにパラメタリゼーション・スキームを必要とする。これらのパラメタリゼーションにおける重要な課題は、係数、具体的には潜在渦（PV）拡散係数（$k$）の適切な値を決定することである。FRモデルからPVフラックスを直接診断することは、PVフラックスに発散成分と回転成分の両方が含まれており、その分解が非一意的であるため困難である。したがって、粗解像度モデルの精度を向上させるために、外部強制（風応力）および海底地形に対して$k$がどのように依存するかを理解する必要がある。

手法
本研究では、逆解法を組み合わせたデュアルモデル・アプローチを採用している：

1. 高解像度（FR）モデル： 10 kmの水平格子を用いた準地衡流・バロトロピック・チャンネルモデルを使用し、個々のメソスケール渦を解像する。このモデルには、帯状の風応力強制、線形底面摩擦、および規定された海底地形（正弦波状の変動およびフーリエモードで定義）が含まれる。実験は、風応力の大きさ（0.1から0.3 N/m²）および地形の振幅（100から300 m）を変えて実施される。
2. 低解像度（CR）モデル： 拡散的PVフラックススキーム（$v'q' = -k\nabla q$）を用いて渦をパラメトリゼーションしたスペクトル解析モデルを使用する。このモデルは、同じ地球物理学的パラメータと強制力を共有しているが、明示的な渦の解像度は持たない。
3. 逆法： 著者らは、フラックスから$k$を直接診断しようとする代わりに、逆のアプローチを用いる。彼らは、FR実験から得られた体積輸送（平均帯状速度、$U$）を「既知」の効果であると仮定する。次に、与えられた風応力と地形の下で、この特定の輸送を生み出すための係数$k_0$（領域平均拡散係数）を求めるために、CRモデルを使用する。
4. 理論的導出： 著者らは、帯状運動量、渦エネルギー、および渦エントロピーの積分制約に基づき、$U$、$k_0$、および風応力（$\tau_0$）を結びつける解析的な関係式を導出する。この導出は、飽和領域において簡略化される三次方程式へとつながる。

主な貢献と結果

• 飽和領域における線形関係： 主要な理論的知見は、渦飽和領域において、PV拡散係数（$k_0$）と帯状風応力（$\tau_0$）の間に線形関係が存在することである。導出された方程式は $k_0 = \Theta \tau_0 \rho_0^{-1} + \mu$ であり、ここで $\Theta$ と $\mu$ は海底地形、幾何学的パラメータ、および基準状態に依存する定数である。
• 逆法の検証： FR輸送データを用いた逆法から得られた$k_0$の値と、理論的予測を比較することにより、著者らは一貫性を実証している。結果は、固定された地形において、風応力が飽和領域へと増加するにつれ、観測された輸送の飽和を維持するために必要な拡散係数が線形に増加することを示している。
• 地形の役割： 本研究は、海底地形が単なる運動量の受動的なシンクではなく、過渡的な渦を生成し、定常的な渦場を形成する能動的な役割を果たすことを確認している。導出された係数（$\Theta$ および $\mu$）は、海底地形の振幅およびスペクトル特性（フーリエモード）に明示的に依存する。
• 積分制約： 本研究は、物理的に妥当なパラメトリゼーションを策定する上で、積分制約（角運動量、エネルギー、およびエントロピーの保存）の重要性を再確認している。また、顕著な地形を持つチャンネルにおいては、一般化されたブレセトンの定理により、全PVフラックスはゼロではなく、地形形態応力によってバランスされることを強調しており、異なる層や領域間における拡散係数の特定の関係性を必要とすることを指摘している。

意義
本論文は、PV拡散係数を決定するための物理に基づいた枠組みを提供することで、その意義を主張している。渦飽和領域内における$k$の風応力への線形依存性を確立することにより、本研究は、高解像度実験から得られた既知の輸送値を用いて、粗解像度モデルにおける「自由な」パラメータである$k$を制約する方法を提示している。このアプローチは、一定または純粋にヒューリスティックな係数を超えたものであり、気候モデルや海洋モデルにおける渦飽和を正確に表現するには、拡散係数が強制力や地形の粗さに応じて動的に変化する必要があることを示唆している。これらの知見は、渦解像ダイナミクスとパラメトリゼーションされた低解像度モデルとの間の溝を埋めるための、逆法の使用を支持するものである。