Ocean Circulation Research Articles

Energy Conversion by Helical Hydrokinetic Turbine in A Pipe

Hydrokinetic turbines are mechanisms designed for the purpose of utilizing the kinetic energy present in the movement of water bodies like rivers, tidal currents, or ocean currents, and transforming it into electrical power. These turbines function based on a principle akin to that of wind turbines; however, they are positioned underwater to harness the energy of the water flow. This study focuses on the fundamentals of hydrokinetic turbines and presents existing research. Additionally, simulations have been conducted to observe how the hydrokinetic turbine responds hydrodynamically inside a pipe. A three-bladed vertical-axis helical hydrokinetic turbine was installed within a circular conduit and subjected to analysis under varying flow conditions. The k-ω SST turbulence model was employed in the analyses. The results indicated that increasing the turbine's angular velocity initially raises the torque and the power coefficient until a peak is reached, after which the power coefficient decreases. The highest power coefficient was observed at a flow velocity of 2 m/s. Additionally, consistent with previous studies, the hydrokinetic turbine within the pipe surpassed the Betz limit.

Order of magnitude increase in power from flow-induced vibrations

Natural hydro environments such as rivers and ocean currents provide abundant and essentially free sources of kinetic energy. Although significant research has been devoted to harnessing this power through traditional turbine systems, there are also studies on the novel use of flow-induced vibration of a bluff body. The large oscillations of an elliptical cylinder resulting from the newly-discovered hyper branch have raised interesting questions about its potential as a source of effective renewable energy. This study investigates the power extraction characteristics of the hyper branch over a range of flow velocities. A factor-of-ten increase in the maximum time-mean power coefficient relative to an established flow-induced vibration system — the vortex induced vibration for aquatic clean energy converter is demonstrated. This suggests the possibility of a new and effective strategy to utilise flow-induced vibrations for energy generation and could help drive the commercial implementation and uptake of oscillating energy converters.

Numerical experiments investigating the influence of drag on trajectory patterns of floating macroalgae

Abstract Ocean currents are a crucial means of dispersing natural and human-made materials on the ocean surface. Macroalgae are among the most conspicuous natural dispersers, often called the ‘tumbleweeds of the ocean.’ Despite numerous studies on the subject, the relative influence of wind and surface currents on the trajectory of macroalgal dispersal remains uncertain. Previous studies have focused on kelp rafts of varying sizes, making it challenging to determine the impact of wind versus currents. These studies have also disregarded the macroalgae’s drag characteristics and surface area, which have been shown to impact the trajectory and accumulation of floating flotsam. This numerical study aims to shed light on the relative influence of wind and currents and the role of drag in determining the course and accumulation of macroalgae. By comparing simulations of virtual kelp ‘particles’ that incorporate drag and those without, this study focused on solitary kelp plants and considered the impact of morphological characteristics, flow-field combinations, and the presence of Stokes drift. Our results show that virtual kelp particles generally followed ocean currents, but the inclusion of drag caused deviations from purely Lagrangian particles’ trajectories and sheds light on the complex interplay of factors affecting macroalgal dispersal in the ocean.

A functional-group-based perspective on the response of marine phytoplankton to mesoscale eddies

This study analyzed the response of marine phytoplankton to environmental changes induced by mesoscale warm eddies through the lens of functional groups, highlighting the complex interactions within the ecosystem. It was found that warm eddies significantly affected phytoplankton distribution, with cell abundance in the center being only 75.60 cells/L, compared to 1095.00 cells/L in the periphery. Vertical transport within warm eddies altered light conditions, affecting photophilic diatoms more, while increased temperatures favored the growth of warm-water dinoflagellates. This study also emphasized that ocean currents were significant factors, showing correlations with various functional groups and playing a key role in material transport and phytoplankton distribution. Additionally, the distinct responses of different functional groups to temperature and salinity underscored their unique adaptations to environmental changes. In periods without warm eddies, phytoplankton primarily congregated in shallower water layers.

Mapping Ocean Gateways from 3D Seismic Data Using Color Processing and Analogues from Seismic Data

Mapping sedimentary features to delineate ancient ocean gateways from 3D seismic data has long been a challenge due to their limited vertical resolution. We present the technique of color processing of seismic data which improves the vertical resolution to very few seismic samples corresponding to 5–10 m. The result is a representation of the seismic impedance contrast correlation in Red-Green-Blue (RGB) colors, which creates images similar to satellite images that assist the geological interpretation because the RGB colors can be associated with mud-rich and sand-rich sediments. Using this approach, we have extracted six different sedimentary bedforms from different bottom flow settings and correlated these with the bedform-velocity matrix by Stow et al. (2009). Then we studied bedform assemblies in their regional context for the area covered by the 3D seismic data set. We conclude with comparing a modern and an ancient analogue – both from 3D seismic data – and infer paleo flow direction and velocity using oceanographic measurements available for the modern analogue. The method provides insights into ocean bottom currents and their impact on the formation of contourite depositional systems in modern and ancient ocean gateways which contribute to the reconstruction of the paleoclimate as an input to present day climate change models.

Evaporite sequences as archives for Mg isotope compositions of seawater - Evidence from a Tethys marginal shelf basin in the Anisian

Evaporites have recently been suggested as a potential archive for recording the Mg isotope compositions (δ26Mg) of coeval seawater. However, episodic dolomitization during the deposition of massive evaporites could cause considerable Mg removal and isotopic fractionation. To constrain the hydrological changes and influence of dolomitization on ambient brine δ26Mg, we present petrographic and mineralogical features, as well as Mg-C-O-Sr isotope data extracted from carbonate phases of a middle Triassic (ca. 247 Myr) marine anhydrite-dolostone sequence from a drill core in eastern China. The drilled lithologies are characterized by massive dolostone layers in the lower part, followed by an upward decline in dolomite contents accompanied by a rise in anhydrite. Multiple lines of evidence consistently point to a syn-depositional origin for the dolostones in a marginal basin of the Tethys Ocean and a lack of diagenetic alteration since deposition. We reconstructed the dynamic changes of δ26Mg values of basin waters based on Mg isotope compositions in dolomite leachates. According to our findings, the δ26Mg values of the basin waters were remarkably high (about 0.38 ± 0.05‰) at the onset of evaporation, indicating a significant Mg sink of the massive dolomitization. The δ26Mg of brine in the basin then changed towards the value of coeval seawater (about −0.32 ± 0.05‰) starting with the deposition of evaporites. Concurrently, 87Sr/86Sr ratios of dolomites shift from radiogenic values towards contemporaneous seawater composition. Our results demonstrate that the evaporite basin was not strictly restricted, and water exchange with the open ocean never ceased. Modeling calculation reveals that, even when the seawater exchange rate is far below the average ocean circulation, δ26Mg of brine in the basin will reach the value of open ocean within 1 Myr, completely removing the influence of early dolomitization. We suggested that massive marine evaporite sequences have the potential to record the long-term evolution of seawater Mg isotopes.

High-resolution UK’37 sea surface temperature reconstruction in the southwestern Atlantic Ocean during the Holocene

Modern sea surface temperatures (SST) were assessed using the alkenone calibration index (UK′37; SST-UK′37) values from twenty-one core-top sediment samples collected along the Southwestern Atlantic Ocean and validated by comparison with satellite-derived SSTs. The SST-UK′37 values revealed an annual signal between 28 and 30°S and a winter signal (83% of the region samples) between 30 and 34° S, with deviations falling within the UK′37 calibration error range. The paleotemperature (SST-UK′37) reconstruction through the Holocene was performed using two high-resolution marine sediment cores from 33 to 34°S latitude band. Millennial-scale variations dominate most of the SST-UK′37 results with changes up to ∼6 °C. Considering the modern ocean circulation dynamics, the SST-UK’37 in our records likely stems from the combined influence of the Subtropical Shelf Front (STSF) position and the interaction between the warm, salty, and nutrient-poor Subtropical Shelf Water (STSW) and the cold, fresh, and nutrient-rich Subantarctic Shelf Water (SASW). Our SST-UK’37 results also suggest that the northernmost position limit of the STSF reached around 33°S during the mid-to-late Holocene, while after ∼5500 cal yr BP the Plata Plume Water became the dominant influence on the regional SST-UK’37 signal. Spectral analyses revealed statistically significant oscillations of ∼8600 and ∼ 2250 years in record MDBT 557, and ∼ 5300 and ∼ 2600 years in record MDBT 561 which suggests the influence of Southern Westerly Winds (SWW) on the SST-UK’37, and therefore, on the STSF dynamic through the Holocene.

Numerical Modeling and Analysis of Pendant Installation Method Dynamics Using Absolute Nodal Coordinate Formulation

Accurately simulating the deployment process of coupled systems in deep-sea environments remains a significant challenge. This study employs the Absolute Nodal Coordinate Formulation (ANCF) to dynamically model and analyze multi-body systems based on the Pendant Installation Method (PIM). Utilizing the principle of energy conversion, this study calculates the stiffness, generalized elastic forces, mass matrices, and Morison equation, formulating a motion equation for the dynamic coupling of nonlinear time-domain forces in cables during pendulum deployment, which is numerically solved using the implicit generalized-α method. By comparing the simulation results of this model with those from the catenary theory model, the advanced modeling capabilities of this model are validated. Lastly, the sensitivity of the multi-body system under various boundary conditions is analyzed. The results indicate that deployment operations are more effective in environments with strong ocean currents. Furthermore, upon comparing the impacts of structural mass and deployment depth on the system, it was found that deployment depth has a more significant effect. Consequently, the findings of this study provide a scientific basis for formulating subsequent optimization strategies.

Singular Phenomenon Analysis of Wind-Driven Circulation System Based on Galerkin Low-Order Model

Ocean circulation plays an important role in the formation and occurrence of extreme climate events. The study shows that the periodic variation of ocean circulation under strong wind stress is closely related to climate oscillation. Ocean circulation is a nonlinear dynamic system, which shows complex nonlinear characteristics, so the essence behind ocean circulation has not been clearly explained. Therefore, the response and evolution of the circulation system to wind stress are studied based on the bifurcation and catastrophe theories in nonlinear dynamics. First, the quasi-geostrophic gyre equation and the normalized gravity model are introduced and developed to study ocean circulation driven by wind stress, and solved using the Galerkin method. Then, the bifurcation and catastrophe behaviors of the system governed by the low-order ocean circulation model during the change in wind stress intensity and the coexistence of multiple equilibria in the circulation system are studied in detail. The results show that saddle and unstable nodes appear in the system after a cusp catastrophe. With the change in model parameters, the unstable node becomes the unstable focus, and then the subcritical Hopf bifurcation occurs. The system forms a bistable interval when the system undergoes a catastrophe twice, and the system shows hysteresis. In addition, multiple equilibrium states are coexisting in the circulating system, and the unstable equilibrium state always changes into a stable equilibrium state through vortex movement. Therefore, there is a route for the system to induce short-term climate oscillation, that is, in the multi-stable equilibrium state of the system, the vortex oscillates after being subject to small disturbances, and then triggers climate oscillation.

Sea state contributions to thermosteric sea-level in high-resolution ocean-wave coupled simulations

This study examines the impact of wave-induced processes (WIPs) in modulating thermosteric sea-level changes, highlighting the need to include these processes in future sea-level rise assessments and climate projections. The impact of wave-induced processes on thermosteric sea-level changes is investigated using coupled ocean-wave simulations. These simulations include the effects of Stokes-Coriolis forcing, sea-state dependent surface stress and energy fluxes, and wave-induced mixing. The experiments use a high-resolution configuration of the Geesthacht COAstal Model SysTem (GCOAST), covering the Northeast Atlantic, the North Sea and the Baltic Sea. The GCOAST system uses the Nucleus for European Modelling of the Ocean (NEMO) ocean model to account for wave-ocean interactions and ocean circulation. It is fully coupled with the WAM spectral wind wave model. The aim is to accurately quantify the sea state contribution to thermosteric sea level variability and trends over a 26-year period (1992–2017). The ability of wave-ocean coupled simulations to reveal the contribution of sea state to sea level variability and surge is demonstrated. It is clear that wave-induced processes (WIPs) play a significant role in sea surface dynamics, ocean mixing (mixed layer thickness) and modulation of air-sea fluxes (e.g. heat flux) in both winter (10–20%) and summer (10%), which in turn affect thermosteric sea level variability. The North Atlantic (in summer) and the Norwegian Trench (in winter) show significant contributions (40%) to the thermosteric sea-level variability due to wave-induced processes. The influence of WIPs on thermosteric sea level trends in the North Atlantic is up to the order of 1 mm yr-1 in both winter and summer, in the open ocean and at the shelf break. Smaller contributions are observed over the shelf areas of the North Sea. This study underscores the crucial role of WIPs in modulating sea-level changes and highlights the importance of including these processes in future sea-level rise assessments and climate projections.

From source to sink: part 1-characterization and Lagrangian tracking of riverine microplastics in the Mediterranean Basin.

The Mediterranean Sea is one of the most critically polluted areas due to its semi-enclosed structure and its highly anthropized shoreline. Rivers are significant vectors for pollutant transfers from the continental to the marine environment. In this context, a 3D Lagrangian simulation of the dispersion of riverine microplastics (MPs) was performed, which included the application of a recently developed model that reassessed the MP fluxes discharged by rivers. MP physical properties from river samples were further investigated to approximate vertical displacement in modeled ocean currents. The use of a high-resolution circulation model, integrating Stokes drift, turbulent diffusion, and MP sinking and rising velocities, enabled us to establish stock balances. Our simulation suggested that 65% of river inputs may be made of floating MPs drifting in the surface layer and 35% of dense MPs sinking to deeper layers. The Eastern Mediterranean tends to accumulate floating MPs, primarily originating from the Western Mediterranean Basin, where major river sources are concentrated. After 2 years of simulation, modeled stranding sequestered 90% of the MP inputs, indicating relatively short average residence times from a few days to months at most for particles at sea. Although spatial distribution patterns stabilized after this period and a steady state may have been approached, the surface concentrations we modeled generally remained below field observations. This suggested either an underestimation of sources (rivers and unaccounted sources), by a factor of 6 at most, or an overestimation of MP withdrawal through stranding, to be reduced from 90 to around 60% or less if unaccounted sinks were considered.

Bathymetry-constrained ocean geostrophic currents play a key role in shaping the sea ice circulation in the Canada Basin, Arctic Ocean

Abstract Satellite-based observations and a pan-Arctic coupled sea ice-ocean model are utilized to study the effect of ocean geostrophic currents on large-scale sea ice circulation in the Canada Basin, Arctic Ocean. We find that surface winds primarily drive sea ice drifts in the west-east direction, while the geostrophic currents in the Beaufort Gyre promote north-south ice drifts. Wind fluctuations can create variable ice drifts, yet geostrophic currents respond more slowly due to their larger vertical scale, serving as a slowly-evolving conveyor belt for maintaining the anticyclonic ice circulation. It is further demonstrated that the bathymetry can regulate the movement of sea ice via constraining the expansion of ocean circulation. This mechanism is indirect in the sense that the ice is far from the seafloor. Our research underscores the necessity of considering the bathymetry-constrained geostrophic currents in understanding Arctic sea ice dynamics. With the rapid retreat of Arctic sea ice, the multi-scale interactions between ice drifts and ocean currents may have significant implications for the Arctic ecosystem, climate, and shipping corridors.

Arctic September Sea Ice Concentration Biases in CMIP6 Models and Their Relationships with Other Model Variables

Abstract The models that participated in the Coupled Model Intercomparison Project (CMIP) exhibit large biases in Arctic sea ice climatology that seem related to biases in seasonal atmospheric and oceanic circulations. Using historical runs of 34 CMIP6 models from 1979 to 2014, we investigate the links between the climatological sea ice concentration (SIC) biases in September and atmospheric and oceanic model climatologies. The main intermodel spread of September SIC is well described by two leading EOFs, which together explain ∼65% of its variance. The first EOF represents an underestimation or overestimation of SIC in the whole Arctic, while the second EOF describes opposite SIC biases in the Atlantic and Pacific sectors. Regression analysis indicates that the two SIC modes are closely related to departures from the multimodel mean of Arctic surface heat fluxes during summer, primarily shortwave and longwave radiation, with incoming Atlantic Water playing a role in the Atlantic sector. Local and global links with summer cloud cover, low-level humidity, upper or lower troposphere temperature/circulation, and oceanic variables are also found. As illustrated for three climate models, the local relationships with the SIC biases are mostly similar in the Arctic across the models but show varying degrees of Atlantic inflow influence. On a global scale, a strong influence of the summer atmospheric circulation on September SIC is suggested for one of the three models, while the atmospheric influence is primarily via thermodynamics in the other two. Clear links to the North Atlantic oceanic circulation are seen in one of the models.

Scaling of the emission of low frequency radiation by transient surface currents in laser–solid interactions

Abstract It has been suggested that ‘transient surface currents’ caused by multi-MeV fast electrons can be responsible for the emission of low frequency radiation (e.g. in THz range) from ultra-intense laser–solid interactions. This mechanism has been analyzed, and analytic upper bounds on the intensity, electric field amplitude, and normalized vector potential have been developed and tested against 1D EM Particle-in-Cell simulations. The ‘transient surface current’ mechanism is effective and sufficiently efficient to fully account for all radiation that has been emitted in experiments so far.

Anti-disturbance control strategy in capture stage for AUV dynamic base docking with optical guided constraints

This article investigates the challenge of potential optical guidance failure at the end of the underwater docking process due to relative attitude deviation between the Autonomous Underwater Vehicle (AUV) and the recovery platform, especially in the context of a moving base. A novel anti-disturbance control strategy is presented, accounting for environmental dynamics and optical guidance limitations. Firstly, we establish constraint models such as the effective field of view and the guidance beam associated with the optical guidance method. We then divide the recovery control strategy area based on the constraint model and propose the Line-of-sight Considering Visual Constraints (LOSCVC) docking control strategy. This strategy utilizes the line-of-sight method and the virtual leader method to make motion decisions based on the relative attitude offset and completes the docking operation by tracking the motion of the virtual leader. Furthermore, we consider the effect of ocean currents on the docking process of the moving base and design a disturbance compensation control law based on the Lyapunov redesign method. This strategy ensures the AUV maintains effective and continuous optical guidance during the capture phase, thereby achieving precise, safe, and rapid docking control. Simulation results validate the effectiveness and robustness of the proposed docking control strategy.

Future sea ice weakening amplifies wind-driven trends in surface stress and Arctic Ocean spin-up

Arctic sea ice mediates atmosphere-ocean momentum transfer, which drives upper ocean circulation. How Arctic Ocean surface stress and velocity respond to sea ice decline and changing winds under global warming is unclear. Here we show that state-of-the-art climate models consistently predict an increase in future (2015–2100) ocean surface stress in response to increased surface wind speed, declining sea ice area, and a weaker ice pack. While wind speeds increase most during fall (+2.2% per decade), surface stress rises most in winter (+5.1% per decade) being amplified by reduced internal ice stress. This is because, as sea ice concentration decreases in a warming climate, less energy is dissipated by the weaker ice pack, resulting in more momentum transfer to the ocean. The increased momentum transfer accelerates Arctic Ocean surface velocity (+31–47% by 2100), leading to elevated ocean kinetic energy and enhanced vertical mixing. The enhanced surface stress also increases the Beaufort Gyre Ekman convergence and freshwater content, impacting Arctic marine ecosystems and the downstream ocean circulation. The impacts of projected changes are profound, but different and simplified model formulations of atmosphere-ice-ocean momentum transfer introduce considerable uncertainty, highlighting the need for improved coupling in climate models.

Sea surface circulation in the Baltic Sea: decomposed components and pattern recognition

By decomposing the total sea surface current into its geostrophic and ageostrophic components, we examined the contribution of each to the long-term variability of the total sea surface current. Our findings demonstrate the importance of geostrophic currents in Baltic Sea gyre formations. Additionally, ageostrophic currents contribute significantly to the flow across the region. Quantifying the difference between total sea surface current fields has revealed two dominant general sea surface circulation patterns in the Baltic Sea, whose characteristics correspond to the monthly mean climatology of sea surface current fields in May and December. Subsequently, a machine learning technique was employed to effectively detect the types of sea surface circulation patterns using wind vectors and sea level anomaly fields. This underscored the combined influence of sea level anomaly-driven and wind-driven components in shaping surface current vectors in the Baltic Sea, consistent with geostrophic and ageostrophic decompositions.

Large nitrogen cycle perturbations during the Early Triassic hyperthermal

Ocean temperature, redox state, circulation, and nutrient levels regulate the marine nitrogen (N) cycle, yet their specific impacts during greenhouse intervals remain poorly understood. Here, we examined the Smithian–Spathian hyperthermal event (∼250.5 Ma) during the Early Triassic greenhouse using stable N isotopes (δ15N) from sedimentary records in the Nanpanjiang Basin and Northern Yangtze Basin of South China. The δ15N profiles in both basins reveal consistent trends that correspond to fluctuations in temperature across the Smithian–Spathian transition. The Smithian hyperthermal interval exhibited low δ15N values (mostly <+2‰), indicating N deficiency and enhanced biological N2 fixation. Bacterial blooms and the release of the potent greenhouse gas N2O, enhanced by high temperatures, may have triggered positive feedback mechanisms that sustained the warming and contributed to the late Smithian extinction. During subsequent cooling across the Smithian–Spathian transition, δ15N increased to a range of +3‰ to +7‰, likely reflecting signals of partial denitrification based on reconstruction of the NO3− inventory associated with oceanic cooling and oxygenation. The prevailing increases in sedimentary TOC/TN ratios signify heightened deamination (removal of amino groups) and N recycling across the Smithian–Spathian boundary. This transition is likely attributed to cooling-driven amelioration of seawater stratification and anoxia from the Smithian to the Spathian, which resulted in increased nitrate availability in the photic zone. Eukaryotic algae thrived while prokaryotes declined during the Spathian, evidenced by elevated δ15N and Δ13Ccarb-org (the difference between carbonate and organic carbon isotopes). The proliferation of eukaryotic algae had a positive impact on environmental conditions and facilitated biotic recovery due to more efficient burial of organic particles. A notable latitudinal gradient in the N cycle response was observed, with low-latitude regions showing a more pronounced response to cooling compared to mid-latitude areas. This significant gradient may suggest that the rapid recovery of nutrient cycles, despite relatively small decline in high temperature, was a key factor in the amelioration of climate and recovery of life in low latitudes. These findings highlight that the marine N cycle is highly sensitive to temperature changes, particularly in low-latitude regions, and that changes in N cycle under high-temperature conditions may be a critical life-limiting factor.

Coupling of potential habitat models with particle tracking experiments to examine larval fish dispersal and connectivity in deep water regions.

Computing Lagrangian trajectories with ocean circulation models is a powerful way to infer larval dispersal pathways and connectivity. Defining release areas and timing of particles to represent larval habitat realistically is critical to obtaining representative dispersal pathways. However, it is challenging due to spatial and temporal variability in larval density. Forward-tracking particle experiments were conducted to study larval connectivity of four species (neritic or mesopelagic) in the Gulf of Mexico's (GoM) deep-water region. A seasonal climatology coupled with predicted potential larval habitat models based on generalized additive models was used to delimit the particle dispersal origin. Two contrasting mesoscale circulation patterns were examined: (1) high Loop Current (LC) intrusion, absence of recently detached LC anticyclonic eddies (LC-ACE), and no interaction between LC-ACEs and the semi-permanent cyclonic eddy (CE) in the Bay of Campeche (BoC), and (2) limited LC intrusion, a recently detached LC-ACE, and interaction between LC-ACEs and the BoC's CE. To simulate larval transport, virtual larvae were randomly released in the potential habitats and advected for 30 days with the velocity fields of the HYbrid Coordinate Ocean Model with hourly-resolution assimilation. Potential habitat location and size played a major role in dispersal and connectivity. A greater percentage of particles were retained in potential habitats restricted to the southern BoC, suggesting lower connectivity with other GoM regions than those encompassing most of the BoC or the central Gulf. Mesoscale feature interactions in the western GoM and BoC led to greater dispersion along the western basin. By contrast, the absence of ACE-CE interaction in the BoC led to greater retention and less connectivity between the southern and northern GoM. Under high LC intrusion, particles seeded north of the Yucatan Shelf were advected through the Florida Straits and dispersed within the GoM. Coupling potential habitat models with particle experiments can help characterize the dispersal and connectivity of fish larvae in oceanic systems.

Multi-AUVs cooperative path planning in 3D underwater terrain and vortex environments based on improved multi-objective particle swarm optimization algorithm

This paper presents an improved multi-objective particle swarm optimization (IMOPSO) algorithm for cooperative path planning of multiple autonomous underwater vehicles (multi-AUVs) in 3D underwater terrain and vortex environment. Traditional PSO algorithms often struggle with local optima in high-dimensional, complex search spaces. To address this, we propose enhancements including Sine double-sequence perturbation interpolation for initial particle population generation, adaptive differential evolution strategies, iterative local search, and traction operations to enhance global search capability and avoid local optima. Additionally, cubic spline interpolation is employed to smooth the planned paths. The objective function is tailored to the ocean currents environment, considering path length, smoothness, deflection angle, and current impact, with a penalty function method to handle multi-AUVs cooperation and environmental constraints optimization for optimal energy consumption. Simulation results demonstrate that the proposed algorithm effectively plans short-range, high-security global paths for multi-AUVs, verifying its effectiveness and practicality in meeting the navigation requirements in ocean currents environments.