CABARET in the ocean gyres

A new high-resolution Eulerian numerical method is proposed for modelling quasigeostrophic ocean dynamics in eddying regimes. The method is based on a novel, second-order non-dissipative and low-dispersive conservative advection scheme called CABARET. The properties of the new method are compared with those of several high-resolution Eulerian methods for linear advection and gas dynamics. Then, the CABARET method is applied to the classical model of the double-gyre ocean circulation and its performance is contrasted against that of the common vorticity-preserving Arakawa method. In turbulent regimes, the new method permits credible numerical simulations on much coarser computational grids.     

Effects of adjusting vertical resolution on the eastern equatorial Pacific cold tongue

The vertical resolution of LICOM1.0 (LASG/IAP Climate System Ocean Model, version 1.0) is adjusted by increasing the level amount within the upper 150 m while keeping the total of levels. It is found that the eastern equatorial Pacific cold tongue is sensitive to the adjustment. Compared with the simulation of the original level scheme, the adjusting yields a more realistic strucature of cold tongue extending from the coast of Peru to the equator, as well as a temperature minimum at Costa Rica coast, north of the cold tongue. In the original scheme experiment, the sharp heating by net surface heat flux at the beginning of spin-up leads to a great warming in the eastern equatorial Pacific Ocean. The weak vertical advection due to a too thick mixed layer in the coarse vertical structure also accounts for the warm bias. The fact that most significant improvements of the upper 50 m temperature appear at the region of the thinnest mixed layer indicates the necessity of fine vertical resolution for the eastern equatorial Pacific Ocean. However, the westward extension of equatorial cold tongue, a defect in the original scheme, gets even more serious in the adjusting scheme due to the intensified vertical velocity and hence vertical advection in the central-eastern equatorial Pacific Ocean.     

采用质点跟踪方法对物质输运方程平流项数值格式的改进

用数值模式对河口海岸地区的物质输运进行计算时,平流项的数值格式必须要能对物质浓度锋面进行正确处理,以避免产生过多的数值耗散或频散。本文中设计了一种在网格内设置一些质点并对质点进行跟踪的格式计算平流项。结果表明,质点跟踪格式在一维情形下无频散和几乎没有耗散,在二维情形下无频散和在水深变化剧烈的地方基本避免了垂向数值耗散。与其他数值格式的耗散性和频散性相比,本文中设计的数值格式明显地提高了物质输运方程中平流项的计算精度,在河口海洋物质输运的计算中具有较大的应用价值。     

Second-order moment advection scheme applied to Arctic Ocean simulation

We apply the second-order moment (SOM) advection scheme of (Prather, M.J. 1986. Numerical advection by conservation of second-order moments. J. Geophys. Res. 91, 6671–6681.) to the simulation of the large-scale circulation of the Arctic Ocean with a coupled ocean–sea-ice model. Compared to three other advection schemes commonly employed in ocean simulations (centred differences, flux corrected transport, and multidimensional positive definite advection transport), the SOM method helps preserve the vertical structure of Arctic water masses. The depth, thickness and hydrographic properties of the Arctic Surface Water and the Arctic Atlantic Layer are better represented with SOM than with any of the other three advection algorithms. We also present a convenient method for calculating the implicit numerical diffusivity of upstream based schemes, such as the SOM method, and discuss three approaches for improving the monotonicity properties of the SOM algorithm.     

1980～1999年热带太平洋次表层盐度年际变化同化数据分析

通过统计方法利用一套海洋同化数据分析了热带太平洋次表层的盐度变化特征.结果表明次表层盐度的年际变化与ENSO相关,且次表层盐度信号区域呈东西方向“跷跷板”的分布.对影响这些次表层的盐度信号区域平均的纬向平流、经向平流、垂直运动和淡水通量异常等因素进行了分析,并且与影响表层盐度年际变化模态的影响因素差异进行了比较,结果表明,纬向平流的异常对表层盐度的异常变化影响较大,而对次表层盐度异常有较大影响的是海水的垂直运动异常.     

斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式

依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍     

两种陆面方案对陆雾与海雾模拟效果的对比研究

利用WRF(Weather Research and Forecasting)模式,针对1次陆雾与1次同时包含近海雾区与开阔海域雾区的海雾个例,开展了2种陆面方案——SLAB(five-layer thermal diffusion scheme)方案与Noah(Noah land surface scheme)方案—...     

Improved representation of topographic effects by a vertical adaptive grid in Vector-Ocean-Model (VOM). Part II: Simulations in unstructured adaptive grids

A static adaptive grid approximates the topography and defines the vertical resolution in Vector-Ocean-Model (VOM). The adaptation to topography creates unstructured grids, which are organised in a one-dimensional vector by column-wise storage of only wet cells. The model’s name reflects this data structure. The intention of VOM is better resolving flow and stratification near topographic boundaries in Z-coordinates. This is the second part of a publication that describes the generation of adaptive grids (part I), and simulations with VOM in unstructured grids (this part). Adaptive grids generated for a synthetic topography in part I include shelf, continental slope, and ocean. Three of those grids are here utilised in upwelling simulations. Under the same forcing increased vertical resolution at seabed and slopes yields a significant increase in flow energy as compared to coarser grids. Results allow explaining the surface intensification of a continental slope jet by vertical displacements of water masses in the seabed Ekman layer. Results in unstructured grids are almost identical to reference simulations in equidistant grids where the respective smallest grid size of unstructured grids was used. Negative effects of grids on predicted flow and stratification are absent also over particularly rough topography, as demonstrated by using vertical velocity as most sensitive indicator. In a further simulation an overflow governed by the advection of water mass properties is presented to demonstrate the conservation properties of the model. After 5 months of simulation the predicted domain average temperature deviated by 10⁻⁸ from the initial temperature field. Compared to equidistant grids the advection/diffusion scheme looses about one order of magnitude in accuracy when used in an unstructured grid. The results of VOM, being defined in Z-coordinates, are void of coordinate transformation errors. In an arbitrary topography unforced zero-flow remains quiescent in a stratification that only varies in the vertical. VOM due to its depth-independent vertical resolution appears particularly suitable for simulations of ocean-shelf exchange.     

长江口浑浊带核心区北槽水动力特征研究

本文根据2015年长江口洪季浑浊带水域座底三脚架全水深范围观测数据,对北槽中下段纵横向流、余流结构、边界层特征等进行了深入研究,结果显示:（1）大潮潮流矢量特征在垂向上的变化差异最为显著;北槽中下段水流出现明显的横向输移,指向航槽北面;（2）横向流速在大潮期间的变化（垂线平均）为-0.10~0.39 m/s,小潮期间的变化为-0.13~0.25 m/s;横向平流对于北槽纵向物质输运有显著调整作用;（3）余流垂向差异大,底部几乎为0,小潮期间余流更强;周期性再悬浮的细颗粒泥沙可作为北槽最大浑浊带的重要物质来源;（4）传统六点对数流速拟合法获得摩阻流速误差较大,近底1 m区域高分辨流速剖面是准确获取底边界层参数的有效途径。     

The application of a Godunov-type shock capturing scheme for the simulation of waves from deep water up to the swash zone

A two-dimensional vertical (2DV) non-hydrostatic boundary fitted model based on a Godunov-type shock-capturing scheme is introduced and applied to the simulation of waves from deep water up to the swash zone. The effects of shoaling, breaking, surf zone dissipation and swash motions are considered. The application of a Godunov-type shock-capturing algorithm together with an implicit solver on a standard staggered grid is proposed as a new approach in the 2DV simulation of large gradient problems such as wave breaking and hydraulic jumps. The complete form of conservative Reynolds averaged Navier–Stokes (RANS) equations are solved using an implicit finite volume method with a pressure correction technique. The horizontal advection of the horizontal velocity is solved by an explicit predictor–corrector method. Fluxes are predicted by an exact Riemann solver and corrected by a downwind scheme. A simple total variation diminishing (TVD) method with a monotonic upstream-centered scheme for conservation laws (MUSCL) limiter function is employed to eliminate undesirable oscillations across discontinuities. Validation of the model is carried out by comparing the results of the simulations with several experimental test cases of wave breaking and run-up and the analytical solution to linear short waves in deep water. Promising performance of the model has been observed.     

Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.     

Effects in a climate model of slope tapering in neutral physics schemes

In many global ocean climate models, mesoscale eddies are parameterized as along isopycnal diffusion and eddy-induced advection (or equivalently skew-diffusion). The eddy-induced advection flattens isopycnals and acts as a sink of available potential energy, whereas the isopycnal diffusion mixes tracers along neutral directions. While much effort has gone into estimating diffusivities associated with this closure, less attention has been paid to the details of how this closure (which tries to flatten isopycnals) interacts with the mixed layer (in which vertical mixing tries to drive the isopycnals vertical). In order to maintain numerical stability, models often stipulate a maximum slope S_max which in combination with the thickness diffusivity A_gm defines a maximum eddy-induced advective transport A_gm1S_max. In this paper, we examine the impact of changing S_max within the GFDL global coupled climate model. We show that this parameter produces significant changes in wintertime mixed layer depth, with implications for wintertime temperatures in key regions, the distribution of precipitation, and the vertical structure of heat uptake. Smaller changes are seen in details of ventilation and currents, and even smaller changes as regards the overall hydrography. The results suggest that not only the value of the coefficient, but the details of the tapering scheme, need to be considered when comparing isopycnal mixing schemes in models.     

黄海、渤海盐度准三维数值预报模式

黄海、渤海盐度的垂直结构具有典型的自模性，而其水平分布又受平流、水平扩散效应及径流等因素的影响。本文根据黄海、渤海实测资料拟合了盐度垂直剖面的自模函数，并结合描述表、底层盐度及上均匀层厚度这３个特征量水平分布的方程，给出盐度三维结构的准三维模式。在模式中，综合考虑了海面风和热输入的强迫作用以及流场的平流、侧向混合及底层混合的影响因素，同时还考虑了径流、蒸发及降水的作用，较客观地反映了盐度的三维分布及其变化的物理过程。试报结果分析表明，模式的功能较好，结果令人满意。     

北太平洋副热带模态水形成区混合层热动力过程诊断分析

利用NCEP海洋数据和COADS海气通量资料,通过诊断分析,揭示了海表热力强迫、垂直夹卷、埃克曼平流和地转平流效应在北太平洋副热带模态水形成过程中的贡献。研究表明,在北太平洋副热带3个模态水形成海域冬季混合层降温过程中,海表热力强迫和垂直夹卷效应是主导因素,二者的相对贡献分别约为67%和19%(西部模态水)、53%和21%(中部模态水)、65%和30%(东部模态水);并且在东部模态水形成海域,埃克曼平流和地转平流皆是暖平流效应,而在西部和中部模态水形成海域,仅有地转平流是暖平流效应。进一步的分析表明,海洋平流(地转平流、埃克曼平流)对北太平洋副热带模态水形成海域秋、冬季混合层温度的年际、年代际异常有显著影响,在西部模态水形成海域,海表热力强迫(62%)和地转平流(32%)是导致混合层温度年际、年代际变化的主要因子;在中部模态水形成海域,混合层温度的年际、年代际变化是埃克曼平流(32%)、地转平流(30%)和海表热力强迫(25%)共同作用的结果;相对而言,东部模态水形成海域混合层温度的年际、年代际异常主要受海表热力强迫(67%)控制。     

A generalized coordinate ocean model and a comparison of the bottom boundary layer dynamics in terrain-following and in z-level grids

Sensitivity studies with a new generalized coordinate ocean model are performed in order to compare the behavior of bottom boundary layers (BBLs) when terrain-following (sigma or combined sigma and z-level) or z-level vertical grids are used, but most other numerical aspects remain unchanged. The model uses a second-order turbulence closure scheme that provides surface and BBL mixing and results in a quite realistic climatology and deep water masses after 100 year simulations with a coarse resolution (1° × 1°) basin-scale terrain-following grid. However, with the same turbulence scheme but using a z-level grid, the model was unable to produce dense water masses in the deep ocean. The latter is a known problem for coarse resolution z-level models, unless they include highly empirical BBL schemes.A set of dense water overflow experiments with high-resolution grids (10 and 2.5 km) are used to investigate the influence of model parameters such as horizontal diffusivity, vertical mixing, horizontal resolution, and vertical resolution on the simulation of bottom layers for the different coordinate systems. Increasing horizontal diffusivity causes a thinner BBL and a bottom plume that extends further downslope in a sigma grid, but causes a thicker BBL and limited downslope plume extension in a z-level grid. A major difference in the behavior of the BBL in the two grids is due to the larger vertical mixing generated by the turbulence scheme over the step-like topography in the z-level grid, compared to a smaller vertical mixing and a more stably stratified BBL in the sigma grid. Therefore, the dense plume is able to maintain its water mass better and penetrates farther downslope in the sigma grid than in the z-level grid. Increasing horizontal and vertical resolution in the z-level grid converges the results toward those obtained by a much coarser resolution sigma coordinate grid, but some differences remain due to the basic differences in the mixing process in the BBL.     

北太平洋海表温度及各贡献因子的变化

采用1958年1月至2007年12月SODA海洋上层温度的月平均资料,基于海温变化方程和统计分析方法,分析了北太平洋海表面温度(SST)异常特征及各局地因子贡献比例的变化。结果表明,伴随着1976/1977风场最强中心位置的南北移动,形成了两个北太平洋SST年际-年代际变化的异常中心:一个是位于30°N附近的副热带海盆内区,SST异常主要受风应力强度的主导;一个是位于40°N附近的副热带和副极地环流交汇区,SST异常主要受风应力旋度的位置即风场位置的影响。在副热带海盆内区,最强降温发生在1978-1982年,SST异常的主要局地贡献因子为海表热通量和经向平流,二者所占比例和约为50%~60%,均为同相增温或降温作用,余项所占比例约为20%~50%。在副热带和副极地环流交汇区,海盆内区和西部边界区的SST异常的跃变时间同为1975年,但是内区的垂直混合项的跃变时间早于西部5年左右。SST异常的主要贡献因子为海表热通量和经向平流,但在1983-1988年海温强降温期间,经向平流项贡献大于海表热通量项的贡献。两个区域的垂直混合项均为降温贡献,虽然量值小却显示出很强的年代际变化信号。平流项中经向平流最大,垂直平流最小。     

中国青岛市近年来灰霾污染数值模拟研究

基于SO₂和PM_2.5的小时观测数据,本研究发现青岛市2014—2018年采暖期间SO2和PM_2.5平均浓度呈现下降的趋势,但是高浓度PM_2.5在不同年份偏移明显,尤其是2015年高浓度PM_2.5显著高于其它年份。为阐明这一时期高浓度PM_2.5的物理化学过程,利用区域气象化学模式Weather Research and Forecasting-Community Multiscale Air Quality Model(WRF-CMAQ)对2015年12月~2016年1月的PM_2.5累积和清除的物理化学过程分析,发现浓度累积主要受源排放、气溶胶过程和垂直对流过程影响,清除的主要控制过程为垂向扩散和水平对流。同时,对2015年12月29~31日和2016年1月15~17日两次高浓度PM2.5污染事件的分析发现,PM_2.5的主要成分为二次无机离子(SNA)和有机气溶胶(OA),其中SNA占比44%,OA占近30%。SNA主要由气溶胶反应过程、垂直和水平对流过程控制,OA主要受垂直和水平对流过程影响。本研究揭示的影响PM2.5组分的不同物理化学过程,一方面利于污染减排政策更有针对性的制定;另一方面从模式开发的角度,改进PM_2.5高比例组分,如SNA和OA的模拟能力,对于PM_2.5生消的机制解析具有重要意义。     

The Influence of Physical Factors on the Variation of Phytoplankton and Nutrients in the Bohai Sea

The cycle of the phytoplankton in a coastal water is controlled by the biological processes, solar radiation, water temperature and physical transport processes. A three-dimensional ecosystem dynamic model is adopted in this study to investigate the influence of different physical factors on the variation of phytoplankton and nutrients in the Bohai Sea. The simulation is carried out for the year 1982. The simulated annual cycle of the primary production and nutrients are in reasonable agreement with the observations in the pattern. Vertical mixing can both affect the vertical transportation of nutrients and horizontal distribution of primary production. In winter the vertical distribution of nutrients is homogeneous because of the intensive mixing, while in summer there is a high value of nutrients in the depth about 15 m due to the stratification. The high primary production plague and the weak mixing center is positional correspondence. The production of phytoplankton is sensitive to the photosynthetically active radiation, which is strongly influenced by the transparency. The increase of the transparency can promote the production in spring and autumn significantly, but has little effect on the production in summer. The change of the transparency can both affect the occurrence time and the amplitude of the phytoplankton bloom dramatically. Horizontal advection does not affect the variation trend of the annual cycle of chlorophyll-a, but does affect the relative magnitude of the phytoplankton bloom, especially in summer. Horizontal advection can dramatically alter the horizontal distribution of chlorophyll-a. The maximum concentration of chlorophyll-a without horizontal advection in summer is twice as high than that with advection and the high chlorophyll-a areas locate along the coast. The river discharge only has regional influence on the ecosystem. The Huanghe River with high nitrate concentration influ-ences the annual cycle of nitrogen of the Laizhou Bay significantly.     

Acoustic and optical variations during rapid downward motion episodes in the deep north-western Mediterranean Sea

An Acoustic Doppler Current Profiler (ADCP) was moored at the deep-sea site of the ANTARES neutrino telescope near Toulon, France, thus providing a unique opportunity to compare high-resolution acoustic and optical observations between 70 and 170 m above the sea bed at 2475 m. The ADCP measured downward vertical currents of magnitudes up to 0.03 m s⁻¹ in late winter and early spring 2006. In the same period, observations were made of enhanced levels of acoustic reflection, interpreted as suspended particles including zooplankton, by a factor of about 10 and of horizontal currents reaching 0.35 m s⁻¹. These observations coincided with high light levels detected by the telescope, interpreted as increased bioluminescence. During winter 2006 deep dense-water formation occurred in the Ligurian subbasin, thus providing a possible explanation for these observations. However, the 10-20 days quasi-periodic episodes of high levels of acoustic reflection, light and large vertical currents continuing into the summer are not direct evidence of this process. It is hypothesized that the main process allowing for suspended material to be moved vertically later in the year is local advection, linked with topographic boundary current instabilities along the rim of the ‘Northern Current’.     

Finite differencing of the advection term

The advection term in the finite difference diffusion equation can give rise to the false computational oscillation unless the ratio of the grid size to the coefficient of eddy diffusivity is smaller than a critical value. The slightly uncentered differencing scheme is useful for the advection term.