Coast-fitting tomographic inversion that is based on function expansion using three types of normal modes (the Dirichlet, Neumann, and open boundary modes) is proposed to reconstruct current fields from the coastal acoustic tomography (CAT) data. The superiority of the method was validated while using CAT data that were obtained in 2015 in the Dalian Bay. The semidiurnal tidal and residual current fields were accurately reconstructed over the entire model domain surrounded by coasts and open boundaries. The proposed method was effective, particularly around the peripheral regions of the tomography domain and the near-coast regions outside the domain, where accurate results are not expected from the conventional inverse method based on function expansion by Fourier function series with no coast fittings. The error velocity for the semidiurnal tidal currents was 2.2 cm s⁻¹, which was calculated from the root-mean-square-difference between the CAT-observed and inverted range-averaged currents that were obtained along the nine peripheral transmission paths. The error velocity for the residual currents estimated from the 12-h mean net residual transport at the bay mouth was 0.9 cm s⁻¹. The errors were significantly smaller than the amplitude of the tidal and residual currents.

The first coastal acoustic tomography (CAT) experiment site of the Neko-Seto Channel was revisited to elucidate the propagation and generation characteristics of the M₂ and M₄ tidal currents with a second CAT experiment, which was conducted from 3–6 April 2018 (local time). Two-dimensional flow fields of the M₂ and M₄ tidal currents and the residual current were reconstructed using a coast-fitting inversion model with the reciprocal travel-time data of five acoustic stations. The M₂ tidal current is a progressive-type wave that propagates eastward at a speed of 0.7 ms⁻¹, much slower than expected for free progressive tides in this region (19 ms⁻¹). The M₄ nonlinear current constructed an out-of-phase relationship between the western and eastern halves of the tomography domain, implying the generation of standing-type waves. Such nonlinear processes led to flood- and ebb-dominant tidal current asymmetries for the western and eastern halves of the model domain, respectively. The two-day mean residual currents constructed a northeastward current with a maximum speed of 0.3 ms⁻¹ in the western half of the model domain and a clockwise rotation in the eastern half. The averaged inversion errors were 0.03 ms⁻¹, significantly smaller than the amplitude of the aforementioned currents.

Herein, the temporal variability of the deep overflow through the Kerama Gap between the East China Sea and the Philippine Sea is investigated based on observational data combined with reanalysis data obtained during 2004–2011. The observations and model results show a strong bottom-intensified flow intruding into the deep Okinawa Trough. The observed deep overflow shows intraseasonal variations that are enhanced from August to November. The variability in the deep overflow via the Kerama Gap is well-correlated with the density changes near its sill depth in the Philippine Sea. Additionally, some portion of the dense water originates from a region east of Miyakojima, which can be related to the northeastward-flowing Ryukyu Current at intermediate depths. In contrast, three extreme deep overflow events indicate that the arriving mesoscale eddies propagated from the east resulted in an increase in the density near the Kerama Gap sill than that on the Okinawa Trough side. The density difference associated with the baroclinic pressure gradient across the Kerama Gap forced the deep overflow into the Okinawa Trough. The volume transport of the deep overflow computed by integrating the cross-sectional velocity and through hydraulic theory are 0.14 and 0.11 Sv (1 Sv = 10⁶ m³/s), respectively.

To accurately reconstruct large-area three-dimensional current fields in coastal regions, simultaneous observations with ten coastal acoustic tomography (CAT) stations and two high-frequency radar (HFR) stations were performed in the Xiangshan Bay (XSB) on 4–5 December 2020. The section-averaged velocity that was observed by CAT and the radial velocity that was observed by HFR were, for the first time, synchronously assimilated into a three-dimensional barotropic ocean model. Compared with acoustic Doppler current profile data, the velocities of the model assimilating both CAT and HFR data had the highest accuracy according to root mean square differences (RMSDs), ranging from 0.05 to 0.08 m/s for all the vertical layers. For the models individually assimilating CAT and HFR, the values in the vertical layers ranged from 0.07 to 0.12 m/s and 0.08 to 0.13 m/s, respectively. A harmonic analysis of the model grid data showed that the spatial mean amplitudes of the M2, M4, and residual currents were 0.66, 0.14, and 0.09 m/s, respectively. Furthermore, the standing wave characteristics of the M2 current and M4 associated-oscillation in the inner XSB, mouth of the Xiangshan fjord, were better captured by the model assimilating both CAT and HFR. Our study demonstrates the advances in three-dimensional tidal current analysis using a model that assimilates both CAT and HFR data, especially in regions with complex coastal geography.