A generalized model of pelagic biogeochemistry for the global ocean ecosystem. Part II: Numerical simulations

This paper presents a global ocean implementation of a multi-component model of marine pelagic biogeochemistry coupled on-line with an ocean general circulation model forced with climatological surface fields (PELAgic biogeochemistry for Global Ocean Simulations, PELAGOS). The final objective is the inclusion of this model as a component in an Earth System model for climate studies. The pelagic model is based on a functional stoichiometric representation of marine biogeochemical cycles and allows simulating the dynamics of C, N, P, Si, O and Fe taking into account the variation of their elemental ratios in the functional groups. The model also includes a parameterization of variable chlorophyll/carbon ratio in phytoplankton, carrying chl as a prognostic variable. The first part of the paper analyzes the contribution of non-local advective–diffusive terms and local vertical processes to the simulated chl distributions. The comparison of the three experiments shows that the mean chl distribution at higher latitudes is largely determined by mixing processes, while vertical advection controls the distribution in the equatorial upwelling regions. Horizontal advective and diffusive processes are necessary mechanisms for the shape of chl distribution in the sub-tropical Pacific. In the second part, the results have been compared with existing datasets of satellite-derived chlorophyll, surface nutrients, estimates of phytoplankton community composition and primary production data. The agreement is reasonable both in terms of the spatial distribution of annual means and of the seasonal variability in different dynamical oceanographic regions. Results indicate that some of the model biases in chl and surface nutrients distributions can be related to deficiencies in the simulation of physical processes such as advection and mixing. Other discrepancies are attributed to inadequate parameterizations of phytoplankton functional groups. The model has skill in reproducing the overall distribution of large and small phytoplankton but tends to underestimate diatoms in the northern higher latitudes and overestimate nanophytoplankton with respect to picoautotrophs in oligotrophic regions. The performance of the model is discussed in the context of its use in climate studies and an approach for improving the parameterization of functional groups in deterministic models is outlined.     

基于Vega的虚拟海洋视景仿真研究

虚拟海洋视景仿真系统足虚拟现实技术在军事领域中的应用.由于海浪运动的随机性,虚拟海洋的计算机生成成为视景仿真领域的研究难点.文中详细讨论了两种海浪模型,重点研究Vega中海浪原理和海洋模块的计算机实现,基于Vega软件生成海浪.从效果图中可以看出,模拟的海浪效果比较好,能够满足实际应用的需求.     

Operational surface drift prediction using linear and non-linear hyper-ensemble statistics on atmospheric and ocean models

Multimodel super-ensemble forecasts, which exploit the power of an optimal local combination of individual models usually show superior forecasting skills when compared to individual models because they allow for local correction and/or bias removal. Deterministic approaches to the problem of surface drift are often limited by strong assumptions on the underlying physics. A new approach based on linear and non-linear optimization is proposed, using hyper-ensemble deduced statistics to forecast at short time scale Lagrangian drifts from combined atmospheric and ocean operational models and local observations that were made available during the MREA04 field experiment along the West coast of Portugal. Optimization methods are based on a training/forecast cycle. The performance and the limitations of the hyper-ensembles and the individual models are discussed. Results suggest that our statistical methods reduce the position errors significantly for 12 to 48 h forecasts and hence compete with pure deterministic approaches.     

Evaluation of two atmospheric models for wind–wave modelling in the NW Mediterranean

The NW Mediterranean experiences, as illustrated by the last decade, strong and rapidly varying storms with severe waves and winds. This has motivated a continuous validation of models and the efforts to improve wave and wind predictions. In this paper we use two atmospherics models, MASS (from SMC-Meteorological Office of Catalunya) and ARPEGE (from Météo-France), to force two third generation wave models: WAM and SWAN. The evaluation and comparison has been carried out for two severe storms registered in November 2001 and March–April 2002.The ARPEGE and MASS models predicted higher 10 m wind speeds than coastal meteorological stations, a fact attributed to local land influences. Regarding the 10 m wind direction, models do not present large differences, although considerable deviations from recorded data were found during some dates. ARPEGE presents less scatter and lower errors than MASS when compared with QuikSCAT data.The 10m wind fields from both atmospheric models were used to force the two selected wave models and analyse the errors and sensitivities when predicting severe wave storms. The wave model simulations show some interesting results; during the storm, the spatial wave pattern using ARPEGE showed a higher maximum, although the values of significant wave height at the buoys were lower than the ones forced by MASS (with both WAM and SWAN). The SWAN simulations show a better agreement in predicting the growing and waning of the storm peaks. The prediction of mean period was improved when using the ARPEGE wind field. However the underestimation by SWAN due to the large energy at high frequencies was evident. Validation of spectral shape predictions showed that it still has considerable error when predicting the full frequency spectra. The storms showed bimodal spectral features which were not always reproduced by wave models and are likely to be responsible for part of the discrepancies.     

海洋声场预报及战术应用

随着科学技术的飞速发展,对海洋声场做出合乎要求的快速精确预报显得尤为重要。射线模型、抛物线模型、简正波等模型适合于不同具体的海洋声场环境,这些传播模型可以实现声场数值预报的实时性。基于上述模型建立的海洋声场预报系统及对其战术应用进行的研究,对声呐的设计、使用及舰艇作战有着重要的意义。本文在上述模型基础上利用Matlab和Fortran Powerstation工具建立了海洋声场预报系统,使用声场预报系统对典型的海洋环境,即海洋锋、涡漩的声场能量传播规律进行了研究,并提出了在反潜作战中对海洋锋、涡漩、深海会聚区和影区等海洋环境的战术应用。     

Ventilation of bottom water in the North Sea–Baltic Sea transition zone

The transition zone between the North Sea and the Baltic Sea is a highly dynamic region where a general estuarine circulation forms a regional scale frontal system from northern Kattegat to the Arkona Sea. This system is characterized by an upper low saline out?owing Baltic water mass from the in?owing saline Skagerrak bottom water to the Kattegat and Belt Sea area. Large and rapid ?uctuations of the frontal system are caused by barotropic transports, forced by changing sea level difference between northern Kattegat and the western Baltic Sea, and this results in high variability of the hydrographic conditions and also in frequent in- and out?ow events to the Baltic. The dynamics in the region are here analyzed by a regional model of the transition zone, covering the area from the northern Kattegat to the Arkona Sea. The model is validated against water level, temperature and salinity measurements from the region, and the transports through the Danish straits are related to previous estimates and empirical relations. A sensitivity study quantify the role of bathymetry, the tidally induced mixing and the in?owing Skagerrak bottom water for ventilating the bottom water with Skagerrak water or surface water.Furthermore, the dynamics in the region is analyzed with tracers representing the age of the water. The distribution of age tracers with different boundary conditions are analyzed, and the role of advection and mixing for ventilating the bottom water is quanti?ed in terms of the water age. It is shown that the Great Belt area is a very dynamical area where bottom water is ventilated with surface water. The interannual variation of the ventilation of bottom water in the period 2001–2003 is analyzed by various age tracers and related to observed oxygen conditions, and it is shown that the extreme hypoxic event in the autumn 2002 in the southern Kattegat, the Great Belt and in the western Baltic Sea coincide with an unusual low vertical ventilation rate in the Great Belt area, but normal advection rates of bottom water from the northern Kattegat. This indicates that during this particular event, and probably in general, ventilation of bottom water in the Great Belt has signi?cant in?uence on oxygen conditions in the southern part of the region and for ventilation of bottom waters in the western Baltic Sea. In contrast, the central Kattegat is primarily ventilated by advection of bottom water from the Skagerrak. An age tracer representing the ventilation rate of bottom water with either Skagerrak water or surface water is shown to be inversely correlated to the observed oxygen distribution in the region.     

海洋船舶的轨迹估计算法

准确估计和预测海洋导航中机动船舶的轨迹是改善海上安全和保安的重要工具。因此,许多常规海洋导航系统和船舶交通管理和报告服务都为此目的配备了雷达设施。然而,船舶操纵轨迹的预测的准确性主要取决于船舶位置,速度和加速度的估计的良好性。因此,本研究提出了一个基于曲线运动模型的机动海洋船只模型,其基于用于相同目的的线性位置模型的测量。此外,还假定与白高斯噪声相关联的系统状态和测量模型。扩展卡尔曼滤波器被提出作为用于估计位置,速度和加速度的自适应滤波器算法,用于预测机动的海洋船舶轨迹。最后,提出的模型被实施,并获得成功的计算结果,在本研究中的海洋导航中的船舶操纵轨迹的预测。     

A non-geometrically similar model for predicting the wake field of full-scale ships

The scale effect leads to large discrepancies between the wake fields of model-scale and actual ships, and causes differences in cavitation performance and exciting forces tests in predicting the performance of actual ships. Therefore, when test data from ship models are directly applied to predict the performance of actual ships, test results must be subjected to empirical corrections. This study proposes a method for the reverse design of the hull model. Compared to a geometrically similar hull model, the wake field generated by the modified model is closer to that of an actual ship. A non- geometrically similar model of a Korean Research Institute of Ship and Ocean Engineering (KRISO)’s container ship (KCS) was designed. Numerical simulations were performed using this model, and its results were compared with full-scale calculation results. The deformation method of getting the wake field of full-scale ships by the non-geometrically similar model is applied to the KCS successfully.     

用于船舶伴流场预报的一种非几何相似的模型（英文）

The scale effect leads to large discrepancies between the wake fields of model-scale and actual ships, and causes differences in cavitation performance and exciting forces tests in predicting the performance of actual ships. Therefore, when test data from ship models are directly applied to predict the performance of actual ships, test results must be subjected to empirical corrections. This study proposes a method for the reverse design of the hull model. Compared to a geometrically similar hull model, the wake field generated by the modified model is closer to that of an actual ship. A non- geometrically similar model of a Korean Research Institute of Ship and Ocean Engineering(KRISO)'s container ship(KCS) was designed. Numerical simulations were performed using this model, and its results were compared with full-scale calculation results. The deformation method of getting the wake field of full-scale ships by the non-geometrically similar model is applied to the KCS successfully.     

潜艇在海洋密度锋中受绕运动建模与仿真研究

根据分析潜艇通过海洋密度锋区的受力特点,建立潜艇在海洋密度锋中的受力和垂直面运动模型。通过潜艇垂直面运动仿真,分析自由运动、PID操舵控制、操舵加均衡控制和综合运用增速、操舵加均衡控制等情形下海洋密度锋对潜艇纵倾角、潜浮角和升沉运动的影响,结果表明综合运用操舵、增速和均衡进行控制可以更为有效地控制潜艇在海洋锋区航行深度和运动姿态。     

Mesoscale subduction at the Almeria–Oran front: Part 1: Ageostrophic flow

This paper presents a detailed diagnostic analysis of hydrographic and current meter data from three, rapidly repeated, fine-scale surveys of the Almeria–Oran front. Instability of the frontal boundary, between surface waters of Atlantic and Mediterranean origin, is shown to provide a mechanism for significant heat transfer from the surface layers to the deep ocean in winter. The data were collected during the second observational phase of the EU funded OMEGA project on RRS Discovery cruise 224 during December 1996. High resolution hydrographic measurements using the towed undulating CTD vehicle, SeaSoar, traced the subduction of Mediterranean Surface Water across the Almeria–Oran front. This subduction is shown to result from a significant baroclinic component to the instability of the frontal jet. The Q-vector formulation of the omega equation is combined with a scale analysis to quantitatively diagnose vertical transport resulting from mesoscale ageostrophic circulation. The analyses are presented and discussed in the presence of satellite and airborne remotely sensed data; which provide the basis for a thorough and novel approach to the determination of observational error.     

A new approach for deriving temperature and salinity fields in the Indian Ocean using artificial neural networks

This work reports a new methodology for deriving monthly averages of temperature (T) and salinity (S) fields for the Indian Ocean based on the use of an artificial neural network (ANN). Investigation and analysis were performed for this region with two distinct datasets: (1) monthly climatological data for T and S fields (in 1° × 1° grid boxes) at standard depth levels of the World Ocean Atlas 1994 (WOA94), and; (2) heterogeneous randomly distributed in situ ARGO, ocean station data (OSD) and profiling (PFL) floats. A further numerical experiment was conducted with these two distinct datasets to train the neural network model. Nonlinear regression mapping utilizing a multilayer perceptron (MLP) is employed to tackle nonlinearity in the data. This study reveals that a feed-forward type of network with a resilient backpropagation algorithm is best suited for deriving T and S fields; this is demonstrated by independently using WOA94 and in situ data, which thus tests the robustness of the ANN model. The suppleness of the T and S fields derived from the ANN model provides the freedom to generate a new grid at any desired level with a high degree of accuracy. Comprehensive training, testing and validation exercises were performed to demonstrate the robustness of the model and the consistency of the derived fields. The study points out that the parameters derived from the ANN model using scattered, inhomogeneous in situ data show very good agreement with state-of-the-art WOA climatological data. Using this approach, improvements in ocean climatology can be expected to occur in a synergistic manner with in situ observations. Our investigation of the Indian Ocean reveals that this approach can be extended to model global oceans.     

Open boundary conditions in stratified ocean models

In most applications of numerical ocean models, artificial boundaries are introduced to limit the domain. Along such a boundary we need to apply what is often referred to as an open boundary condition (OBC). In this paper a number of local methods used in barotropic ocean models are applied and discussed for the stratified case using a normal mode approach. The OBCs are the simple conditions: clamped, prescribed and zero gradient; the radiation conditions: Camerlengo–O'Brien, Orlanski and a method of characteristics based on linear equations; and a sponge type condition: the flow relaxation scheme. The OBCs have been implemented in a 3-layer ocean model and examples of how the various OBCs perform for three simple flow situations are investigated. The cases are: internal wave radiation, a quasi-steady coastal jet and the response to a storm moving across a strait. It is found that the flow relaxation scheme and the method based on characteristics perform well for the test cases in general, although some of the simpler methods give better results in individual cases.     

Gaussian anamorphosis extension of the DEnKF for combined state parameter estimation: Application to a 1D ocean ecosystem model

We consider the problem of combined state-parameter estimations in biased nonlinear models with non-Gaussian extensions of the Deterministic Ensemble Kalman Filter (DEnKF). We focus on the particular framework of ocean ecosystem models. Such models present important obstacles to the use of data assimilation methods based on Kalman filtering due to the non-linearity of the models, the constraints of positiveness that apply to the variables and parameters, and the non-Gaussian distribution of the variables in which they result.We present extensions of the DEnKF dealing with these difficulties by introducing a nonlinear change of variables (anamorphosis function) in order to execute the analysis step with Gaussian transformed variables and parameters. Several strategies to build the anamorphosis functions are investigated and compared within the framework of twin experiments realized in a simple 1D ocean ecosystem model. A solution to the problem of the specification of the observation error for transformed observations is suggested. The study highlights the inability of the plain DEnKF with a simple post-processing of the negative values to properly estimate parameters when constraints of positiveness apply to the variables. It goes on to show that the introduction of the Gaussian anamorphosis can remedy these assimilation biases.     

A numerical model of the long term flow along the Malin-Hebrides shelf

A three dimensional hydrodynamic model of the Malin-Hebrides shelf region is used to investigate the spatial variability of the wind and tidally induced residual flow in the region and the influence of flow from the Irish Sea and along the shelf edge. By this means it is possible to understand the spatial variability in the long term observed flow fields in the region and the range of driving forces producing this flow. The model uses a sigma coordinate grid in the vertical with a finer grid in the near surface and near bed shear layers. The vertical diffusion of momentum in the model is parameterised using an eddy viscosity coefficient which is derived from turbulence energy closure models. Two different turbulence models are used to compute the eddy viscosity, namely a two-equation (itq²−q²ℓ) model which has prognostic equations for both turbulence energy and mixing length and a simpler model in which the mixing length is a specified algebraic function of the water depth.The wind induced response to spatially and temporally constant orthogonal wind stresses, namely westerly and southerly winds of 1 N m⁻², are derived from the model. By using orthogonal winds and assuming linearity, then to first order the response to any wind direction can be derived. Computed flows show a uniform wind driven surface layer of magnitude about 3% of the wind speed and direction 15 ° to the right of the wind, in deep water. Currents at depth particularly in the shelf edge and near coastal region show significant spatial variability which is related to variations in bottom topography and the coastline.Calculations show that tidal residual flows are only significant in the near coastal regions where the tidal current is strong and exhibits spatial variability. Flow into the region from the Irish Sea through the North Channel although having its greatest influence in the near coastal region, does affect currents near the shelf edge region. Again the spatial variability of the flow is influenced by topographic effects.A detailed examination of wind induced current profiles together with turbulence, mixing length and viscosity, at a number of locations in the model from deep ocean to shallow near coastal, shows that both turbulence models yield comparable results, with the mixing length in the two equation model showing a similar dependence to that specified in the simpler turbulence model.Calculations clearly show that flow along the shelf edge area to the west of Ireland and from the Irish Sea entering the region, together with local wind forcing can have a major effect upon currents along the Malin-Hebrides shelf. The flow fields show significant spatial variability in the region, comparable to those deduced from long term tracer measurements. The spatial variability found in the calculations suggests that a very intense measurement programme together with inflow measurements into the area is required to understand the circulation in the region, and provide data sets suitable for a rigorous model validation.     

High-resolution synthetic monitoring by a 4-dimensional variational data assimilation system in the northwestern North Pacific

A data assimilation system is applied to the integrated monitoring of oceanic states in the northwestern North Pacific by combining a high resolution ocean general circulation model with an adjoint method. A comparison of assimilation results with observations shows that the system is better able to represent synoptic features of ocean circulation than do models or data alone. Furthermore, meso-scale features associated with frontal structures and eddies, which are often seen in the Kuroshio and Oyashio extension regions and the Sea of Japan, are better defined in the assimilation results. These features suggest that our 4D-VAR high-resolution data assimilation system is capable of providing time series data which satisfy the model physics and fit the observations, and hence the ocean state derived from our system has greater information content than that obtained from earlier methods.     

海洋平台的振动试验方案及其模型设计

针对海洋平台的振动具有低频、高幅的特点，提出了振动主动控制的试验方案。振动主动控制可改善台上人员的工作条件、提高生产效率、延长海洋平台及台载机械使用寿命，已成为提升海洋平台综合性能的潜在技术手段。该试验装置包括振动试验模型、激振器、作动器、主机及DSP板、低通滤波器、信号放大器等。根据对海洋平台的类型、特点及振动控制研究的分析，设计了用于海洋平台振动主动控制试验的两个模型，一个用于研究水平面两个方向的振动，另一个用于研究纵向振动。海洋平台振动试验模型采用柔性较强的材料，支架模型采用刚性较强的槽钢。对于第一个模型进行了必要的强度和刚度计算。设计的模型可用于海洋平台的概念性振动主动控制试验，对海洋平台的振动主动控制提供理论指导。     

Development of ocean carriers' behaviour model focusing on their 'cost and tariff' based on the spatial general equilibrium

Spatial computable general equilibrium (SCGE) theory has been applied to an international trade model to evaluate tariff and fiscal policies of a country. Those models can not be applied for the evaluation of transport policies such as port development and tax/subsidy policy against transportation sectors, since they do not deal with ocean freight and ocean carriers explicitly in the model. Ocean freight often varies with the changes of competitive conditions and/or demand/ supply balances in the short run, while it should reflect the actual expenditure of the carriers in the long run. The model proposed here considers the profit maximization behaviour of ocean carriers, and deals with ocean freight explicitly. The model is applied to four major economic regions; Japan, USA, EU and Asia. A multi-level function composed of the Cobb-Douglas function is adopted to produce reliable parameters of the production function for many industries.     

Development of ocean carriers' behaviour model focusing on their 'cost and tariff' based on the spatial general equilibrium

Spatial computable general equilibrium (SCGE) theory has been applied to an international trade model to evaluate tariff and fiscal policies of a country. Those models can not be applied for the evaluation of transport policies such as port development and tax/subsidy policy against transportation sectors, since they do not deal with ocean freight and ocean carriers explicitly in the model. Ocean freight often varies with the changes of competitive conditions and/or demand/ supply balances in the short run, while it should reflect the actual expenditure of the carriers in the long run. The model proposed here considers the profit maximization behaviour of ocean carriers, and deals with ocean freight explicitly. The model is applied to four major economic regions; Japan, USA, EU and Asia. A multi-level function composed of the Cobb-Douglas function is adopted to produce reliable parameters of the production function for many industries.     

Decreasing oxygen concentrations in the Gullmar Fjord, Sweden, as confirmed by benthic foraminifera, and the possible association with NAO

The deep Gullmar Fjord on the Swedish west coast may be subject to a trend of decreasing oxygen content of the bottom-water. Due to a non-uniform frequency of oxygen measurements before and after 1980, and small differences in the minimum concentrations, it is difficult to assess with certainty if this is true or not. We have used benthic foraminifera from three ultra-high resolution sediment records to try to find out if such a trend exists. The sediment records comprise the time between 1930 and 1996, with an almost year-based temporal resolution. The foraminiferal records show a significant change in faunal composition during the mid-1970's when it altered from a common Skagerrak fauna to an assemblage characterized by an opportunistic, low oxygen tolerant species; Stainforthia fusiformis. This species is an indicator of low oxygen environments of many Scandinavian fjords. The timing of this faunal change coincides with one of the first recorded severe low oxygen events in the Gullmar Fjord in the winter of 1973/1974. The change suggests that the Gullmar Fjord has actually been subject to a decrease in bottom-water oxygen concentrations. The reason for this change is not clear but here, for the first time, we present sediment records that suggest that there is a causal relation between variations in the North Atlantic Oscillation (NAO) indices and the marine environment. In the early 1970s the NAO-indices altered from a mainly negative phase towards a mainly positive period which is still ongoing. Positive NAO-indices (winter values) in the Skagerrak region are generally marked by westerly winds which normally prevent the exchange of bottom-water in the fjords, i.e., the oxygenation of the sea floor and the supply of oxygen to benthic life.