Mesoscale eddies play vital roles in ocean processes. Although previous studies focused on eddy surface features and individual three-dimensional (3D) eddy cases in the northwestern Pacific Ocean, the analysis of unique eddy 3D regional characteristics is still lacking. A 3D eddy detection scheme is applied to 9 years (2000–2008) of eddy-resolving Regional Ocean Modeling System (ROMS) output to obtain a 3D eddy dataset from the surface to a depth of 1 000 m in the northwestern Pacific Ocean (15°–35°N, 120°–145°E). The 3D characteristics of mesoscale eddies are analyzed in two regions, namely, Box1 (Subtropical Countercurrent, 15°–25°N, 120°–145°E) and Box2 (Southern Kuroshio Extension, 25°–35°N, 120°–145°E). In Box1, the current is characterized by strong vertical shear and weak horizontal shear. In Box2, the current is characterized by the strong Kuroshio, topographic effect, and the westward propagation of Rossby waves. The results indicate the importance of baroclinic instability in Box1, whereas in Box2, both the barotropic and baroclinic instability are important. Moreover, the mesoscale eddies’ properties in Box1 and Box2 are distinct. The eddies in Box1 have larger number and radius but a shorter lifetime. By contrast, Box2 has fewer eddies, which have smaller radius but longer lifetime. Vertically, more eddies are detected at the subsurface than at the surface in both regions; the depth of 650 m is the turning point in Box1. Above this depth, the number of cyclonic eddies (CEs) is larger than that of anticyclonic eddies (AEs). In Box2, the number of CEs is dominant vertically. Eddy kinetic energy (EKE) and mean normalized relative vorticity in Box2 are significantly higher than those in Box1. With increasing depth, the attenuation trend of EKE and relative vorticity of Box1 become greater than those of Box2. Furthermore, the upper ocean (about 300 m in depth) contains 68.6% of the eddies (instantaneous eddy). Only 16.6% of the eddies extend to 1 000 m. In addition, about 87% of the eddies are bowl-shaped eddies in the two regions. Only about 3% are cone-shaped eddies. With increasing depth of the eddies, the proportion of bowl-shaped eddies gradually decreases. Conversely, the cone- and lens-shaped eddies are equal in number at 700–1 000 m, accounting for about 30% each. Studying the 3D characteristics of eddies in two different regions of the northwestern Pacific Ocean is an important stepping stone for discussing the different eddy generation mechanisms.

We introduced the Coupled Model Intercomparison Project Phase 6 (CMIP6) Ocean Model Intercomparison Project CORE2-forced (OMIP-1) experiment by using the First Institute of Oceanography Earth System Model version 2.0 (FIO-ESM v2.0), and comprehensively evaluated the simulation results. Unlike other OMIP models, FIO-ESM v2.0 includes a coupled ocean surface wave component model that takes into account non-breaking surface wave-induced vertical mixing in the ocean and effect of surface wave Stokes drift on air-sea momentum and heat fluxes in the climate system. A sub-layer sea surface temperature (SST) diurnal cycle parameterization was also employed to take into account effect of SST diurnal cycle on air-sea heat fluxes to improve simulations of air-sea interactions. Evaluations show that mean values and long-term trends of significant wave height were adequately reproduced in the FIO-ESM v2.0 OMIP-1 simulations, and there is a reasonable fit between the SST diurnal cycle obtained from in situ observations and that parameterized by FIO-ESM v2.0. Evaluations of model drift, temperature, salinity, mixed layer depth, and the Atlantic Meridional Overturning Circulation show that the model performs well in the FIO-ESM v2.0 OMIP-1 simulation. However, the summer sea ice extent of the Arctic and Antarctic is underestimated.

We report field measurements of vertical profiles of the turbulent diffusivity and temperature at different stations in the South China Sea (SCS). Our study shows that the measured turbulent diffusivity follows a power-law distribution with a varying exponent in water layers. Similar multiple-layer scaling regimes were also observed from the temperature fluctuations. Combining turbulent diffusivity and temperature fluctuations, the vertical structure of temperature was revealed. Furthermore, we discussed the temperature profiles in each layer. A constant function of a dimensionless temperature profile was found in water layers that have identical turbulence conditions. Our results reveal the multiple-layer structure of temperature in the SCS. This study contributes to the understanding of the vertical structure of multiple layers in the SCS and provides clues for exploring the physical mechanism for maintaining the temperature structure.

In this paper, we propose a hybrid forecasting model to improve the forecasting accuracy for depth-averaged current velocities (DACVs) of underwater gliders. The hybrid model is based on a discrete wavelet transform (DWT), a deep belief network (DBN), and a least squares support vector machine (LSSVM). The original DACV series are first decomposed into several high- and one low-frequency subseries by DWT. Then, DBN is used for high-frequency component forecasting, and the LSSVM model is adopted for low-frequency subseries. The effectiveness of the proposed model is verified by two groups of DACV data from sea trials in the South China Sea. Based on four general error criteria, the forecast performance of the proposed model is demonstrated. The comparison models include some well-recognized single models and some related hybrid models. The performance of the proposed model outperformed those of the other methods indicated above.

The velocity and direction of internal waves (IWs) are important parameters of the ocean, however, traditional observation methods can only obtain the average parameters of IWs for a single location or large area. Herein, a new method based on optical flow is proposed to derive the phase velocity vectors of IWs from X-band marine radar images. First, the X-band marine radar image sequence is averaged, and ramp correction is used to reduce the attenuation of gray values with increasing radial range. Second, the average propagation direction of the IWs is determined using the two-dimensional Fourier transform of the radar images; two radial profiles along this direction are selected from two adjacent radar images; and then, the average phase velocity of the IWs is estimated from these radial profiles. Third, the averaged radar images are processed via histogram equalization and binarization to reduce the influence of noise on the radar images. Fourth, a weighting factor is determined using the average phase velocity of a reference point; the phase velocities on the wave crest of the IWs are subsequently estimated via the optical flow method. Finally, the proposed method is validated using X-band marine radar image sequences observed on an oil platform in the South China Sea, and the error of the phase velocity is calculated to be 0.000 3–0.073 8 m/s. The application conditions of the proposed method are also discussed using two different types of IW packets.

Coastal wetlands are located in the ecotone of interaction between the land surface and sea, and anthropogenic activities extensively interfere with these wetlands through the reclamation of large tidal wetlands and destruction of the function of the ecosystems. In this study, we investigated the dynamic evolutionary characteristics of the Bohai Rim coastal area over the past 40 years using the Modified Normalized Difference Water Index, the fractal dimension, object-oriented classification, the land-use transfer trajectory, and regression analysis. Additionally, we quantified and monitored the evolution of reclamation and analyzed the correlation between reclamation and coastal wetlands based on 99 Landsat-2, -5, and -8 images (at 60 m and 30 m spatial resolution) over the period 1980–2019. The results are as follows. (1) The coastline of the Bohai Rim increased by 1 631.2 km from 1980 to 2019 with a zigzag variation. The artificial coastline increased by 2 946.1 km, whereas the natural coastline decreased by 90%. (2) The area of man-made wetlands increased by 3 736.9 km², the area of construction land increased by 1 008.4 km², and the natural wetland area decreased by 66%. The decrease of tidal flats is the main contributor to the decrease of natural wetland area (takes account for 91.1%). Coastal areas are affected by intense human disturbance, which was taken place across a large area of tidal flats and caused the landscape to fragment and be more heterogeneous. The coastal zone development activities were primarily concentrated in the southern Laizhou Bay, the Yellow River Delta, the Bohai Bay, the northern Liaodong Bay, and the Pulandian Bay. The solidified shorelines and increase in sea level have resulted in intertidal wetlands decreasing and impaired wetland ecology. (3) There is a good agreement between reclamation and the size of the coastal wetlands. Both land reclamation and the reduction in coastal wetland areas are significantly related to the population size, fishery output value, and urbanization rate. In summary, human activities, such as the construction of aquaculture ponds and salt pans, industrialization, and urbanization, are the primary forces that influence the environmental changes in the coastal region. This study is beneficial for establishing and improving the systems for the rational development and utilization of natural resources, and provides theoretical references for restoring wetland ecology and managing future reclamation activities in other coastal zone-related areas.

The role of wave breaking (WB) in the ocean dynamics in the Bohai Sea, China under typhoon condition is systematically investigated utilizing a coupled wave-current model. The influences of WB on ocean dynamics and processes (mixing coefficient, temperature, mixed layer depth, and current) during the entire typhoon period (including the pre-typhoon, during-typhoon and after-typhoon stages) are comprehensively detected and discussed. Experimental results show that WB greatly enhances the turbulent mixing at about top 10 m depth under typhoon condition, the increase can be up to 10 times that of the normal weather. At the same time, WB generally strengthens the sea surface cooling by ~1.2°C at the during-typhoon stage, about 3 times that in normal weather. The mixed layer depth, is rapidly increased by ~1.6–3.6 m during typhoon due to WB, particularly, the deepening is stronger in the region from 120.5°E to 121.0°E on account of close to the typhoon eye. In addition, WB renders the current speed more uniformly within the entire depth in the Bohai Sea, the change in speed is ~0.2 m/s, whereas the alternation in current vector is generally opposite to the wind direction except for the typhoon eye region, reflecting that WB has an inhibitory effect on the typhoon-forced current change. The effects of WB on vertical mixing coefficient response to the typhoon rapidly, while the impacts of WB on temperature, and mixed layer depth present hysteretic responses to typhoon. Finally, the mechanisms and distribution characteristics of WB-induced mixing and tidal mixing are compared under typhoon condition.

This study presents an efficient Boussinesq-type wave model accelerated by a single Graphics Processing Unit (GPU). The model uses the hybrid finite volume and finite difference method to solve weakly dispersive and nonlinear Boussinesq equations in the horizontal plane, enabling the model to have the shock-capturing ability to deal with breaking waves and moving shoreline properly. The code is written in CUDA C. To achieve better performance, the model uses cyclic reduction technique to solve massive tridiagonal linear systems and overlapped tiling/shared memory to reduce global memory access and enhance data reuse. Four numerical tests are conducted to validate the GPU implementation. The performance of the GPU model is evaluated by running a series of numerical simulations on two GPU platforms with different hardware configurations. Compared with the CPU version, the maximum speedup ratios for single-precision and double-precision calculations are 55.56 and 32.57, respectively.

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.

Smaller mesoscale eddies (SMEs) have an important effect on the transmission of ocean temperatures, salinity, energy, and marine biochemical processes. However, traditional altimeters, the dominant sensors used to identify and track eddies, have made it challenging to observe SMEs accurately due to resolution limitations. Eddies drive local upwelling or downwelling, leaving signatures on sea surface temperatures (SSTs) and chlorophyll concentrations (Chls). SST can be observed by spaceborne infrared sensors, and Chl can be measured by ocean color remote sensing. Therefore, multisatellite observations provide an opportunity to obtain information to characterize SMEs. In this paper, an eddy detection algorithm based on SST and Chl images is proposed, which identifies eddies by characterizing the spatial and temporal distribution of SST and Chl data. The algorithm is applied to characterize and analyze SMEs in the Kuroshio Extension. Statistical results on their distribution and seasonal variability are shown, and the formation processes are preliminarily discussed. SMEs generation may be contributed by horizontal strain instability, the interaction of topographic obstacles and currents, and wind stress curl.