海洋の混合層深度（MLD）とは？求め方と季節変化の読み方

水温や塩分のプロファイルを読むと、表層付近には値がほとんど変わらない浅い層があり、その下で値が急に変わり始めることに気づきます。この表層のほぼ一様な層が混合層で、その深さを表す混合層深度（MLD）は、海洋物理学で最もよく使われる指標の一つです。

MLDを使うと、風による混合や表層の冷却が上層をどれだけ深くかき混ぜたかを、一つの数値で表せます。また、亜表層酸素極大やモード水の層とも密接に関係しています。

この記事では、混合層深度が物理的に何を意味するのか、Argoプロファイルからどのように計算されるのか、同じ水柱でも方法によってMLDの値が変わりうる理由、そしてOceanGraphの実際のArgoデータを使って季節変化を読む方法を説明します。

まだプロファイルを読んだことがない場合は、先に 海洋の水温・塩分プロファイルを読み解く を参照してください。MLDは、そこで扱う表層の読み方の延長線上にあります。

混合層深度が重要な理由

混合層は、海が大気と熱・気体・運動量をもっとも直接的にやり取りする部分です。その深さは、次のようなことに影響します。

• 同じだけ表層が暖まったり冷えたりしたときに、どれだけの熱量が動くか
• 表層付近の栄養塩や酸素などが、どれだけ深くまでかき混ぜられるか
• 季節躍層が海面からどのくらいの深さに位置するか
• 二つのプロファイルがどれだけ比較可能か（混合層が20 mの場合と150 mの場合とでは、「表層」の値が代表する水の量がまったく異なる）

このため、MLDは海面と大気の熱収支、生物地球化学的な循環、モード水の形成、季節予報で使われる混合層の熱量など、海洋学のさまざまな場面に登場します。

混合層深度とは何か

物理的には、混合層とは、主に風によるかき混ぜと、表層の冷却や蒸発によって生じる対流によって均質化された表層付近の層です。この層では水温・塩分・密度がほぼ一様である一方、その下では勾配が強まり、水柱はより成層した状態になります。MLDは、このほぼ一様な層が終わる深さを表します。

ArgoプロファイルからMLDを求める方法

OceanGraphは、個々のArgoプロファイルから、探索用の指標として、関連する三つの量にしきい値基準を適用してMLDを導出します。

使用する三つの量

• ポテンシャル水温（θ）：圧縮による水温の上昇分を取り除き、鉛直方向の変化が圧力の違いではなく水そのものの性質の違いを反映するようにしたもの
• 絶対塩分（SA）：実用塩分、圧力、位置からTEOS-10規格に基づいて求めたもの
• 0 dbar基準のポテンシャル密度偏差（σ₀）：0 dbarを基準に求めたポテンシャル密度から1000 kg/m³を引いた量で、水温と塩分の効果を一つの密度指標として表したもの

水温と塩分は密度への効果を互いに打ち消し合うことがあるため、片方だけを見る方法では、密度としては混合しているのに水温では混合していない（あるいはその逆の）ケースを見逃すことがあります。そのため、三つをあわせて使うことで、「混合している」という物理的な意味により近い定義になります。

三つのしきい値基準

MLDは、浅い基準深度（10 dbar）を起点とした、三つの独立したしきい値判定によって推定されます。

• 水温しきい値（Δθ）：ポテンシャル水温が基準深度での値から初めて0.5°C以上変わる深さ
• 塩分しきい値（ΔSA）：絶対塩分が基準深度での値から初めて0.05 g/kg以上変わる深さ
• 密度しきい値（Δσ₀）：ポテンシャル密度偏差が基準深度での値から初めて0.125 kg/m³以上変わる深さ

各基準について、OceanGraphはしきい値を挟む最初の二つの圧力点を探し、その間を線形補間して、しきい値と交差する圧力を求めます。したがって、二点のうち深い側の圧力をそのままMLDにするわけではありません。こうして三つの基準を独立に計算し、補間した交差圧力のうち、もっとも浅い値をMLDとして採用します。ただし、プロファイル中のどこでも三つのしきい値に達しない場合は、強引に値を決めず、未定義として扱います。

OceanGraphの値を読むのではなく、自分で計算するために基準を選ぶ場合は、混合層深度（MLD）の求め方：基準の選び方と結果への影響で、こうした数値の背後にある4つの判断を一つずつ扱っています。

データ品質の要件

しきい値による方法が意味を持つのは、プロファイルの表層付近に十分な数の測定点がある場合に限られます。MLDを計算する前に、次の二つを確認します。

• 基準深度の範囲：ポテンシャル水温（θ）と絶対塩分（SA）の両方に欠損していない数値がある圧力行だけを使います。最浅行の圧力は30 dbar以下、最深行の圧力は10 dbar以上でなければなりません。最浅行が10〜30 dbarにある場合は、その値を基準値として使います。
• 対象行の数：0〜2,000 dbarに3点以上、そのうち50 dbar以浅に2点以上必要です。

圧力から深さへの変換

Argoが直接測る鉛直座標は深さではなく圧力なので、OceanGraphは推定したMLDをプロファイルと同じ圧力座標のdbarで記録・表示します。地図のツールチップと時系列鉛直断面図の重ね合わせにも同じ値を使い、断面図の圧力軸と単位がずれないようにしています。一方、幾何学的な深さへ定量的に換算する場合は、この表示上の取り決めとは別に緯度を考慮した処理が必要です。

MLDの季節変化を読む

多くの中高緯度の海域では、表層の熱の出入りのバランスによって、MLDははっきりとした季節変化を示します。

夏と冬でプロファイルが違う理由から知りたい場合は、計算方法へ進む前に 海の季節は何mの深さまで届く？ を参照してください。

冬に深くなる

冬は、表層の冷却と強い風によって乱流混合が強まり、表層付近の水の浮力が弱まるため、混合がより深くまで届き、混合層は厚くなります。また、地域によっては、この冬の深い混合そのものがモード水を形成する仕組みでもあります。モード水とは？北太平洋亜熱帯モード水を例に形成と見分け方を解説 では、冬の混合層が亜表層の層として残るまでを説明しています。

夏に浅くなる

夏は、太陽による加熱が、混合によってその熱を下方へ分配する速さを上回るため、表層付近には浅く強く成層した層ができ、MLDは浅くなって、場合によっては数十メートル程度になります。

地域による違い

季節変化の大きさは海域によって異なります。黒潮続流の南側にあたる北太平洋西部のように、冬の熱損失が大きい亜熱帯・中緯度域では、季節によるMLDの変動幅が大きくなる傾向があります。一方、季節による強制力が弱い海域や、淡水の流入によって恒常的に強い成層がある海域では、変動幅は小さくなる傾向があります。

例：季節の異なる二つのプロファイルを比べる

似たような場所で観測された、晩冬と晩夏の二つのプロファイルを比べてみましょう。

典型的には、次のような違いが見られます。

• 晩冬のプロファイルでは、水温・塩分がほぼ一様な深い層が100 mをはるかに超えて広がり、密度しきい値と水温しきい値がほぼ同じ深さに達する
• 晩夏のプロファイルでは、一様な層はわずか数十メートルしかなく、その下に鋭い水温躍層が続く
• 両方のプロファイルの、混合層より深いところの水は比較的よく似ている

このように実際のデータを比べると、季節変化を直接確認できます。また、多数のプロファイルを一つずつ全体の形から読むのではなく、上層の状態を一つの数値で手早く比較できるという、MLDの使いやすさもよく分かります。

MLDと他の構造とのつながり

MLDはそれ単体で役立つだけでなく、OceanGraphでは、ほかの派生指標を計算する基準にもなっています。

• 亜表層酸素極大（SOM）：探索範囲は計算されたMLDのすぐ下から始まるため、MLDが変わるとSOMが見つかる場所も変わりうる
• モード水：冬の深い混合によって形成されるため、MLDが深い海域・季節は、モード水の形成が期待される場所でもある

MLDを読むときによくある誤り

MLDを厳密な物理的境界として扱う

しきい値法によるMLDはあくまで推定値であり、正確な物理的な壁ではありません。実際の混合層は緩やかに移り変わることがあり、正確な深さは選んだしきい値や基準深度に左右されます。

MLDが計算できたかを確認せずに比較する

MLDが欠損しているからといって、値がゼロや非常に深いことを意味するわけではなく、表層付近のデータ点が十分でなかったか、どのしきい値にも達しなかったことを意味します。欠損値をそのまま平均計算に混ぜてしまうと、結果が気づかないうちに偏ります。

基準深度を意識しない

すべてのしきい値は、本当の海面ではなく浅い基準深度を基準に測られているため、その基準深度付近で通常と異なる挙動を示すプロファイルは、MLDをそのまま信用する前にもう一度確認する価値があります。

季節や場所を考慮せずにMLDを読む

MLDが20 mというのは、多くの緯度で夏なら珍しくありませんが、冬であれば注目すべき浅さかもしれません。MLDは必ず、そのプロファイルの日付と場所とあわせて解釈しましょう。

従来の流れ：Pythonでプロファイルから計算する

MLDを求める一般的な方法は、ArgoのNetCDFファイルをダウンロードし、GSWなどのツールボックスで水温・塩分をポテンシャル水温・絶対塩分に変換し、選んだしきい値の定義を適用し、多数のプロファイルについてループ処理して時系列や地図を作る、というものです。

これは標準的で柔軟な研究手法ですが、実際には次のような準備が必要です。

• 特定のMLD定義を選び、一貫した形で実装する
• 表層付近のデータが不十分なプロファイルへの対応を考える
• 圧力ではなく物理的な深さが必要な場合は、緯度を考慮した圧力・深度変換を用意する
• 季節や地域をまたいで多数のプロファイルを比較するには、これらすべてを大規模に管理する

ある海域の季節的なMLDのパターンを把握したいだけであれば、全体像を確認する前の準備としては作業が多すぎることがあります。

OceanGraphでMLDの分布を確認する

MLDを計算するコードを書く前に、まずパターンを直接見てみると効果的です。

• 海域と期間で検索する
• MLDの地図表示を有効にして、多数のプロファイルにわたるMLDの変化を一目で確認する
• 同じ海域について、冬の検索結果と夏の検索結果でマーカーの色を見比べる
• 個々のプロファイルを開き、あるMLD値の背景にある水温・塩分構造を確認する

参考ページ：

• Mixed Layer Depth (MLD)（英語）
• Analysis Lab: Vertical Profiles（JSONアップロード、英語）
• 海洋の水温・塩分プロファイルを読み解く
• BGC Argoの溶存酸素プロファイルを読み解く

OceanGraphで混合層深度を探索する

MLDの定義を理解したら、季節や海域ごとにMLDで色分けした検索結果を比べると、実際のArgoデータに現れる季節変化を確認できます。OceanGraphでは、自分でしきい値計算のコードを書かなくても、多数のプロファイルにわたるMLDを一度に見比べられます。

よくある質問

混合層深度には唯一正しい定義がありますか

いいえ。MLDは選んだしきい値基準と基準深度によって変わるため、異なる妥当な定義を使えば、同じプロファイルでも多少異なる深さが得られることがあります。したがって、もっとも重要なのは、プロファイルを比較するときに一貫した定義を使うことです。

なぜOceanGraphは水温・塩分・密度をあわせて使うのですか

OceanGraphは、三つの量を一つの密度に基づく基準へまとめているわけではありません。水温・塩分・密度のしきい値を独立に判定し、三つの交差深度のうち最も浅いものを採用します。そのため、MLDは保守的な、探索用の指標になります。どの基準が最初に反応するかによって、結果が水温だけ、塩分だけ、あるいは密度だけで決まることもあります。それでも三つをあわせて使うことで、どれか一つだけを見る方法では見逃すケースを拾えるようになりますが、その結果は厳密な意味での密度に基づいた定義とは異なります。

プロファイルによってはなぜMLDが表示されないのですか

MLD計算は、プロファイルの表層付近に十分なデータ点がない場合はスキップされ、プロファイル中のどこでもしきい値基準に達しない場合は未定義として扱われます。どちらの場合も、近似値ではなく欠損として表示されます。

MLDは亜表層酸素極大とどう関係していますか

OceanGraphは元の酸素観測範囲全体でピーク形状を評価した後、計算されたMLDより深いピークだけを候補にします。そのため、MLDが変わると、選ばれるピークや、候補が残るかどうかが変わります。さらに、酸素の観測範囲も独立して影響し、同じピークなら、より深くまで、またはより浅い圧力から観測されたプロファイルほどprominenceとhalf-prominence幅が維持または増加して、形状のしきい値を通過しやすくなります。ただし、競合するピークの順位も変わり得るため、最終的に選ばれるとは限りません。

MLDのパターンを探索するのにPythonは必要ですか

必要ありません。独自のMLD定義や大規模な統計解析にはPythonが役立ちますが、まずは季節や地域によるMLDのパターンを目で見て比較し、独自に計算する必要があるかを判断できます。

まとめ

混合層深度は、表層付近の水柱の状態を一つの数値で表すため、さまざまな海洋解析に使われています。MLDを読むときは、しきい値法による推定値であることを理解し、未定義の値を無理に補わず、季節と場所をあわせて確認することが大切です。

MLDの季節変化を理解するには、実際のプロファイルのグラフと地図上のMLD分布を見比べる方法が分かりやすいでしょう。OceanGraphでは、地図上の分布と個々のプロファイルを続けて確認できます。