Fujiwhara Effect: How Binary Cyclones Challenge IMD Forecasts

Context

The Fujiwhara effect occurs when two tropical cyclones approach within approximately 1,400 kilometres of each other. They initiate a mutual orbital dance around a shared centre of mass known as a barycentre. This binary interaction causes erratic track deviations, rapid intensity fluctuations, and potential vortex mergers that disrupt standard single-vortex numerical weather prediction models.

In the North Indian Ocean basin, simultaneous storm systems across the Arabian Sea and the Bay of Bengal can exert non-linear teleconnections. Upper-tropospheric outflow from one vortex distorts the ambient steering flow of the other. As rising sea surface temperatures expand the window for concurrent cyclogenesis, operational agencies like the India Meteorological Department face compounding forecasting challenges that directly cut into coastal evacuation lead times.

The India Meteorological Department closely monitors binary cyclonic systems due to forecasting anomalies during simultaneous cyclogenesis across the North Indian Ocean basin. In October 2023, the basin generated concurrent twin disturbances: Extremely Severe Cyclonic Storm Tej over the Arabian Sea and Very Severe Cyclonic Storm Hamoon over the Bay of Bengal.

According to reports from the India Meteorological Department, upper-tropospheric outflow from Cyclone Tej altered the steering and wind shear environment around Cyclone Hamoon. This dynamical interaction created two operational challenges:

• Abrupt Track Recurvature: Cyclone Hamoon veered sharply northeastwards toward Bangladesh instead of following the climatologically expected west-northwest vector toward the Indian coast.
• Rapid Category Shifts: Cyclone Hamoon changed intensity categories four times within 24 hours prior to landfall, rapidly escalating to a Very Severe Cyclonic Storm before weakening prior to crossing the coast.

These sudden operational shifts highlight how binary atmospheric interactions strain traditional early warning timelines.

The Fujiwhara effect describes the physical process where two sufficiently close cyclonic vortices orbit around a common barycentre through mutual advective steering. Japanese meteorologist Sakuhei Fujiwhara first documented this phenomenon in 1921 through water tank experiments showing that co-rotating fluid vortices naturally attract and rotate about a shared centroid.

Under guidelines established by the World Meteorological Organization, binary tropical cyclone interaction initiates when storm centres fall within approximately 1,400 kilometres, or roughly 10 to 14 degrees of latitude. When crossing this threshold, outer circulating wind fields begin steering each vortex, producing a cyclonic mutual rotation.

Planetary rotation dictates the orbital direction through the Coriolis force. Interacting systems rotate counter-clockwise around their shared barycenter in the Northern Hemisphere and clockwise in the Southern Hemisphere. The interaction strength depends on the relative size, core intensity, and initial separation distance of the two tropical cyclones.

Binary cyclonic interactions resolve into five structurally distinct operational outcomes depending on storm size, intensity disparity, and core separation distance. Atmospheric scientists L. E. Carr and R. L. Elsberry classified these interactions to objectively diagnose structural evolutions between two tropical systems:

1. Elastic Interaction: The two cyclones approach the interaction threshold and execute mutual cyclonic rotation around their barycenter. Both systems alter their tracks via mutual cyclonic steering but retain their core structures and eventually separate without merging.
2. Partial Straining-Out: When one cyclone possesses significantly greater intensity or size, horizontal wind shear from the dominant storm strips away outer convective rainbands from the smaller partner. The weaker system suffers structural degradation but maintains an identifiable vortex core.
3. Complete Straining-Out (One-Sided Ejection): The shearing forces exerted by the dominant system completely tear apart the weaker companion. The smaller cyclone is dismantled into an elongated vorticity filament or sheared band, functionally destroying its organized core.
4. Partial Merger: The two cyclonic circulations approach close enough for convective mixing, where outer cores blend while secondary vorticity maxima temporarily persist before being absorbed.
5. Complete Merger: Two tropical cyclones approach within a critical separation distance, typically under 750 to 800 kilometres. Their central cores coalesce directly into a single, broader cyclonic circulation.

The North Indian Ocean basin rarely witnesses direct core-merger events because the Indian subcontinent physically bifurcates the basin into the Arabian Sea and the Bay of Bengal. However, simultaneous cyclogenesis and dynamical teleconnections occur regularly across these separate sub-basins.

Two regional historical events illustrate how these dynamics unfold in practice:

• Inter-Basin Teleconnections (North Indian Ocean): The India Meteorological Department recorded simultaneous twin cyclogenesis in October 2018 when Very Severe Cyclonic Storm Titli formed in the Bay of Bengal while Very Severe Cyclonic Storm Luban developed concurrently over the Arabian Sea. While the landmass prevented physical core coalescence, large-scale moisture convergence and upper-tropospheric outflow created subtle inter-basin atmospheric linkages.
• Direct Open-Ocean Merger (South Indian Ocean): In April 2021, Severe Tropical Cyclone Seroja engaged in a pronounced Fujiwhara interaction with Tropical Cyclone Odette off Western Australia. The two storms executed a rapid mutual cyclonic orbit that culminated in Cyclone Seroja absorbing Cyclone Odette before making landfall.

Predictive skill drops noticeably during binary events because standard numerical weather prediction frameworks struggle with multi-vortex non-linear dynamics. Research published by the American Meteorological Society demonstrates that conventional single-vortex tracking models fail to account for non-linear horizontal advection and diabatic heating asymmetries induced by a companion storm.

The India Meteorological Department maintains strict verification metrics for cyclone tracking. For the five-year period from 2019 to 2023, average operational track forecast errors across the North Indian Ocean stood at 72 kilometres at 24 hours, 112 kilometres at 48 hours, and 156 kilometres at 72 hours.

During 2023, the operational Multi-Model Ensemble track forecast errors across all North Indian Ocean storms registered at 59 kilometres at 24 hours, 101 kilometres at 48 hours, 168 kilometres at 72 hours, 178 kilometres at 96 hours, and 199 kilometres at 120 hours.

Operational verification findings from the IMD confirmed that disturbances undergoing sharp track recurvature and rapid intensity shifts—such as Cyclone Tej—produced significantly elevated forecast errors compared to straight-moving baseline storms. These deviations impair coastal preparedness in three primary ways:

• Distorted Steering Currents: Outflow from one cyclone introduces unexpected vertical wind shear over the other, displacing the convective core.
• Rapid Intensity Volatility: The companion vortex can suppress inflow channels or abruptly inject moisture, causing unpredicted intensification or rapid decay.
• Landfall Dislocation: Sudden orbital loops create misplaced coastal evacuations, deploying emergency resources to administrative districts that avoid the ultimate storm strike.

Rising ocean temperatures alter the thermodynamic baseline of the North Indian Ocean, creating conditions conducive to simultaneous cyclogenesis. According to assessment reports by the Ministry of Earth Sciences, sea surface temperatures across the Arabian Sea and equatorial Indian Ocean have warmed at rates exceeding 0.12°C to 0.15°C per decade.

Two interdependent factors drive this increasing frequency of dual systems:

• Elevated Ocean Heat Content: Sustained warming elevates the Tropical Cyclone Heat Potential, which measures thermal energy stored in the upper ocean down to the 26°C isotherm. Higher heat potential ensures that even when a leading cyclonic system upwells cooler sub-surface water, surrounding waters remain warm enough to sustain a second concurrent vortex.
• Coaligned Dynamic Triggers: Technical studies by the National Centre for Medium Range Weather Forecasting show that Equatorial Rossby waves and the active convective phase of the Madden-Julian Oscillation deliver the planetary vorticity and moisture convergence necessary for twin storms to develop simultaneously across the Arabian Sea and Bay of Bengal.

Managing binary cyclonic threats requires upgrading deterministic forecasting into dynamic ensemble frameworks backed by high-resolution observation networks. To capture track divergence in multi-vortex scenarios, the National Centre for Medium Range Weather Forecasting operates the NCMRWF Unified Model global ensemble prediction system (NEPS-G) at 12-kilometre horizontal resolution and a dedicated regional system (NEPS-R) at 4-kilometre resolution.

High-resolution regional grids resolve localized wind shears that govern whether twin vortices merge or repel. Complementing these ensemble predictions, the India Meteorological Department operationally executes the Hurricane Weather Research and Forecasting model. This system utilises movable vortex-following nests updated four times daily to project intensity and tracking pathways across 120-hour lead periods.

Observational platforms supply the baseline data required to initialise these numerical equations:

• Ocean Moored Buoy Networks: In-situ data assimilation from the Ocean Moored Buoy Network (OMNI) deployed by the National Institute of Ocean Technology and RAMA moored buoys managed by the Indian National Centre for Ocean Information Services captures real-time upper-ocean thermal profiles.
• Space-Based Scatterometry: Satellite scatterometer wind retrievals from Oceansat-3 (SCAT-3 Ku-band) allow forecasters to accurately identify distinct surface circulation centres and map asymmetric wind radii during complex binary events.

Combining high-resolution numerical models with real-time ocean and satellite observations provides disaster managers with the predictive reliability needed to coordinate safe evacuations along vulnerable coastlines.

• Japanese meteorologist Sakuhei Fujiwhara first identified the Fujiwhara effect in 1921, demonstrating how co-rotating fluid vortices orbit around a shared centroid.
• Binary tropical cyclone interaction typically initiates when two storms approach within approximately 1,400 kilometres (10 to 14 degrees of latitude) of each other.
• The Coriolis force causes interacting cyclones to rotate counter-clockwise in the Northern Hemisphere and clockwise in the Southern Hemisphere around their common barycenter.
• Binary interactions produce five structural outcomes: elastic interaction, partial straining-out, complete straining-out, partial merger, and complete merger (which occurs under 750 to 800 kilometres separation).
• Simultaneous storms like Cyclone Tej and Cyclone Hamoon in October 2023 illustrate how inter-basin teleconnections induce erratic track recurvatures and rapid intensity category changes.
• Operational forecasting relies on the NCMRWF Unified Model ensemble prediction systems (NEPS-G at 12 km and NEPS-R at 4 km) along with Oceansat-3 scatterometer data to resolve multi-vortex tracks.