Recent studies have revealed that the predominant tidal constituents have seasonal variations at some locations. However, how to accurately obtain these variations remains a problem for the traditional harmonic analysis (HA) due to the tradeoff between length of time window and resolution of constituents. Therefore, a method named as “two-step HA” is developed in this study, which consists of both long- and short-time-window HA. Through a series of ideal experiments, practical application at two tidal gauges and comparison with the traditional HA, the feasibility and accuracy of the two-step HA are verified: The two-step HA performs better than the traditional HA in estimating monthly amplitudes and phases for the predominant constituents, whether they have seasonal variability or not. In addition to capturing variations of the predominant constituents at tidal gauges, the two-step HA would be useful in investigation of the coherence and incoherence of internal tides.

Between June 2015 and June 2017, two pressure-recording inverted echo sounders (PIESs) and five current and pressure-recording inverted echo sounders (CPIESs) deployed along a section across the Kerama Gap acquired a dataset of ocean bottom pressure records in which there was significant 21-day variability (P_bot21). The P_bot21, which was particularly strong from July–December 2016, was coherent with wind stress curl (WSC) on the continental shelf of the East China Sea (ECS) with a squared coherence of 0.65 for a 3-day time lag. A barotropic ocean model demonstrated the generation, propagation, and dissipation of P_bot21. The modeled results show that the P_bot21 driven by coastal ocean WSC in the ECS propagated toward the Ryukyu Island Chain (RIC), while deep ocean WSC could not induce such variability. On the continental shelf, the P_bot21 was generated nearly synchronously with the WSC from the coastline to the southeast but dissipated within a few days due to the effect of bottom friction. The detection of P_bot21 by the moored array was dependent on the 21-day WSC patterns on the continental shelf. The P_bot21 driven southeast of the Changjiang Estuary by the WSC was detected while the P_bot21 generated northeast of the Changjiang Estuary was not.

Many typhoons pass through the East China Sea (ECS) and the oceanic responses to typhoons on the ECS shelf are very energetic. However, these responses are not well studied because of the complicated background oceanic environment. The sea surface temperature (SST) response to a severe Typhoon Rananim in August 2004 on the ECS shelf was observed by the merged cloud-penetrating microwave and infrared SST data. The observed SST response shows an extensive SST cooling with a maximum cooling of 3°C on the ECS shelf and the SST cooling lags the typhoon by about one day. A numerical model is designed to simulate the oceanic responses to Rananim. The numerical model reasonably simulates the observed SST response and thereby provides a more comprehensive investigation on the oceanic temperature and current responses. The simulation shows that Rananim deepens the ocean mix layer by more than 10 m on the ECS shelf and causes a cooling in the whole mixed layer. Both upwelling and entrainment are responsible for the cooling. Rananim significantly deforms the background Taiwan Warm Current on the ECS shelf and generates strong Ekman current at the surface. After the typhoon disappears, the surface current rotates clockwise and vertically, the current is featured by near inertial oscillation with upward propagating phase.

A strong spring Wyrtki jet (WJ) presents in May 2013 in the eastern equatorial Indian Ocean. The entire buildup and retreat processes of the spring WJ were well captured by two adjacent Acoustic Doppler Current Profilers mounted on the mooring systems. The observed zonal jet behaved as one intraseasonal event with the significant features of abrupt emergence as well as slow disappearance. Further research illustrate that the pronounced surface westerly wind burst during late-April to mid-May, associated with the active phase of a robust eastward-propagating Madden–Julian oscillation in the tropical Indian Ocean, was the dominant reason for the rapid acceleration of surface WJ. In contrasting, the governing mechanism for the jet termination was equatorial wave dynamics rather than wind forcing. The decomposition analysis of equatorial waves and the corresponding changes in the ocean thermocline demonstrated that strong WJ was produced rapidly by the wind-generated oceanic downwelling equatorial Kelvin wave and was terminated subsequently by the westward-propagating equatorial Rossby wave reflecting from eastern boundaries of the Indian Ocean.

Interannual variability (IAV) in the barrier layer thickness (BLT) and forcing mechanisms in the eastern equatorial Indian Ocean (EEIO) and Bay of Bengal (BoB) are examined using monthly Argo data sets during 2002–2017. The BLT during November–January (NDJ) in the EEIO shows strong IAV, which is associated with the Indian Ocean dipole mode (IOD), with the IOD leading the BLT by two months. During the negative IOD phase, the westerly wind anomalies driving the downwelling Kelvin waves increase the isothermal layer depth (ILD). Moreover, the variability in the mixed layer depth (MLD) is complex. Affected by the Wyrtki jet, the MLD presents negative anomalies west of 85°E and strong positive anomalies between 85°E and 93°E. Therefore, the BLT shows positive anomalies except between 86°E and 92°E in the EEIO. Additionally, the IAV in the BLT during December–February (DJF) in the BoB is also investigated. In the eastern and northeastern BoB, the IAV in the BLT is remotely forced by equatorial zonal wind stress anomalies associated with the El Niño-Southern Oscillation (ENSO). In the western BoB, the regional surface wind forcing-related ENSO modulates the BLT variations.

The comprehensive three-dimensional structures of an anti-cyclonic mesoscale eddy (AE) in the subtropical northwestern Pacific Ocean were investigated by combining the Argo floats profiles with enhanced vertical and temporal sampling and satellite altimetry data. The AE originated near the Kuroshio Extension and then propagated westward with mean velocity of 8.9 cm/s. Significant changes and evolutions during the AE’s growing stage (T1) and further growing stage (T2) were revealed through composite analysis. In the composite eddy core, maximum temperature (T) and salinity (S) anomalies were of 1.7 (1.9)°C and 0.04 (0.07) psu in T1 (T2) period, respectively. The composite T anomalies showed positive in almost whole depth, but the S anomalies exhibited a sandwich-like pattern. The eddy’s intensification and its influence on the intermediate ocean became more significant during its growth. The trapping depth increased from 400×10⁴ Pa to 580×10⁴ Pa while it was growing up, which means more water volume, heat and salt content in deeper layers can be transported. The AE was strongly nonlinear in upper oceans and can yield a typical mean volume transport of 0.17×10⁶ m³/s and a mean heat and salt transport anomaly of 3.6×10¹¹ W and –2.1×10³ kg/s during the observation period. The Energy analysis showed that eddy potential and kinetic energy increased notably as it propagated westward and the baroclinic instability is the major energy source of the eddy growth. The variation of the remained Argo float trapped within the eddy indicated significant water advection during the eddy’s propagation.

In this study, oil spill experiments were performed in a water tank to determine changes in the surface scattering characteristics during the emulsification of oil spills. A C-band fully-polarimetric microwave scatterometer and a vector network analyzer were used to observe films of the following oils: crude oil with an asphalt content below 3% that is prone to emulsification (type A), fresh crude oil extracted from an oilfield (type B), and industrial crude oil that was dehydrated and purified (type C). The difference in the backscatter results between the emulsified oil film and the calm water surface under C-band microwaves and the influence of the emulsification of the oil film on the backscatter were analyzed in detail. The results demonstrate that under a low-wind and no-waves condition (the maximum wave height was below than 3 mm), the emulsification of crude oil could modulated the backscatter through changes in the surface roughness and the dielectric constant, where the surface roughness had the dominant effect. The surface backscatters of the type B oil were greater than that of the type C oil in both the emulsified and non-emulsified states. In the non-emulsified state, the average differences in the backscatter between the type B and C oils were 2.19 dB, 2.63 dB, and 2.21 dB for the polarization modes of VV, HH, and HV/VH, respectively. Smaller corresponding average differences of 0.98 dB, 1.49 dB, and 1.5 dB were found for the emulsified state with a 20% moisture constant for the oil film. The results demonstrated that the surface roughness of the different oil films could vary due to the differences in the oil compositions and the oil film properties, which in turn affect the backscatter of the oil film surface.

In previous studies, Lagrangian analyses were used to assess large-scale ocean circulation, and the Lagrangian coherent structure could also reveal the evolution of the two-dimensional structure of the mesoscale eddies. However, few studies have demonstrated the three-dimensional structure of the mesoscale eddies via Lagrangian analysis. Compared with previous studies, which investigated the eddy structure via a Eulerian view, we used a Lagrangian view to provide a different perspective to study the eddy structure. An idealized cyclonic mesoscale eddy is built up over a seamount, and it presents downwelling inside the eddy and upwelling alongside the eddy formed within a closed circulation system. This structure is difficult to display via a Eulerian analysis. However, the trajectories of particles can well demonstrate the full cycle: the fluid sank and rotated inside the eddies, converged to the upwelling zone of the bottom layer and returned to the surface through upwelling. We also applied a Lagrangian analysis to a realistic simulation. As a significant phenomenon in the South China Sea, the dipole structure of the anticyclonic eddy (AE)/cyclonic eddy (CE) pair off of central Vietnam has been well studied but mainly at the sea surface. With a Lagrangian analysis, we illustrate the three-dimensional structure of the eddy pair: the fluid sank (rose) and rotated inside the AE (CE). More importantly, the trajectories of the particles suggested that there was no fluid exchange between the two eddies since the strong boundary jet separates them from each other. All the conclusions above have been verified and are supported by the computational error estimate. With a selected time step and integral period, the computational errors always present small values, although they increase with strong divergent and vertical diffusive flow.

The sea surface temperature (SST) seasonal cycle in the eastern equatorial Pacific (EEP) plays an important role in the El Niño–Southern Oscillation (ENSO) phenomenon. However, the reasonable simulation of SST seasonal cycle in the EEP is still a challenge for climate models. In this paper, we evaluated the performance of 17 CMIP6 climate models in simulating the seasonal cycle in the EEP and compared them with 43 CMIP5 climate models. In general, only CESM2 and SAM0-UNICON are able to successfully capture the annual mean SST characteristics, and the results showed that CMIP6 models have no fundamental improvement in the model annual mean bias. For the seasonal cycle, 14 out of 17 climate models are able to represent the major characteristics of the observed SST annual evolution. In spring, 12 models capture the 1–2 months leading the eastern equatorial Pacific region 1 (EP1; 5°S–5°N, 110°–85°W) against the eastern equatorial Pacific region 2 (EP2; 5°S–5°N, 140°–110°W). In autumn, only two models, GISS-E2-G and SAM0-UNICON, correctly show that the EP1 and EP2 SSTs vary in phase. For the CMIP6 MME SST simulation in EP1, both the cold bias along the equator in the warm phase and the warm bias in the cold phase lead to a weaker annual SST cycle in the CGCMs, which is similar to the CMIP5 results. However, both the seasonal cold bias and warm bias are considerably decreased for CMIP6, which leads the annual SST cycle to more closely reflect the observation. For the CMIP6 MME SST simulation in EP2, the amplitude is similar to the observed value due to the quasi-constant cold bias throughout the year, although the cold bias is clearly improved after August compared with CMIP5 models. Overall, although SAM0-UNICON successfully captured the seasonal cycle characteristics in the EEP and the improvement from CMIP5 to CMIP6 in simulating EEP SST is clear, the fundamental climate models simulated biases still exist.

This study examines wave reflection by a multi-chamber partially perforated caisson breakwater based on potential theory. A quadratic pressure drop boundary condition at perforated walls is adopted, which can well consider the effect of wave height on the wave dissipation by perforated walls. The matched eigenfunction expansions with iterative calculations are applied to develop an analytical solution for the present problem. The convergences of both the iterative calculations and the series solution itself are confirmed to be satisfactory. The calculation results of the present analytical solution are in excellent agreement with the numerical results of a multi-domain boundary element solution. Also, the predictions by the present solution are in reasonable agreement with experimental data in literature. Major factors that affect the reflection coefficient of the perforated caisson breakwater are examined by calculation examples. The analysis results show that the multi-chamber perforated caisson breakwater has a better wave energy dissipation function (lower reflection coefficient) than the single-chamber type over a broad range of wave frequency and may perform better if the perforated walls have larger porosities. When the porosities of the perforated walls decrease along the incident wave direction, the perforated caisson breakwater can achieve a lower reflection coefficient. The present analytical solution is simple and reliable, and it can be used as an efficient tool for analyzing the hydrodynamic performance of perforated breakwaters in preliminary engineering design.