Wind Stress Anomalies Research Articles

Role of Reemergence in the Central North Pacific Revealed by a Mixed Layer Heat Budget Analysis

AbstractWintertime sea surface temperature (SST) anomalies that disappear in summer and recur in the following winter are known as “reemergence.” The role of reemergence in recurrent SST anomalies in an area in the central North Pacific is investigated quantitatively for the first time by conducting an online mixed layer heat budget analysis with a realistic ocean model simulation. In contrast to past studies that have emphasized the importance of vertical entrainment of subsurface temperature anomalies, it is shown that several mechanisms are operating in recurrence of SST anomalies during 1977–1999. In particular, anomalous Ekman meridional advection of the mean meridional temperature gradient induced by zonal wind stress anomalies plays an important role in many years. Also, coincidence in the sign of SST anomalies in winter and surface heat flux anomalies in the following autumn leads to recurrence of SST anomalies.

On the Generation and Salinity Impacts of Intraseasonal Westward Jets in the Equatorial Indian Ocean

AbstractWhile westerly winds dominate the equatorial Indian Ocean and generate the well‐known eastward flowing Wyrtki Jets during boreal spring and fall, there is evidence of a strong westward surface jet during winter that is swifter than eastward currents during that season. A weaker westward jet is found in summer. In this study, we report the occurrence, characteristics, and intraseasonal variability of this westward jet and its impact on mixed layer salinity in the equatorial Indian Ocean using the HYbrid Coordinate Ocean Model (HYCOM) reanalysis with the Navy Coupled Ocean Data Assimilation (NCODA). The westward jet typically occurs in the upper 50 m, above an eastward flowing equatorial undercurrent, with peak westward volume transport of approximately −8 Sv. The westward jet builds up gradually, decays rapidly, and is primarily forced by local intraseasonal wind stress anomalies generated by atmospheric intraseasonal convection. Westward acceleration of the jet occurs when the dominant intraseasonal westward wind anomaly is not balanced by the zonal pressure gradient (ZPG) force. The intraseasonal westward jet generates strong horizontal advection and is the leading cause of mixed layer freshening in the western equatorial Indian Ocean. Without it, a saltier mixed layer would persist and weaken any barrier layers. Existing barrier layers are strengthened following the passage of freshwater‐laden westward jets. Deceleration of the westward jet occurs when the eastward ZPG becomes increasingly important and the westward intraseasonal wind anomalies weaken. A rapid reversal of atmospheric intraseasonal convection‐driven surface winds eventually terminates the westward jet.

Can Tropical Pacific Winds Enhance the Footprint of the Interdecadal Pacific Oscillation on the Upper-Ocean Heat Content in the South China Sea?

AbstractIn this study, an enhanced footprint of the interdecadal Pacific oscillation (IPO) on the upper-ocean heat content (OHC) in the South China Sea (SCS) since the 1990s is revealed. The negative OHC–IPO correlation is significant (r= −0.71) during 1990–2010 [period 2 (P2)], whereas it is statistically insignificant during 1960–80 [period 1 (P1)]. Analyses show that the scope of the equatorial Pacific wind anomalies is wider during P2 compared with that during P1 due to a larger east–west SST gradient and enhanced tropical warming in the Indian Ocean. When the IPO is negative during P2, a wider scope of the wind stress anomalies associated with the IPO could lead to 1) the southward migration of the North Equatorial Current bifurcation latitude (NECBL) by affecting the wind stress curl over the key region where it is near the climatological NECBL and 2) an increase in the interbasin pressure gradient (sea surface height difference) between the western Pacific and the SCS; these two processes strengthen the Kuroshio and weaken the Luzon Strait transport (LST) or SCS throughflow into the SCS. Also, 3) the equatorial Pacific wind anomalies are wide enough to directly weaken the LST in the SCS through the “island rule.” These three pathways finally change the oceanic gyre in the SCS and increase the OHC. Our results suggest that the scope of the tropical wind stress is the crucial factor when we consider the relationship between the upper ocean thermal conditions in the SCS and the Pacific variability.

Phase synchronisation in the Kuroshio Current System

Abstract. The Kuroshio Current System in the North Pacific displays path transitions on a decadal timescale. It is known that both internal variability involving barotropic and baroclinic instabilities and remote Rossby waves induced by North Pacific wind stress anomalies are involved in these path transitions. However, the precise coupling of both processes and its consequences for the dominant decadal transition timescale are still under discussion. Here, we analyse the output of a multi-centennial high-resolution global climate model simulation and study phase synchronisation between Pacific zonal wind stress anomalies and Kuroshio Current System path variability. We apply the Hilbert transform technique to determine the phase and find epochs where such phase synchronisation appears. The physics of this synchronisation are shown to occur through the effect of the vertical motion of isopycnals, as induced by the propagating Rossby waves, on the instabilities of the Kuroshio Current System.

The influence of wintertime SST variability in the Western North Pacific on ENSO diversity

Sea surface temperature anomalies (SSTa) in the Western North Pacific (WNP) have been linked to the development of El Niño Southern Oscillation (ENSO) events a full year in advance. However, the contribution of the WNP precursor to the temporal evolution and spatial complexity of ENSO remains unclear. Using the preindustrial experiment of the Community Earth System Model as the control climate, a partially coupled experiment is conducted in which WNP SSTa are restored to the model climatology. By comparing the perturbed experiment to the control, we are able to clearly characterize ENSO’s response to WNP SST variability. We find that the WNP is predominantly linked to eastern Pacific (EP) ENSO events. Without SST variability in the WNP, central Pacific (CP) ENSO events are more likely to develop. This variation in ENSO flavor is controlled by how the WNP projects onto wind stress anomalies in the western equatorial Pacific, which in turn impacts the discharge and recharge of ocean heat during the ENSO cycle. Specifically, the removal of SSTa in the WNP weakens the buildup of ocean heat in the western equatorial Pacific, which then hinders the development of EP-type events.

Interannual Variability of the Global Meridional Overturning Circulation Dominated by Pacific Variability

AbstractThe most prominent feature of the time-mean global meridional overturning circulation (MOC) is the Atlantic MOC (AMOC). However, interannual variability of the global MOC is shown here to be dominated by Pacific MOC (PMOC) variability over the full depth of the ocean at most latitudes. This dominance of interannual PMOC variability is robust across modern climate models and an observational state estimate. PMOC interannual variability has large-scale organization, its most prominent feature being a cross-equatorial cell spanning the tropics. Idealized experiments show that this variability is almost entirely wind driven. Interannual anomalies of zonal mean zonal wind stress produce zonally integrated Ekman transport anomalies that are larger in the Pacific Ocean than in the Atlantic Ocean, simply because the Pacific is wider than the Atlantic at most latitudes. This contrast in Ekman transport variability implies greater variability in the near-surface branch of the PMOC when compared with the near-surface branch of the AMOC. These near-surface variations in turn drive compensating flow anomalies below the Ekman layer. Because the baroclinic adjustment time is longer than a year at most latitudes, these compensating flow anomalies have baroclinic structure spanning the full depth of the ocean. Additional analysis reveals that interannual PMOC variations are the dominant contribution to interannual variations of the global meridional heat transport. There is also evidence of interaction between interannual PMOC variability and El Niño–Southern Oscillation.

Wind-Forced Variability of the Remote Meridional Overturning Circulation

AbstractThe mechanisms by which time-dependent wind stress anomalies at midlatitudes can force variability in the meridional overturning circulation at low latitudes are explored. It is shown that winds are effective at forcing remote variability in the overturning circulation when forcing periods are near the midlatitude baroclinic Rossby wave basin-crossing time. Remote overturning is required by an imbalance in the midlatitude mass storage and release resulting from the dependence of the Rossby wave phase speed on latitude. A heuristic theory is developed that predicts the strength and frequency dependence of the remote overturning well when compared to a two-layer numerical model. The theory indicates that the variable overturning strength, relative to the anomalous Ekman transport, depends primarily on the ratio of the meridional spatial scale of the anomalous wind stress curl to its latitude. For strongly forced systems, a mean deep western boundary current can also significantly enhance the overturning variability at all latitudes. For sufficiently large thermocline displacements, the deep western boundary current alternates between interior and near-boundary pathways in response to fluctuations in the wind, leading to large anomalies in the volume of North Atlantic Deep Water stored at midlatitudes and in the downstream deep western boundary current transport.

The role of the South Pacific in modulating Tropical Pacific variability

Tropical Pacific variability (TPV) heavily influences global climate, but much is still unknown about its drivers. We examine the impact of South Pacific variability on the modes of TPV: the El Niño-Southern Oscillation (ENSO) and the Interdecadal Pacific Oscillation (IPO). We conduct idealised coupled experiments in which we suppress temperature and salinity variability at all oceanic levels in the South Pacific. This reduces decadal variability in the equatorial Pacific by ~30% and distorts the spatial pattern of the IPO. There is little change to overall interannual variability, however there is a decrease in the magnitude of the largest 5% of both El Niño and La Niña sea-surface temperature (SST) anomalies. Possible reasons for this include: (i) reduced decadal variability means that interannual SST variability is superposed onto a ‘flatter’ background signal, (ii) suppressing South Pacific variability leads to the alteration of coupled processes linking the South and equatorial Pacific. A small but significant mean state change arising from the imposed suppression may also contribute to the weakened extreme ENSO SST anomalies. The magnitude of both extreme El Niño and La Niña SST anomalies are reduced, and the associated spatial patterns of change of upper ocean heat content and wind stress anomalies are markedly different for both types of events.

Heat Distribution in the Southeast Pacific Is Only Weakly Sensitive to High‐Latitude Heat Flux and Wind Stress

AbstractThe Southern Ocean features regionally varying ventilation pathways that transport heat and carbon from the surface ocean to the interior thermocline on timescales of decades to centuries, but the factors that control the distribution of heat along these pathways are not well understood. In this study, we use a global ocean state estimate (ECCOv4) to (1) define the recently ventilated interior Pacific (RVP) using numerical passive tracer experiments over a 10‐year period and (2) use an adjoint approach to calculate the sensitivities of the RVP heat content (RVPh) to changes in net heat flux and wind stress. We find that RVPh is most sensitive to local heat flux and wind stress anomalies north of the sea surface height contours that delineate the Antarctic Circumpolar Current, with especially high sensitivities over the South Pacific Gyre. Surprisingly, RVPh is not especially sensitive to changes at higher latitudes. We perform a set of step response experiments over the South Pacific Gyre, the subduction region, and the high‐latitude Southern Ocean. In consistency with the adjoint sensitivity fields, RVPh is most sensitive to wind stress curl over the subtropical gyre, which alter isopycnal heave, and it is only weakly sensitive to changes at higher latitudes. Our results suggest that despite the localized nature of mode water subduction hot spots, changes in basin‐scale pressure gradients are an important controlling factor on RVPh. Because basin‐scale wind stress is expected to change in the coming decades to centuries, our results may have implications for climate, via the atmosphere/ocean partitioning of heat.

The Influence of Seasonal and Interannual Variability on Surface Chlorophyll-a Off the Western Lesser Sunda Islands

Novel and long-term satellite data (2003–2017) are analyzed to investigate the variability of ocean surface chlorophyll-a (Chl-a) concentration off the Western Lesser Sunda Islands (WLSI) under the influence of the Indonesian Australian monsoon, the El Nino-Southern Oscillation, and the Indian Ocean Dipole (IOD). In this article, we first analyzed the seasonal variability of Chl-a, and then describe the relationship among the sea surface Chl-a, sea surface temperature (SST), and sea surface wind stress in the region. Our results demonstrate that prevailing southeasterly wind stress plays a pivotal role in generating the Chl-a maxima off the WLSI. Particularly on seasonal time scale, the strengthening of southeasterly wind stress (up to ∼0.01 N·m−2) during the southeast monsoon season produces enhanced Chl-a concentrations (0.59 mg·m−3) associated with sea surface cooling (∼28.8 °C) in the area of study. In contrast, the Chl-a maxima completely vanished during the northwest monsoon season. On interannual time scale, the largest positive Chl-a and wind stress anomalies and the coolest SST anomaly are observed in 2006 when El Nino and positive IOD events occur at the same time. Meanwhile, the greatest negative Chl-a anomaly is prevailed during the 2016 negative IOD event. This article demonstrates that wind variability is the essential factor in determining the magnitude of the Chl-a maxima off the WLSI.

Distinctive characteristics of upwelling along the Peninsular Malaysia's east coast during 2009/10 and 2015/16 El Niños

We investigate the southwest monsoon average structure and interannual variability of upwelling along the Peninsular Malaysia's east coast (PMEC) based on satellite-derived sea surface temperature (SST), reanalysis wind data and cruise observations. The upwelling area prompted by southwest monsoon winds is represented by a patch of cooler water along the coast. This cooler water is also an outcome of the cooler water flooded from the Karimata Strait, nonetheless, localised upwelling conditions demonstrating their substantial impact to further cool down this already cooled water. Results from empirical orthogonal function (EOF) analysis revealed that variation in SST was dominated by the first EOF mode, which accountable of 80.99% of total variance and its loading map demonstrated a persistence of SST front, represented the signature of upwelling. The time series of the first EOF mode was highly correlated with the alongshore wind stress anomalies, indicating the Ekman transport driven nature of upwelling interannual variability. We demonstrate that a delayed El Niño Southern Oscillation (ENSO) effect as an important factor for the interannual variations of upwelling along the PMEC. Our analysis focusing on the two recent strongest El Niños, 2009/10 and 2015/16 revealed that the intense northwest Pacific anticyclone anomaly (NWPAA) and the occurrence of the northerly wind anomalies from the subtropical western Pacific (SWP) in 2010 weakened the upwelling along the PMEC. While in 2016, the occurrence of the ‘C-shape’ wind anomalies generated by the warming of the southeast tropical Indian Ocean (SETIO) weakened the upwelling along the PMEC.

Influence of Intraseasonal Oscillation on the Asymmetric Decays of El Niño and La Niña

Warm and cold phases of El Nino–Southern Oscillation (ENSO) exhibit a significant asymmetry in their decay speed. To explore the physical mechanism responsible for this asymmetric decay speed, the asymmetric features of anomalous sea surface temperature (SST) and atmospheric circulation over the tropical Western Pacific (WP) in El Nino and La Nina mature-to-decay phases are analyzed. It is found that the interannual standard deviations of outgoing longwave radiation and 850 hPa zonal wind anomalies over the equatorial WP during El Nino (La Nina) mature-to-decay phases are much stronger (weaker) than the intraseasonal standard deviations. It seems that the weakened (enhanced) intraseasonal oscillation during El Nino (La Nina) tends to favor a stronger (weaker) interannual variation of the atmospheric wind, resulting in asymmetric equatorial WP zonal wind anomalies in El Nino and La Nina decay phases. Numerical experiments demonstrate that such asymmetric zonal wind stress anomalies during El Nino and La Nina decay phases can lead to an asymmetric decay speed of SST anomalies in the central-eastern equatorial Pacific through stimulating different equatorial Kelvin waves. The largest negative anomaly over the Nino3 region caused by the zonal wind stress anomalies during El Nino can be threefold greater than the positive Nino3 SSTA anomalies during La Nina, indicating that the stronger zonal wind stress anomalies over the equatorial WP play an important role in the faster decay speed during El Nino.

Interannual Variability of the North Pacific Mixed Layer Associated with the Spring Tibetan Plateau Thermal Forcing

AbstractBy using multisourced data and two sets of sensitivity runs from the coupled general circulation model CESM1.2.0, we investigated the effects of the spring [March, April, and May (MAM)] surface sensible heating over the Tibetan Plateau (SHTP) on the interannual variability of the North Pacific Ocean sea surface temperature (SST) and mixed layer. The results indicated that an above-normal MAM SHTP can generate a Rossby wave downstream and form an anomalous equivalent barotropic anticyclone over the North Pacific, inducing anticyclonic wind stress anomalies. As a result of Ekman transport and Ekman pumping, sea currents converge near 40°N, accompanied by weak downwelling motion. The mixed layer heat budget diagnosis indicates that the net heat fluxes, together with meridional advection anomalies, contributed significantly to changes in the mixed layer temperature (MLT). As a result, the SST anomalies (SSTAs) and MLT anomalies both present a horseshoelike pattern. In addition, the significant warm SSTAs show a maximum in the late spring, but the significant warm MLT anomalies centered under the sea surface (25-m depth) could be sustained until summer, acting like a signal storage for the anomalous spring SHTP. Moreover, the midlatitude ocean–atmosphere interaction provides a positive feedback on the development of the anomalous anticyclone over the North Pacific, since the SSTA pattern could strengthen the oceanic front and induce more active transient eddy activities. The eddy vorticity forcing that is dominant among the total atmospheric forcings tends to produce an equivalent barotropic atmospheric high pressure, which in turn intensifies the initial anomalous anticyclone.

Contributions of Atmosphere–Ocean Interaction and Low-Frequency Variation to Intensity of Strong El Niño Events since 1979

Evolutions of oceanic and atmospheric anomalies in the equatorial Pacific during four strong El Niños (1982/83, 1991/92, 1997/98, and 2015/16) since 1979 are compared. The contributions of the atmosphere–ocean coupling to El Niño–associated sea surface temperature anomalies (SSTA) are identified and their association with low-level winds as well as different time-scale variations is examined. Although overall SSTA in the central and eastern equatorial Pacific is strongest and comparable in the 1997/98 and 2015/16 El Niños, the associated subsurface ocean temperature as well as deep convection and surface wind stress anomalies in the central and eastern equatorial Pacific are weaker during 2015/16 than that during 1997/98. That may be associated with a variation of the wind–SST and wind–thermocline interactions. Both the wind–SST and wind–thermocline interactions play a less important role during 2015/16 than during 1997/98. Such differences are associated with the differences of the low-level westerly wind as well as the contribution of different time-scale variations in different events. Similar to the interannual time-scale variation, the intraseasonal–interseasonal time-scale component always has positive contributions to the intensity of all four strong El Niños. Interestingly, the role of the interdecadal-trend time-scale component varies with event. The contribution is negligible during the 1982/83 El Niño, negative during the 1997/98 El Niño, and positive during the 1991/92 and 2015/16 El Niños. Thus, in addition to the atmosphere–ocean coupling at intraseasonal to interannual time scales, interdecadal and longer time-scale variations may play an important and sometimes crucial role in determining the intensity of El Niño.

The Effect of Wind Stress Anomalies and Location in Driving Pacific Subtropical Cells and Tropical Climate

Abstract The importance of subtropical and extratropical zonal wind stress anomalies on Pacific subtropical cell (STC) strength is assessed through several idealized and realistic numerical experiments with a global ocean model. Different zonal wind stress anomalies are employed, and their intensity is strengthened or weakened with respect to the climatological value throughout a suite of simulations. Subtropical strengthened (weakened) zonal wind stress anomalies result in increased (decreased) STC meridional mass and energy transport. When upwelling of subsurface water into the tropics is intensified (reduced), a distinct cold (warm) anomaly appears in the equatorial thermocline and up to the surface, resulting in significant tropical sea surface temperature (SST) anomalies. The use of realistic wind stress anomalies also suggests a potential impact of midlatitude atmospheric modes of variability on tropical climate through STC dynamics. The remotely driven response is compared with a set of simulations where an equatorial zonal wind stress anomaly is imposed. A dynamically distinct response is achieved, whereby the equatorial thermocline adjusts to the wind stress anomaly, resulting in significant equatorial SST anomalies as in the remotely forced simulations but with no role for STCs. Significant anomalies in Indonesian Throughflow transport are generated only when equatorial wind stress anomalies are applied, leading to remarkable heat content anomalies in the Indian Ocean. Equatorial wind stress anomalies do not involve modifications of STC transport but could set up the appropriate initial conditions for a tropical–extratropical teleconnection involving Hadley cells, exciting an STC anomalous transport, which ultimately feeds back on the tropics.

The roles of tropical and subtropical wind stress anomalies in the El Ni\xf1o Modoki onset

The El Niño Modoki is separated into two types, El Niño Modoki I and El Niño Modoki II, which are characterized by substantially different mechanism. According to comparisons between the evolution of air-sea coupled processes and heat budget analyses, this study discovered the distinct dynamical mechanism of the onset of El Niño Modoki II compared to that of El Niño Modoki I. For El Niño Modoki I, westerly wind anomalies in the western Pacific appear earlier than warm sea surface temperature (SST) anomalies in the central tropical Pacific; in addition, zonal advection is responsible for changes in mixed-layer temperature in the Niño-4 region during the developing phase in spring and summer. However, warm SST anomalies accompany easterly anomalies in the central tropical Pacific, and there are no westerly anomalies in the western tropical Pacific during the developing phase in spring and summer for El Niño Modoki II. For El Niño Modoki II, Ekman feedback is the major contributor to changes in mixed-layer temperature. According to reanalysis data, the 2.5-layer model results and coupled model simulations, anomalous westerly wind stress in the subtropical North Pacific and easterly wind stress in the central–eastern tropical Pacific can lead to equatorward and westward oceanic current anomalies, thus enhancing meridional and zonal convergences at upper layer. Because anomalous convergence at upper layer can inhibit upwelling motion in the central tropical Pacific, SST anomalies in the central tropical Pacific tend to be warm, even without anomalous westerlies in the western Pacific.

The Observed Impacts of the Two Types of El Niño on the North Equatorial Countercurrent in the Pacific Ocean

AbstractThis study investigates the interannual variations in the North Equatorial Counter Current (NECC) associated with the eastern Pacific and the central Pacific types of El Niño. Using observational analysis and ocean simulations, we show that the wind stress anomalies during the two El Niño types are of comparable amplitude but have different spatial structures, which results in significant and distinct variations in the NECC. The NECC shifts southward and intensifies during the developing phase of El Niño, but the variations are confined in the central eastern Pacific for the eastern Pacific type and the western central Pacific for the central Pacific type. These differences can be attributed to modulations in equatorial Kelvin wave and tropical Rossby wave propagation as well as Ekman pumping.

The Importance of a Properly Represented Stratosphere for Northern Hemisphere Surface Variability in the Atmosphere and the Ocean

Major sudden stratospheric warmings (SSWs) are extreme events during boreal winter, which not only impact tropospheric weather up to three months but also can influence oceanic variability through wind stress and heat flux anomalies. In the North Atlantic region, SSWs have the potential to modulate deep convection in the Labrador Sea and thereby the strength of the Atlantic meridional overturning circulation. The impact of SSWs on the Northern Hemisphere surface climate is investigated in two coupled climate models: a stratosphere-resolving (high top) and a non-stratosphere-resolving (low top) model. In both configurations, a robust link between SSWs and a negative NAO is detected, which leads to shallower-than-normal North Atlantic mixed layer depth. The frequency of SSWs and the persistence of this link is better captured in the high-top model. Significant differences occur over the Pacific region, where an unrealistically persistent Aleutian low is observed in the low-top configuration. An overrepresentation of SSWs during El Niño conditions in the low-top model is the main cause for this artifact. Our results underline the importance of a proper representation of the stratosphere in a coupled climate model for a consistent surface response in both the atmosphere and the ocean, which, among others, may have implications for oceanic deep convection in the subpolar North Atlantic.

Western Pacific Oceanic Heat Content: A Better Predictor of La Niña Than of El Niño

AbstractThe western equatorial Pacific oceanic heat content (warm water volume in the west or WWVw) is the best El Niño–Southern Oscillation (ENSO) predictor beyond 1‐year lead. Using observations and selected Coupled Model Intercomparison Project Phase 5 simulations, we show that a discharged WWVw in boreal fall is a better predictor of La Niña than a recharged WWVw for El Niño 13 months later, both in terms of occurrence and amplitude. These results are robust when considering the heat content across the entire equatorial Pacific (WWV) at shorter lead times, including all Coupled Model Intercomparison Project Phase 5 models or excluding Niño‐Niña and Niña‐Niño phase transitions. Suggested mechanisms for this asymmetry include (1) the negatively skewed WWVw distribution with stronger discharges related to stronger wind stress anomalies during El Niño and (2) the stronger positive Bjerknes feedback loop during El Niño. The possible role of stronger subseasonal wind variations during El Niño is also discussed.

A New Hybrid Coupled Model of Atmosphere, Ocean Physics, and Ocean Biogeochemistry to Represent Biogeophysical Feedback Effects in the Tropical Pacific

AbstractMultiple processes are involved in modulating the El Niño‐Southern Oscillation (ENSO) in the tropical Pacific, and these processes are neither well‐represented nor well‐understood in climate models. A new hybrid coupled model (HCM) of atmosphere, ocean physics, and ocean biogeochemistry (AOPB) is developed to represent the feedback from ocean biogeochemistry onto ocean physics via modulating the penetration of shortwave radiation in the upper ocean. An ocean biogeochemistry model is coupled with a simplified ocean‐atmosphere system consisting of an ocean general circulation model (OGCM) and a statistical atmospheric model for interannual anomalies of wind stress (τ). The HCM AOPB serves as a simple Earth system for the tropical Pacific to represent the coupling among the atmospheric and physical and biogeochemical ocean components. Model experiments are performed to illustrate this new model's ability to depict the mean ocean state and interannual variability associated with the ENSO. The relationships among anomaly fields are analyzed to illustrate the ocean biogeochemistry‐induced heating feedback and its modulating effects on the ENSO, which is characterized by a negative feedback. The underlying processes and mechanisms are analyzed and can be attributed to dominant modulation of the penetrative solar radiation through the base of the mixed layer (ML). It is demonstrated that the ocean biogeochemistry‐induced negative feedback is mainly driven by more solar radiation penetrating out of the ML during El Nino and less penetrating during La Nina. Further model applications to studies on these processes and biogeochemical‐physical interactions are discussed.