No abstract

International audiencePrevious studies using coupled general circulation models (GCMs) suggest that the atmosphere model plays a dominant role in the modeled El Niño–Southern Oscillation (ENSO), and that intermodel differences in the thermodynamical damping of sea surface temperatures (SSTs) are a dominant contributor to the ENSO amplitude diversity. This study presents a detailed analysis of the shortwave flux feedback (αSW) in 12 Coupled Model Intercomparison Project phase 3 (CMIP3) simulations, motivated by findings that αSW is the primary contributor to model thermodynamical damping errors.A “feedback decomposition method,” developed to elucidate the αSW biases, shows that all models underestimate the dynamical atmospheric response to SSTs in the eastern equatorial Pacific, leading to underestimated αSW values. Biases in the cloud response to dynamics and the shortwave interception by clouds also contribute to errors in αSW. Changes in the αSW feedback between the coupled and corresponding atmosphere-only simulations are related to changes in the mean dynamics.A large nonlinearity is found in the observed and modeled SW flux feedback, hidden when linearly calculating αSW. In the observations, two physical mechanisms are proposed to explain this nonlinearity: 1) a weaker subsidence response to cold SST anomalies than the ascent response to warm SST anomalies and 2) a nonlinear high-level cloud cover response to SST. The shortwave flux feedback nonlinearity tends to be underestimated by the models, linked to an underestimated nonlinearity in the dynamical response to SST. The process-based methodology presented in this study may help to correct model ENSO atmospheric biases, ultimately leading to an improved simulation of ENSO in GCMs

Several studies using ocean-atmosphere GCMs suggest that the atmospheric component plays a dominant role in the modelled ENSO. To help elucidate these findings, the two main atmosphere feedbacks relevant to ENSO, the Bjerknes positive feedback (µ) and the heat flux negative feedback (α), are analysed here in 12 coupled GCMs.We find that the models generally underestimate both feedbacks, leading to an error compensation. The strength of α is inversely related to the ENSO amplitude in the models and the latent heat and shortwave flux components of this feedback dominate. Furthermore, the shortwave component could help explain the model diversity in both overall α and ENSO amplitude.

Currently, most operational forecasting models use latitude-longitude grids, whose convergence of meridians toward the poles limits parallel scaling. Quasi-uniform grids might avoid this limitation. Thuburn et al. and Ringler et al. have developed a method for arbitrarily structured, orthogonal C grids called TRiSK, which has many of the desirable properties of the C grid on latitude-longitude grids but which works on a variety of quasi-uniform grids. Here, five quasi-uniform, orthogonal grids of the sphere are investigated using TRiSK to solve the shallow-water equations.Some of the advantages and disadvantages of the hexagonal and triangular icosahedra, a ''Voronoi-ized'' cubed sphere, a Voronoi-ized skipped latitude-longitude grid, and a grid of kites in comparison to a full latitude-longitude grid are demonstrated. It is shown that the hexagonal icosahedron gives the most accurate results (for least computational cost). All of the grids suffer from spurious computational modes; this is especially true of the kite grid, despite it having exactly twice as many velocity degrees of freedom as height degrees of freedom. However, the computational modes are easiest to control on the hexagonal icosahedron since they consist of vorticity oscillations on the dual grid that can be controlled using a diffusive advection scheme for potential vorticity.

A new theoretical framework is derived for parameterization of subgrid physical processes in atmospheric models; the application to parameterization of convection and boundary layer fluxes is a particular focus. The derivation is based on conditional filtering, which uses a set of quasi-Lagrangian labels to pick out different regions of the fluid, such as convective updrafts and environment, before applying a spatial filter. This results in a set of coupled prognostic equations for the different fluid components, including subfilter-scale flux terms and entrainment/detrainment terms. The framework can accommodate different types of approaches to parameterization, such as local turbulence approaches and mass flux approaches. It provides a natural way to distinguish between local and nonlocal transport processes and makes a clearer conceptual link to schemes based on coherent structures such as convective plumes or thermals than the straightforward application of a filter without the quasi-Lagrangian labels. The framework should facilitate the unification of different approaches to parameterization by highlighting the different approximations made and by helping to ensure that budgets of energy, entropy, and momentum are handled consistently and without double counting. The framework also points to various ways in which traditional parameterizations might be extended, for example, by including additional prognostic variables. One possibility is to allow the large-scale dynamics of all the fluid components to be handled by the dynamical core. This has the potential to improve several aspects of convection-dynamics coupling, such as dynamical memory, the location of compensating subsidence, and the propagation of convection to neighboring grid columns.