A multi-model study of atmosphere predictability in coupled ocean–atmosphere systems

Of great importance for guiding numerical weather and climate predictions, understanding predictability of the atmosphere in the ocean − atmosphere coupled system is the first and critical step to understand predictability of the Earth system. However, previous predictability studies based on prefect model assumption usually depend on a certain model. Here we apply the predictability study with the Nonlinear Local Lyapunov Exponent and Attractor Radius to the products of multiple re-analyses and forecast models in several operational centers to realize general predictability of the atmosphere in the Earth system. We first investigated the predictability characteristics of the atmosphere in NCEP, ECMWF and UKMO coupled systems and some of their uncoupled counterparts and other uncoupled systems. Although the ECMWF Integrated Forecast System shows higher skills in geopotential height over the tropics, there is no certain model providing the most precise forecast for all variables on all levels and the multi-model ensemble not always outperforms a single model. Improved low-frequency signals from the air − sea and stratosphere − troposphere interactions that extend predictability of the atmosphere in coupled system suggests the significance of air − sea coupling and stratosphere simulation in practical forecast development, although uncertainties exist in the model representation for physical processes in air − sea interactions and upper troposphere. These inspire further exploration on predictability of ocean and stratosphere as well as sea − ice and land processes to advance our understanding of interactions of Earth system components, thus enhancing weather − climate prediction skills.

     

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.     

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.     

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.     

Interannual variations of the boreal summer intraseasonal variability predicted by ten atmosphere–ocean coupled models

The reproducibility of boreal summer intraseasonal variability (ISV) and its interannual variation by dynamical models are assessed through diagnosing 21-year retrospective forecasts from ten state-of-the-art ocean–atmosphere coupled prediction models. To facilitate the assessment, we have defined the strength of ISV activity by the standard deviation of 20–90 days filtered precipitation during the boreal summer of each year. The observed climatological ISV activity exhibits its largest values over the western North Pacific and Indian monsoon regions. The notable interannual variation of ISV activity is found primarily over the western North Pacific in observation while most models have the largest variability over the central tropical Pacific and exhibit a wide range of variability in spatial patterns that are different from observation. Although the models have large systematic biases in spatial pattern of dominant variability, the leading EOF modes of the ISV activity in the models are closely linked to the models’ El Nino-Southern Oscillation (ENSO), which is a feature that resembles the observed ISV and ENSO relationship. The ENSO-induced easterly vertical shear anomalies in the western and central tropical Pacific, where the summer mean vertical wind shear is weak, result in ENSO-related changes of ISV activity in both observation and models. It is found that the principal components of the predicted dominant modes of ISV activity fluctuate in a very similar way with observed ones. The model biases in the dominant modes are systematic and related to the external SST forcing. Thus the statistical correction method of this study based on singular value decomposition is capable of removing a large portion of the systematic errors in the predicted spatial patterns. The 21-year-averaged pattern correlation skill increases from 0.25 to 0.65 over the entire Asian monsoon region after applying the bias correction method to the multi-model ensemble mean prediction.     

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.     

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.     

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).     

Uncertainty in the ocean–atmosphere feedbacks associated with ENSO in the reanalysis products

The evolution of El Ni?o-Southern Oscillation (ENSO) variability can be characterized by various ocean–atmosphere feedbacks, for example, the influence of ENSO related sea surface temperature (SST) variability on the low-level wind and surface heat fluxes in the equatorial tropical Pacific, which in turn affects the evolution of the SST. An analysis of these feedbacks requires physically consistent observational data sets. Availability of various reanalysis data sets produced during the last 15?years provides such an opportunity. A consolidated estimate of ocean surface fluxes based on multiple reanalyses also helps understand biases in ENSO predictions and simulations from climate models. In this paper, the intensity and the spatial structure of ocean–atmosphere feedback terms (precipitation, surface wind stress, and ocean surface heat flux) associated with ENSO are evaluated for six different reanalysis products. The analysis provides an estimate for the feedback terms that could be used for model validation studies. The analysis includes the robustness of the estimate across different reanalyses. Results show that one of the “coupled” reanalysis among the six investigated is closer to the ensemble mean of the results, suggesting that the coupled data assimilation may have the potential to better capture the overall atmosphere–ocean feedback processes associated with ENSO than the uncoupled ones.     

How well do coupled models replicate ocean energetics relevant to ENSO?

Accurate replication of the processes associated with the energetics of the tropical ocean is necessary if coupled GCMs are to simulate the physics of ENSO correctly, including the transfer of energy from the winds to the ocean thermocline and energy dissipation during the ENSO cycle. Here, we analyze ocean energetics in coupled GCMs in terms of two integral parameters describing net energy loss in the system using the approach recently proposed by Brown and Fedorov (J Clim 23:1563?C1580, 2010a) and Fedorov (J Clim 20:1108?C1117, 2007). These parameters are (1) the efficiency ?? of the conversion of wind power into the buoyancy power that controls the rate of change of the available potential energy (APE) in the ocean and (2) the e-folding rate ?? that characterizes the damping of APE by turbulent diffusion and other processes. Estimating these two parameters for coupled models reveals potential deficiencies (and large differences) in how state-of-the-art coupled GCMs reproduce the ocean energetics as compared to ocean-only models and data assimilating models. The majority of the coupled models we analyzed show a lower efficiency (values of ?? in the range of 10?C50% versus 50?C60% for ocean-only simulations or reanalysis) and a relatively strong energy damping (values of ??^?1 in the range 0.4?C1?years versus 0.9?C1.2?years). These differences in the model energetics appear to reflect differences in the simulated thermal structure of the tropical ocean, the structure of ocean equatorial currents, and deficiencies in the way coupled models simulate ENSO.     

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.     

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.     

Future projections of the surface heat and water budgets of the Mediterranean Sea in an ensemble of coupled atmosphere–ocean regional climate models

Within the CIRCE project “Climate change and Impact Research: the Mediterranean Environment”, an ensemble of high resolution coupled atmosphere–ocean regional climate models (AORCMs) are used to simulate the Mediterranean climate for the period 1950–2050. For the first time, realistic net surface air-sea fluxes are obtained. The sea surface temperature (SST) variability is consistent with the atmospheric forcing above it and oceanic constraints. The surface fluxes respond to external forcing under a warming climate and show an equivalent trend in all models. This study focuses on the present day and on the evolution of the heat and water budget over the Mediterranean Sea under the SRES-A1B scenario. On the contrary to previous studies, the net total heat budget is negative over the present period in all AORCMs and satisfies the heat closure budget controlled by a net positive heat gain at the strait of Gibraltar in the present climate. Under climate change scenario, some models predict a warming of the Mediterranean Sea from the ocean surface (positive net heat flux) in addition to the positive flux at the strait of Gibraltar for the 2021–2050 period. The shortwave and latent flux are increasing and the longwave and sensible fluxes are decreasing compared to the 1961–1990 period due to a reduction of the cloud cover and an increase in greenhouse gases (GHGs) and SSTs over the 2021–2050 period. The AORCMs provide a good estimates of the water budget with a drying of the region during the twenty-first century. For the ensemble mean, he decrease in precipitation and runoff is about 10 and 15% respectively and the increase in evaporation is much weaker, about 2% compared to the 1961–1990 period which confirm results obtained in recent studies. Despite a clear consistency in the trends and results between the models, this study also underlines important differences in the model set-ups, methodology and choices of some physical parameters inducing some difference in the various air-sea fluxes. An evaluation of the uncertainty sources and possible improvement for future generation of AORCMs highlights the importance of the parameterisation of the ocean albedo, rivers and cloud cover.     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

Seasonal prediction of global sea level anomalies using an ocean–atmosphere dynamical model

Advanced warning of extreme sea level events is an invaluable tool for coastal communities, allowing the implementation of management policies and strategies to minimise loss of life and infrastructure damage. This study is an initial attempt to apply a dynamical coupled ocean–atmosphere model to the prediction of seasonal sea level anomalies (SLA) globally for up to 7 months in advance. We assess the ability of the Australian Bureau of Meteorology’s operational seasonal dynamical forecast system, the Predictive Ocean Atmosphere Model for Australia (POAMA), to predict seasonal SLA, using gridded satellite altimeter observation-based analyses over the period 1993–2010 and model reanalysis over 1981–2010. Hindcasts from POAMA are based on a 33-member ensemble of seasonal forecasts that are initialised once per month for the period 1981–2010. Our results show POAMA demonstrates high skill in the equatorial Pacific basin and consistently exhibits more skill globally than a forecast based on persistence. Model predictability estimates indicate there is scope for improvement in the higher latitudes and in the Atlantic and Southern Oceans. Most characteristics of the asymmetric SLA fields generated by El-Nino/La Nina events are well represented by POAMA, although the forecast amplitude weakens with increasing lead-time.     

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.     

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.     

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.