Natural variability of summer rainfall over China in HadCM3

Summer rainfall over China has shown decadal variability in the past half century, which has resulted in major north–south shifts in rainfall with important implications for flooding and water resource management. This study has demonstrated how multi-century climate model simulations can be used to explore interdecadal natural variability in the climate system in order to address important questions around recent changes in Chinese summer rainfall, and whether or not anthropogenic climate change is playing a role. Using a 1,000-year simulation of HadCM3 with constant pre-industrial external forcing, the dominant modes of total and interdecadal natural variability in Chinese summer rainfall have been analysed. It has been shown that these modes are comparable in magnitude and in temporal and spatial characteristics to those observed in the latter part of the twentieth century. However, despite 1,000 years of model simulation it has not been possible to demonstrate that these modes are related to similar variations in the global circulation and surface temperature forcing occurring during the latter half of the twentieth century. This may be in part due to model biases. Consequently, recent changes in the spatial distribution of Chinese summer rainfall cannot be attributed solely to natural variability, nor has it been possible to eliminate the likelihood that anthropogenic climate change has been the driving factor. It is more likely that both play a role.     

The impact of changes in atmospheric and land surface physics on the thermohaline circulation response to anthropogenic forcing in HadCM2 and HadCM3

We have undertaken a comparative study of the mechanisms which drive the response of the Atlantic thermohaline circulation (THC) to a fourfold increase in CO₂ over 70 years with stabilisation thereafter in HadCM2 and HadCM3. In both models, the THC changes are driven by surface flux forcing, with advection (and diffusion in HadCM2) acting in the opposite sense to limit the circulation change. In both cases, heat fluxes are more important than those of freshwater. We find that different patterns of heat flux forcing in HadCM2 and HadCM3 are the prime determinants of the differing response in the two models. The increased northerly component to the near surface winds (associated with an increase in reflective low level cloud), leads to enhanced heat loss in the west-central North Atlantic, which in turn tends to steepen the steric gradient and strengthen the THC. By contrast, in HadCM3 the winds become more westerly rather than northerly, there is no dynamically-forced enhancement of surface heat loss, and the heat flux in the North Atlantic continues to be strongly positive, relative to the control, leading to a reduction in the meridional steric gradient, and a weaker overturning circulation. Differences in atmospheric response patterns appear to be caused by improvements to atmospheric and land surface physics, and suggest that the THC response in HadCM2 is less credible than in HadCM3.     

Downscaling of South America present climate driven by 4-member HadCM3 runs

The objective of this work is to evaluate climate simulations over South America using the regional Eta Model driven by four members of an ensemble of the UK Met Office Hadley Centre HadCM3 global model. The Eta Model has been modified with the purpose of performing long-term decadal integrations and has shown to reproduce “present climate”—the period 1961–1990—reasonably well when forced by HadCM3. The global model lateral conditions with a resolution of 2.5° latitude?×?3.75° longitude were provided at a frequency of 6?h. Each member of the global model ensemble has a different climate sensitivity, and the four members were selected to span the range of uncertainty encompassed by the ensemble. The Eta Model nested in the HadCM3 global model was configured with 40-km horizontal resolution and 38 layers in the vertical. No large-scale internal nudging was applied. Results are shown for austral summer and winter at present climate defined as 1961–90. The upper and low-level circulation patterns produced by the Eta-CPTEC/HadCM3 experiment set-up show good agreement with reanalysis data and the mean precipitation and temperature with CRU observation data. The spread in the downscaled mean precipitation and temperature is small when compared against model errors. On the other hand, the benefits in using an ensemble is clear in the improved representation of the seasonal cycle by the ensemble mean over any one realization. El Ni?o and La Ni?a years were identified in the HadCM3 member runs based on the NOAA Climate Prediction Center criterion of sea surface temperature anomalies in the Ni?o 3.4 area. The frequency of the El Ni?o and La Ni?a events in the studied period is underestimated by HadCM3. The precipitation and temperature anomalies typical of these events are reproduced by most of the Eta-CPTEC/HadCM3 ensemble, although small displacements of the positions of the anomalies occur. This experiment configuration is the first step on the implementation of Eta-CPTEC/HadCM3 upcoming experiments on climate change studies that are discussed in a companion paper.     

Systematic optimisation and climate simulation of FAMOUS, a fast version of HadCM3

FAMOUS is an unfluxadjusted coupled atmosphere-ocean general circulation model (AOGCM) based on the Met Office Hadley Centre AOGCM HadCM3. Its parametrisations of physical and dynamical processes are almost identical to those of HadCM3, but by virtue of reduced horizontal and vertical resolution and increased timestep it runs about ten times faster. The speed of FAMOUS means that parameter sensitivities can be investigated more thoroughly than in slower higher-resolution models, with the result that it can be tuned closer to its target climatology. We demonstrate a simple method for systematic tuning of parameters, resulting in a configuration of FAMOUS whose climatology is significantly more realistic than would be expected for a model of its resolution and speed. FAMOUS has been tuned to reproduce the behaviour of HadCM3 as nearly as possible, in order that experiments with each model are of maximum relevance to the physical interpretation of the other. Analysis of the control climate and climate change simulation of FAMOUS show that it possesses sufficient skill for its intended purposes in Earth system science as a tool for long-timescale integrations and for large ensembles of integrations, when HadCM3 cannot be afforded. Thus, it can help to bridge the gap between models of intermediate complexity and the higher-resolution AOGCMs used for policy-relevant climate prediction.     

Influence of Low-frequency Solar Forcing on the East Asian Winter Monsoon Based on HadCM3 and Observations

In this study, we investigate the influence of low-frequency solar forcing on the East Asian winter monsoon(EAWM)by analyzing a four-member ensemble of 600-year simulations performed with Had CM3(Hadley Centre Coupled Model,version 3). We find that the EAWM is strengthened when total solar irradiance(TSI) increases on the multidecadal time scale. The model results indicate that positive TSI anomalies can result in the weakening of Atlantic meridional overturning circulation, causing negative sea surface temperature(SST) anomalies in the North Atlantic. Especially for the subtropical North Atlantic, the negative SST anomalies can excite an anomalous Rossby wave train that moves from the subtropical North Atlantic to the Greenland Sea and finally to Siberia. In this process, the positive sea-ice feedback over the Greenland Sea further enhances the Rossby wave. The wave train can reach the Siberian region, and strengthen the Siberian high. As a result, low-level East Asian winter circulation is strengthened and the surface air temperature in East Asia decreases. Overall,when solar forcing is stronger on the multidecadal time scale, the EAWM is typically stronger than normal. Finally, a similar linkage can be observed between the EAWM and solar forcing during the period 1850–1970.     

Climate response to the physiological impact of carbon dioxide on plants in the Met Office Unified Model HadCM3

The concentration of carbon dioxide in the atmosphere acts to control the stomatal conductance of plants. There is observational and modelling evidence that an increase in the atmospheric concentration of CO₂ would suppress the evapotranspiration (ET) rate over land. This process is known as CO₂ physiological forcing and has been shown to induce changes in surface temperature and continental runoff. We analyse two transient climate simulations for the twenty-first century to isolate the climate response to the CO₂ physiological forcing. The land surface warming associated with the decreased ET rate is accompanied by an increase in the atmospheric lapse rate, an increase in specific humidity, but a decrease in relative humidity and stratiform cloud over land. We find that the water vapour feedback more than compensates for the decrease in latent heat flux over land as far as the budget of atmospheric water vapour is concerned. There is evidence that surface snow, water vapour and cloudiness respond to the CO₂ physiological forcing and all contribute to further warm the climate system. The climate response to the CO₂ physiological forcing has a quite different signature to that from the CO₂ radiative forcing, especially in terms of the changes in the temperature vertical profile and surface energy budget over land.     

The internal climate variability of HadCM3, a version of the Hadley Centre coupled model without flux adjustments

We examine the internal climate variability of a 1000?year long integration of the third version of the Hadley Centre coupled model (HadCM3). The model requires no flux adjustment, needs no spin up procedure prior to coupling and has a stable climate in the global mean. The principal aims are (1) to validate the internal climate variability against observed climate variability, (2) to examine the model for any periodic modes of variability, (3) to use the model estimate of internal climate variability to asses the probability of occurrence of observed trends in climate variables, and (4) to compare HadCM3 with the previous version of the Hadley Centre model, HadCM2. The magnitude and frequency characteristics of the variability of the global mean surface temperature of HadCM3 on annual to decadal time scales is in good agreement with the observations. Observed upward trends in temperature over the last 20?years and longer are inconsistent with the internal variability of the model. The simulated spatial pattern of surface temperature variability is qualitatively similar to that observed, although there is an overestimation of the land temperature variability and regional errors in ocean temperature variability. The model simulates an El Niño Southern Oscillation with an irregular 3–4?year cycle, and with a teleconnection pattern which is much more like the observations than was found in HadCM2. The interdecadal variability of the model ocean in the tropical Pacific, North Pacific and North Atlantic is broadly similar to that in the real world with none of the simulated patterns having any periodic behaviour. HadCM3 simulates an Arctic Oscillation/North Atlantic Oscillation (NAO) in Northern Hemisphere winter which has a spatial pattern consistent with the observations in the Atlantic region, but has too much teleconnection with the North Pacific. The recent observed upward trend in the NAO index is inconsistent with the model internal variability. The variability of the simulated zonal mean atmospheric temperature shows some marked differences to the observed zonal mean temperature variability, although the comparison is confounded by the sparse observational network and its possible contamination by a climate change signal.     

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     

On the consideration of mesoscale transports in climate modelling

Summary The dynamical effect of land surface heterogeneity on heat fluxes in the atmospheric boundary layer (ABL) is investigated using numerical simulations with a non-hydrostatic model over a wide range of grid resolutions. It is commonly assumed that mesoscale or dynamical fluxes associated with mesoscale and convective circulations simulated by a high-resolution model (subgrid (SG) model) on the subgrid scale of a climate model (large-scale (LS) model) represent additional processes in the ABL, which are not considered by the turbulence scheme of the LS-model, and which can be parameterized using the SG-model. The present study investigates the usefulness of this methodology for small-scale and large-scale idealized heterogeneities using a SG-model resolving mesoscale or even microscale circulations to compute the mesoscale fluxes on the scale of the LS-model. It is shown that the dynamical transports as derived from the SG-model should not be used to correct the parameterized turbulent fluxes of the LS-model. The reason is that the subgrid circulations simulated by the SG-model interact with the fields of wind and scalars in the ABL, which results in reduced turbulent fluxes in the ABL. Thus the methodology of previous studies to use mesoscale/dynamical fluxes for the correction of flux profiles simulated by climate models seems to be questionable.     

Evaluation of HadCM2 and Direct Use of Daily GCM Data in Impact Assessment Studies

This paper investigates two important aspects of methods used to explore possible effects of climatic changes on agricultural productivity on regional spatial scales. First, an evaluation of precipitation and near surface air temperature in two successive versions of the Hadley Centre General Circulation Model (GCM) has been performed to consider to what extent GCMs are capable of simulating the mean and variability of local climates. This is explored by comparing the output of an individual GCM grid box with three station observations. Several ancillary issues associated with the comparisons of observations of daily precipitation and model output that affect the statistical results are also discussed. Finally, daily data from the control and sulphate runs of the latest Hadley Centre GCM (HadCM2) have been used directly as input to the CERES-Wheat model, and the modelled yield distribution is compared to that produced with the historical data series. Our results imply that for this particular grid box covering the study region in central France, the daily raw data from HadCM2 experiment can be used directly to assess the potential impact of the greenhouse gas and sulphate aerosol radiative induced forcings and the associated climatic change on average regional winter wheat production. On the other hand, less confidence should be placed on their use regarding the estimation of future agricultural risk and variability assessment. Furthermore, a possibly more severe methodological problem that has arisen from our study is the inability of CERES-Wheat to simulate the waterlogging effects of excessive soil water on crop growth and development. Finally, we assess the potential impact of changing climate on regional winter wheat production by using the daily data from the sulphate integration up to the end of the 21st century.     

On the role of eddy mechanisms in the meridional energy transports

Summary Meridional transports of sensible and latent heat associated with standing eddies were computed from climatic mean data, and are compared with information available in the literature. The standing eddy flux of latent heat has a main maximum in the latitude of the subtropical anticyclones, where longitudinal contrasts in the latent heat content of the atmosphere are also pronounced. The standing eddy flux of sensible heat has, in winter, a main maximum in the latitude of the subpolar lows, where marked land-sea contrasts in temperature occur. Longitudinal variations in the energy content of the atmosphere account for constrasts in the latitudinal and seasonal pattern of the standing eddy fluxes of sensible and latent heat.

---

Zusammenfassung Mit ortsfesten Störmechanismen verbundene Meridionaltransporte von fühlbarer und latenter Wärme werden auf Grund klimatologischer Daten berechnet und mit den in der Literatur verfügbaren Werten verglichen. Die Störungsbewegung von latenter Wärme hat ein Hauptmaximum in der Breitenlage der subtropischen Hochdruckzellen, wo auch die zonalen Unterschiede im latenten Wärmegehalt der Luft ausgeprägt sind. Die entsprechende Störungsbewegung der fühlbaren Wärme hat im Winter ein Hauptmaximum in der Breitenlage der subpolaren Tiefdruckgebiete, wo auch starke ozeanisch-kontinentale Temperaturgegensätze auftreten. Zonale Unterschiede im Energiegehalt der Atmosphäre bedingen breitenmäßige und jahreszeitliche Gegensätze in der Störungsbewegung von fühlbarer und latenter Wärme.

---

Résumé On calcule ici au moyen de données climatologiques les transports méridiens de chaleurs latente et sensible liés à des mécanismes perturbateurs immobiles. Les valeurs ainsi obtenues sont comparées à celles que l'on trouve dans la littérature. Le mouvement perturbateur de la chaleur latente a un maximum principal dans les latitudes des cellules anticycloniques subtropicales. C'est là aussi que l'on rencontre de grandes différences zonales du contenu de l'air en chaleur latente. Le mouvement perturbateur correspondant de la chaleur sensible a, en hiver, un maximum principal dans les latitudes des dépressions subpolaires où l'on rencontre aussi de fortes oppositions de températures entre les continents et les océans. Des différences zonales du contenu de l'atmosphère en énergie provoquent des contradictions saisonnières et en latitude du mouvement perturbateur des chaleurs latente et sensible.

---

---

With 4 Figures     

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.     

Annual oceanic heat transports computed from an atmospheric model

Mean annual estimates of the oceanic poleward energy transport are obtained using a global atmospheric general circulation model. The computations are carried out by using the atmospheric model to determine the net annual heat flux into the ocean on an 8° × 10° grid. Assuming no net annual heat storage, the annual surface heat fluxes into any zonal band must be accompanied by a corresponding meridional heat transport in the ocean. Heat is transported northward at all latitudes in the Atlantic Ocean and is transported poleward in both hemispheres in the Pacific Ocean. To account for the net northward transport throughout the Atlantic, heat is transported into the Atlantic from the Indian and Pacific basins. The results are compared with several recent direct and indirect calculations of oceanic meridional heat transports.     

Mean annual poleward energy transports by the oceans in the southern hemisphere

A survey is made of the published estimates of the components of the poleward flux of energy by the atmosphere in the Southern Hemisphere in order to determine the total atmospheric transport. Together with recent measurements by satellite of the Earth's radiation budget this allows a new estimate of the required poleward energy transport by the oceans in the Southern Hemisphere for mean annual conditions. Results show that the ocean and atmosphere each contribute similar amounts for 0–30°S and that the ocean probably also transports about one third of the total at 60°S. The latter is in contrast to similar latitudes in the Northern Hemisphere where the ocean transport is negligible, but consistent with the different distribution of land and sea in the two hemispheres.     

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.     

Water mass properties and transports in the Arabian Sea from Argo observations

The information acquired from Argo floats such as temperature and salinity profiles is used to study water mass properties in the Arabian Sea from 2002 to 2004. An examination of water mass structure at different locations reveals the presence of high salinity water of marginal seas in the Arabian Sea. During the southwest monsoon season, the impact of the early onset of southwesterlies is noticed in the upper ocean temperature and salinity structure over the Western Arabian Sea (WAS) during 2002. Surface density variations are found to be more during the southwest monsoon season due to strong wind forcing. Argo temperature and salinity profiles showed that the winter cooling and the formation of Arabian Sea High Salinity Water (ASHSW) over the Northern Arabian Sea (NAS) began during the second half of November within the upper 100 m depth. In the NAS, the Persian Gulf Water (PGW) salinity is above 36, as PGW moves towards the south along isopycnal layer of 26.6σ_θ (σ_θ is potential density) salinity decreases. It is observed that the PGW high salinity water is not continuously prominent over the WAS in 2002 and in 2003. In the WAS the 27.2σ_θ isopycnal layer depth, corresponding to Red Sea Water (RSW), did not exactly follow the pattern of isotherms as is seen in the northern and eastern Arabian Sea. The variability related to RSW salinity is due to the underwater currents. The present study also confirms that RSW is prominent in the southeast Arabian Sea at the potential density of 27.2 with a maximum in summer monsoon compared to other seasons. The observed peak in the salinity at 27.2 density level during the spring intermonsoon is due to the influence of winter time spreading of RSW to the south of Socotra in 2002. Westward movement of Argo floats in the region east of Socotra during the winter is evident in both the observations and model studies. Water mass properties change when they move away from their source region due to the consistent horizontal advection. The changes in the water mass properties along the Argo float trajectory are confirmed by comparing with the climatological mean monthly values from the World Ocean Atlas 2001 data set.     

Constraints on dynamical transports of energy on a spherical planet

The factors which control the total flux of energy across a latitude belt in the atmosphere—ocean system are determined by comparing the flux resulting from various approximations with the observed flux, and by using a one-dimensional heat-balance climate model to calculate the sensitivity of the flux to the efficiency of the dynamical transports. The results show that, as long as a hemisphere is in equilibrium and as long as the structure of the atmosphere—ocean system is dominated by the planetary scale, the total flux is constrained to peak near 35° latitude, the flux per unit area to peak near 45° latitude, and the magnitude of the flux is determined primarily by the solar constant, the size of the earth, the tilt of the earth's axis, and the hemispheric mean albedo. The magnitude of the flux is insensitive to the structure and dynamics of the atmosphere—ocean system, in part because of the high efficiency of the dynamical transport mechanisms and in part because of the negative correlation between local planetary albedo and local thermal emissions to space. These results explain why the total flux and its latitudinal distribution calculated in various experiments with the GFDL general circulation model experienced relatively little modification when the hydrological cycle, mountains, or oceans were removed from the system.