On the mechanisms in a tropical ocean–global atmosphere coupled general circulation model. Part II: interannual variability and its relation to the seasonal cycle

 The thirty year simulation of the coupled global atmosphere-tropical Pacific Ocean general circulation model of the Laboratoire de Métérologie Dynamique and the Laboratoire d’Océanographie Dynamique et de Climatologie presented in Part I is further investigated in order to understand the mechanisms of interannual variability. The model does simulate interannual events with ENSO characteristics; the dominant periodicity is quasi-biennial, though strong events are separated by four year intervals. The mechanism that is responsible for seasonal oscillations, identified in Part I, is also active in interannual variability with the difference that now the Western Pacific is dynamically involved. A warm interannual phase is associated with an equatorward shift of the ITCZ in the Western and Central Pacific. The coupling between the ITCZ and the ocean circulation is then responsible for the cooling of the equatorial subsurface by the draining mechanism. Cold subsurface temperature anomalies then propagate eastward along the mean equatorial thermocline. Upon reaching the Eastern Pacific where the mean thermocline is shallow, cold subsurface anomalies affect surface temperatures and reverse the phase of the oscillation. The preferred season for efficient eastward propagation of thermocline depth temperature anomalies is boreal autumn, when draining of equatorial waters towards higher latitudes is weaker than in spring by a factor of six. In that way, the annual cycle acts as a dam that synchronizes lower frequency oscillations. Received: 7 April 1997 / Accepted: 15 July 1998     

Upper-bound general circulation of coupled ocean–atmosphere: Part 1. Atmosphere

We consider the general atmospheric circulation within the deductive framework of our climate theory. The preceding three parts of this theory have reduced the troposphere to the tropical and polar air masses and determined their temperature and the surface latitude of their dividing boundary, which provide the prior thermal constraint for the present dynamical derivation. Drawing upon its similar material conservation as the thermal property, the (columnar) potential vorticity (PV) is assumed homogenized as well in air masses, which moreover has a zero tropical value owing to the hemispheric symmetry. Inverting this PV field produces an upper-bound zonal wind that resembles the prevailing wind, suggesting that the latter may be explained as the maximum macroscopic motion extractable by random eddies – within the confine of the thermal differentiation.With the polar front determined in conjunction with the zonal wind, the approximate leveling of the isobars at the surface and high aloft specifies the tropopause, which is colder and higher in the tropics than in the polar region. The zonal wind drives the meridional circulation via the Ekman dynamics, and the preeminence of the Hadley cell stems from the singular Ekman convergence at the equator that allows it to supply the upward mass flux in the ITCZ demanded by the global energy balance.     

NAO–ocean circulation interactions in a coupled general circulation model

The interplay between the North Atlantic Oscillation (NAO) and the large scale ocean circulation is inspected in a twentieth century simulation conducted with a state-of-the-art coupled general circulation model. Significant lead–lag covariance between oceanic and tropospheric variables suggests that the system supports a damped oscillatory mode involving an active ocean–atmosphere coupling, with a typical NAO-like space structure and a 5 years timescale, qualitatively consistent with a mid-latitude delayed oscillator paradigm. The two essential processes governing the oscillation are (1) a negative feedback between ocean gyre circulation and the high latitude SST meridional gradient and (2) a positive feedback between SST and the NAO. The atmospheric NAO pattern appears to have a weaker projection on the ocean meridional overturning, compared to the gyre circulation, which leads to a secondary role for the thermohaline circulation in driving the meridional heat transport, and thus the oscillatory mode.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere–ocean global coupled model. Part I: a mechanism governing the variability of ocean convection in a preindustrial experiment

A preindustrial climate experiment was conducted with the third version of the CNRM global atmosphere–ocean–sea ice coupled model (CNRM-CM3) for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). This experiment is used to investigate the main physical processes involved in the variability of the North Atlantic ocean convection and the induced variability of the Atlantic meridional overturning circulation (MOC). Three ocean convection sites are simulated, in the Labrador, Irminger and Greenland–Iceland–Norwegian (GIN) Seas in agreement with observations. A mechanism linking the variability of the Arctic sea ice cover and convection in the GIN Seas is highlighted. Contrary to previous suggested mechanisms, in CNRM-CM3 the latter is not modulated by the variability of freshwater export through Fram Strait. Instead, the variability of convection is mainly driven by the variability of the sea ice edge position in the Greenland Sea. In this area, the surface freshwater balance is dominated by the freshwater input due to the melting of sea ice. The ice edge position is modulated either by northwestward geostrophic current anomalies or by an intensification of northerly winds. In the model, stronger than average northerly winds force simultaneous intense convective events in the Irminger and GIN Seas. Convection interacts with the thermohaline circulation on timescales of 5–10 years, which translates into MOC anomalies propagating southward from the convection sites.     

Seasonal climate hindcasts with Eta model nested in CPTEC coupled ocean–atmosphere general circulation model

This work evaluates the added value of the downscaling technique employed with the Eta model nested in the CPTEC atmospheric general circulation model and in the CPTEC coupled ocean?Catmosphere general circulation model (CGCM). The focus is on the austral summer season, December?CJanuary?CFebruary, with three members each year. Precipitation, latent heat flux, and shortwave radiation flux at the surface hindcast by the models are compared with observational data and model analyses. The global models generally overestimate the precipitation over South America and tropical Atlantic. The CGCM and the nested Eta (Eta + C) both produce a split in the ITCZ precipitation band. The Eta + C produces better precipitation pattern for the studied season. The Eta model reduces the excessive latent heat flux generated by these global models, in particular the Eta + C. Comparison against PIRATA buoys data shows that the Eta + C results in the smallest precipitation and shortwave radiation forecast errors. The Eta + C comparatively best results are though as a consequence of both: the regional model resolution/physics and smaller errors on the lateral boundary conditions provided by the CGCM.     

Predictability of the subtropical dipole modes in a coupled ocean–atmosphere model

Predictability of the subtropical dipole modes is assessed using the SINTEX-F coupled model. Despite the known difficulty in predicting subtropical climate due to large internal variability of the atmosphere and weak ocean–atmosphere coupling, it is shown for the first time that the coupled model can successfully predict the South Atlantic Subtropical Dipole (SASD) 1 season ahead, and the prediction skill is better than the persistence in all the 1–12 month lead hindcast experiments. There is a prediction barrier in austral winter due to the seasonal phase locking of the SASD to austral summer. The prediction skill is lower for the Indian Ocean Subtropical Dipole (IOSD) than for the SASD, and only slightly better than the persistence till 6-month lead because of the low predictability of the sea surface temperature anomaly in its southwestern pole. However, for some strong IOSD events in the last three decades, the model can predict them 1 season ahead. The co-occurrence of the negative SASD and IOSD in 1997/1998 austral summer can be predicted from July 1st of 1997. This is because the negative sea level pressure anomalies over the South Atlantic and the southern Indian Ocean in September–October (November–December) that trigger the occurrence of the negative SASD and IOSD are related to the well predicted tropical Indian Ocean Dipole (El Niño/Southern Oscillation). Owing to the overall good performances of the SINTEX-F model in predicting the SASD, some strong IOSD, and El Niño/Southern Oscillation, the prediction skill of the southern African summer precipitation is high in the SINTEX-F model.     

Impact of resolving the diurnal cycle in an ocean–atmosphere GCM. Part 2: A diurnally coupled CGCM

Coupled ocean atmosphere general circulation models (GCM) are typically coupled once every 24 h, excluding the diurnal cycle from the upper ocean. Previous studies attempting to examine the role of the diurnal cycle of the upper ocean and particularly of diurnal SST variability have used models unable to resolve the processes of interest. In part 1 of this study a high vertical resolution ocean GCM configuration with modified physics was developed that could resolve the diurnal cycle in the upper ocean. In this study it is coupled every 3 h to atmospheric GCM to examine the sensitivity of the mean climate simulation and aspects of its variability to the inclusion of diurnal ocean-atmosphere coupling. The inclusion of the diurnal cycle leads to a tropics wide increase in mean sea surface temperature (SST), with the strongest signal being across the equatorial Pacific where the warming increases from 0.2°C in the central and western Pacific to over 0.3°C in the eastern equatorial Pacific. Much of this warming is shown to be a direct consequence of the rectification of daily mean SST by the diurnal variability of SST. The warming of the equatorial Pacific leads to a redistribution of precipitation from the Inter tropical convergence zone (ITCZ) toward the equator. In the western Pacific there is an increase in precipitation between Papa new guinea and 170°E of up to 1.2 mm/day, improving the simulation compared to climatology. Pacific sub tropical cells are increased in strength by about 10%, in line with results of part 1 of this study, due to the modification of the exchange of momentum between the equatorially divergent Ekman currents and the geostropic convergence at depth, effectively increasing the dynamical response of the tropical Pacific to zonal wind stresses. During the spring relaxation of the Pacific trade winds, a large diurnal cycle of SST increases the seasonal warming of the equatorial Pacific. When the trade winds then re-intensify, the increase in the dynamical response of the ocean leads to a stronger equatorial upwelling. These two processes both lead to stronger seasonal basin scale feedbacks in the coupled system, increasing the strength of the seasonal cycle of the tropical Pacific sector by around 10%. This means that the diurnal cycle in the upper ocean plays a part in the coupled feedbacks between ocean and atmosphere that maintain the basic state and the timing of the seasonal cycle of SST and trade winds in the tropical Pacific. The Madden–Julian Oscillation (MJO) is examined by use of a large scale MJO index, lag correlations and composites of events. The inclusion of the diurnal cycle leads to a reduction in overall MJO activity. Precipitation composites show that the MJO is stronger and more coherent when the diurnal cycle of coupling is resolved, with the propagation and different phases being far more distinct both locally and to larger lead times across the tropical Indo-Pacific. Part one of this study showed that that diurnal variability of SST is modulated by the MJO and therefore increases the intraseasonal SST response to the different phases of the MJO. Precipitation-based composites of SST variability confirm this increase in the coupled simulations. It is argued that including this has increased the thermodynamical coupling of the ocean and atmosphere on the timescale of the MJO (20–100 days), accounting for the improvement in the MJO strength and coherency seen in composites of precipitation and SST. These results show that the diurnal cycle of ocean–atmosphere interaction has profound impact on a range of up-scale variability in the tropical climate and as such, it is an important feature of the modelled climate system which is currently either neglected or poorly resolved in state of the art coupled models.     

Dynamics and predictability of a low-order wind-driven ocean–atmosphere coupled model

The dynamics of a low-order coupled wind-driven ocean–atmosphere system is investigated with emphasis on its predictability properties. The low-order coupled deterministic system is composed of a baroclinic atmosphere for which 12 dominant dynamical modes are only retained (Charney and Straus in J Atmos Sci 37:1157–1176, 1980) and a wind-driven, quasi-geostrophic and reduced-gravity shallow ocean whose field is truncated to four dominant modes able to reproduce the large scale oceanic gyres (Pierini in J Phys Oceanogr 41:1585–1604, 2011). The two models are coupled through mechanical forcings only. The analysis of its dynamics reveals first that under aperiodic atmospheric forcings only dominant single gyres (clockwise or counterclockwise) appear, while for periodic atmospheric solutions the double gyres emerge. In the present model domain setting context, this feature is related to the level of truncation of the atmospheric fields, as indicated by a preliminary analysis of the impact of higher wavenumber (“synoptic” scale) modes on the development of oceanic gyres. In the latter case, double gyres appear in the presence of a chaotic atmosphere. Second the dynamical quantities characterizing the short-term predictability (Lyapunov exponents, Lyapunov dimension, Kolmogorov–Sinaï (KS) entropy) displays a complex dependence as a function of the key parameters of the system, namely the coupling strength and the external thermal forcing. In particular, the KS-entropy is increasing as a function of the coupling in most of the experiments, implying an increase of the rate of loss of information about the localization of the system on its attractor. Finally the dynamics of the error is explored and indicates, in particular, a rich variety of short term behaviors of the error in the atmosphere depending on the (relative) amplitude of the initial error affecting the ocean, from polynomial (at ² + bt ³ + ct ⁴) up to exponential-like evolutions. These features are explained and analyzed in the light of the recent findings on error growth (Nicolis et al. in J Atmos Sci 66:766–778, 2009).     

A study on air–sea interaction on the simulated seasonal climate in an ocean–atmosphere coupled model

This study investigates the effects of air–sea interaction upon simulated tropical climatology, focusing on the boreal summer mean precipitation and the embedded intra-seasonal oscillation (ISO) signal. Both the daily coupling of ocean–atmosphere and the diurnal variation of sea surface temperature (SST) at every time step by accounting for the ocean mixed layer and surface-energy budget at the ocean surface are considered. The ocean–atmosphere coupled model component of the global/regional integrated model system has been utilized. Results from the coupled model show better precipitation climatology than those from the atmosphere-only model, through the inclusion of SST–cloudiness–precipitation feedback in the coupled system. Cooling the ocean surface in the coupled model is mainly responsible for the improved precipitation climatology, whereas neither the coupling itself nor the diurnal variation in the SST influences the simulated climatology. However, the inclusion of the diurnal cycle in the SST shows a distinct improvement of the simulated ISO signal, by either decreasing or increasing the magnitude of spectral powers, as compared to the simulation results that exclude the diurnal variation of the SST in coupled models.     

Bi-decadal variability excited in the coupled ocean–atmosphere system by strong tropical volcanic eruptions

Decadal and bi-decadal climate responses to tropical strong volcanic eruptions (SVEs) are inspected in an ensemble simulation covering the last millennium based on the Max Planck Institute—Earth system model. An unprecedentedly large collection of pre-industrial SVEs (up to 45) producing a peak annual-average top-of-atmosphere radiative perturbation larger than ?1.5 Wm^?2 is investigated by composite analysis. Post-eruption oceanic and atmospheric anomalies coherently describe a fluctuation in the coupled ocean–atmosphere system with an average length of 20–25 years. The study provides a new physically consistent theoretical framework to interpret decadal Northern Hemisphere (NH) regional winter climates variability during the last millennium. The fluctuation particularly involves interactions between the Atlantic meridional overturning circulation and the North Atlantic gyre circulation closely linked to the state of the winter North Atlantic Oscillation. It is characterized by major distinctive details. Among them, the most prominent are: (a) a strong signal amplification in the Arctic region which allows for a sustained strengthened teleconnection between the North Pacific and the North Atlantic during the first post-eruption decade and which entails important implications from oceanic heat transport and from post-eruption sea ice dynamics, and (b) an anomalous surface winter warming emerging over the Scandinavian/Western Russian region around 10–12 years after a major eruption. The simulated long-term climate response to SVEs depends, to some extent, on background conditions. Consequently, ensemble simulations spanning different phases of background multidecadal and longer climate variability are necessary to constrain the range of possible post-eruption decadal evolution of NH regional winter climates.     

Thirty-two-year ocean–atmosphere coupled downscaling of global reanalysis over the Intra-American Seas

This study examines the oceanic and atmospheric variability over the Intra-American Seas (IAS) from a 32-year integration of a 15-km coupled regional climate model consisting of the Regional Spectral Model (RSM) for the atmosphere and the Regional Ocean Modeling System (ROMS) for the ocean. It is forced at the lateral boundaries by National Centers for Environmental Prediction-Department of Energy (NCEP-DOE R-2) atmospheric global reanalysis and Simplified Ocean Data Assimilation global oceanic reanalysis. This coupled downscaling integration is a free run without any heat flux correction and is referred as the Regional Ocean–Atmosphere coupled downscaling of global Reanalysis over the Intra-American Seas (ROARS). The paper examines the fidelity of ROARS with respect to independent observations that are both satellite based and in situ. In order to provide a perspective on the fidelity of the ROARS simulation, we also compare it with the Climate Forecast System Reanalysis (CFSR), a modern global ocean–atmosphere reanalysis product. Our analysis reveals that ROARS exhibits reasonable climatology and interannual variability over the IAS region, with climatological SST errors less than 1 °C except along the coastlines. The anomaly correlation of the monthly SST and precipitation anomalies in ROARS are well over 0.5 over the Gulf of Mexico, Caribbean Sea, Western Atlantic and Eastern Pacific Oceans. A highlight of the ROARS simulation is its resolution of the loop current and the episodic eddy events off of it. This is rather poorly simulated in the CFSR. This is also reflected in the simulated, albeit, higher variance of the sea surface height in ROARS and the lack of any variability in the sea surface height of the CFSR over the IAS. However the anomaly correlations of the monthly heat content anomalies of ROARS are comparatively lower, especially over the Gulf of Mexico and the Caribbean Sea. This is a result of ROARS exhibiting a bias of underestimation (overestimation) of high (low) clouds. ROARS like CFSR is also able to capture the Caribbean Low Level Jet and its seasonal variability reasonably well.     

Impact of initialization procedures on the predictive skill of a coupled ocean–atmosphere model

The sensitivity of the predictive skill of a decadal climate prediction system is investigated with respect to details of the initialization procedure. For this purpose, the coupled ocean–atmosphere UCLA/MITgcm climate model is initialized using the following three different initialization approaches: full state initialization (FSI), anomaly initialization (AI) and FSI employing heat flux and freshwater flux corrections (FC). The ocean initial conditions are provided by the German contribution to Estimating the Circulation and Climate of the Ocean state estimate (GECCO project), from which ensembles of decadal hindcasts are initialized every 5 years from 1961 to 2001. The predictive skill for sea surface temperature (SST), sea surface height (SSH) and the Atlantic meridional overturning circulation (AMOC) is assessed against the GECCO synthesis. In regions with a deep mixed layer the predictive skill for SST anomalies remains significant for up to a decade in the FC experiment. By contrast, FSI shows less persistent skill in the North Atlantic and AI does not show high skill in the extratropical Southern Hemisphere, but appears to be more skillful in the tropics. In the extratropics, the improved skill is related to the ability of the FC initialization method to better represent the mixed layer depth, and the highest skill occurs during wintertime. The correlation skill for the spatially averaged North Atlantic SSH hindcasts remains significant up to a decade only for FC. The North Atlantic MOC initialized hindcasts show high correlation values in the first pentad while correlation remains significant in the following pentad too for FSI and FC. Overall, for the current setup, the FC approach appears to lead to the best results, followed by the FSI and AI procedures.     

Climate entropy budget of the HadCM3 atmosphere–ocean general circulation model and of FAMOUS, its low-resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ～50 mW m^?2 K^?1, a value intermediate in the range 30–70 mW m^?2 K^?1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (～38 mW m^?2 K^?1). Another 13 mW m^?2 K^?1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m^?2 K^?1. The material entropy production within the ocean due to turbulent mixing is ～1 mW m^?2 K^?1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m^?2 K^?1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.     

A highly nonlinear coupled mode of decadal variability in a mid-latitude ocean–atmosphere model

This study examines mid-latitude climate variability in a model that couples turbulent oceanic and atmospheric flows through an active oceanic mixed layer. Intrinsic ocean dynamics of the inertial recirculation regions combines with nonlinear atmospheric sensitivity to sea-surface temperature (SST) anomalies to play a dominant role in the variability of the coupled system.Intrinsic low-frequency variability arises in the model atmosphere; when run in a stand-alone mode, it is characterized by irregular transitions between preferred high-latitude and less frequent low-latitude zonal-flow states. When the atmosphere is coupled to the ocean, the low-latitude state occurrences exhibit a statistically significant signal in a broad 5–15-year band. A similar signal is found in the time series of the model ocean's energy in this coupled simulation. Accompanying uncoupled ocean-only and atmosphere-only integrations are characterized by a decrease in the decadal-band variability, relative to the coupled integration; their spectra are indistinguishable from a red spectrum.The time scale of the coupled interdecadal oscillation is set by the nonlinear adjustment of the ocean's inertial recirculations to the high-latitude and low-latitude atmospheric forcing regimes. This adjustment involves, in turn, SST changes resulting in long-term ocean–atmosphere heat-flux anomalies that induce the atmospheric regime transitions.     

Coupling ocean–atmosphere intensity determines ocean chlorophyll-induced SST change in the tropical Pacific

Phytoplankton pigments (e.g., chlorophyll-a) absorb solar radiation in the upper ocean and induce a pronounced radiant heating effect (chlorophyll effect) on the climate. However, the ocean chlorophyll-induced heating effect on the mean climate state in the tropical Pacific has not been understood well. Here, a hybrid coupled model (HCM) of the atmosphere, ocean physics and biogeochemistry is used to investigate the chlorophyll effect on sea surface temperature (SST) in the eastern equatorial Pacific; a tunable coefficient, α, is introduced to represent the coupling intensity between the atmosphere and ocean in the HCM. The modeling results show that the chlorophyll effect on the mean-state SST is sensitively dependent on α (the coupling intensity). At weakly represented coupling intensity (0 ≤ α < 1.01), the chlorophyll effect tends to induce an SST cooling in the eastern equatorial Pacific, whereas an SST warming emerges at the strongly represented coupling intensity (α ≥ 1.01). Thus, a threshold exists for the coupling intensity (about α = 1.01) at which the sign of SST responses can change. Mechanisms and processes are illustrated to understand the different SST responses. In the weak coupling cases, indirect dynamical cooling processes (the adjustment of ocean circulation, enhanced vertical mixing, and upwelling) tend to dominate the SST cooling. In the strong coupling cases, the persistent warming induced by chlorophyll in the southern subtropical Pacific tends to induce cross-equatorial northerly winds, which shifts to anomalous westerly winds in the eastern equatorial Pacific, consequently reducing the evaporative cooling and weakening indirect dynamical cooling; eventually, SST warming maintains in the eastern equatorial Pacific. These results provide new insights into the biogeochemical feedback on the climate and bio-physical interactions in the tropical Pacific.

     

Impact of diurnal atmosphere–ocean coupling on tropical climate simulations using a coupled GCM

The impacts of diurnal atmosphere–ocean (air–sea) coupling on tropical climate simulations are investigated using the SNU coupled GCM. To investigate the effect of the atmospheric and oceanic diurnal cycles on a climate simulation, a 1-day air–sea coupling interval experiment is compared to a 2-h coupling experiment. As previous studies have suggested, cold temperature biases over equatorial western Pacific regions are significantly reduced when diurnal air–sea coupling strategy is implemented. This warming is initiated by diurnal rectification and amplified further by the air–sea coupled feedbacks. In addition to its effect on the mean climatology, the diurnal coupling has also a distinctive impact on the amplitude of the El Nino-Southern Oscillation (ENSO). It is demonstrated that a weakening of the ENSO magnitude is caused by reduced (increased) surface net heat fluxes into the ocean during El Nino (La Nina) events. Primarily, decreased (increased) incoming shortwave radiation during El Nino (La Nina) due to cloud shading is responsible for the net heat fluxes associated with ENSO.     

Isolating mesoscale coupled ocean–atmosphere interactions in the Kuroshio Extension region

The Kuroshio Extension region is characterized by energetic oceanic mesoscale and frontal variability that alters the air–sea fluxes that can influence large-scale climate variability in the North Pacific. We investigate this mesoscale air-sea coupling using a regional eddy-resolving coupled ocean–atmosphere (OA) model that downscales the observed large-scale climate variability from 2001 to 2007. The model simulates many aspects of the observed seasonal cycle of OA coupling strength for both momentum and turbulent heat fluxes. We introduce a new modeling approach to study the scale-dependence of two well-known mechanisms for the surface wind response to mesoscale sea surface temperatures (SSTs), namely, the ‘vertical mixing mechanism’ (VMM) and the ‘pressure adjustment mechanism’ (PAM). We compare the fully coupled model to the same model with an online, 2-D spatial smoother applied to remove the mesoscale SST field felt by the atmosphere. Both VMM and PAM are found to be active during the strong wintertime peak seen in the coupling strength in both the model and observations. For VMM, large-scale SST gradients surprisingly generate coupling between downwind SST gradient and wind stress divergence that is often stronger than the coupling on the mesoscale, indicating their joint importance in OA interaction in this region. In contrast, VMM coupling between crosswind SST gradient and wind stress curl occurs only on the mesoscale, and not over large-scale SST gradients, indicating the essential role of the ocean mesocale. For PAM, the model results indicate that coupling between the Laplacian of sea level pressure and surface wind convergence occurs for both mesoscale and large-scale processes, but inclusion of the mesoscale roughly doubles the coupling strength. Coupling between latent heat flux and SST is found to be significant throughout the entire seasonal cycle in both fully coupled mode and large-scale coupled mode, with peak coupling during winter months. The atmospheric response to the oceanic mesoscale SST is also studied by comparing the fully coupled run to an uncoupled atmospheric model forced with smoothed SST prescribed from the coupled run. Precipitation anomalies are found to be forced by surface wind convergence patterns that are driven by mesoscale SST gradients, indicating the importance of the ocean forcing the atmosphere at this scale.     

ENSO at 6ka and 21ka from ocean–atmosphere coupled model simulations

We analyze how the characteristics of El Niño-Southern Oscillation (ENSO) are changed in coupled ocean–atmosphere simulations of the mid-Holocene (MH) and the Last Glacial Maximum (LGM) performed as part of the Paleoclimate Modeling Intercomparison Project phase 2 (PMIP2). Comparison of the model results with present day observations show that most of the models reproduce the large scale features of the tropical Pacific like the SST gradient, the mean SST and the mean seasonal cycles. All models simulate the ENSO variability, although with different skill. Our analyses show that several relationships between El Niño amplitude and the mean state across the different control simulations are still valid for simulations of the MH and the LGM. Results for the MH show a consistent El Niño amplitude decrease. It can be related to the large scale atmospheric circulation changes. While the Northern Hemisphere receives more insolation during the summer time, the Asian summer monsoon system is strengthened which leads to the enhancement of the Walker circulation. Easterlies prevailing over the central eastern Pacific induce an equatorial upwelling that damps the El Niño development. Results are less conclusive for 21ka. Large scale dynamic competes with changes in local heat fluxes, so that model shows a wide range of responses, as it is the case in future climate projections.