a b s t r a c tLarge-scale planetary waves are diagnosed from an analysis of profiles retrieved from the Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft during its scientific mapping phase. The analysis is conducted by assimilating thermal profiles and total dust opacity retrievals into a Mars global circulation model. Transient waves are largest throughout the northern hemisphere autumn, winter and spring period and almost absent during the summer. The southern hemisphere exhibits generally weaker transient wave behaviour. A striking feature of the low-altitude transient waves in the analysis is that they show a broad subsidiary minimum in amplitude centred on the winter solstice, a period when the thermal contrast between the summer hemisphere and the winter pole is strongest and baroclinic wave activity might be expected to be strong. This behaviour, here called the 'solsticial pause,' is present in every year of the analysis. This strong pause is under-represented in many independent model experiments, which tend to produce relatively uniform baroclinic wave activity throughout the winter. This paper documents and diagnoses the transient wave solsticial pause found in the analysis; a companion paper investigates the origin of the phenomenon in a series of model experiments.

The planet Mars hosts an atmosphere that is perhaps the closest in terms of its meteorology and climate to that of the Earth. But Mars differs from Earth in its greater distance from the Sun, its smaller size, its lack of liquid oceans and its thinner atmosphere, composed mainly of CO 2 . These factors give Mars a rather different climate to that of the Earth. In this article we review various aspects of the martian climate system from a physicist's viewpoint, focusing on the processes that control the martian environment and comparing these with corresponding processes on Earth. These include the radiative and thermodynamical processes that determine the surface temperature and vertical structure of the atmosphere, the fluid dynamics of its atmospheric motions, and the key cycles of mineral dust and volatile transport. In many ways, the climate of Mars is as complicated and diverse as that of the Earth, with complex nonlinear feedbacks that affect its response to variations in external forcing. Recent work has shown that the martian climate is anything but static, but is almost certainly in a continual state of transient response to slowly varying insolation associated with cyclic variations in its orbit and rotation. We conclude with a discussion of the physical processes underlying these longterm climate variations on Mars, and an overview of some of the most intriguing outstanding problems that should be a focus for future observational and theoretical studies.

Current methods for initializing coupled atmosphere–ocean forecasts often rely on the use of separate atmosphere and ocean analyses, the combination of which can leave the coupled system imbalanced at the beginning of the forecast, potentially accelerating the development of errors. Using a series of experiments with the European Centre for Medium-Range Weather Forecasts coupled system, the magnitude and extent of these so-called initialization shocks is quantified, and their impact on forecast skill measured. It is found that forecasts initialized by separate oceanic and atmospheric analyses do exhibit initialization shocks in lower atmospheric temperature, when compared to forecasts initialized using a coupled data assimilation method. These shocks result in as much as a doubling of root-mean-square error on the first day of the forecast in some regions, and in increases that are sustained for the duration of the 10-day forecasts performed here. However, the impacts of this choice of initialization on forecast skill, assessed using independent datasets, were found to be negligible, at least over the limited period studied. Larger initialization shocks are found to follow a change in either the atmosphere or ocean model component between the analysis and forecast phases: changes in the ocean component can lead to sea surface temperature shocks of more than 0.5 K in some equatorial regions during the first day of the forecast. Implications for the development of coupled forecast systems, particularly with respect to coupled data assimilation methods, are discussed.

The redistribution of a finite amount of Martian surface dust during global dust storms and in the intervening periods has been modelled in a dust lifting version of the UK Mars General Circulation Model. When using a constant, uniform threshold in the model's wind stress lifting parameterisation and assuming an unlimited supply of surface dust, multiannual simulations displayed some variability in dust lifting activity from year to year, arising from internal variability manifested through surface wind stress, but dust storms were limited in size and formed within a relatively short seasonal window. Lifting thresholds were then allowed to vary at each model gridpoint, dependent on the rates of emission or deposition of dust. This enhanced interannual variability in dust storm magnitude and timing, such that model storms could cover most of the observed ranges in size and initiation date within a single multiannual simulation. Peak storm magnitude in a given year was primarily determined by the availability of surface dust at a number of key sites in the southern hemisphere. The observed global dust storm (GDS) frequency of roughly one in every three years was approximately reproduced, but the model failed to generate these GDSs spontaneously in the southern hemisphere, where they have typically been observed to initiate. After several years of simulation, the surface dust density field showed good qualitative agreement with the observed pattern of Martian surface dust cover. The model produced a net northward cross-equatorial dust mass flux, which necessitated the addition of an artificial threshold decrease rate in order to allow the continued generation of dust storms over the course of a multiannual simulation. At standard model resolution, for the southward mass flux due to cross-equatorial flushing storms to offset the northward flux due to GDSs on a timescale of ∼ 3 3 years would require an increase in the former by a factor of 3-4. Results at higher model resolution and uncertainties in dust vertical profiles mean that this nevertheless appears to be a plausible explanation for the observed GDS frequency.