基于HYCOM的海洋环流仿真建模研究

HYCOM海洋模式是垂向混合坐标模式的代表。基于HYCOM海洋模式,采用三层嵌套技术对东海南部海域的海洋环流进行仿真建模研究。通过对仿真结果的验证表明,HYCOM模式的垂向混合坐标方案能够很好的模拟地形复杂海域的环流。     

海洋模式中潮致混合调和分析参数化方案的建立

潮致混合对海洋环流的调整起着重要作用。陆架环流的数值模拟中如果不考虑潮汐作用,往往不能得到与观测相符的垂向温盐结构。本文基于调和分析方法,建立了一套潮致混合参数化方案。该方案通过对垂向混合系数进行调和分析,从而得到随时间变化的潮致混合系数。将该方案用于黄海冷水团数值模拟的结果显示,其能够得到与在数值模式开边界直接加入潮汐强迫相当的冷水团温盐结构。和直接引入潮汐强迫相比,这一潮致混合参数化方案的优势在于,它能够大大节省数值模拟计算机时,因此有望显著提高大规模高分辨率的海洋环流及气候模式的模拟能力和计算效率。     

A modeling study of the characteristics and mechanism of the westward coastal current during summer in the northwestern South China Sea

The characteristics and dynamical mechanism of summer-time coastal current over the North South China Sea shelf have been investigated based on a high resolution unstructuredgrid finite volume community ocean model (FVCOM). Modeldata comparison demonstrates that the model describes and explains well the coastal dynamics over the North South China Sea shelf. The coastal current on the North South China Sea shelf is greatly influenced by monsoon and the freshwater discharge of the Pearl River. Strong southwesterly wind drives the coastal current northeastward. However, under weak southwest monsoon, the coastal current west of the Pearl River estuary (PRE) advects toward the southwest, and splits into two parts when reaching east of the Qiongzhou Strait, with one branch entering the Gulf of Tonkin through the Qiongzhou Strait, transporting low salinity water into the Gulf of Tonkin, and the other part flows cyclonic and interacts with the northeastward current around southeast of Hainan Island, forming a cyclonic eddy east of the Qiongzhou Strait. A variety of model experiments focused on freshwater discharge, wind forcing, tidal rectification, and stratification are performed to study the physical mechanism of the southwestward coastal current which is usually against the summer wind. Process-oriented experiment results indicate that the southwest monsoon and freshwater discharge are important factors influencing the formation of southwestward coastal current during summer.     

Baroclinic wind adjustment processes in the Mediterranean Sea

The wind-driven general circulation of the Mediterranean Sea is studied using a primitive equation model. The model uses a 0.25° horizontal resolution and eight or 16 levels in the vertical. The model uses the Mediterranean basin geometry, and the Strait of Gibraltar is closed. The vertical density structure is initialized with annual average data, and the temperature and salinity values are fixed at the surface to simulate perpetual annual mean conditions. The wind forcing consists of monthly mean climatological stresses.The results show that the general circulation of the Mediterranean Sea has a multiple time-scale character (seasonal excursions and steady state amplitudes are comparable) and it is composed by sub-basin scale gyres corresponding to the scale of the wind stress curl centers. The steady state circulation (annual mean average) is determined by a Sverdrup balacne modified by viscous effects.The unsteady vertically integrated transport circulation consists of sub-basin scale gyres similar to the steady state transport components, which amplify seasonally and the partial or total reversal of the currents in many subportions of the basin. The gyres can be stationary in position or propagating. This seasonal ocean response is partly constituted by Rossby modes due to the wind stress curl annual harmonic. The baroclinic circulation shows the seasonal shift of the North African Current from a position along the African coasts during winter to the center of the Balearic and Ionian basin during summer.     

Evaluation of a dynamically downscaled atmospheric reanalyse in the prospect of forcing long term simulations of the ocean circulation in the Gulf of Lions

The paper evaluates atmospheric reanalysis as possible forcing of model simulations of the ocean circulation inter-annual variability in the Gulf of Lions in the Western Mediterranean Sea between 1990 and 2000. The sensitivity of the coastal atmospheric patterns to the model resolution is investigated using the REMO regional climate model (18 km, 1 h), and the recent global atmospheric reanalysis ERA40 (125 km, 6 h). At scales from a few years to a few days, both atmospheric data sets exhibit a very similar weather, and agreement between REMO and ERA40 is especially good on the seasonal cycle and at the daily variability scale. At smaller scales, REMO reproduces more realistic spatio-temporal patterns in the ocean forcing: specific wind systems, particular atmospheric behaviour on the shelf, diurnal cycle, sea-breeze. Ocean twin experiments (1990–1993) clearly underline REMO skills to drive dominant oceanic processes in this microtidal area. Finer wind patterns induce a more realistic circulation and hydrology of the shelf water: unique shelf circulation, upwelling, temperature and salinity exchanges at the shelf break. The hourly sampling of REMO introduces a diurnal forcing which enhances the behaviour of the ocean mixed layer. In addition, the more numerous wind extremes modify the exchanges at the shelf break: favouring the export of dense shelf water, enhancing the mesoscale variability and the interactions of the along slope current with the bathymetry.     

Cascades of dense water around the world ocean

Dense water overflow off continental shelves (cascading) is one of the contributing processes of shelf-deep ocean exchange, and of topical interest to climate studies and nutrient fluxes. Dense water originating from cooling, evaporation, freezing and salinization on a shallow shelf spills over the shelf edge and may develop as near-bottom gravity current or an intermediate-depth intrusion. It is difficult to observe in nature due to its intermittent character.This paper provides an extensive inventory of observed cases of water cascades around the World Ocean, summarises their locations and individual properties, and provides statistics of the identified cases. The search for cascading was carried out using oceanographic databases and a literature review. This study identified 61 confirmed cases world-wide, including 25 cases in the Arctic seas, 12 at mid-latitudes, seven in sub-tropical and tropical regions, and 17 off the Antarctic shelves. Eighteen cascades had not been reported before. We analyze a set of numerical parameters of dense water cascades, allowing us to quantify, compare and contrast the properties of water cascades. The overall average density contrast between the confirmed cascades and ambient water is 0.37 (kg/m³); it can be as much as 2 (kg/m³) on some Arctic shelves. Frequently initiated by strong cooling at the surface, cascades often remain colder through the descent, thus supplying the deep ocean with colder and fresher water. In non-dimensional variables, the data from all climate zones fit well to a unique curve, which represents a relationship between a cascade’s internal structure and the parameters describing its forcing. On average, the down-slope volumetric flux provided by dense water cascades is estimated as 0.05 to 0.08 Sv per 100 km of shelf edge.Regional terms: Arctic, Antarctic, North Atlantic Ocean, Barents Sea, Mediterranean Sea, Skagerrak, Tasman Sea, Sea of Okhotsk     

Advantages of vertically adaptive coordinates in numerical models of stratified shelf seas

Shelf seas such as the North Sea and the Baltic Sea are characterised by spatially and temporally varying stratification that is highly relevant for their physical dynamics and the evolution of their ecosystems. Stratification may vary from unstably stratified (e.g., due to convective surface cooling) to strongly stratified with density jumps of up to 10 kg/m³ per m (e.g., in overflows into the Baltic Sea). Stratification has a direct impact on vertical turbulent transports (e.g., of nutrients) and influences the entrainment rate of ambient water into dense bottom currents which in turn determine the stratification of and oxygen supply to, e.g., the central Baltic Sea. Moreover, the suppression of the vertical diffusivity at the summer thermocline is one of the limiting factors for the vertical exchange of nutrients in the North Sea. Due to limitations of computational resources and since the locations of such density jumps (either by salinity or temperature) are predicted by the model simulation itself, predefined vertical coordinates cannot always reliably resolve these features. Thus, all shelf sea models with a predefined vertical coordinate distribution are inherently subject to under-resolution of the density structure.To solve this problem, Burchard and Beckers (2004) and Hofmeister et al. (2010) developed the concept of vertically adaptive coordinates for ocean models, where zooming of vertical coordinates at locations of strong stratification (and shear) is imposed. This is achieved by solving a diffusion equation for the position of the coordinates (with the diffusivity being proportional to the stratification or shear frequencies). We will show for a coupled model system of the North Sea and the Baltic Sea (resolution ∼ 1.8 km) how numerical mixing is substantially reduced and model results become significantly more realistic when vertically adaptive coordinates are applied. We additionally demonstrate that vertically adaptive coordinates perform well in simulating the two dynamically different regions North Sea and Baltic Sea with a single parameter set.An analysis of the computational overhead of the adaptive coordinates indicates an increase of 5–8% in runtime. This is still less expensive than adding more sigma-layers to reduce spurious numerical mixing.     

Parameterization of vertical mixing in the Weddell Sea

A series of vertical mixing schemes implemented in a circumpolar coupled ice–ocean model of the BRIOS family is validated against observations of hydrography and sea ice coverage in the Weddell Sea. Assessed parameterizations include the Richardson number-dependent Pacanowski–Philander scheme, the Mellor–Yamada turbulent closure scheme, the K-profile parameterization, a bulk mixed layer model and the ocean penetrative plume scheme (OPPS). Combinations of the Pacanowski–Philander parameterization or the OPPS with a simple diagnostic model depending on the Monin–Obukhov length yield particularly good results. In contrast, experiments using a constant diffusivity and the traditional convective adjustment cannot reproduce the observations. An underestimation of wind-driven mixing in summer leads to an accumulation of salt in the winter water layer, inducing deep convection in the central Weddell Sea and a homogenization of the water column. Large upward heat fluxes in these simulations lead to the formation of unrealistic, large polynyas in the central Weddell Sea after only a few years of integration. Furthermore, spurious open-ocean convection affects the basin-scale circulation and leads to a significant overestimation of meridional overturning rates. We conclude that an adequate parameterization of both wind-induced mixing and buoyancy-driven convection is crucial for realistic simulations of processes in seasonally ice-covered seas.     

南海环流的三维数值模拟

本文采用Backhaus的三维海流模式,运用半隐式及C-网格方法求解基本方程,对南海各季的平均海流流场进行了数值模拟。将模拟结果与已有的研究结果进行比较,其主要流系基本上是相符的,夏季表层基本上为一反气旋型环流;冬季则转变成气旋型;在冬季,从50m层起开始显露出“南海暖流”的存在。这些表明南海海流的某些主要特征基本上已经被此模型成功地再现出来,同时此模型又给出了垂直方向各层的海流情况,在目前尚缺乏深层实测海流资料的情况下,上述深层海流的模拟结果有一定参考意义。     

Comparison of gravity current mixing parameterizations and calibration using a high-resolution 3D nonhydrostatic spectral element model

In light of the pressing need for development and testing of reliable parameterizations of gravity current entrainment in ocean general circulation models, two existing entrainment parameterization schemes, K-profile parameterization (KPP) and one based on Turner’s work (TP), are compared using idealized experiments of dense water flow over a constant-slope wedge using the HYbrid Coordinate Ocean Model (HYCOM). It is found that the gravity current entrainment resulting from KPP and TP differ significantly from one another. Parameters of KPP and TP are then calibrated using results from the high-order nonhydrostatic spectral element model Nek5000. It is shown that a very good agreement can be reached between the HYCOM simulations with KPP and TP, even though these schemes are quite different from each other.     

Trans-Atlantic responses of Calanus finmarchicus populations to basin-scale forcing associated with the North Atlantic Oscillation

Populations of the copepod species Calanus finmarchicus often dominate the springtime biomass and secondary production of shelf ecosystems throughout the North Atlantic Ocean. Recently, it has been hypothesised that interannual to interdecadal fluctuations observed in such populations are driven primarily by climate-associated changes in ocean circulation. Here, we compare evidence from the North Sea and Gulf of Maine/Western Scotian Shelf (GoM/WSS) linking fluctuations in C. finmarchicus abundance to changes in ocean circulation associated with the North Atlantic Oscillation (NAO). A particularly striking contrast emerges from this Trans-Atlantic comparison: whereas the North Sea C. finmarchicus population exhibits a negative correlation with the NAO index, the GoM/WSS population exhibits a more complex, positive association with the index. The physical processes underlying these contrasting population responses are discussed in the context of regional- to basin-scale circulation changes associated with the NAO.     

A carbon budget in Tokyo Bay

Organic carbon flux from eutrophicated Tokyo Bay to the Pacific Ocean is estimated as 260 ton C day^–1 based on the horizontal gradient of COD and the dispersion coefficient at the bay mouth. Also, carbon flux from the air or from the open ocean to Tokyo Bay is estimated as 156 ton C day^–1. If we suppose that five percent of the coastal seas in the world might be eutrophicated as Tokyo Bay and the organic carbon flux from the shelf to the open ocean in other coastal seas might be one third of that in Tokyo Bay, 1.12 G tons year^–1 would be transported from the eutrophicated coastal seas to the open ocean and such carbon flux may account for the missing sink in the global carbon budget.     

Turbulence modeling in ocean circulation problems

A physical formulation of the problem is considered. A mathematical model and a numerical algorithm of the turbulence model as part of the ocean circulation model for simulations for decades are formulated. The model is based on the evolution equations for turbulent kinetic energy (TKE) and the frequency of its viscous dissipation. A numerical solution algorithm for both the circulation model and the turbulence model is based on implicit schemes of splitting with respect to physical processes and geometric coordinates. For the turbulence model, this provided analytical solutions at a splitting step related to TKE generation and dissipation. Numerical experiments have been performed with a model of the joint circulation of the North Atlantic, the Arctic Ocean, and the Bering Sea to reproduce the annual cycle and synoptic disturbances of ocean characteristics. The model has a resolution of 0.25° in latitude and longitude and 40 levels in the vertical, which are compressed toward the surface to reproduce the process of developed turbulence better. The results are compared with observations and with the results of simulations using traditional parameterizations of the upper ocean mixing. It is shown that the model reproduces ocean characteristics correctly, only slightly increasing the computation time in comparison with simple parameterizations. Spatial and temporal characteristics of turbulence are analyzed.     

Formation and transport of intermediate water masses in a model of the Pacific Ocean

A basin-wide ocean general circulation model of the Pacific Ocean was used to investigate how the interior restoration in the Okhotsk Sea and the isopycnal diffusion affect the circulation and intermediate water masses. Four numerical experiments were conducted, including a run with the same isopycnal and thickness diffusivity of 1.0×103 m2/s, a run employing the interior restoration of temperature and salinity in the Okhotsk Sea with a time scale of 3 months, a run that is the same as the first run except for the enhanced isopycnal mixing, and a final run with the combination of the restoration in the Okhotsk Sea and large isopycnal diffusivity. Simulated results show that the intermediate water masses reproduced in the first run are relatively weak. An increase in isopycnal diffusivity can improve the simulation of both Antarctic and North Pacific intermediate waters, mainly increasing the transport in the interior ocean, but inhibiting the outflow from the Okhotsk Sea. The interior restoration generates the reverse current from the observation in the Okhotsk Sea, whereas the simulation of the temperature and salinity is improved in the high latitude region of the Northern Hemisphere because of the reasonable source of the North Pacific Intermediate Water. A comparison of vertical profiles of temperature and salinity along 50°N between the simulation and observations demonstrates that the vertical mixing in the source region of intermediate water masses is very important.     

Assessment of a NEMO-based downscaling experiment for the North-Western Mediterranean region: Impacts on the Northern Current and comparison with ADCP data and altimetry products

The aim of this paper is to assess the impact of numerical resolution on the simulation of the circulation in the North-Western Mediterranean Sea by comparing two realistic model configurations of 1/64° and 1/16°, the low resolution configuration (LRC) being also used at the high resolution configuration (HRC) open boundaries. The study mainly focuses on the Northern Current (NC), the major feature of the regional circulation in the area, influencing important processes such as meanders giving birth to eddies and NC intrusions on the shelf, strongly conditioning the interactions between coastal and off-shore waters. Both configurations use the NEMO code with the LRC covering the entire Mediterranean basin and the HRC being restricted to the North-Western Mediterranean Sea. The two numerical configurations are assessed by using hull-mounted ADCP data and geostrophic velocities derived from AVISO altimetry at key locations along the NC path. The two configurations develop a different solution regarding the current dynamics but both indicate an important variability of the NC position along its general path. The LRC simulations systematically exhibit a NC shifted too far off the coast and also spuriously displaced seaward along the Gulf of Lions shelf edge while high resolution simulations better fit the remote and in situ observations. In addition, the HRC is able to simulate more realistic coastal features such as a confined coastal jet in good agreement with the ADCP measurements. As satisfactory agreement is reached with the various observations used at different space and time scales, the development of this NEMO high resolution configuration stands as a promising experiment for various process studies or operational oriented applications.     

中国边缘海氮循环过程和源汇格局——以南海为例

边缘海氮循环过程研究是全球海洋氮循环研究的重要组成部分,对全球氮源汇格局有显著影响,进而对全球气候变化产生反馈.人为活动和气候变化又是影响边缘海关键氮循环过程速率的重要因素.南海作为中国和西北太平洋最大的边缘海,是边缘海氮循环研究的理想场所.本文详细总结了南海近岸和海盆区的氮源汇过程及其内循环过程的最新研究进展,结果显...     

Effects of winds,tides and storm surges on ocean surface waves in the Sea of Japan

Ocean surface waves are strongly forced by high wind conditions associated with winter storms in the Sea of Japan. They are also modulated by tides and storm surges. The effects of the variability in surface wind forcing, tides and storm surges on the waves are investigated using a wave model, a high-resolution atmospheric mesoscale model and a hydrodynamic ocean circulation model. Five month-long wave model simulations are inducted to examine the sensitivity of ocean waves to various wind forcing fields, tides and storm surges during January 1997. Compared with observed mean wave parameters, results indicate that the high frequency variability in the surface wind filed has very great effect on wave simulation. Tides and storm surges have a significant impact on the waves in nearshores of the Tsushima-kaihyō, but not for other regions in the Sea of Japan. High spatial and temporal resolution and good quality surface wind products will be crucial for the prediction of surface waves in the JES and other marginal seas, especially near the coastal regions.     

Cross-shelf exchange in a model of the Ross Sea circulation and biogeochemistry

Transport of warm, nutrient-rich Circumpolar Deep Water (CDW) onto Antarctic continental shelves and coastal seas has important effects on physical and biological processes. The present study investigates the locations of this transport and its dynamics in the Ross Sea with a high-resolution three-dimensional numerical model. The model circulation is forced by daily wind stress along with heat and salt fluxes calculated from atmospheric climatologies by bulk formulae. All surface fluxes are modified by an imposed climatological ice cover. Waters under the Ross Ice Shelf are not included explicitly, but their effect on temperature and salinity is imposed in a buffer zone at the southern end of the model domain. A simple nutrient uptake is calculated based on the climatological chlorophyll distribution and Monod uptake kinetics.Model circulation is strongly affected by bottom topography, due to weak stratification, and agrees with schematics of the general flow and long-term current measurements except near the southern boundary. The sea-surface temperature is similar to satellite estimates except that the warmest simulated temperatures are slightly higher than observations. There is a significant correlation between the curvature of the shelf break and the transport across the shelf break. A momentum term balance shows that momentum advection helps to force flow across the shelf break in specific locations due to the curvature of the bathymetry (that is, the isobaths curve in front of the flow). For the model to create a strong intrusion of CDW onto the shelf, it appears two mechanisms are necessary. First, CDW is driven onto the shelf at least partially due to momentum advection and the curvature of the shelf break; then, the general circulation on the shelf takes the CDW into the interior.     

中国近海温跃层判定方法的研究

根据东海以及南海东北部多组资料,探讨了拟阶梯函数法和垂向梯度法在浅海区（水深〈200m）、陆架坡折海域（水深在200m左右）和深水开阔海区（水深〉200m）的应用情况,分析了温跃层在陆架坡折海域的水平衔接问题。对比分析表明：在浅海区,两种方法结果比较一致,特别对某些垂直温度梯度较弱的曲线,拟阶梯函数法和实际水文廓线符合更好。陆架坡折海域,垂向梯度法对于温度梯度相对较弱的跃层刻画不好,而拟阶梯函数法可以对较弱温跃层进行判定,与实际垂直水文结构一致。深水开阔区域,垂向梯度法所得温跃层特征量与实际符合较好,而拟阶梯函数法获得的温跃层强度偏大。因此,对本研究海域来说,当水深≤200m时,采用拟阶梯函数法,当水深〉200m时采用垂向梯度法,最小判定标准统一采用0．05℃／m,可以有效避免陆架坡折海域的水平不连续问题,与实际水文廓线符合较好。     

The connection between Labrador Sea buoyancy loss, deep western boundary current strength, and Gulf Stream path in an ocean circulation model

The sensitivity of the North Atlantic gyre circulation to high latitude buoyancy forcing is explored in a global, non-eddy resolving ocean general circulation model. Increased buoyancy forcing strengthens the deep western boundary current, the northern recirculation gyre, and the North Atlantic Current, which leads to a more realistic Gulf Stream path. High latitude density fluxes and surface water mass transformation are strongly dependent on the choice of sea ice and salinity restoring boundary conditions. Coupling the ocean model to a prognostic sea ice model results in much greater buoyancy loss in the Labrador Sea compared to simulations in which the ocean is forced by prescribed sea ice boundary conditions. A comparison of bulk flux forced hindcast simulations which differ only in their sea ice and salinity restoring forcings reveals the effects of a mixed thermohaline boundary condition transport feedback whereby small, positive temperature and salinity anomalies in subpolar regions are amplified when the gyre spins up as a result of increased buoyancy loss and convection. The primary buoyancy flux effects of the sea ice which cause the simulations to diverge are ice melt, which is less physical in the diagnostic sea ice model, and insulation of the ocean, which is less physical with the prognostic sea ice model. Increased salinity restoring ensures a more realistic net winter buoyancy loss in the Labrador Sea, but it is found that improvements in the Gulf Stream simulation can only be achieved with the excessive buoyancy loss associated with weak salinity restoring.