EC-Earth V2.2: description and validation of a new seamless earth system prediction model

EC-Earth, a new Earth system model based on the operational seasonal forecast system of the European Centre for Medium-Range Weather Forecasts (ECMWF), is presented. The performance of version 2.2 (V2.2) of the model is compared to observations, reanalysis data and other coupled atmosphere–ocean-sea ice models. The large-scale physical characteristics of the atmosphere, ocean and sea ice are well simulated. When compared to other coupled models with similar complexity, the model performs well in simulating tropospheric fields and dynamic variables, and performs less in simulating surface temperature and fluxes. The surface temperatures are too cold, with the exception of the Southern Ocean region and parts of the Northern Hemisphere extratropics. The main patterns of interannual climate variability are well represented. Experiments with enhanced CO₂ concentrations show well-known responses of Arctic amplification, land-sea contrasts, tropospheric warming and stratospheric cooling. The global climate sensitivity of the current version of EC-Earth is slightly less than 1?K/(W?m^?2). An intensification of the hydrological cycle is found and strong regional changes in precipitation, affecting monsoon characteristics. The results show that a coupled model based on an operational seasonal prediction system can be used for climate studies, supporting emerging seamless prediction strategies.     

Major modes of short-term climate variability in the newly developed NUIST Earth System Model (NESM)

A coupled earth system model(ESM) has been developed at the Nanjing University of Information Science and Technology(NUIST) by using version 5.3 of the European Centre Hamburg Model(ECHAM), version 3.4 of the Nucleus for European Modelling of the Ocean(NEMO), and version 4.1 of the Los Alamos sea ice model(CICE). The model is referred to as NUIST ESM1(NESM1). Comprehensive and quantitative metrics are used to assess the model's major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model's assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring–fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific(CP)-ENSO and eastern Pacific(EP)-ENSO; however, the equatorial SST variability,biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden–Julian Oscillation(MJO), and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version(T159) demonstrates improved simulation with respect to the climatology, interannual variance, monsoon–ENSO lead–lag correlation, spatial structures of the leading mode of the Asian–Australian monsoon rainfall variability, and the eastward propagation of the MJO.     

The simulation of SST, sea ice extents and ocean heat transports in a version of the Hadley Centre coupled model without flux adjustments

Results are presented from a new version of the Hadley Centre coupled model (HadCM3) that does not require flux adjustments to prevent large climate drifts in the simulation. The model has both an improved atmosphere and ocean component. In particular, the ocean has a 1.25° × 1.25° degree horizontal resolution and leads to a considerably improved simulation of ocean heat transports compared to earlier versions with a coarser resolution ocean component. The model does not have any spin up procedure prior to coupling and the simulation has been run for over 400 years starting from observed initial conditions. The sea surface temperature (SST) and sea ice simulation are shown to be stable and realistic. The trend in global mean SST is less than 0.009 °C per century. In part, the improved simulation is a consequence of a greater compatibility of the atmosphere and ocean model heat budgets. The atmospheric model surface heat and momentum budget are evaluated by comparing with climatological ship-based estimates. Similarly the ocean model simulation of poleward heat transports is compared with direct ship-based observations for a number of sections across the globe. Despite the limitations of the observed datasets, it is shown that the coupled model is able to reproduce many aspects of the observed heat budget. Received: 1 October 1998 / Accepted: 20 July 1999     

Description and evaluation of the bergen climate model: ARPEGE coupled with MICOM

A new coupled atmosphere–ocean–sea ice model has been developed, named the Bergen Climate Model (BCM). It consists of the atmospheric model ARPEGE/IFS, together with a global version of the ocean model MICOM including a dynamic–thermodynamic sea ice model. The coupling between the two models uses the OASIS software package. The new model concept is described, and results from a 300-year control integration is evaluated against observational data. In BCM, both the atmosphere and the ocean components use grids which can be irregular and have non-matching coastlines. Much effort has been put into the development of optimal interpolation schemes between the models, in particular the non-trivial problem of flux conservation in the coastal areas. A flux adjustment technique has been applied to the heat and fresh-water fluxes. There is, however, a weak drift in global mean sea-surface temperature (SST) and sea-surface salinity (SSS) of respectively 0.1 °C and 0.02 psu per century. The model gives a realistic simulation of the radiation balance at the top-of-the-atmosphere, and the net surface fluxes of longwave, shortwave, and turbulent heat fluxes are within observed values. Both global and total zonal means of cloud cover and precipitation are fairly close to observations, and errors are mainly related to the strength and positioning of the Hadley cell. The mean sea-level pressure (SLP) is well simulated, and both the mean state and the interannual standard deviation show realistic features. The SST field is several degrees too cold in the equatorial upwelling area in the Pacific, and about 1 °C too warm along the eastern margins of the oceans, and in the polar regions. The deviation from Levitus salinity is typically 0.1 psu – 0.4 psu, with a tendency for positive anomalies in the Northern Hemisphere, and negative in the Southern Hemisphere. The sea-ice distribution is realistic, but with too thin ice in the Arctic Ocean and too small ice coverage in the Southern Ocean. These model deficiencies have a strong influence on the surface air temperatures in these regions. Horizontal oceanic mass transports are in the lower range of those observed. The strength of the meridional overturning in the Atlantic is 18 Sv. An analysis of the large-scale variability in the model climate reveals realistic El Niño – Southern Oscillation (ENSO) and North Atlantic–Arctic Oscillation (NAO/AO) characteristics in the SLP and surface temperatures, including spatial patterns, frequencies, and strength. While the NAO/AO spectrum is white in SLP and red in temperature, the ENSO spectrum shows an energy maximum near 3 years.     

Arctic climate change in 21st century CMIP5 simulations with EC-Earth

The Arctic climate change is analyzed in an ensemble of future projection simulations performed with the global coupled climate model EC-Earth2.3. EC-Earth simulates the twentieth century Arctic climate relatively well but the Arctic is about 2 K too cold and the sea ice thickness and extent are overestimated. In the twenty-first century, the results show a continuation and strengthening of the Arctic trends observed over the recent decades, which leads to a dramatically changed Arctic climate, especially in the high emission scenario RCP8.5. The annually averaged Arctic mean near-surface temperature increases by 12 K in RCP8.5, with largest warming in the Barents Sea region. The warming is most pronounced in winter and autumn and in the lower atmosphere. The Arctic winter temperature inversion is reduced in all scenarios and disappears in RCP8.5. The Arctic becomes ice free in September in all RCP8.5 simulations after a rapid reduction event without recovery around year 2060. Taking into account the overestimation of ice in the twentieth century, our model results indicate a likely ice-free Arctic in September around 2040. Sea ice reductions are most pronounced in the Barents Sea in all RCPs, which lead to the most dramatic changes in this region. Here, surface heat fluxes are strongly enhanced and the cloudiness is substantially decreased. The meridional heat flux into the Arctic is reduced in the atmosphere but increases in the ocean. This oceanic increase is dominated by an enhanced heat flux into the Barents Sea, which strongly contributes to the large sea ice reduction and surface-air warming in this region. Increased precipitation and river runoff lead to more freshwater input into the Arctic Ocean. However, most of the additional freshwater is stored in the Arctic Ocean while the total Arctic freshwater export only slightly increases.     

An assessment of the surface climate in the NCEP climate forecast system reanalysis

This paper analyzes surface climate variability in the climate forecast system reanalysis (CFSR) recently completed at the National Centers for Environmental Prediction (NCEP). The CFSR represents a new generation of reanalysis effort with first guess from a coupled atmosphere?Cocean?Csea ice?Cland forecast system. This study focuses on the analysis of climate variability for a set of surface variables including precipitation, surface air 2-m temperature (T2m), and surface heat fluxes. None of these quantities are assimilated directly and thus an assessment of their variability provides an independent measure of the accuracy. The CFSR is compared with observational estimates and three previous reanalyses (the NCEP/NCAR reanalysis or R1, the NCEP/DOE reanalysis or R2, and the ERA40 produced by the European Centre for Medium-Range Weather Forecasts). The CFSR has improved time-mean precipitation distribution over various regions compared to the three previous reanalyses, leading to a better representation of freshwater flux (evaporation minus precipitation). For interannual variability, the CFSR shows improved precipitation correlation with observations over the Indian Ocean, Maritime Continent, and western Pacific. The T2m of the CFSR is superior to R1 and R2 with more realistic interannual variability and long-term trend. On the other hand, the CFSR overestimates downward solar radiation flux over the tropical Western Hemisphere warm pool, consistent with a negative cloudiness bias and a positive sea surface temperature bias. Meanwhile, the evaporative latent heat flux in CFSR appears to be larger than other observational estimates over most of the globe. A few deficiencies in the long-term variations are identified in the CFSR. Firstly, dramatic changes are found around 1998?C2001 in the global average of a number of variables, possibly related to the changes in the assimilated satellite observations. Secondly, the use of multiple streams for the CFSR induces spurious jumps in soil moisture between adjacent streams. Thirdly, there is an inconsistency in long-term sea ice extent variations over the Arctic regions between the CFSR and other observations with the CFSR showing smaller sea ice extent before 1997 and larger extent starting in 1997. These deficiencies may have impacts on the application of the CFSR for climate diagnoses and predictions. Relationships between surface heat fluxes and SST tendency and between SST and precipitation are analyzed and compared with observational estimates and other reanalyses. Global mean fields of surface heat and water fluxes together with radiation fluxes at the top of the atmosphere are documented and presented over the entire globe, and for the ocean and land separately.     

Energy transports by ocean and atmosphere based on an entropy extremum principle. Part II: Two-dimensional transports

Summary The maximum entropy production (MEP) principle used in Part J has been extended to separate the two-dimensional required energy transports determined from Nimbus 7 satellite net radiation measurements into atmospheric and oceanic components. In terms of the meridional component of the ocean transport vectors, results show northward ocean transports throughout the entire Atlantic ocean from southern hemisphere high latitudes to northern hemisphere polar regions, southward transports throughout the entire Indian Ocean, and poleward transports separated at approximately 10°S in the Pacific Ocean. The ocean transport patterns are consistent with well-known features concerning heat transport within the three ocean basins. However, uncertainty remains in the magnitudes of the transports. Because of the large remaining discrepancies between published estimates based on direct measurements and indirect estimates derived from energy budget methods, assessing the accuracy of the magnitudes is difficult, although there is evidence that the limited model resolution leads to synergistic biases in the North Atlantic and North Pacific. In terms of the crossmeridional energy transport component, results suggest that most of the net energy transfer in the tropics takes place within the ocean. In the southern hemisphere high latitudes, the Pacific and Indian Oceans export heat cross-meridionally to the Atlantic Ocean through the passages below Cape Horn and the Cape of Good Hope, although the magnitudes of these inter-ocean heat exchanges are small. Another important aspect of the southern hemisphere results is that poleward transports are dominated by the atmospheric component with strong zonal asymmetry. By contrast, in the northern hemisphere, atmospheric transports over the ocean are generally weaker than the corresponding southern hemisphere terms, indicating that the northern hemisphere oceans are relatively more effective in transferring heat poleward. Finally, poleward atmospheric transports over the continental areas exceed those over the ocean at equivalent latitudes as a result of the generally greater energy deficits over the land areas.With 7 Figures     

Towards a parametrization of river discharges into ocean general circulation models: a closure through energy conservation

Diagnostic methods are defined in order to compare two numerical simulations of ocean dynamics in a region of freshwater influence. The first one is a river plume simulation based on a high resolution numerical configuration of the POM coastal ocean model in which mixing parametrizations have been previously defined. The second one is a simulation based on the NEMO Global Ocean Model used for climate simulations in its half-a-degree configuration in which a river inflow is represented as precipitation on two coastal grid cells. Both simulations are forced with the same freshwater inflows and wind stresses. The divergence of volumetric fluxes above and below the halocline are compared. Results show that when an upwelling wind blows, the two models display similar behavior although the impact of lack of precision can be observed in the NEMO configuration. When a downwelling wind blows, the NEMO Global Ocean configuration can not reproduce the coastally trapped baroclinic dynamics because its grid resolution is too coarse. To find a parametrization to help represent these dynamics in ocean general circulation models, a method based on energy conservation is investigated. This method shows that it is possible to link the energy fluxes provided by river inflows to the divergence of energy fluxes integrated over the grid cells of ocean general circulation models. A parametrization of the dynamics created by freshwater inflows is deduced from this method. This enabled creation of a box model that proved to have the same behavior as the fluxes previously computed from the high resolution configuration.     

Review of simulations of climate variability and change with the GFDL R30 coupled climate model

A review is presented of the development and simulation characteristics of the most recent version of a global coupled model for climate variability and change studies at the Geophysical Fluid Dynamics Laboratory, as well as a review of the climate change experiments performed with the model. The atmospheric portion of the coupled model uses a spectral technique with rhomboidal 30 truncation, which corresponds to a transform grid with a resolution of approximately 3.75° longitude by 2.25° latitude. The ocean component has a resolution of approximately 1.875° longitude by 2.25° latitude. Relatively simple formulations of river routing, sea ice, and land surface processes are included. Two primary versions of the coupled model are described, differing in their initialization techniques and in the specification of sub-grid scale oceanic mixing of heat and salt. For each model a stable control integration of near millennial scale duration has been conducted, and the characteristics of both the time-mean and variability are described and compared to observations. A review is presented of a suite of climate change experiments conducted with these models using both idealized and realistic estimates of time-varying radiative forcing. Some experiments include estimates of forcing from past changes in volcanic aerosols and solar irradiance. The experiments performed are described, and some of the central findings are highlighted. In particular, the observed increase in global mean surface temperature is largely contained within the spread of simulated global mean temperatures from an ensemble of experiments using observationally-derived estimates of the changes in radiative forcing from increasing greenhouse gases and sulfate aerosols.     

A fast version of LASG/IAP climate system model and its 1000-year control integration

A fast version of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geo- physical Fluid Dynamics(LASG)/Institute of Atmospheric Physics(IAP)climate system model is briefly documented.The fast coupled model employs a low resolution version of the atmospheric component Grid Atmospheric Model of IAP/LASG(GAMIL),with the other parts of the model,namely an oceanic component LASG/IAP Climate Ocean Model(LICOM),land component Common Land Model(CLM),and sea ice component from National Center for Atmospheric Research Community Climate System Model (NCAR CCSM2),as the same as in the standard version of LASG/IAP Flexible Global Ocean Atmosphere Land System model(FGOALS g).The parameterizations of physical and dynamical processes of the at- mospheric component in the fast version are identical to the standard version,although some parameter values are different.However,by virtue of reduced horizontal resolution and increased time-step of the most time-consuming atmospheric component,it runs faster by a factor of 3 and can serve as a useful tool for long- term and large-ensemble integrations.A 1000-year control simulation of the present-day climate has been completed without flux adjustments.The final 600 years of this simulation has virtually no trends in global mean sea surface temperatures and is recommended for internal variability studies.Several aspects of the control simulation’s mean climate and variability are evaluated against the observational or reanalysis data. The strengths and weaknesses of the control simulation are evaluated.The mean atmospheric circulation is well simulated,except in high latitudes.The Asian-Australian monsoonal meridional cell shows realistic features,however,an artificial rainfall center is located to the eastern periphery of the Tibetan Plateau persists throughout the year.The mean bias of SST resembles that of the standard version,appearing as a"double ITCZ"(Inter-Tropical Convergence Zone)associated with a westward extension of the equatorial eastern Pacific cold tongue.The sea ice extent is acceptable but has a higher concentration.The strength of Atlantic meridional overturning is 27.5 Sv.Evidence from the 600-year simulation suggests a modulation of internal variability on ENSO frequency,since both regular and irregular oscillations of ENSO are found during the different time periods of the long-term simulation.     

Simulation of sea ice in FGOALS-g2: Climatology and late 20th century changes

Sea ice is an important component in the Earth’s climate system. Coupled climate system models are indispensable tools for the study of sea ice, its internal processes, interaction with other components, and projection of future changes. This paper evaluates the simulation of sea ice by the Flexible Global Ocean-Atmosphere-Land System model Grid-point Version 2 (FGOALS-g2), in the fifth phase of the Coupled Model Inter-comparison Project (CMIP5), with a focus on historical experiments and late 20th century simulation. Through analysis, we find that FGOALS-g2 produces reasonable Arctic and Antarctic sea ice climatology and variability. Sea ice spatial distribution and seasonal change characteristics are well captured. The decrease of Arctic sea ice extent in the late 20th century is reproduced in simulations, although the decrease trend is lower compared with observations. Simulated Antarctic sea ice shows a reasonable distribution and seasonal cycle with high accordance to the amplitude of winter-summer changes. Large improvement is achieved as compared with FGOALS-g1.0 in CMIP3. Diagnosis of atmospheric and oceanic forcing on sea ice reveals several shortcomings and major aspects to improve upon in the future: (1) ocean model improvements to remove the artificial island at the North Pole; (2) higher resolution of the atmosphere model for better simulation of important features such as, among others, the Icelandic Low and westerly wind over the Southern Ocean; and (3) ocean model improvements to accurately receive freshwater input from land, and higher resolution for resolving major water channels in the Canadian Arctic Archipelago.     

High-resolution numerical simulation of the Pacific Ocean

The present study developed Pacific Ocean models from the Research Institute for Applied Mechanics Ocean Model (RIAMOM) with very high horizontal (1/6° and 1/12°) and vertical (70 levels) resolutions. The hydrographic features of the simulations show good agreement with observed climatological features. Solution differences between the 1/6° and 1/12° models are small for general features of various physical components, but large for eddy fields and the strengths of western boundary currents and their extensions. However, the two high-resolution models show realistic climatological features of Pacific Ocean circulation patterns. Volume transports through major straits in the northwestern Pacific Ocean were also simulated and compared with previous observational results.     

Freshwater transports in HadCM3

The hydrological cycle can influence climate through a great variety of processes. A good representation of the hydrological cycle in climate models is therefore crucial. Attempts to analyse the global hydrological cycle are hampered by a deficiency of suitable observations, particularly over the oceans. Fully coupled general circulation models are potentially powerful tools in interpreting the limited observational data in the context of large-scale freshwater exchanges. We have looked at large-scale aspects of the global freshwater budget in a simulation, of over 1000 years, by the Hadley Centre coupled climate model (HadCM3). Many aspects of the global hydrological cycle are well represented, but the model hydrological cycle appears to be too strong, with overly large precipitation and evaporation components in comparison with the observational datasets we have used. We show that the ocean basin-scale meridional transports of freshwater come into near balance with the surface freshwater fluxes on a time scale of about 400 years, with the major change being a relative increase of freshwater transport from the Southern Ocean into the Atlantic Ocean. Comparison with observations, supported by sensitivity tests, suggests that the major cause of a drift to more saline condition in the model Atlantic is an overestimate of evaporation, although other freshwater budget components may also play a role. The increase in ocean freshwater transport into the Atlantic during the simulation, primarily coming from the overturning circulation component, which changes from divergent to convergent, acts to balance this freshwater budget deficit. The stability of the thermohaline circulation in HadCM3 may be affected by these freshwater transport changes and this question is examined in the context of an existing conceptual model.     

Seasonal oceanic heat transports computed from an atmospheric model

Seasonal estimates of the oceanic poleward heat transport are obtained using a climate model that is a global atmospheric general circulation model on an 8° × 10° grid. The climate model is used to calculate the surface heat flux into each ocean grid point for each day of the year. The rate of ocean heat storage is calculated using climatological surface temperatures, mixed layer depths, and ice amounts. By assuming that the rate of change of heat storage in the deep ocean is spatially constant, the horizontal transports are calculated from the vertical fluxes and the upper ocean storage rates. The oceanic meridional transport for each latitude and for each ocean basin are derived, and results are compared with other calculations of the seasonal transports. In the Northern Hemisphere, comparisons between the simulated seasonal transports indicate that the annual variation is much greater in the Pacific than in the Atlantic.     

Sensitivity studies with the North Atlantic sigma coordinate Princeton Ocean Model

The sigma coordinate, Princeton Ocean Model (POM) has been configured for the North Atlantic Ocean between 5°N and 50°N as part of data assimilation, model predictability and intercomparison studies. The model uses a curvilinear orthogonal grid with higher resolution in the western North Atlantic and lower resolution in the eastern North Atlantic. A series of experiments, each one of a 10-year duration, are performed to evaluate the sensitivity of the ocean mean state and variability to model parameters and model configuration; these experiments include open vs. closed boundary conditions, low vs. high resolution grids, and different choices of diffusion and viscosity. The results show that the use of closed boundaries together with near-boundary buffer zones where temperature and salinity are relaxed towards the observed values give less realistic flows, weaker recirculation gyres and less realistic Gulf Stream separation than do open boundary conditions. The experiments show that the sensitivity of the ocean variability in the model to the choice of the Smagorinsky diffusion and viscosity coefficients significantly differs from one region to another and largely depends on other attributes such as the mean position of the Gulf Stream in each simulation. A 50% change in model resolution in the Gulf Stream region has a larger effect on ocean variability than a change of diffusivity by a factor of 10. In areas where either the high or the low resolution models have sufficient resolution, as in the Gulf of Mexico, they are able to produce variability comparable to that observed from altimeter data; elsewhere, model variability is underestimated.     

Factors that affect the amplitude of El Nino in global coupled climate models

Historically, El Nino-like events simulated in global coupled climate models have had reduced amplitude compared to observations. Here, El Nino-like phenomena are compared in ten sensitivity experiments using two recent global coupled models. These models have various combinations of horizontal and vertical ocean resolution, ocean physics, and atmospheric model resolution. It is demonstrated that the lower the value of the ocean background vertical diffusivity, the greater the amplitude of El Nino variability which is related primarily to a sharper equatorial thermocline. Among models with low background vertical diffusivity, stronger equatorial zonal wind stress is associated with relatively higher amplitude El Nino variability along with more realistic east–west sea surface temperature (SST) gradient along the equator. The SST seasonal cycle in the eastern tropical Pacific has too much of a semiannual component with a double intertropical convergence zone (ITCZ) in all experiments, and thus does not affect, nor is it affected by, the amplitude of El Nino variability. Systematic errors affecting the spatial variability of El Nino in the experiments are characterized by the eastern equatorial Pacific cold tongue regime extending too far westward into the warm pool. The time scales of interannual variability (as represented by time series of Nino3 SSTs) show significant power in the 3–4 year ENSO band and 2–2.5 year tropospheric biennial oscillation (TBO) band in the model experiments. The TBO periods in the models agree well with the observations, while the ENSO periods are near the short end of the range of 3–6 years observed during the period 1950–94. The close association between interannual variability of equatorial eastern Pacific SSTs and large-scale SST patterns is represented by significant correlations between Nino3 time series and the PC time series of the first EOFs of near-global SSTs in the models and observations. Received: 17 April 2000 / Accepted: 17 August 2000     

Impact of ocean model resolution on CCSM climate simulations

The current literature provides compelling evidence suggesting that an eddy-resolving (as opposed to eddy-permitting or eddy-parameterized) ocean component model will significantly impact the simulation of the large-scale climate, although this has not been fully tested to date in multi-decadal global coupled climate simulations. The purpose of this paper is to examine how resolved ocean fronts and eddies impact the simulation of large-scale climate. The model used for this study is the NCAR Community Climate System Model version 3.5 (CCSM3.5)—the forerunner to CCSM4. Two experiments are reported here. The control experiment is a 155-year present-day climate simulation using a 0.5° atmosphere component (zonal resolution 0.625 meridional resolution 0.5°; land surface component at the same resolution) coupled to ocean and sea-ice components with zonal resolution of 1.2° and meridional resolution varying from 0.27° at the equator to 0.54° in the mid-latitudes. The second simulation uses the same atmospheric and land-surface models coupled to eddy-resolving 0.1° ocean and sea-ice component models. The simulations are compared in terms of how the representation of smaller scale features in the time mean ocean circulation and ocean eddies impact the mean and variable climate. In terms of the global mean surface temperature, the enhanced ocean resolution leads to a ubiquitous surface warming with a global mean surface temperature increase of about 0.2?°C relative to the control. The warming is largest in the Arctic and regions of strong ocean fronts and ocean eddy activity (i.e., Southern Ocean, western boundary currents). The Arctic warming is associated with significant losses of sea-ice in the high-resolution simulation. The sea surface temperature gradients in the North Atlantic, in particular, are better resolved in the high-resolution model leading to significantly sharper temperature gradients and associated large-scale shifts in the rainfall. In the extra-tropics, the interannual temperature variability is increased with the resolved eddies, and a notable increases in the amplitude of the El Ni?o and the Southern Oscillation is also detected. Changes in global temperature anomaly teleconnections and local air-sea feedbacks are also documented and show large changes in ocean–atmosphere coupling. In particular, local air-sea feedbacks are significantly modified by the increased ocean resolution. In the high-resolution simulation in the extra-tropics there is compelling evidence of stronger forcing of the atmosphere by SST variability arising from ocean dynamics. This coupling is very weak or absent in the low-resolution model.     

Modeling the tropical Pacific Ocean using a regional coupled climate model

A high-resolution tropical Pacific general circulation model (GCM) coupled to a global atmospheric GCM is described in this paper. The atmosphere component is the 5°×4°global general circulation model of the Institute of Atmospheric Physics (IAP) with 9 levels in the vertical direction. The ocean component with a horizontal resolution of 0.5°, is based on a low-resolution model (2°×1°in longitude-latitude).Simulations of the ocean component are first compared with its previous version. Results show that the enhanced ocean horizontal resolution allows an improved ocean state to be simulated; this involves (1) an apparent decrease in errors in the tropical Pacific cold tongue region, which exists in many ocean models,(2) more realistic large-scale flows, and (3) an improved ability to simulate the interannual variability and a reduced root mean square error (RMSE) in a long time integration. In coupling these component models, a monthly "linear-regression" method is employed to correct the model's exchanged flux between the sea and the atmosphere. A 100-year integration conducted with the coupled GCM (CGCM) shows the effectiveness of such a method in reducing climate drift. Results from years 70 to 100 are described.The model produces a reasonably realistic annual cycle of equatorial SST. The large SSTA is confined to the eastern equatorial Pacific with little propagation. Irregular warm and cold events alternate with a broad spectrum of periods between 24 and 50 months, which is very realistic. But the simulated variability is weaker than the observed and is also asymmetric in the sense of the amplitude of the warm and cold events.     

Modelling Indonesian rainfall with a coupled regional model

Long-term high-resolution coupled climate model simulations using the Max Planck Institute Regional Climate Model and the Max Planck Institute Ocean Model have been performed with boundary forcings from two reanalyses: firstly from the European Centre for Medium-Range Weather Forecasts, and secondly from the joint reanalysis of the National Centers for Environmental Prediction and the National Center for Atmospheric Research. This study employs a special coupling setup using a regional atmospheric model and a global ocean model. The latter model applies a special conformal grid from a bipolar orthogonal spherical coordinate system, which allows irregular positions of the poles and focuses on the detail over the Maritime Continent. The coupled model was able to simulate stable and realistic rainfall variabilities without flux correction and at two different ocean resolutions. The coupled system is integrated for a period between 1979 and 1993 and the results are then compared to those from uncoupled runs and from observation. The results show improved performance after coupling: a remarkable reduction of overestimated rainfall over the sea for the atmospheric model and of warm SST biases for the ocean model. There is no significant change in rainfall variability at higher ocean model resolution, but the ocean circulation shows less transport variability within the Makassar Strait in comparison to observations. This paper has not been published or considered by any other journal in any language.     

The global heat balance: heat transports in the atmosphere and ocean

The heat budget has been computed locally over the entire globe for each month of 1988 using compatible top-of-the-atmosphere radiation from the Earth Radiation Budget Experiment combined with European Centre for Medium Range Weather Forecasts atmospheric data. The effective heat sources and sinks (diabatic heating) and effective moisture sources and sinks for the atmosphere are computed and combined to produce overall estimates of the atmospheric energy divergence and the net flux through the Earth's surface. On an annual mean basis, this is directly related to the divergence of the ocean heat transport, and new computations of the ocean heat transport are made for the ocean basins. Results are presented for January and July, and the annual mean for 1988, along with a comprehensive discussion of errors. While the current results are believed to be the best available at present, there are substantial shortcomings remaining in the estimates of the atmospheric heat and moisture budgets. The issues, which are also present in all previous studies, arise from the diurnal cycle, problems with atmospheric divergence, vertical resolution, spurious mass imbalances, initialized versus uninitialized atmospheric analyses, and postprocessing to produce the atmospheric archive on pressure surfaces. Over land, additional problems arise from the complex surface topography, so that computed surface fluxes are more reliable over the oceans. The use of zonal means to compute ocean transports is shown to produce misleading results because a considerable part of the implied ocean transports is through the land. The need to compute the heat budget locally is demonstrated and results indicate lower ocean transports than in previous residual calculations which are therefore more compatible with direct ocean estimates. A Poisson equation is solved with appropriate boundary conditions of zero normal heat flux through the continental boundaries to obtain the ocean heat transport. Because of the poor observational data base, adjustments to the surface fluxes are necessary over the southern oceans. Error bars are estimated based on the large-scale spurious residuals over land of 30 W m^–2 over 1000 km scales (10¹² m²). In the Atlantic Ocean, a northward transport emerges at all latitudes with peak values of 1.1±0.2 PW (1 standard error) at 20 to 30°N. Comparable values are achieved in the Pacific at 20°N, so that the total is 2.1±0.3 PW. The peak southward transport is at 15 to 20°S of 1.9±0.3 PW made up of strong components from both the Pacific and Indian Oceans and with a heat flux from the Pacific into the Indian Ocean in the Indonesian throughflow. The pattern of poleward heat fluxes is suggestive of a strong role for Ekman transports in the tropical regions.