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.     

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.     

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.     

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     

On the mechanisms in a tropical ocean–global atmosphere coupled general circulation model. Part I: mean state and the seasonal cycle

 The mechanisms responsible for the mean state and the seasonal and interannual variations of the coupled tropical Pacific-global atmosphere system are investigated by analyzing a thirty year simulation, where the LMD global atmospheric model and the LODYC tropical Pacific model are coupled using the delocalized physics method. No flux correction is needed over the tropical region. The coupled model reaches its regime state roughly after one year of integration in spite of the fact that the ocean is initialized from rest. Departures from the mean state are characterized by oscillations with dominant periodicites at annual, biennial and quadriennial time scales. In our model, equatorial sea surface temperature and wind stress fluctuations evolved in phase. In the Central Pacific during boreal autumn, the sea surface temperature is cold, the wind stress is strong, and the Inter Tropical Convergence Zone (ITCZ) is shifted northwards. The northward shift of the ITCZ enhances atmospheric and oceanic subsidence between the equator and the latitude of organized convention. In turn, the stronger oceanic subsidence reinforces equatorward convergence of water masses at the thermocline depth which, being not balanced by equatorial upwelling, deepens the equatorial thermocline. An equivalent view is that the deepening of the thermocline proceeds from the weakening of the meridional draining of near-surface equatorial waters. The inverse picture prevails during spring, when the equatorial sea surface temperatures are warm. Thus temperature anomalies tend to appear at the thermocline level, in phase opposition to the surface conditions. These subsurface temperature fluctuations propagate from the Central Pacific eastwards along the thermocline; when reaching the surface in the Eastern Pacific, they trigger the reversal of sea surface temperature anomalies. The whole oscillation is synchronized by the apparent meridional motion of the sun, through the seasonal oscillation of the ITCZ. This possible mechanism is partly supported by the observed seasonal reversal of vorticity between the equator and the ITCZ, and by observational evidence of eastward propagating subsurface temperature anomalies at the thermocline level. Received: 7 April 1997 / Accepted: 15 July 1998     

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

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

An efficient climate model with a 3D ocean and statistical–dynamical atmosphere*

We describe a coupled climate model of intermediate complexity designed for use in global warming experiments. The atmospheric component is a two-dimensional (zonally averaged) statistical-dynamical model based on the Goddard Institute for Space Study's atmospheric general circulation model (GCM). In contrast with energy-balance models used in some climate models of intermediate complexity, this model includes full representation of the hydrological and momentum cycles. It also has parameterizations of the main physical processes, including a sophisticated radiation code. The ocean component is a coarse resolution ocean GCM with simplified global geometry based on the Geophysical Fluid Dynamics Laboratory modular ocean model. Because of the simplified geometry the resolution in the western boundary layers can be readily increased compared to conventional coarse resolution models, without increasing the model's computational requirements in a significant way. The ocean model's efficiency is also greatly increased by using a moderate degree of asynchronous coupling between the oceanic momentum and tracer fields. We demonstrate that this still allows an accurate simulation of transient behavior, including the seasonal cycle. A 100 years simulation with the model requires less than 8 hours on a state-of the art workstation. The main novelty of the model is therefore a combination of computational efficiency, statistical-dynamical atmosphere and 3D ocean. Long-term present-day climate simulations are carried out using the coupled model with and without flux adjustments, and with either the Gent-McWilliams (GM) parametrization scheme or horizontal diffusion (HD) in the ocean. Deep ocean temperatures systematically decrease in the runs without flux adjustment. We demonstrate that the mismatch between heat transports in the uncoupled states of two models is the main cause for the systematic drift. In addition, changes in the circulation and sea-ice formation also contribute to the drift. Flux adjustments in the freshwater fluxes are shown to have a stabilizing effect on the thermohaline circulation in the model, whereas the adjustments in the heat fluxes tend to weaken the global "conveyor". To evaluate the model's response to transient external forcing global warming simulations are also carried out with the flux-adjusted version of the coupled model. The coupled model reproduces reasonably well the behavior of more sophisticated coupled GCMs for both current climate and for the global warming scenarios.     

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.     

The climatology of the Canadian Climate Centre general circulation model as obtained from a five‐year simulation

Abstract

The climatology of the Canadian Climate Centre atmospheric general circulation model is presented and compared with the observed climatology of the atmosphere. The model climatology is obtained from a simulation over five annual cycles and the results are presented in terms of averages for the four seasons.

The climatology of the model is discussed in terms of zonally and time averaged values of the primary atmospheric variables as well as in terms of the spatial distributions of the important surface parameters and of the rotational and divergent components of the tropospheric flow. Some measure of model variability is also presented.

The model is generally quite successful in reproducing the mean observed climatology of the atmosphere.     

Coupled atmosphere–ocean data assimilation experiments with a low-order model and CMIP5 model data

Climate Dynamics - Coupled atmosphere–ocean data assimilation (DA) experiments are performed for estimating the Atlantic meridional overturning circulation (AMOC). Recovery of the AMOC with...     

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.     

Time-dependent response of a zonally averaged ocean–atmosphere–sea ice model to Milankovitch forcing

An ocean–atmosphere–sea ice model is developed to explore the time-dependent response of climate to Milankovitch forcing for the time interval 5–3 Myr BP. The ocean component is a zonally averaged model of the circulation in five basins (Arctic, Atlantic, Indian, Pacific, and Southern Oceans). The atmospheric component is a one-dimensional (latitudinal) energy balance model, and the sea-ice component is a thermodynamic model. Two numerical experiments are conducted. The first experiment does not include sea ice and the Arctic Ocean; the second experiment does. Results from the two experiments are used to investigate (1) the response of annual mean surface air and ocean temperatures to Milankovitch forcing, and (2) the role of sea ice in this response. In both experiments, the response of air temperature is dominated by obliquity cycles at most latitudes. On the other hand, the response of ocean temperature varies with latitude and depth. Deep water formed between 45°N and 65°N in the Atlantic Ocean mainly responds to precession. In contrast, deep water formed south of 60°S responds to obliquity when sea ice is not included. Sea ice acts as a time-integrator of summer insolation changes such that annual mean sea-ice conditions mainly respond to obliquity. Thus, in the presence of sea ice, air temperature changes over the sea ice are amplified, and temperature changes in deep water of southern origin are suppressed since water below sea ice is kept near the freezing point.     

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.     

Wind–wave stabilization by a foam layer between the atmosphere and the ocean

The study is motivated by recent findings of the decrease in the momentum transfer from strong winds to sea. The Kelvin–Helmholtz instability (KHI) of a three-fluid system of air, foam and water is examined within the range of intermediately short surface waves. The foam-layer thickness necessary for effective separation of the atmosphere and the ocean is estimated. Due to high density contrasts in the three-fluid system, even a relatively thin foam layer between the atmosphere and the ocean can provide a significant stabilization of the water surface by the wavelength shift of the instability towards smaller scales. It is conjectured that such stabilization qualitatively explains the observed reduction of roughness and drag.     

An atmosphere–ocean regional climate model for the Mediterranean area: assessment of a present climate simulation

We present an atmosphere–ocean regional climate model for the Mediterranean basin, called the PROTHEUS system, composed by the regional climate model RegCM3 as the atmospheric component and by a regional configuration of the MITgcm model as the oceanic component. The model is applied to an area encompassing the Mediterranean Sea and compared to a stand-alone version of its atmospheric component. An assessment of the model performances is done by using available observational datasets. Despite a persistent bias, the PROTHEUS system is able to capture the inter-annual variability of seasonal sea surface temperature (SST) and also the fine scale spatio-temporal evolution of observed SST anomalies, with spatial correlation as high as 0.7 during summer. The close inspection of a 10-day strong wind event during the summer of 2000 proves the capability of the PROTHEUS system to correctly describe the daily evolution of SST under strong air–sea interaction conditions. As a consequence of the model’s skill in reproducing observed SST and wind fields, we expect a reliable estimation of air–sea fluxes. The model skill in reproducing climatological land surface fields is in line with that of state of the art regional climate 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.     

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.