Extreme value statistics of the total energy in an intermediate-complexity model of the midlatitude atmospheric jet. part 1: Stationary case

A baroclinic model of intermediate complexity for the atmospheric jet at middle latitudes is used as a stochastic generator of atmosphere-like time series. In this case, time series of the total energy of the system are considered. Statistical inference of extreme values is applied to sequences of yearly maxima extracted from the time series in the rigorous setting provided by extreme value theory. The generalized extreme value (GEV) family of distributions is used here as a basic model, both for its qualities of simplicity and its generality. Several physically plausible values of the parameter T-E which represents the forced equator-to-pole temperature gradient and is responsible for setting the average baroclinicity in the atmospheric model, are used to generate stationary time series of the total energy. Estimates of the three GEV parameters-location, scale, and shape-are inferred by maximum likelihood methods. Standard statistical diagnostics, such as return level and quantile-quantile plots, are systematically applied to assess goodness-of-fit. The GEV parameters of location and scale are found to have a piecewise smooth, monotonically increasing dependence on T-E. The shape parameter also increases with T-E but is always negative, as is required a priori by the boundedness of the total energy. The sensitivity of the statistical inferences is studied with respect to the selection procedure of the maxima: the roles occupied by the length of the sequences of maxima and by the length of data blocks over which the maxima are computed are critically analyzed. Issues related to model sensitivity are also explored by varying the resolution of the system. The method used in this paper is put forward as a rigorous framework for the statistical analysis of extremes of observed data, to study the past and present climate and to characterize its variations.