表面海水中の二酸化炭素分圧の推定値および長期変化傾向の見積もり方法

表面海水中の二酸化炭素分圧は、水温や塩分の変化、下層の海水との混合、生物活動などの影響を受けるため時空間的に大きく変動し、そのすべての変動を観測結果だけで把握することは困難です。そのため、Ishii et al.(2011)の手法に基づき、これまで蓄積された気象庁の海洋観測データを用いて、次の手順で東経137度および東経165度における月ごとの二酸化炭素分圧の推定値と長期変化傾向を求めました。なお、この手法では、海洋の二酸化炭素の観測項目として広く測定されている全炭酸濃度、全アルカリ度、pH、二酸化炭素分圧のうち2つの要素の値が分かれば、化学平衡の関係からそのほかの要素の値が計算で求めることができることを利用しています(二酸化炭素分圧と二酸化炭素濃度の関係については海洋の二酸化炭素の観測項目と方法を参照)。

以下の手順で東経137度および東経165度それぞれの緯度1度ごと、月ごとの表面海水中の二酸化炭素分圧の推定値と、その長期変化傾向を求めます。

1. 表面海水中の全炭酸濃度は、海域や季節によって変化し、水温や塩分との間で特徴的な相関関係を示します。その相関関係を利用し、表面海水中の全炭酸濃度を時間・水温・塩分の経験式で表します。
2. 1の経験式に海面水温と海面塩分のデータ（下表の解析値）を与えることで月ごとの表面海水中の全炭酸濃度を求めます。
3. この海域では、塩分で規格化した表面海水中の全アルカリ度に長期的な変化傾向が見られないというTakatani et al.(2014)の報告を基に、診断の対象としている海域の解析期間における全アルカリ度を一定と仮定します。
4. 2で得られた月ごとの全炭酸濃度と3の全アルカリ度から表面海水中の二酸化炭素分圧を算出し、表面海水中の二酸化炭素分圧の長期変化傾向を求めます。

表面海水中の二酸化炭素分圧の推定値算出と長期変化傾向の見積もりに使用したデータは下表のとおりです。

上記の手法で得られた二酸化炭素分圧の推定値の不確かさは、観測で得られた二酸化炭素分圧との誤差を基に評価しています。下図は東経137度および東経165度における表面海水中の二酸化炭素分圧の推定値の、観測値に対する誤差から求めた平均二乗誤差（RMSE）と誤差のヒストグラムを表しています。ここでRMSEとは、個々の推定値の誤差を二乗し期間平均してから平方根を取ったもので、推定誤差の標準的な大きさを示す指数です。RMSEは次の式で表すことができます。