Statistical inversion of South Atlantic circulation in an abyssal neutral density layer
Statistical inversion of South Atlantic circulation in an abyssal neutral density layer
复制标题
DOI：
10.1357/0022240054663240
复制
发表时间：
2005-07-01
影响因子：
0.5
通讯作者：
Herbei, R
中科院分区：
文献类型：
--
作者：
McKeague, IW;Nicholls, G;Herbei, R
本文介绍了一种估计稳态海洋环流和示踪剂场的贝叶斯反演方法。它是基于一个准水平流模型和一个偏微分方程求解器的示踪剂场对流扩散方程的计算解决方案的正问题。现有的海洋环流反演方法的一个典型特征是预处理阶段，其中示踪剂数据在规则网格上插值，并且在随后的反演中忽略插值误差。我们的方法只使用插值数据在这些网格点，有相邻的水文站。通过以综合方式利用基于物理的模型，该方法提供了统计上统一的反演和示踪剂场重建，具有最小的数据平滑。解决这个问题包括在存在一些假设（先验信息）的情况下找到有关环流和示踪剂场的信息;由此产生的后验概率分布总结了我们对这些场的了解。我们开发了一个马尔可夫链蒙特卡罗模拟过程中提取的信息（分析棘手）后验分布的模型中的所有参数的“解决方案”的不确定性表示的变化，在模拟运行的输出。我们的方法的目的是寻找时间平均准水平流和示踪剂领域的深海中性密度层在南大西洋。
This paper introduces a Bayesian inversion approach to estimating steady state ocean circulation and tracer fields. It is based on a quasi-horizontal flow model and a PDE solver for the forward problem of computing solutions to the tracer field advection-diffusion equations. A typical feature of existing ocean circulation inverse methods is a preprocessing stage in which the tracer data are interpolated over a regular grid and the interpolation error is ignored in the subsequent inversion. Our approach only uses interpolated data at those grid points that have neighboring hydrographic stations. By exploiting physically-based models in an integrated fashion, the method provides a statistically unified inversion and tracer field reconstruction with minimal data smoothing. Solving the problem consists of finding information about the circulation and tracer fields in the presence of a number of assumptions (prior information); the resulting posterior probability distribution summarizes what we can know about these fields. We develop a Markov chain Monte Carlo simulation procedure to extract information from the (analytically intractable) posterior distribution of all the parameters in the model; uncertainty about the "solution" is represented by variation in the output of the simulation runs. Our approach is aimed at finding the time-averaged quasi-horizontal flow and tracer fields for an abyssal neutral density layer in the South Atlantic.