On the necessity of creating a series of oceanographic climatic research atlases (SOCRAT)

In this paper, an integration method is proposed for a system of ocean dynamics equations (3DPEM) for creating a series of atlases with the joint title SOCRAT so that circulation in lower ocean layers can be adjusted to circulation in upper layers, avoiding smoothing of the intensity in upper layers caused by long integration or poor resolution. Specific simulations are presented, which were carried out by the author with a group of colleagues. These results conserve the realistic features of circulation thanks to the short integration period and high resolution. It is shown that due to long integration and/or simulation with rough resolution, simulations lose the main feature of ocean dynamics: intense alongshore currents; as a result, a nonexistant stable physical state uncharacteristic of an ocean is imposed.     

A global barotropic ocean model driven by synoptic atmospheric disturbances for detecting seafloor vertical displacements from in situ ocean bottom pressure measurements

A global barotropic ocean model forced by atmospheric disturbances is developed for the detection of seafloor vertical displacements from in situ ocean bottom pressure (OBP) data. The model accuracy is validated by deep-sea OBP data at more than 100 sites obtained over the global ocean. Parameters and boundary conditions including the horizontal resolution incorporated in the ocean model are tested in order to accurately simulate the nontidal (>2 days) OBP variations. The horizontal resolution is found to the factor that most significantly affects the simulated result. The finer the horizontal resolution applied, the smaller the model variability is. The model accuracy is highest when the horizontal resolution is 1/12°, but deteriorates when the horizontal resolution is finer than 1/12°. This may indicate a failure of the energy dissipation parameterization in the barotropic ocean model. Using the developed 1/12° model, the root-mean-square of the observed nontidal OBP component can be reduced by 18 % as an average of all the OBP data used. It is found that the 1/12° model is useful for the detection of a slow seafloor vertical displacement of centimeters related to the 2011 Tohoku-Oki earthquake from in situ OBP records near the hypocenter of the earthquake.     

Marine geodesy— a direction for the future

The sensitivity of the ocean mixed layer response to different parameters such as model horizontal resolution, vertical temperature gradient and eye size is investigated in response to moving Indian Ocean cyclone. For this, a one and one-half layer wind-driven reduced-gravity ocean model is forced with synthetic cyclonic vortex. The sensitivity studies are carried out for a cyclone moving along idealized tracks, initially and then for three observed cyclones (TC 01A), (TC02B) and (TC04A) during 2004. Increasing model resolution resulted in stronger ocean response. The role of initial vertical temperature gradient in modulating the ocean response is found to be important.     

A comparison of two global ocean-ice coupled models with different horizontal resolutions

A global eddy-permitting ocean-ice coupled model with a horizontal resolution of 0.25 by 0.25 is established on the basis of Modular Ocean Model version 4 (MOM4) and Sea Ice Simulator (SIS). Simulation results are compared with those of an intermediate resolution ocean-ice coupled model with a horizontal resolution of about 1 by 1 . The results show that the simulated ocean temperature, ocean current and sea ice concentration from the eddy-permitting model are better than those from the intermediate resolution model. However, both the two models have the common problem of ocean general circulation models (OGCMs) that the majority of the simulated summer sea surface temperature (SST) is too warm while the majority of the simulated subsurface summer temperature is too cold. Further numerical experiments show that this problem can be alleviated by incorporating the non-breaking surface wave-induced vertical mixing into the vertical mixing scheme for both eddy-permitting and intermediate resolution models.     

北极海域海面风场和海浪遥感观测能力分析

卫星遥感是开展北极海域海面风场和海浪分布特征与变化规律研究的重要手段。本文基于在轨多源卫星遥感数据,从遥感观测空间覆盖、时间覆盖和多源卫星遥感数据融合等方面开展北极海域海面风场与海浪遥感观测能力分析,研究主要结果为:基于ASCAT和HY-2A散射计可实现北极海域海面风场遥感观测,通过多星联合观测可获取北极海域时空分辨率优于12 h和0.1°的海面风场遥感融合数据;基于HY-2A、CryoSat-2、SARAL和Sentinel-3高度计可实现北极海域海浪遥感观测,同样通过多星联合观测可获取北极海域时空分辨率优于1 d和0.25°的海浪有效波高遥感融合数据;基于2016年北极海面风场和海浪遥感融合数据,分析得出北极海域海面风场和海浪在2月处于极大值,然后逐渐减小,7月最小,随后开始逐渐增大。本研究表明,基于多源散射计和高度计遥感观测可实现北极海域海面风场和海浪的高时空分辨率遥感业务化监测。     

Simulation of the spatiotemporal variability of the World Ocean sea surface hight by the INM climate models

The results of simulations of the World Ocean sea surface hight (SSH) in by various versions of the Climate Model of the Institute of Numerical Mathematics, Russian Academy of Sciences, are compared with the CNES-CLS09 fields of the mean dynamic topography (deviation of the ocean level from the geoid). Three models with different ocean blocks are considered which slightly differ in numerical schemes and have various horizontal spatial resolution, i.e., the INMCM4 model, which participated in the Climate Model Intercomparison Project (CMIP Phase 5, resolution of 1° × 1/2°); the INMCM5 model, which participates in the next project, CMIP6 (resolution of 1/2° × 1/4°); and the advanced INMCM-ER eddy-resolving model (resolution of 1/6° × 1/8°). It is shown that an increase in the spatial resolution improves the reproduction of ocean currents (with Agulhas and Kuroshio currents as examples) and their variability. A probable cause of relatively high errors in the reproduction of the SSH of Southern and Indian oceans is discussed.     

HF radar ocean current algorithm based on MUSIC and the validation experiments

Wuhan University's ocean state measuring and analyzing radar (OSMAR2000), working at around 7.5 MHz in the low region of the HF band with a 120-m-long linear receiving antenna array, can measure ocean surface current at ranges of up to 200 km. An ocean surface current algorithm based on direction finding (DF) using the multiple signal classification (MUSIC) method is developed for the OSMAR2000 radar. This paper describes the OSMAR2000 ocean surface current algorithm based on MUSIC and the validation experiments in the East China Sea. The results of the ocean surface current measurements demonstrate that the OSMAR2000 ocean surface current algorithm based on MUSIC is feasible for the long range of ocean surface current mapping with a sufficient bearing resolution.     

Improved representation of topographic effects by a vertical adaptive grid in vector-ocean-model (VOM). Part I: Generation of adaptive grids

This is the first part of a publication that describes the generation of adaptive grids (this part), and simulations with vector-ocean-model (VOM) in unstructured grids resulting from the adaptation (part II). A static vertical adaptive grid in z-coordinates allows improving the approximation of topography and vertical resolution at slopes. Adaptive grids use elements from a set of grid sizes by multiplying a basic smallest cell size with powers of two, as in cell division. Grids with locally isotropic vertical resolution at surface, seabed, and slopes can be generated whereby resolution decreases towards the ocean interior. The adaptation to topography yields unstructured grids that are organised in a one-dimensional vector by column-wise storage of cells, discarding land cells. The vector storage suggested the model’s name. Grids are generated by an iterative procedure that relies on rules, i.e. criteria and directives to control the grid structure in favour of a good representation of physics and smooth numerical operations. The directives govern vertical resolution at sea surface and seabed, and at slopes. For the latter vertical resolution is extended in the horizontal. In the ocean interior horizontal distances between changes in grid size can be controlled for the sake of smooth numerics. The use of a z-grid that avoids transformation errors, the depth-independence of vertical resolution, and the lateral extension of vertical resolution at slopes towards the ocean interior are the most significant differences of adaptive grids in comparison to vertical coordinate transformations. Unstructured grids do not rely on a smoothing of topography and can be used within any of the horizontal Arakawa-grids. For the same topography directives allow creating various grids as demonstrated for a shelf-ocean topography. The number of cells per column in two unstructured grids generated for the North Atlantic may locally well exceed typical layer numbers in conventional model matrices. But the domain average is similar to layer numbers of today’s ocean models. Thus, with the same investment of cells per domain a higher resolution in slope regions can be achieved by unstructured grids as compared to conventional z-grids.     

Representation of topography by porous barriers and objective interpolation of topographic data

We present a porous medium approach to representing topography, and a new algorithm for the objective interpolation of topography, for use in ocean circulation models of fixed resolution. The representation and algorithm makes use of two concepts; impermeable thin walls and porous barriers. Impermeable thin walls allow the representation of knife-edge sub-grid-scale barriers that block lateral flow between model grid cells. Porous barriers permit the sub-grid scale geometry to modulate lateral transport as a function of elevation. We find that the porous representation and the resulting interpolated topography retains key features, such as overflow sill depths, without compromising other dynamically relevant aspects, such as mean ocean depth for a cell. The accurate representation of the ocean depth is illustrated in a simple model of a tsunami that has a cross-basin travel time very much less dependent on horizontal resolution than when using conventional topographic interpolation and representation.     

平潭近岸海域浮子漂移轨迹及其数值模拟

平潭近岸海域岸线曲折,周边岛礁众多,海底地形复杂,是福建省海难事故的高发区.本研究建立了平潭近岸海域海上目标物漂移轨迹的预测系统,该系统通过风场和流场的数值模型获取海面动力环境信息,采用拉格朗日算法实现对海上目标物漂移轨迹的预测追踪.其中海流模型采用ROMS（regional ocean modeling system）模型构建,模型水平方向上最高分辨率为100m,垂向上分为16层,并考虑干湿边界,以体现复杂海岸线和水深地形.通过验证分析,潮位、流速和流向的模拟平均绝对误差分别为0.20 m、0.12 m/s和26°.通过平潭近岸2个浮子实验,结果表明,浮子漂移过程中受潮流和局地地形的影响明显,对漂移模型在平潭近岸海域的适用性进行初步验证,浮子模拟轨迹与实际漂移过程基本吻合,模拟时段内最大偏差距离为2.8km,系统可以为平潭近岸海域海上突发事故应急决策提供参考.     

海洋水文气象多时空尺度资料来源分析

目前,海洋数据获取量日益增多,根据时空分辨率、时空范围、垂向分层情况等的不同,可将海洋数据分为多种类型。了解和掌握现有的海洋数据资料,并加以利用,对国家经济发展、国防安全、军事活动等具有重要意义。经过长时间的研究和积累,本文总结了国内外海洋数据资料发布的几个主要数据源,同时归纳了几类常用的海洋数据资料,包括地形水深数据、卫星遥感数据、海洋调查观测数据、再分析产品数据,并对各类数据适用的时间尺度、空间尺度进行分析,旨在为海洋领域的科研工作者提供信息支持,为我国海洋大数据的建设提供参考。     

Parametric estimation of ocean surface currents with HF radar

Ocean surface currents can be estimated, over a large coastal area, by utilizing the backscatter of high frequency (HF) radar waves from ocean gravity waves. Although the overall backscatter mechanism is complicated, the surface current information is contained within the spectral characteristics of two dominant Bragg components. The accuracy of the current estimate, following the usual FFT-based spectral estimate, is limited by the frequency resolution of the FFT and the time-varying characteristics of the Bragg components. This paper describes a high resolution parametric estimation of the ocean currents based on a recently proposed technique for analyzing time-varying signals. This technique, together with a time-domain ocean clutter model, allows all the Bragg signal information to be extracted from the two dominant eigenvalues and eigenvectors of a matrix constructed from the radar data. Using signals from an operational coastal surveillance radar, current estimates made using this technique are compared with those estimated by the conventional FFT-based method     

Validating and assessing the sensitivity of the climate model with an ocean general circulation model developed at the Institute of Atmospheric Physics,Russian Academy of Sciences

A new version of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), climate model (CM) has been developed using an ocean general circulation model instead of the statistical-dynamical ocean model applied in the previous version. The spatial resolution of the new ocean model is 3° in latitude and 5° in longitude, with 25 unevenly spaced vertical levels. In the previous version of the oceanic model, as in the atmospheric model, the horizontal resolution was 4.5° in latitude and 6° in longitude, with four vertical levels (the upper quasi-homogeneous layer, seasonal thermocline, abyssal ocean, and bottom friction layer). There is no correction for the heat and momentum fluxes between the atmosphere and ocean in the new version of the IAP RAS CM. Numerical experiments with the IAP RAS CM have been performed under current initial and boundary conditions, as well as with an increasing concentration of atmospheric carbon dioxide. The main simulated atmospheric and oceanic fields agree quite well with observational data. The new version’s equilibrium temperature sensitivity to atmospheric CO₂ doubling was found to be 2.9 K. This value lies in the mid-range of estimates (2–4.5 K) obtained from simulations with state-of-the-art models of different complexities.     

利用志愿观测船提升表层海洋观测能力

对表层海洋动力学及其调控机制的认识对准确认识海洋自然变化过程及其未来对气候变化的响应具有重要意义。志愿观测船作为自动观测系统的载体和样品采集平台的优势使其在提升人类对于表层海洋的观测能力方面具有极大优势。本文系统性地回顾了志愿观测船的概念、发展历史、现状和未来发展趋势。基于志愿观测船的诸多海洋观测项目的成功显示了,随着快速发展的传感器技术,志愿观测船能够以经济有效的方式显著地提升表层海洋现场观测的时间和空间分辨率。未来的海洋观测系统应该综合利用包括志愿观测船在内的多种海洋观测平台的优势。     

Synthetic aperture radar imaging of ocean waves during the marineland experiment

X- andL-band simultaneously obtained synthetic aperture radar (SAR) data of ocean gravity waves collected during the Marineland Experiment were analyzed using wave contrast measurements. The Marineland data collected in 1975 represents a unique historical data set for testing still-evolving theoretical models of the SAR ocean wave imaging process. The wave contrast measurements referred to are direct measurements of the backscatter variation between wave crests and troughs. These modulation depth measurements, which are indicators of wave detectability, were made as a function of: a) the settings used in processing the SAR signal histories to partially account for wave motion; b) wave propagation direction with respect to radar look direction for bothX- andL-band SAR data; c) SAR resolution; and d) number of coherent looks. The contrast measurements indicated that ocean waves imaged by a SAR are most discernible whenX-band frequency is used (as compared toL-band), and when the ocean waves are traveling in the range direction. Ocean waves can be detected by bothX- andL-band SAR, provided that the radar surface resolution is small compared to the ocean wavelength (at least 1/4 of the ocean wavelength is indicated by this work). Finally, wave detection withL-band SAR can be improved by adjusting the focal distance and rotation of the cylindrical telescope in the SAR optical processor to account for wave motion. The latter adjustments are found to be proportional to a value that is near the wave phase velocity.     

全海洋浅地层剖面仪及其应用

全海洋浅地层剖面仪(Topographic Parametric Sonar,TOPAS)PS 018系统是目前世界上最先进的浅地层剖面仪之一。该系统是全海洋宽带非线性差频浅地层剖面仪,可对海底地层进行全方位测量,同时还兼有测量水深的功能,最大地层穿透深度为150 m,最小分辨率为0.3 m。系统多种发射信号(Ricker波、Burst波和Chirp波)的选取方便了操作者使用,从理论上实现了全海洋测量功能。从实测剖面分析,该系统是中、深水地层测量的理想测量系统。     

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.     

High-resolution sea ice in long-term global ocean GCM integrations

The resolution of the sea-ice component of a coarse-resolution global ocean general circulation model (GCM) has been enhanced to about 22 km in the Southern Ocean. The ocean GCM is designed for long-term integrations suitable for investigations of the deep-ocean equilibrium response to changes in southern hemisphere high-latitude processes. The space and time scales of the high-resolution sea-ice component are commensurate with those of the resolution of satellite passive-microwave sea-ice data. This provides the opportunity for a rigorous evaluation of simulated sea-ice characteristics. It is found that the satellite-derived continuous high ice concentration of the interior winter ice pack can only be captured when vertical oceanic mixing is modified in a way that less local, intermittent convection occurs. Furthermore, the width and the variability of the coastal polynyas around the Antarctic continent and its ice shelves are best captured when some form of ice-shelf melting is accounted for. The width of the wintertime ice edge is reasonably reproduced, while its variability remains underestimated, closely following the coarse-grid pattern of the ocean model due to its high dependence on ocean temperature. Additional variability besides daily winds, e.g. in form of idealized tidal currents, improves the temporal and spatial ice-edge variability, while leads in the interior ice pack become more abundant, more in line with the fine-scale satellite-derived texture. The coast- or ice-shelf line is described on the fine grid based on satellite passive-microwave data. This method requires parts of a coarse coastal ocean grid cell to be covered by an inert layer of “fast ice” or “ice shelf”. Reasonable long-term global deep-ocean properties can only be achieved when these areas are not inert, i.e. are exposed to heat flux and ice growth, or when the vertical mixing parameterization allows for excessive open-ocean convection. The model area exposed to cold high-latitude atmospheric conditions thus being most decisive for a realistic representation of the long-term deep-ocean properties, suggests that high-latitude coastlines are definitely in need of being represented at high resolution, including ice sheets and their effects on the heat and freshwater flux for the ocean.     

HY-1C卫星CZI载荷的黄海绿潮提取研究

海洋一号C（HY-1C）卫星是中国首颗海洋水色业务卫星,其搭载的海岸带成像仪（CZI）在近海海洋环境监测中正发挥越来越重要的作用;随着搭载有相同传感器的HY-1D卫星发射,双星组网观测,可形成3天2次的高频次、大范围对海观测能力,在海洋漂浮藻类、海洋溢油等目标探测方面具备优异的效能。高空间分辨率光学数据中包含了丰富的海洋环境信息,给特定目标的识别提取带来一定干扰。本研究面向HY-1C卫星CZI载荷开展中国近海漂浮藻类识别提取的业务化应用需求,发展基于藻类缩放指数与虚拟基线高度融合的海洋漂浮藻类识别提取算法,算法优选适用于无短波红外波段国产数据的虚拟基线高度指数来增强藻类信号,通过藻类缩放指数滑动窗口运算,有效剔除高空间分辨率光学数据中的复杂干扰信息,实现了基于CZI数据的海洋漂浮藻类高精度提取,且具有较好的计算运行效率。此外,结合准同步高分卫星16 m多光谱数据,开展CZI数据含藻像元的不确定性分析,发现CZI数据反演结果对近海小斑块漂浮藻类存在不可忽视的高估现象。研究还指出,光学数据用于漂浮藻类监测,其不确定性不仅来源于传感器的空间分辨率差异,还与海洋漂浮藻类形态特征的空间分异性有关。明晰海洋漂浮藻类的形态学空间分异特征,将有助于提高光学数据反演结果的精度,并阐明不确定性。     

浅层物探技术在近海灾害地质与工程地质调查中的应用

随着海洋资源的开发与滨岸工程的建设,越来越多地遇到海洋地质灾害的影响与危害。为了避免并降低海洋地质灾害的影响程度,应开展海洋灾害地质与工程地质调查。二十多年来的实践经验表明,若采用以高分辨率地震为主,配合浅地层剖面、旁侧声纳、测深及磁力测量等综合浅层物探技术,对于了解海洋地质灾害因素与海底浅部工程地质特性等方面都能取得良好的效果,这是一种既经济又较快速地提供成果资料的好方法。