The New Hadley Centre Climate Model (Hadgem1): Evaluation of Coupled Simulations
Total Page：16
File Type：pdf, Size：1020Kb
1APRIL 2006 J OHNS ET AL. 1327 The New Hadley Centre Climate Model (HadGEM1): Evaluation of Coupled Simulations T. C. JOHNS,* C. F. DURMAN,* H. T. BANKS,* M. J. ROBERTS,* A. J. MCLAREN,* J. K. RIDLEY,* ϩ C. A. SENIOR,* K. D. WILLIAMS,* A. JONES,* G. J. RICKARD, S. CUSACK,* W. J. INGRAM,# M. CRUCIFIX,* D. M. H. SEXTON,* M. M. JOSHI,* B.-W. DONG,* H. SPENCER,@ R. S. R. HILL,* J. M. GREGORY,& A. B. KEEN,* A. K. PARDAENS,* J. A. LOWE,* A. BODAS-SALCEDO,* S. STARK,* AND Y. SEARL* *Met Office, Exeter, United Kingdom ϩNIWA, Wellington, New Zealand #Met Office, Exeter, and Department of Physics, University of Oxford, Oxford, United Kingdom @CGAM, University of Reading, Reading, United Kingdom &Met Office, Exeter, and CGAM, University of Reading, Reading, United Kingdom (Manuscript received 30 May 2005, in final form 9 December 2005) ABSTRACT A new coupled general circulation climate model developed at the Met Office’s Hadley Centre is pre- sented, and aspects of its performance in climate simulations run for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4) documented with reference to previous models. The Hadley Centre Global Environmental Model version 1 (HadGEM1) is built around a new atmospheric dynamical core; uses higher resolution than the previous Hadley Centre model, HadCM3; and contains several improvements in its formulation including interactive atmospheric aerosols (sulphate, black carbon, biomass burning, and sea salt) plus their direct and indirect effects. The ocean component also has higher resolution and incorporates a sea ice component more advanced than HadCM3 in terms of both dynamics and thermodynamics. HadGEM1 thus permits experiments including some interactive processes not fea- sible with HadCM3. The simulation of present-day mean climate in HadGEM1 is significantly better overall in comparison to HadCM3, although some deficiencies exist in the simulation of tropical climate and El Niño variability. We quantify the overall improvement using a quasi-objective climate index encompassing a range of atmospheric, oceanic, and sea ice variables. It arises partly from higher resolution but also from greater fidelity in modeling dynamical and physical processes, for example, in the representation of clouds and sea ice. HadGEM1 has a similar effective climate sensitivity (2.8 K) to aCO2 doubling as HadCM3 (3.1 K), although there are significant regional differences in their response patterns, especially in the Tropics. HadGEM1 is anticipated to be used as the basis both for higher-resolution and higher-complexity Earth System studies in the near future. 1. Introduction The initial motivation for the development of HadGEM1 was to achieve better scientific performance In this paper, we present the Hadley Centre Global than the previous Hadley Centre model, HadCM3 Environmental Model version 1 (HadGEM1), a new (Pope et al. 2000; Gordon et al. 2000), delivering higher climate model developed at the Hadley Centre, and physical resolution and greater experimental flexibility. describe some aspects of its performance. HadGEM1 is HadGEM1 has at its center a new nonhydrostatic at- intended as a platform for incorporating components of mospheric dynamical core [“New Dynamics” (ND), the environmental system other than just physical cli- Davies et al. (2005)], employing a semi-implicit, semi- mate. Lagrangian time integration scheme that is also used in the operational Numerical Weather Prediction (NWP) System of the United Kingdom Meteorological Office Corresponding author address: Timothy C. Johns, Hadley Cen- tre, Met Office, FitzRoy Road, Exeter EX1 3PB, United King- (Met Office). The atmospheric physical parameteriza- dom. tion schemes in HadGEM1 have been improved con- E-mail: [email protected] siderably and remain closely in step with the NWP Unauthenticated | Downloaded 10/06/21 12:04 PM UTC JCLI3712 1328 JOURNAL OF CLIMATE VOLUME 19 model, thereby continuing the Unified Model (UM) component, the Hadley Centre Atmospheric Model strategy for NWP and climate modeling at the Met Of- version 4 (HadAM4) (Webb et al. 2001), and the fice. HadGEM1 is also designed from the outset as a coupled model, the Hadley Centre Coupled Eddy- potential community climate model for U.K. university permitting Model (HadCEM) (Roberts et al. 2004). scientists and to be suitable for further development The atmospheric physics improvements in HadAM4 toward an Earth System modeling capability, that is, were developed within HadCM3’s old dynamical core enabling coupling of physical climate with chemistry however, so building the atmospheric component of and ecosystem subcomponents. However, HadGEM1 HadGEM1 around ND necessitated recoding of the as described herein remains a conventional atmo- HadAM4 physics schemes to interface them correctly. sphere–ocean general circulation model (AOGCM) but a. Computational cost and grid resolution incorporating significant scientific changes to HadCM3, particularly with regard to its atmospheric and sea ice HadCM3 (Pope et al. 2000; Gordon et al. 2000) be- components. came operational late in 1997. Projecting a normal rate We built on established project and team-based of increase in computer power up to 2003/04 (the target model development approaches and, similarly to other for operational delivery of HadGEM1) implied that a centers (e.g., GFDL Global Atmospheric Model De- model having approximately 15–20 times the computa- velopment Team 2004), drew on scientific expertise tional cost of HadCM3 would be affordable by that spread across a large team—both within the Met Office time. This was factored in as roughly a 6–8 times in- and via external scientific contacts. Scientific perfor- crease from higher resolution, and 2–3 times from ad- mance issues identified in development were tackled by ditional physical complexity. expert teams, working over periods of a year or so and The standard atmospheric component of HadGEM1 considering both climate and NWP model configura- (HadGAM1 N96L38) discussed here uses a horizontal tions, with successive improvements being consolidated resolution of 1.25° ϫ 1.875° in latitude and longitude into the HadGEM1 prototype. The final stage of de- with 38 layers in the vertical extending to over 39 km in velopment explored parameter-based tuning to opti- height, given in Table 2 of Martin et al. (2006). (“N96” mize the skill of the model while keeping any imbal- is our shorthand denoting a resolution of 96 two-grid- ances and implied drifts in the coupled model within an length waves, that is, 192 grid points in longitude.) The acceptable tolerance. oceanic component of HadGEM1 (HadGOM1) uses a It is beyond the scope of this paper to provide a full latitude–longitude grid with a zonal resolution of 1° scientific description of HadGEM1, but in section 2 we everywhere and meridional resolution of 1° between briefly describe the main model components and some the poles and 30° latitude, from which it increases factors influencing the model design (other papers and smoothly to 1⁄3° at the equator, giving 360 ϫ 216 grid technical reports referenced therein provide more de- points in total. It has 40 unevenly spaced levels in the tails). In section 3 we outline the experiments per- vertical reducing to 10-m thickness near the surface formed with HadGEM1. In section 4 we discuss some (Table 1). key criteria used to determine acceptable model per- HadGEM1 has 8 times as many atmosphere grid formance and quantify HadGEM1 “skill” compared to points as HadCM3 (2 ϫ 2 ϫ 2 in latitude, longitude, and HadCM3 using a quasi-objective statistical method. In vertical dimensions) but only 3.75 times as many oce- sections 5–7 we evaluate aspects of HadGEM1’s con- anic grid points (1.5 ϫ 1.25 ϫ 2 in latitude, longitude, trol simulation in more detail, drawing comparisons and depth). As neither the atmosphere nor ocean time both with observed climatological fields and HadCM3 steps differ from HadCM3, the total cost increase is (discussing areas of both improvement and degradation around 6 times due to resolution factors alone. Other in model skill). In section 8 we attempt to explain changes to physics and dynamics account for an addi- changes in climate sensitivity and local feedbacks be- tional factor estimated at 2.5 times, giving an overall tween HadCM3 and HadGEM1. Section 9 is our con- cost increase in HadGEM1 of approximately 15 times cluding discussion. compared to HadCM3, consistent with the original de- sign. (This factor applies to our NEC SX6 supercom- 2. Model description puter but may depend on the technical implementa- tions and optimizations applied.) The starting points for building HadGEM1 were the The fact that ocean and atmosphere resolutions are ND dynamical core (Davies et al. 2005) plus the physi- more similar in HadGEM1 and higher than in HadCM3 cal and technical improvements made with respect to is likely to improve coupling in that more spatial detail HadCM3 in two intermediate models: the atmospheric is contained in the fields forcing the ocean. Unauthenticated | Downloaded 10/06/21 12:04 PM UTC 1APRIL 2006 J OHNS ET AL. 1329 TABLE 1. Depth and thickness of ocean model levels (centered physical schemes at such scales are much more likely to on tracer points). be necessary). Depth Thickness Depth Thickness Another key advance in HadGAM1 is the inclusion Level (m) (m) Level (m) (m) as a standard feature of the interactive modeling of 1 5.0 10.0 21 398.60 80.8 atmospheric aerosols, driven by surface and elevated 2 15.0 10.0 22 488.25 98.5 emissions (from both natural and anthropogenic 3 25.0 10.0 23 596.40 117.8 sources), including tropospheric chemical processes as 4 35.0 10.0 24 724.50 138.4 well as physical removal processes such as washout. 5 45.0 10.0 25 873.85 160.3 This scheme is coupled to the atmospheric model dy- 6 55.0 10.0 26 1045.60 183.2 7 65.0 10.0 27 1240.70 207.0 namics and physics, including radiation, cloud micro- 8 75.0 10.0 28 1459.90 231.4 physics, precipitation, and boundary layer, the aerosol 9 85.0 10.0 29 1703.70 256.2 species represented being sulphate, black carbon, bio- 10 95.0 10.0 30 1972.40 281.2 mass smoke, and sea salt (Jones et al.