State Estimation of the North Pacific Ocean by a Four-Dimensional Variational Data Assimilation Experiment

A four-dimensional variational data assimilation system has been applied to an experiment to describe the dynamic state of the North Pacific Ocean. A synthesis of available observational records and a sophisticated ocean general circulation model produces a dynamically consistent dataset, which, in contrast to the nudging approach, provides realistic features of the seasonally-varying ocean circulation with no artificial sources/sinks for temperature and salinity fields. This new dataset enables us to estimate heat and water mass transports in addition to the qualification of water mass formation and movement processes. A sensitivity experiment on our assimilation system reveals that the origin of the North Pacific Intermediate Water can be traced back to the Sea of Okhotsk and the Bering Sea in the subarctic region and to the subtropical Kuroshio region further south. These results demonstrate that our data assimilation system is a very powerful tool for the identification and characterization of ocean variabilities and for our understanding of the dynamic state of ocean circulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Coupling a two-way nested primitive equation model and a statistical SST predictor of the Ligurian Sea via data assimilation

A primitive equation model and a statistical predictor are coupled by data assimilation in order to combine the strength of both approaches. In this work, the system of two-way nested models centred in the Ligurian Sea and the satellite-based ocean forecasting (SOFT) system predicting the sea surface temperature (SST) are used. The data assimilation scheme is a simplified reduced order Kalman filter based on a constant error space. The assimilation of predicted SST improves the forecast of the hydrodynamic model compared to the forecast obtained by assimilating past SST observations used by the statistical predictor. This study shows that the SST of the SOFT predictor can be used to correct atmospheric heat fluxes. Traditionally this is done by relaxing the model SST towards the climatological SST. Therefore, the assimilation of SOFT SST and climatological SST are also compared.     

Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

伴随模式同化系统在修正模式地形中的应用

任何模式都只是实际大气的一种近似模拟，存在误差。鉴于传统的四维资料伴随模式同化系统都是假设模式完全正确仅对初始场进行修正。将伴随模式同化系统应用于修正模式地形误差，通过对不同初始地形的修正试验表明，MM5伴随模式同化系统能很好地对地形误差进行修正，能够反演出与初始气象要素场分辨率相匹配、与模式更协调的地形场，为伴随模式同化系统的广泛应用提供了一种新的思路。     

变分同化方法在Lorenz系统中的简单应用研究

利用Lorenz模式作变分同化数值试验，通过对一个简单系统的讨论，介绍四维变分同化方法。对初值敏感性和观测点的个数及观测值作了对比试验，发现随着模式对初值敏感性的增加，同化效果会越来越差；观测点越少，观测值误差越大，这些都会影响同化效果，甚至导致同化失败。     

改进的切线性模式对一个边界层模式变分资料同化的改善

大气、海洋中边界层物理过程的强非线性会对切线性模式近似（尤其是当积分时间长时）有较大的影响，从而给相应的4维变分资料同化问题的求解造成困难。本构造了Mellor-Yamada湍流封闭模式（level 2.5)的一个改进的切线性模式，相比通常的切线性模式和简化的切线性模式可以提高对非线性扰动的逼近。利用这个改进的切线性模式的伴随模式分别进行了1到7天的变分资料同化试验，得到了满意的结果。而用通常的伴随模式和简化的伴随模式都无法得到令人满意的结果。     

On the oceanic impact of a data-assimilation method in a coupled ocean-land-atmosphere model

 Based on the Kalman filter theory, a new data-assimilation method has been used to improve the 3-D oceanic temperature field of the Center for Ocean-Land-Atmosphere Studies (COLA) coupled general circulation model. This method is applied to assimilate surface and subsurface temperature of in situ measurements from the Pilot Research Moored Array in the Tropical Atlantic project (PIRATA). The assimilation of the PIRATA data produces an improved representation of the thermal state of the ocean and allows a better estimation of other oceanographic quantities, like meridional heat fluxes and zonal currents. The present paper focuses on the tropical Atlantic and, in particular, it contains new reconstructed temperature profiles. One-month forecast experiments during 1999 were performed and the impact of the assimilation is discussed. Received: 24 April 2001 / Accepted: 8 March 2002     

求解变分不等式的Newton迭代的半局部收敛性分析

分析了求解变分不等式Newton方法的半局部收敛性，建立了类似于Kantorovich定理的收敛性结果。该结果不仅为判断Newton方法的收敛性提供了可计算的充分条件，也给出了Newton方法的收敛域以及问题解的存在区域。同时，文章也得到了Newton方法的若干收敛性质，包含收敛阶以及可计算的误差估计式等。