We describe the development, tuning and climate of Planet Simulator (PLASIM)–Grid-ENabled Integrated Earth system model (GENIE), a new intermediate complexity Atmosphere–Ocean General Circulation Model (AOGCM), built by coupling the Planet Simulator to the ocean, sea-ice and land-surface components of the GENIE Earth system model. PLASIM–GENIE supersedes GENIE2, a coupling of GENIE to the Reading Intermediate General Circulation Model (IGCM). The primitive-equation atmosphere includes chaotic, three-dimensional (3-D) motion and interactive radiation and clouds, and dominates the computational load compared to the relatively simpler frictionalgeostrophic ocean, which neglects momentum advection. The model is most appropriate for long-timescale or large ensemble studies where numerical efficiency is prioritised, but lack of data necessitates an internally consistent, coupled calculation of both oceanic and atmospheric fields. A 1000-year simulation with PLASIM–GENIE requires approximately 2 weeks on a single node of a 2.1 GHz AMD 6172 CPU. We demonstrate the tractability of PLASIM–GENIE ensembles by deriving a subjective tuning of the model with a 50member ensemble of 1000-year simulations. The simulated climate is presented considering (i) global fields of seasonal surface air temperature, precipitation, wind, solar and thermal radiation, with comparisons to reanalysis data; (ii) vegetation carbon, soil moisture and aridity index; and (iii) sea surface temperature, salinity and ocean circulation. Considering its resolution, PLASIM–GENIE reproduces the main features of the climate system well and demonstrates usefulness for a wide range of applications.

Fundamental to the redistribution of energy in a planetary atmosphere is the general circulation and its meridional structure. We use a general circulation model of the atmosphere in an aquaplanet configuration with prescribed sea surface temperature and investigate the influence of the gravitational acceleration g on the structure of the circulation. For g=g0=9.81ms−2, three meridional cells exist in each hemisphere. Up to about g/g0=1.4 all cells increase in strength. Further increasing this ratio results in a weakening of the thermally indirect cell, such that a two- and finally a one-cell structure of the meridional circulation develops in each hemisphere. This transition is explained by the primary driver of the thermally direct Hadley cell: the diabatic heating at the equator which is proportional to g. The analysis of the energetics of the atmospheric circulation based on the Lorenz energy cycle supports this finding. For Earth-like gravitational accelerations transient eddies are primarily responsible for the meridional heat flux. For large gravitational accelerations, the direct zonal mean conversion of energy dominates the meridional heat flux.

We describe the development, tuning and climate of Planet Simulator (PLASIM)–Grid-ENabled Integrated Earth system model (GENIE), a new intermediate complexity Atmosphere–Ocean General Circulation Model (AOGCM), built by coupling the Planet Simulator to the ocean, sea-ice and land-surface components of the GENIE Earth system model. PLASIM–GENIE supersedes GENIE-2, a coupling of GENIE to the Reading Intermediate General Circulation Model (IGCM). The primitive-equation atmosphere includes chaotic, three-dimensional (3-D) motion and interactive radiation and clouds, and dominates the computational load compared to the relatively simpler frictional-geostrophic ocean, which neglects momentum advection. The model is most appropriate for long-timescale or large ensemble studies where numerical efficiency is prioritised, but lack of data necessitates an internally consistent, coupled calculation of both oceanic and atmospheric fields. A 1000-year simulation with PLASIM–GENIE requires approximately 2 weeks on a single node of a 2.1 GHz AMD 6172 CPU. We demonstrate the tractability of PLASIM–GENIE ensembles by deriving a subjective tuning of the model with a 50-member ensemble of 1000-year simulations. The simulated climate is presented considering (i) global fields of seasonal surface air temperature, precipitation, wind, solar and thermal radiation, with comparisons to reanalysis data; (ii) vegetation carbon, soil moisture and aridity index; and (iii) sea surface temperature, salinity and ocean circulation. Considering its resolution, PLASIM–GENIE reproduces the main features of the climate system well and demonstrates usefulness for a wide range of applications.

Introduction Conclusions References Tables Figures

Many applications in the evaluation of climate impacts and environmental policy require detailed spatio-temporal projections of future climate. To capture feedbacks from impacted natural or socio-economic systems requires interactive two-way coupling, but this is generally computationally infeasible with even moderately complex general circulation models (GCMs). Dimension reduction using emulation is one solution to this problem, demonstrated here with the GCM PLASIM-ENTS (Planet Simulator coupled with the efficient numerical terrestrial scheme). Our approach generates temporally evolving spatial patterns of climate variables, considering multiple modes of variability in order to capture non-linear feedbacks. The emulator provides a 188-member ensemble of decadally and spatially resolved (~ 5° resolution) seasonal climate data in response to an arbitrary future CO2 concentration and non-CO2 radiative forcing scenario. We present the PLASIM-ENTS coupled model, the construction of its emulator from an ensemble of transient future simulations, an application of the emulator methodology to produce heating and cooling degree-day projections, the validation of the simulator (with respect to empirical data) and the validation of the emulator (with respect to high-complexity models). We also demonstrate the application to estimates of sea-level rise and associated uncertainty.

We explore the potential multistability of the climate for a planet around the habitable zone. We focus on conditions reminiscent to those of the Earth system, but our investigation aims at presenting a general methodology for dealing with exoplanets. We provide a thorough analysis of the non-equilibrium thermodynamical properties of the climate system and explore, using a a flexible climate model, how such properties depend on the energy input of the parent star, on the infrared atmospheric opacity, and on the rotation rate. It is possible to reproduce the multi-stability properties reminiscent of the paleoclimatologically relevant snowball (SB) - warm (W) conditions. We then study the thermodynamics of the W and SB states, clarifying the role of the hydrological cycle in shaping the irreversibility and the efficiency of the W states, and emphasizing the extreme diversity of the SB states, where dry conditions are realized. Thermodynamics provides the clue for studying the tipping points of the system and leads us to constructing parametrizations where the main thermodynamic properties are expressed as functions of the emission temperature of the planet only. Such functions are rather robust with respect to changing the rotation rate of the planet from the current terrestrial one to half of it. We then explore the dynamical range of slowy rotating and phase locked planets. There is a critical rotation rate below which the multi-stability properties are lost. Such critical rotation rate corresponds roughly to the phase lock 2:1 condition. Therefore, if an Earth-like planet is 1:1 phase locked with respect to the parent star, only one climatic state would be compatible with a given set of astronomical and astrophysical parameters. These results have relevance for the general theory of planetary circulation and for the definition of necessary and sufficient conditions for habitability.

We have applied the relation for the mean of the expected values of the maximum excursion in a bounded random walk to estimate the random walk length from time series of eight independent global mean quantities (temperature maximum, summer lag, temperature minimum and winter lag in land and over the ocean) derived from the NOAA-CIRES twentieth century reanalysis (V2) for 1871-2008 and the ECHAM5 IPCC AR4 twentieth century run for 1860-2100, and also the Millenium 3100 yr control run mil01, which was segmented into records of specified period. The results for NOAA-CIRES, ECHAM5, and mil01 (mean of thirty 100 yr segments) are very similar and indicate a random walk length on land of 24 yr and over the ocean of 20 yr. Using three 1000 yr segments from mil01, the random walk lengths increased to 37 yr on land and 33 yr over the ocean. This result indicates that the shorter records may not totally capture the random variability of climate relevant on the time scale of civilizations, for which the random walk length is likely to be about 30 years. For this random walk length, the observed standard deviations of maximum temperature and minimum temperature yield respective expected maximum excursions on land of 1.4 and 2.3 degrees C and over the ocean of 0.5 and 0.7 degrees C, which are substantial fractions of the global warming signal. Citation: Bye, J., K. Fraedrich, E. Kirk, S. Schubert, and X. Zhu (2011), Random walk lengths of about 30 years in global climate, Geophys. Res. Lett., 38, L05806, doi: 10.1029/2010GL046333.

The impact of mountains and ice sheets on the large-scale circulation of the world's oceans is investigated in a series of simulations with a new coupled ocean atmosphere model [Oregon State University University of Victoria model (OSUVic)], in which the height of orography is scaled from 1.5 times the actual height (at T42 resolution) to 0 (no mountains). The results suggest that the effects of mountains and ice sheets on the buoyancy and momentum transfer from the atmosphere to the surface ocean determine the present pattern of deep ocean circulation. Higher mountains reduce water vapor transport from the Pacific and Indian Oceans into the Atlantic Ocean and contribute to increased (decreased) salinities and enhanced (reduced) deepwater formation and meridional overturning circulation in the Atlantic (Pacific). Orographic effects also lead to the observed interhemispheric asymmetry of midlatitude zonal wind stress. The presence of the Antarctic ice sheet cools winter air temperatures by more than 20 degrees C directly above the ice sheet and sets up a polar meridional overturning cell in the atmosphere. The resulting increased meridional temperature gradient strengthens midlatitude westerlies by similar to 25% and shifts them poleward by similar to 10 degrees. This leads to enhanced and poleward-shifted upwelling of deep waters in the Southern Ocean, a stronger Antarctic Circumpolar Current, increased poleward atmospheric moisture transport, and more advection of high-salinity Indian Ocean water into the South Atlantic. Thus, it is the current configuration of mountains and ice sheets on earth that determines the difference in deep-water formation between the Atlantic and the Pacific.

The Planet Simulator is a Model of Intermediate Complexity (MIC), suitable for climate and paleo-climate simulations and time scales up to ten thousand years or more. It runs on a wide range of hardware and operating systems including massive parallel cluster and workstations with UNIX, Solaris, Linux, MAC OS or any other UNIX style operating system. The priorities in development are set to speed, easy handling and portability. Its modular structure allows a problem dependent configuration. A graphical Model Starter (MoSt) can be used to select a model configuration from the repository, set its parameters, compile and run the model. The model can be run either in production mode for maximum performance or in interactive mode using a Graphical User Interface (GUI). This GUI is both a real time visualization interface for all model variables and a tool for tuning and experimenting with the model.

The Planet Simulator is a Model of Intermediate Complexity (MIC), suitable for climate and paleo-climate simulations and time scales up to ten thousand years or more. It runs on a wide range of hardware and operating systems including massive parallel cluster and workstations with UNIX, Solaris, Linux, MAC OS or any other UNIX style operating system. The priorities in development are set to speed, easy handling and portability. Its modular structure allows a problem dependent configuration. A graphical Model Starter (MoSt) can be used to select a model configuration from the repository, set its parameters, compile and run the model. The model can be run either in production mode for maximum performance or in interactive mode using a Graphical User Interface (GUI). This GUI is both a real time visualization interface for all model variables and a tool for tuning and experimenting with the model.

Two kinds of objective functions for parameter optimisation in simplified general circulation models (SGCMs) are introduced and tested with an SGCM employing linear parameterisations for diabatic heating, surface friction and horizontal diffusion. (a) A set of circulation indices is introduced to characterise the zonal mean primary and secondary circulation and the global energetics. The objective function is then given by the distance between the modelled and a reference (e.g. observed) circulation in a state space spanned by these indices. (b) The global and time mean entropy production and kinetic energy dissipation are introduced as additional objective functions, following the maximum entropy production principle. It is found that both methods lead to optimal parameter values close to the standard configuration of the model, though the method of the second kind is restricted to those model parameters associated with internal processes such as heat and momentum fluxes.

Boundary layer turbulence plays a central role in determining the strength of the overall atmospheric circulation by affecting the intensity of exchange of heat, mass, and momentum at the Earth's surface. It is often parameterized using the bulk formula, in which the von‐Karman parameter plays a critical role. Here we conducted a range of sensitivity simulations with an atmospheric general circulation model in which we modified the strength of boundary layer turbulence by varying the von‐Karman parameter. These simulations show that the maximum of entropy production associated with boundary layer dissipation is consistent with the observed value of the von‐Karman parameter of 0.4 and maximizes the planetary rate of entropy production with the global radiative temperature being close to its minimum value. Additional sensitivity simulations were conducted with an increased concentration of atmospheric carbon dioxide, which affects the relative radiative forcing of tropical vs. polar regions. We find that the global climate sensitivity is more‐or‐less independent of the assumed strength of boundary layer turbulence in our idealized setup. The difference in climate sensitivities of tropical and polar regions is at a minimum at a climatic state of MEP.

An application of the Planet Simulator is presented to estimate the maximum effect of vegetation on the Earth’s climate. Four sets of sensitivity experiments are performed: (1) All vegetation related land surface parameters are changed simultaneously. (2) Only one effect of vegetation on climate is considered: albedo, surface roughness and soil hydrology. To identify the nature of vegetation-climate interaction, linear superposition and non-linear interaction of these three effects are compared. (3) The first experimental set-up is repeated but with mixed-layer ocean and thermodynamic sea-ice. (4) The effect of enhanced greenhouse gas concentrations on extreme vegetation climates is analysed repeating the preceding experimental setup with twice the CO2 concentration (compared to set-3). Zusammenfassung Der Planet Simulator wird hier benutzt, um den maximalen Effekt der Vegetation auf das Erdklima

The ability of analysing atmospheric dynamics by idealized experiments using a simplified circulation model is illustrated in three related studies. The investigations focus on the organization of localized strom tracks, their impact on low-frequency variability, and on the response to external thermal forcing. A localized storm track in agreement with observations is forced by a heating dipole embedded in a zonally symmetric field if the dipole orientation corresponds to the Northern Hemisphere winter case. The interaction of two storm tracks leads to low-frequency variability. Spatial resonance between a low-frequency large scale retrograde travelling wave and the storm track eddies is identified causing the fluctuations. Teleconnection pattern remarkably similar to the observed North Atlantic Oscillation (NAO) and Pacific North American (PNA) pattern emerge when the distance of the two storm tracks is set to the observed value of about 150 degrees. While the spatial resonance mechanism forces the NAO-like pattern, baroclinic processes are related to the PNA-like teleconnection. Anomalies induced by large scale thermal forcing strongly depend on the background flow. A non-linear response is observed in the model depending on the sign of the forcing and its position relative to the storm track. A baroclinic and an equivalent barotropic component defines the response. In addition, the change of the space-time variability is affected by eddy feedbacks.

The Planet Simulator is a Model of Intermediate Complexity (MIC) which can be used to run climate and paleo-climate simulations for time scales up to 10 thousand years or more in an acceptable real time. The priorities in development are set to speed, easy handling and portability. Its modular structure allows a problem dependent configuration. Adaptions exist for the atmospheres of Mars and of Saturn's moon Titan. Common coupling interfaces enable the addition of ocean models, ice models, vegetation and more. An interactive mode with a Model Starter (MoSt) and a Graphical User Interface (GUI) can be used to select a model configuration from the available hierarchy, set its parameters and inspect atmospheric fields while changing model parameters on the fly: This is especially useful for teaching, debugging and tuning of parameterizations. This paper gives an overview of the model's features. The complete model including sources and documentation is available at (www.mi.uni-hamburg.de/plasim).