Ocean Circulation Research Articles

Distribution patterns, risk assessment and potential sources of heavy metals in sediment in the Qiongzhou Strait, China

Risk assessment and pollutant source analysis are crucial tools for the management and protection of coastal ecosystems. The distribution patterns, risk assessment, and potential sources of heavy metals (Cd, Cr, Fe, Ni, Cu, Zn, As, Hg, and Pb) in surface sediment were analyzed in the Qiongzhou Strait, China, in summer and autumn of 2022. Heavy metals in autumn showed higher ecological risk than that in summer. Seasonal shifts in ocean currents may result in variations in heavy metal accumulation and dispersion. Cd and Hg were the priority heavy metals found, and according to the Positive Matrix Factorization results, the study area contains five sources of pollution, with natural sources, shipping-related activities, and industrial activities being the primary contributors. This study indicated that pollutants from adjacent areas should be considered for managing the environmental quality of Qiongzhou Strait.

Regional Sea Level Response to External Forcings from the Twentieth to the Twenty-First Century

Abstract Using a range of Detection and Attribution Model Intercomparison Project (DAMIP) simulations from phase 6 of Coupled Model Intercomparison Project (CMIP6), we study the response of dynamic sea level (DSL) to external anthropogenic climate forcing [greenhouse gases (GHGs), aerosols, and stratospheric ozone] with a focus on the differences over the twentieth and twenty-first century. In the second half of the twentieth century, the DSL nonuniformity in the Northern Hemisphere (NH) was relatively small due to a cancellation between the effects of increasing GHGs and aerosols. In contrast, the DSL signal in the Southern Hemisphere (SH) over this period was large because stratospheric ozone depletion reinforced the effects of increasing GHGs. In the twenty-first century, the DSL response has been intensified in the NH because the warming effects of diminishing aerosols have acted to reinforce the effects of increasing GHGs. Meanwhile, the distribution of SH DSL has also become uneven although stratospheric ozone recovery has partially offset the effects of rising GHGs. Using a global ocean circulation model, we decompose the changes in the twenty-first century DSL into distinct responses to surface forcings including sea surface temperature, salinity, and wind stress. Our results show that the dipole-like pattern of DSL in the North Pacific can be attributed largely to sea surface warming, while the dipole-like pattern in the North Atlantic is attributed to subpolar surface salinity freshening. The belted pattern of DSL changes in the Southern Ocean is induced by both surface warming and intensifying/poleward-shifting westerly winds.

Joshi–Malafarina–Narayan singularity in weak magnetic field

The importance and significance of magnetic fields in the astrophysical scenario is well known. Many domains of astrophysical black hole physics such as polarized shadow image, high energy emitting processes and jet formation are dependent on the behavior of the magnetic fields in the vicinity of the compact objects. In light of this, we determine the master equation and master differential equation that determine the spatial behavior of the magnetic field inside a matter distribution or vacuum region, of general spherically symmetric metric, which is immersed in a test magnetic field. We also investigate here the case of JMN-1 singularity immersed in a uniform weak magnetic field and determine the behavior of magnetic fields by defining electromagnetic four potential vector. We find that the tangential component of the magnetic field is discontinuous at the matching surface of the JMN-1 singularity with the external Schwarzschild metric, resulting in surface currents. We define the covariant expression of surface current density in this scenario. We also analyze the behavior of center-of-mass energy of two oppositely charged particles in the geometry of the magnetized JMN-1 singularity. We briefly discuss the possible scenarios which would possess a discontinuous magnetic field and implications of the same and future possibilities in the realm of astrophysics are indicated.

A tipping point in stable isotope composition of Antarctic meteoric waters during Cenozoic glaciation.

Triple oxygen isotopes of Cenozoic intrusive rocks emplaced along the Ross Sea coastline in Antarctica, reveal that meteoric-hydrothermal waters imprinted their stable isotope composition on mineral phases, leaving a clear record of oxygen and hydrogen isotope variations during the establishment of the polar cap. Calculated O- and H-isotope compositions of meteoric waters vary from -9 ± 2‰ and -92 ± 5‰ at 40 ± 0.6 Ma, to -30 and -234 ± 5‰ at 34 ± 1.9 Ma, and intersect the modern Global Meteoric Water Line. These isotopic variations likely depict the combined variations in temperature, humidity, and moisture source regions, resulting from rearrangement of oceanic currents and atmospheric cooling during the onset of continental ice cap. Here, we report a paleo-climatic proxy based on triple oxygen geochemistry of crystalline rocks that reveals changes in the hydrological cycle. We discuss the magnitude of temperature changes at high latitudes during the Eocene-Oligocene climatic transition.

VO2 based multi-functional ultra-wideband terahertz meta-absorber for EMI shielding application

ABSTRACT This paper presents a multifunctional tuneable wideband metasurface absorber based on phase-changing Vanadium dioxide (VO2). The ‘quad isotropic VO2 circular sectors’ are constructed on amorphous silicon dioxide (SiO2), backed with 0.2 µm-thick gold layer to design the proposed absorber. This novel wideband tunable structure can independently be used as a reflector and absorber as an EMI shield; hence, this multi-functionality in a single structure proves novel and unique in the terahertz frequency range. This unit cell offers a bandwidth of 3.5 THz with more than 90% of absorptivity under normal wave incidence. The equivalent circuit model, surface currents, and field distributions are further investigated to provide a better understanding of structural evolution. Moreover, a wideband absorption has been observed up to the 60° incident angle for both the transverse modes, showing a polarization-insensitive behavior due to the symmetricity in the structure. The parametric study explores, optimizes, and demonstrates various other main characteristics and parameters of the proposed absorber. The novel structure offers shielding effectiveness (SE) higher than 40 dB throughout the frequency band under both normal and oblique incidence. This outstanding performance makes it a strong contender for EMI-EMC applications. Furthermore, it shows great potential in the field of terahertz technology, such as terahertz imagers, detectors, sensors, and other wideband applications.

Mathematical modeling of water waves interaction with trapezoidal-shaped breakwater in the presence of current

In this study, an investigation is conducted into the interaction between water waves and a bottom-standing, thick, porous, and trapezoidal-shaped breakwater under the influence of an uniform current. This investigation has been conducted using several effect factors, such as the friction coefficient and porosity of the breakwater and ocean currents, to comprehend various phenomena related to wave energy dissipation, wave load, wave transmission, and wave reflection, etc. The boundary element method-based numerical technique is used to solve the associated boundary value problem. The modeling of the porosity of the breakwater structure is done using the Sollit and cross model [CK Sollitt, RH Cross. Wave transmission through permeable breakwaters. Coast Eng. 1972; 1972: 1827–1846.] of water wave passes through thick porous structures. In the presence of the following and opposite currents, the ‘Doppler-shift’ effect is observed in wave transformation characteristics due to the shift in wave frequencies. The wave-blocking phenomenon in open water and porous regions has been demonstrated. The results show that the following and opposing currents shift the frequencies for which optima in wave reflection occurs.

Joint inversion of direct waves, primaries, and multiples for vertical cable seismic data

By anchoring vertical cable in seawater and receiving the data excited by the towing cable seismic source, the reflection data from subsurface structures with a high signal-to-noise ratio, wide frequency band, and wide azimuth angle can be obtained, thereby achieving high-quality imaging of the target area. However, due to drifting receivers driven by ocean currents and the nonuniform coverage caused by the vertical orientation of the receivers, the traditional imaging method cannot be directly used in the vertical cable seismic data. Furthermore, crosstalk between the primaries and multiples commonly contaminates the final imaging. A new inversion method, based on the Bayesian inference and the assumption that all errors are Gaussian distributed, is presented to address these problems. The direct waves are used to invert the source/receiver locations. The primary and multiple wavefields are accurately reconstructed based on the reconstructed direct wavefield, and the corresponding reflectivity updates are obtained. The model-domain Hessian matrices consisting of different data are approximated by their diagonal elements to mitigate the nonuniform coverage issue, and the related covariance matrices and data-domain Hessian matrices are defined and approximated, respectively, to alleviate the calculation problem. In addition, the zero-lag correlated wavefields in the gradient are consistent with the actual wavefield propagation process, e.g., a primary wave is generated when the direct wavefield encounters reflectivity, avoiding the crosstalk. Numerical examples for synthetic models and field data illustrate that our method can accurately recover the source/receiver locations and provide high-resolution balanced-amplitude imaging results without evident crosstalk.

Maze-enclosed quad symmetric kite shaped SRR based metamaterial absorber for doppler navigation aids and earth exploration satellite remote sensing services

This paper represents a wide incident angle-polarization insensitive, stable in results, and compact-sized metamaterial absorber that can be used for Doppler navigation aids and remote sensing applications in Earth Exploration Satellite Services. The unit cell dimension was subwavelength-sized of 0.112λ X 0.112λ, with a good effective medium ratio (8.9), and achieved polarization insensitivity at a maximum incident angle up to 180°, making it more effective for proposed applications. The absorber resonated at 3.855 GHz, 5.13 GHz, 9.855 GHz, 13.335 GHz, and 16.02 GHz with single negative metamaterial properties due to negative permeability and absorptions of 99.7 %, 99.2 %, 98.6 %, 98.1 %, and 97.9 %, respectively. An RLC equivalent circuit on ADS was analyzed to validate the simulated results using the dipoles created by the surface currents. The absorber is cost-effective as no lumped elements or multilayer substrate materials were used, making fabricating less time-consuming and simple. The simulated results were validated by measuring data from a vector network analyzer after fabricating it. The proposed metamaterial absorber is suitable for use in Doppler navigation and Earth Exploration Satellite Services, and this research can further be extended to beamforming and sensing applications in the future.

VISIR-2: ship weather routing in Python

Abstract. Ship weather routing, which involves suggesting low-emission routes, holds potential for contributing to the decarbonisation of maritime transport. However, including because of a lack of readily deployable open-source and open-language computational models, its quantitative impact has been explored only to a limited extent. As a response, the graph-search VISIR (discoVerIng Safe and effIcient Routes) model has been refactored in Python, incorporating novel features. For motor vessels, the angle of attack of waves has been considered, while for sailboats the combined effects of wind and sea currents are now accounted for. The velocity composition with currents has been refined, now encompassing leeway as well. Provided that the performance curve is available, no restrictions are imposed on the vessel type. A cartographic projection has been introduced. The graph edges are quickly screened for coast intersection via a K-dimensional tree. A least-CO2 algorithm in the presence of dynamic graph edge weights has been implemented and validated, proving a quasi-linear computational performance. The software suite's modularity has been significantly improved, alongside a thorough validation against various benchmarks. For the visualisation of the dynamic environmental fields along the route, isochrone-bounded sectors have been introduced. The resulting VISIR-2 model has been employed in numerical experiments within the Mediterranean Sea for the entirety of 2022, utilising meteo-oceanographic analysis fields. For a 125 m long ferry, the percentage saving of overall CO2 expenditure follows a bi-exponential distribution. Routes with a carbon dioxide saving of at least 2 % with respect to the least-distance route were found for prevailing beam or head seas. Two-digit savings, up to 49 %, were possible for about 10 d in a year. In the case of an 11 m sailboat, time savings increased with the extent of path elongation, particularly during upwind sailing. The sailboat's routes were made approximately 2.4 % faster due to optimisation, with the potential for an additional 0.8 % in savings by factoring in currents. VISIR-2 serves as an integrative model, uniting expertise from meteorology, oceanography, ocean engineering, and computer science, to evaluate the influence of ship routing on decarbonisation efforts within the shipping industry.

Causes and multiyear predictability of the rapid acceleration of U.S. Southeast Sea level rise after 2010

The rate of sea level rise (SLR) along the Southeast Coast of the U.S. increased significantly after 2010. While anthropogenic radiative forcing causes an acceleration of global mean SLR, regional changes in the rate of SLR are strongly influenced by internal variability. Here we use observations and climate models to show that the rapid increase in the rate of SLR along the U.S. Southeast Coast after 2010 is due in part to multidecadal buoyancy-driven Atlantic meridional overturning circulation (AMOC) variations, along with heat transport convergence from wind-driven ocean circulation changes. We show that an initialized decadal prediction system can provide skillful regional SLR predictions induced by AMOC variations 5 years in advance, while wind-driven sea level variations are predictable 2 years in advance. Our results suggest that the rate of coastal SLR and its associated flooding risk along the U.S. southeastern seaboard are potentially predictable on multiyear timescales.

Nonstationary thermodynamic and dynamic contributions to the interannual variability of winter sea ice growth in the Kara-Laptev Seas

Abstract The Kara-Laptev Seas (KLS), known as the “Ice Factory of the Arctic”, witnesses rising instead of falling winter sea ice growth (WSIG) under the shrinkage of Arctic ice. However, knowledge of the large year-to-year variation is still unclear. Combining a seasonal ice concentration budget, a composite analysis, and a typical case study, we study both the interannual variability of WSIG in the KLS and the associated air-sea forcings during 1985–2021. Results quantitatively reveal that, during 1985–2021, thermodynamic melt in the melting season (Apr–Aug) contributed 80.3% to the interannual ice loss difference and promoted the subsequent WSIG by the recovery mechanism in the KLS. This consistent thermodynamic melt is caused by the strengthened summer Beaufort High, transporting heat and introducing a locally positive ice-albedo feedback. However, since 2010, the dynamic growth during the freezing season (Oct–Feb) has increasingly stimulated the WSIG. Typical cases in 2013 and 2017 indicate that the overlying anticyclonic atmospheric regime restricts the ice drift from the KLS and contributes to the dynamic growth of 41.6% of the WSIG difference, while the turbulent-heat-induced thermodynamic growth in winter is down to 58.4%. In short, we reveal an unstable relationship between the summer ice loss and the subsequent WSIG under the background of Arctic warming. Our study points out that the distinct dynamic ice growth driven by surface winds or ocean currents during the freezing season is likely to increase in the near future, with thinner and more mobile seasonal ice predominating in the Arctic.

Northwestern Pacific Oceanic circulation shaped by ENSO

The intricate currents of the Northwest Pacific Ocean, with strong manifestations along the westside rim, connect tropical and subtropical gyres and significantly influence East Asian and global climates. The El Niño/Southern Oscillation (ENSO) originates in the tropical Pacific Ocean and disrupts this ocean circulation system. However, the spatiotemporal dependence of the impact of ENSO events has yet to be elucidated because of the complexities of both ENSO events and circulation systems, as well as the increased availability of observational data. We thus combined altimeter and drifter observations to demonstrate the distinct tropical and subtropical influences of the circulation system on ENSO diversity. During El Niño years, the North Equatorial Current, North Equatorial Countercurrent, Mindanao Current, Indonesian Throughflow, and the subtropical Kuroshio Current and its Extension region exhibit strengthening, while the tropical Kuroshio Current weakens. The tropical impact is characterized by sea level changes in the warm pool, whereas the subtropical influence is driven by variations in the wind stress curl. The tropical and subtropical influences are amplified during the Centra Pacific El Niño years compared to the Eastern Pacific El Niño years. As the globe warms, these impacts are anticipated to intensify. Thus, strengthening observation systems and refining climate models are essential for understanding and projecting the enhancing influences of ENSO on the Northwest Pacific Oceanic circulation.

On the Driving Factors of the Future Changes in the Wintertime Northern‐Hemisphere Atmospheric Waviness

AbstractDespite the significant socioeconomic implications in the link between atmospheric waviness and extreme weather events, future atmospheric waviness trends remain elusive due to uncertainties arising from diverse definitions and insufficient dynamical formalism in existing metrics. This study employs a local wave activity (LWA) metric, whose prognostic equation links wave activity changes to forcing mechanisms, to assess wintertime Northern Hemisphere waviness in ERA5 and HighResMIP data sets. The models generally exhibit high fidelity in reproducing observed waviness, while biases stem primarily from biases in the LWA source, low‐level meridional heat flux, which tend to improve with higher resolutions. Future projections exhibit reduction in LWA, primarily due to suppressed LWA generation, which is mitigated by higher‐resolution models. We found that both biases and reduction of the LWA source are closely associated with sensible heat fluxes from the ocean to the atmosphere, highlighting the potential impacts of resolving ocean currents.

Characterizing the internal wave wakes and synthetic aperture radar image features of underwater moving objects

In this study, a refined and improved theoretical model was proposed for calculating the synthetic aperture radar (SAR) images of internal wave wakes. This model was based on an equivalent source model of internal wave wakes and the velocity-bunching principle. Wake field, scattering field, and SAR image characteristics were systematically calculated for multiple scenarios. The model described the transfer and solution processes by inputting the background flow field temperature and salinity profiles to map the radar image intensity distribution. The physical mechanisms of the internal wave wake generation and its capture by SAR were comprehensively explained. The effects of various parameters on the internal wave wake SAR images were analyzed comprehensively. The SAR image features containing internal wave wakes were determined based on the gray-level co-occurrence matrix and image power spectrum, which revealed that the maximum surface current induced by the internal wave modes of the same amplitude was considerably larger than that induced by surface-wave mode. When traveling far from the water surface, the rate of change in the total scattering coefficient because of the Bragg effect was approximately 20 times that of the slope effect. In SAR images, the recognizability of wakes depends primarily on the maximum intensity rather than the average disturbance on the sea surface. Therefore, SAR images exhibit considerable differences during downwind radar observations, with the most prominent contrast differences occurring at the maximum spreading angle of wake scattering and 1500 m behind underwater objects. In the image power spectrum, the frequency range of the first modal internal waves was approximately 10–30 Hz, whereas the higher-order modal internal waves had frequencies lower than 10 Hz. Furthermore, the frequency range of the surface-wave mode was larger than that of internal-wave modes, approximately 50–60 Hz.

Mechanism of formation and evolution of bound states in the continuum of split-hole all-dielectric metasurface under asymmetric displacement perturbation

Bound states in the continuum (BICs) with ultra-high Q properties have attracted much attention for their perfect localization in the continuous spectral range coexisting with extended waves. In this study, breaking the traditional excitation form of structure breakage or excitation field asymmetry, a monolithic silicon nanodisk array with relative displacement generated by the complete splitting of square nanopores is proposed based on the unique electromagnetic properties of all-dielectric metamaterials. During the introduction of perturbations by asymmetric displacements of splitting holes, it is shown by numerical simulations that two BICs at different wavelengths can be realized. Combined with eigenmodes of group theory, the symmetric matching relationship between the symmetry-protected BICs and the free-space radiation during the evolution process is analytically demonstrated, and the formation mechanism and the evolution law of the BICs excited by this metasurface are deeply investigated. meanwhile, it also provides a theoretical basis for the polarization dependence of quasi-BICs excitation and the ultra-high Q factor expression of BICs. Furthermore, near-field distribution and multipole decomposition show that the field distribution and surface currents support the excitation of BIC-driven toroidal dipole and magnetic quadrupole dual modes. This study not only provides an effective reference for the stability of high-Q resonance wavelengths, but also solves the problem of the lack of universality in analyzing the resonance mechanism based on resonance phenomena, and provides solid theoretical support for the study of displacement-mediated BICs resonance excitation and evolution.

Coupled Meteo–Hydrodynamic Approach in Semi-Enclosed Basins and Sensitivity Assessment of Wind-Driven Current

A coupled numerical approach that combines the WRF model and the Mike 3 (DHI) hydrodynamic model was developed and applied in two semi-enclosed basins in the Ionian Sea (Italy) to assess the wind-driven current. To gain a better understanding of how the sea current field can vary depending on meteorological data forcing, three different scenario were set up. The sensitivity of the sea current pattern was investigated as a function of the type of meteorological forcing and appreciating the differences in the results. The aims of this study are threefold. Firstly, we wish to define an ad hoc procedure to join the model-computed meteorological parameters in the hydrodynamic model. Secondly, we will investigate the feedback from the Mar Piccolo and Mar Grande basins in the Ionian Sea using fully coupled simulations and an uncoupled system where the atmospheric parameters are derived from a ground station. Finally, we will evaluate the results achieved by applying two scenarios of typical meteorological conditions to the study site. The model results highlighted the variability of sea currents depending on meteorological forcing.

Mesoscale eddy-strengthened deep-sea topographic Rossby waves in the southwestern South China Sea.

Topographic Rossby waves (TRWs) dominate the low-frequency variability of deep ocean currents and play a crucial role in energy exchange and material mixing. On the continental slope of the southwestern South China Sea, a deep-water mooring was deployed to observe TRWs for a period of ~ 40days. The TRWs, with a wavelength of 109km, account for 41.3% of the subinertial variations. A ray-tracing model was applied to investigate the propagation and energy source. The results showed that the TRWs propagated from the northeast of the mooring location and were most likely caused by the mesoscale eddy disturbances off the Vietnam coast. This study provides a new perspective on examining the impact of mesoscale eddies off Vietnam on abyssal currents.

Asymmetry of hemispheric climate during MIS 16 versus MIS 12

Marine isotope stage (MIS) 16 (676-621 ka) and MIS 12 (478-424 ka) stand out as two significant glacial periods during the late Quaternary, demonstrating consistent continental ice volume or sea level, as discerned from benthic δ18O records. Intriguingly, the average sea surface temperature (SST) during MIS 16 surpasses that of MIS 12 by approximately 0.54–0.79 °C, while the atmospheric CO2 concentration during MIS 16 is approximately 16 ppmv lower than that of MIS 12. The underlying mechanism driving these discrepancies remains enigmatic. In this study, we present new high-resolution X-ray fluorescence (XRF) and grain size records covering the interval from 700 to 400 ka from Ocean Drilling Program (ODP) Site 1090 in the South Atlantic Ocean. We also integrated benthic δ13C, ice-rafted debris (IRD) and carbonate ion concentration ([CO32−]) records globally, to elucidate the distinct climate evolution during the MIS 16 versus MIS 12. Our findings suggest that MIS 12 is distinguished by a more substantial Antarctic ice-sheet, coupled with a weakened and/or equatorward displacement of southern westerly winds, and expansion of southern source water (SSW) alongside contractive northern source water (NSW). Conversely, the opposite climate state can be observed in MIS 16. The elevated atmospheric CO2 in MIS 12 and diminished atmospheric CO2 in MIS 16 may be attributed to reconfigurations in the ocean carbon reservoir regulated by the reorganization of deep-water masses and meridional ocean circulation, as a result of the divergent evolution of bipolar ice-sheets. Our study underscores the critical role of the interplay between the two polar regions in modulating deep-sea carbon storage on the orbital timescale and atmospheric CO2. This interplay is essential for a comprehensive understanding of past and future carbon cycle adjustments in response to climate change.

Impact of assimilation of absolute dynamic topography on Arctic Ocean circulation

The ocean circulation is typically constrained in operational analysis and forecasting systems through the assimilation of sea level anomaly (SLA) retrievals from satellite altimetry. This approach has limited benefits in the Arctic Ocean and surrounding seas due to data gaps caused by sea ice coverage. Moreover, assimilation of SLA in seasonally ice-free regions may be negatively affected by the quality of the Mean Sea Surface (MSS) used to derive the SLA. Here, we use the Regional Ice Ocean Prediction System (RIOPS) to investigate the impact of assimilating Absolute Dynamic Topography (ADT) fields on the circulation in the Arctic Ocean. This approach avoids the use of a MSS and additionally provides information on sea level in ice covered regions using measurements across leads (openings) in the sea ice. RIOPS uses a coupled ice-ocean model on a 3-4 km grid-resolution pan-Arctic domain together with a multi-variate reduced-order Kalman Filter. The system assimilates satellite altimetry and sea surface temperature together with in situ profile observations. The background error is modified to match the spectral characteristics of the ADT fields, which contain less energy at small scales than traditional SLA due to filtering applied to reduce noise originating in the geoid product used. A series of four-year reanalyses demonstrate significant reductions in innovation statistics with important impacts across the Arctic Ocean. Results suggest that the assimilation of ADT can improve circulation and sea ice drift in the Arctic Ocean, and intensify volume transports through key Arctic gateways and resulting exchanges with the Atlantic Ocean. A reanalysis with a modified Mean Dynamic Topography (MDT) is able to reproduce many of the benefits of the ADT but does not capture the enhanced transports. Assimilation of SLA observations from leads in the sea ice appears to degrade several circulation features; however, these results may be sensitive to errors in MDT. This study highlights the large uncertainties that exist in present operational ocean forecasting systems for the Arctic Ocean due to the relative paucity and reduced quality of observations compared to ice-free areas of the Global Ocean. Moreover, this underscores the need for dedicated and focused efforts to address this critical gap in the Global Ocean Observing System.

Study on Fault Diagnosis Technology for Efficient Swarm Control Operation of Unmanned Surface Vehicles

The purpose of this study is to design a Swarm Control algorithm for the effective mission performance of multiple unmanned surface vehicles (USVs) used for marine research purposes at sea. For this purpose, external force information was utilized for the control of multiple USV swarms using a lead–follow-formation technique. At this time, to efficiently control multiple USVs, the LSTM algorithm was used to learn ocean currents. Then, the predicted ocean currents were used to control USVs, and a study was conducted on behavioral-based control to manage USV formation. In this study, a control system model for several USVs, each equipped with two rear thrusters and a front lateral thruster, was designed. The LSTM algorithm was trained using historical ocean current data to predict the velocity of subsequent ocean currents. These predictions were subsequently utilized as system disturbances to adjust the controller’s thrust. To measure ocean currents at sea as each USV moves, velocity, azimuth, and position data (latitude, longitude) from the GPS units mounted on the USVs were utilized to determine the speed and direction of the hull’s movement. Furthermore, the flow rate was measured using a flow rate sensor on a small USV. The movement and position of the USV were regulated using an Artificial Neural Network-PID (ANN-PID) controller. Subsequently, this study involved a comparative analysis between the results obtained from the designed USV model and those simulated, encompassing the behavioral control rules of the USV swarm and the path traced by the actual USV swarm at sea. The effectiveness of the USV mathematical model and behavior control rules were verified. Through a comparison of the movement paths of the swarm USV with and without the disturbance learning algorithm and the ANN-PID control algorithm applied to the designed simulator, we analyzed the position error and maintenance performance of the swarm formation. Subsequently, we compared the application results.