The ASSET intercomparison of ozone analyses: method and first results

This paper aims to summarise the current performance of ozone data assimilation (DA) systems, to show where they can be improved, and to quantify their errors. It examines 11 sets of ozone analyses from 7 different DA systems. Two are numerical weather prediction (NWP) systems based on general circulation models (GCMs); the other five use chemistry transport models (CTMs). The systems examined contain either linearised or detailed ozone chemistry, or no chemistry at all. In most analyses, MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) ozone data are assimilated; two assimilate SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Chartography) observations instead. Analyses are compared to independent ozone observations covering the troposphere, stratosphere and lower mesosphere during the period July to November 2003. Biases and standard deviations are largest, and show the largest divergence between systems, in the troposphere, in the upper-troposphere/lower-stratosphere, in the upper-stratosphere and mesosphere, and the Antarctic ozone hole region. However, in any particular area, apart from the troposphere, at least one system can be found that agrees well with independent data. In general, none of the differences can be linked to the assimilation technique (Kalman filter, three or four dimensional variational methods, direct inversion) or the system (CTM or NWP system). Where results diverge, a main explanation is the way ozone is modelled. It is important to correctly model transport at the tropical tropopause, to avoid positive biases and excessive structure in the ozone field. In the southern hemisphere ozone hole, only the analyses which correctly model heterogeneous ozone depletion are able to reproduce the near-complete ozone destruction over the pole. In the upper-stratosphere and mesosphere (above 5 hPa), some ozone photochemistry schemes caused large but easily remedied biases. The diurnal cycle of ozone in the mesosphere is not captured, except by the one system that includes a detailed treatment of mesospheric chemistry. These results indicate that when good observations are available for assimilation, the first priority for improving ozone DA systems is to improve the models. The analyses benefit strongly from the good quality of the MIPAS ozone observations. Using the analyses as a transfer standard, it is seen that MIPAS is similar to 5% higher than HALOE (Halogen Occultation Experiment) in the mid and upper stratosphere and mesosphere (above 30 hPa), and of order 10% higher than ozonesonde and HALOE in the lower stratosphere (100 hPa to 30 hPa). Analyses based on SCIAMACHY total column are almost as good as the MIPAS analyses; analyses based on SCIAMACHY limb profiles are worse in some areas, due to problems in the SCIAMACHY retrievals.

Predictability and skill of boreal winter forecasts made with the ECMWF Seasonal Forecasting System II

The performance of boreal winter forecasts made with the European Centre for Medium-Range Weather Forecasts (ECMWF) System 11 Seasonal Forecasting System is investigated through analyses of ensemble hindcasts for the period 1987-2001. The predictability, or signal-to-noise ratio, associated with the forecasts, and the forecast skill are examined. On average, forecasts of 500 hPa geopotential height (GPH) have skill in most of the Tropics and in a few regions of the extratropics. There is broad, but not perfect, agreement between regions of high predictability and regions of high skill. However, model errors are also identified, in particular regions where the forecast ensemble spread appears too small. For individual winters the information provided by t-values, a simple measure of the forecast signal-to-noise ratio, is investigated. For 2 m surface air temperature (T2m), highest t-values are found in the Tropics but there is considerable interannual variability, and in the tropical Atlantic and Indian basins this variability is not directly tied to the El Nino Southern Oscillation. For GPH there is also large interannual variability in t-values, but these variations cannot easily be predicted from the strength of the tropical sea-surface-temperature anomalies. It is argued that the t-values for 500 hPa GPH can give valuable insight into the oceanic forcing of the atmosphere that generates predictable signals in the model. Consequently, t-values may be a useful tool for understanding, at a mechanistic level, forecast successes and failures. Lastly, the extent to which t-values are useful as a predictor of forecast skill is investigated. For T2m, t-values provide a useful predictor of forecast skill in both the Tropics and extratropics. Except in the equatorial east Pacific, most of the information in t-values is associated with interannual variability of the ensemble-mean forecast rather than interannual variability of the ensemble spread. For GPH, however, t-values provide a useful predictor of forecast skill only in the tropical Pacific region.

Dynamical evolution of the 2003 southern hemisphere stratospheric winter using Envisat trace-gas observations

Data from the MIPAS instrument on Envisat, supplemented by meteorological analyses from ECMWF and the Met Office, are used to study the meteorological and trace-gas evolution of the stratosphere in the southern hemisphere during winter and spring 2003. A pole-centred approach is used to interpret the data in the physically meaningful context of the evolving stratospheric polar vortex. The following salient dynamical and transport features are documented and analysed: the merger of anticyclones in the stratosphere; the development of an intense, quasi-stationary anticyclone in spring; the associated top-down breakdown of the polar vortex; the systematic descent of air into the polar vortex; and the formation of a three-dimensional structure of a tracer filament on a planetary scale. The paper confirms and extends existing paradigms of the southern hemisphere vortex evolution. The quality of the MIPAS observations is seen to be generally good. though the water vapour retrievals are unrealistic above 10 hPa in the high-latitude winter.

North Atlantic Oscillation response to transient greenhouse gas forcing and the impact on European winter climate: a CMIP2 multi-model assessment

This study investigates the response of wintertime North Atlantic Oscillation (NAO) to increasing concentrations of atmospheric carbon dioxide (CO2) as simulated by 18 global coupled general circulation models that participated in phase 2 of the Coupled Model Intercomparison Project (CMIP2). NAO has been assessed in control and transient 80-year simulations produced by each model under constant forcing, and 1% per year increasing concentrations of CO2, respectively. Although generally able to simulate the main features of NAO, the majority of models overestimate the observed mean wintertime NAO index of 8 hPa by 5-10 hPa. Furthermore, none of the models, in either the control or perturbed simulations, are able to reproduce decadal trends as strong as that seen in the observed NAO index from 1970-1995. Of the 15 models able to simulate the NAO pressure dipole, 13 predict a positive increase in NAO with increasing CO2 concentrations. The magnitude of the response is generally small and highly model-dependent, which leads to large uncertainty in multi-model estimates such as the median estimate of 0.0061 +/- 0.0036 hPa per %CO2. Although an increase of 0.61 hPa in NAO for a doubling in CO2 represents only a relatively small shift of 0.18 standard deviations in the probability distribution of winter mean NAO, this can cause large relative increases in the probabilities of extreme values of NAO associated with damaging impacts. Despite the large differences in NAO responses, the models robustly predict similar statistically significant changes in winter mean temperature (warmer over most of Europe) and precipitation (an increase over Northern Europe). Although these changes present a pattern similar to that expected due to an increase in the NAO index, linear regression is used to show that the response is much greater than can be attributed to small increases in NAO. NAO trends are not the key contributor to model-predicted climate change in wintertime mean temperature and precipitation over Europe and the Mediterranean region. However, the models' inability to capture the observed decadal variability in NAO might also signify a major deficiency in their ability to simulate the NAO-related responses to climate change.

Assimilation of stratospheric ozone from MIPAS into a global general-circulation model: The September 2002 vortex split

Ozone and temperature profiles from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) have been assimilated, using three-dimensional variational assimilation, into a stratosphere troposphere version of the Met Office numerical weather-prediction system. Analyses are made for the month of September 2002, when there was an unprecedented split in the southern hemisphere polar vortex. The analyses are validated against independent ozone observations from sondes, limb-occultation and total column ozone satellite instruments. Through most of the stratosphere, precision varies from 5 to 15%, and biases are 15% or less of the analysed field. Problems remain in the vortex and below the 60 hPa. level, especially at the tropopause where the analyses have too much ozone and poor agreement with independent data. Analysis problems are largely a result of the model rather than the data, giving confidence in the MIPAS ozone retrievals, though there may be a small high bias in MIPAS ozone in the lower stratosphere. Model issues include an excessive Brewer-Dobson circulation, which results both from known problems with the tracer transport scheme and from the data assimilation of dynamical variables. The extreme conditions of the vortex split reveal large differences between existing linear ozone photochemistry schemes. Despite these issues, the ozone analyses are able to successfully describe the ozone hole split and compare well to other studies of this event. Recommendations are made for the further development of the ozone assimilation system.

Equatorial waves in the lower stratosphere. II: Annual and interannual variability

The annual and interannual variability of idealized, linear, equatorial waves in the lower stratosphere is investigated using the temperature and velocity fields from the ECMWF 15-year re-analysis dataset. Peak Kelvin wave activity occurs during solstice seasons at 100 hPa, during December-February at 70 hPa and in the easterly to westerly quasi-biennial oscillation (QBO) phase transition at 50 hPa. Peak Rossby-gravity wave activity occurs during equinox seasons at 100 hPa, during June-August/September-November at 70 hPa and in the westerly to easterly QBO phase transition at 50 hPa. Although neglect of wind shear means that the results for inertio-gravity waves are likely to be less accurate, they are still qualitatively reasonable and an annual cycle is observed in these waves at 100 hPa and 70 hPa. Inertio-gravity waves with n = 1 are correlated with the QBO at 50 hPa, but the eastward inertio-gravity n = 0 wave is not, due to its very fast vertical group velocity in all background winds. The relative importance of different wave types in driving the QBO at 50 hPa is also discussed. The strongest acceleration appears to be provided by the Kelvin wave while the acceleration provided by the Rossby-gravity wave is negligible. Of the higher-frequency waves, the westward inertio-gravity n = 1 wave appears able to contribute more to the acceleration of the 50 hPa mean zonal wind than the eastward inertio-gravity n = 1 wave.

The influence of physical processes on extratropical singular vectors

An investigation is made of the impact of a full linearized physical (moist) parameterization package on extratropical singular vectors (SVs) using the ECMWF integrated forecasting system (IFS). Comparison is made for one particular period with a dry physical package including only vertical diffusion and surface drag. The crucial extra ingredient in the full package is found to be the large-scale latent heat release. Consistent with basic theory, its inclusion results in a shift to smaller horizontal scales and enhanced growth for the SVs. Whereas, for the dry SVs, T42 resolution is sufficient, the moist SVs require T63 to resolve their structure and growth. A 24-h optimization time appears to be appropriate for the moist SVs because of the larger growth of moist SVs compared with dry SVs. Like dry SVs, moist SVs tend to occur in regions of high baroclinicity, but their location is also influenced by the availability of moisture. The most rapidly growing SVs appear to enhance or reduce large-scale rain in regions ahead of major cold fronts. The enhancement occurs in and ahead of a cyclonic perturbation and the reduction in and ahead of an anticyclonic perturbation. Most of the moist SVs for this situation are slightly modified versions of the dry SVs. However, some occur in new locations and have particularly confined structures. The most rapidly growing SV is shown to exhibit quite linear behavior in the nonlinear model as it grows from 0.5 to 12 hPa in 1 day. For 5 times this amplitude the structure is similar but the growth is about half as the perturbation damps a potential vorticity (PV) trough or produces a cutoff, depending on its sign.

Equatorial waves in the lower stratosphere. I: A novel detection method

Wavenumber-frequency spectral analysis and linear wave theory are combined in a novel method to quantitatively estimate equatorial wave activity in the tropical lower stratosphere. The method requires temperature and velocity observations that are regularly spaced in latitude, longitude and time; it is therefore applied to the ECMWF 15-year re-analysis dataset (ERA-15). Signals consistent with idealized Kelvin and Rossby-gravity waves are found at wavenumbers and frequencies in agreement with previous studies. When averaged over 1981-93, the Kelvin wave explains approximately 1 K-2 of temperature variance on the equator at 100 hPa, while the Rossby-gravity wave explains approximately 1 m(2)s(-2) of meridional wind variance. Some inertio-gravity wave and equatorial Rossby wave signals are also found; however the resolution of ERA-15 is not sufficient for the method to provide an accurate climatology of waves with high meridional structure.

The formation mechanism of the Bonin high in August

The Bonin high is a subtropical anticyclone that is predominant near Japan in the summer. This anticyclone is associated with an equivalent-barotropic structure, often extending throughout the entire troposphere. Although the equivalent-barotropic structure of the Bonin high has been known for years among synopticians because of its importance to the summer climate in east Asia, there are few dynamical explanations for such a structure. The present paper attempts to provide a formation mechanism for the deep ridge near Japan. We propose a new hypothesis that this equivalent-barotropic ridge near Japan is formed as a result of the propagation of stationary Rossby waves along the Asian jet in the upper troposphere (‘the Silk Road pattern’). First, the monthly mean climatology is examined in order to demonstrate this hypothesis. It is shown that the enhanced Asian jet in August is favourable for the propagation of stationary Rossby waves and that the regions of descent over the eastern Mediterranean Sea and the Aral Sea act as two major wave sources. Second, a primitive-equation model is used to simulate the climatology of August. The model successfully simulates the Bonin high with an equivalent-barotropic structure. The upper-tropospheric ridge is found to be enhanced by a height anomaly of more than 80 m at 200 hPa, when a wave packet arrives. Sensitivity experiments are conducted to show that the removal of the diabatic cooling over the Asian jet suppresses the Silk Road pattern and formation of an equivalent-barotropic ridge near Japan, while the removal of the diabatic heating in the western Pacific does not. Copyright © 2003 Royal Meteorological Society

Evaluation of ozone and water vapor fields from the ECMWF reanalysis ERA-40 during 1991-1999 in comparison with UARS satellite and MOZAIC aircraft observations

[ 1] The European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year Reanalysis (ERA-40) ozone and water vapor reanalysis fields during the 1990s have been compared with independent satellite data from the Halogen Occultation Experiment (HALOE) and Microwave Limb Sounder (MLS) instruments on board the Upper Atmosphere Research Satellite (UARS). In addition, ERA-40 has been compared with aircraft data from the Measurements of Ozone and Water Vapour by Airbus In-Service Aircraft (MOZAIC) program. Overall, in comparison with the values derived from the independent observations, the upper stratosphere in ERA-40 has about 5 - 10% more ozone and 15 - 20% less water vapor. This dry bias in the reanalysis appears to be global and extends into the middle stratosphere down to 40 hPa. Most of the discrepancies and seasonal variations between ERA-40 and the independent observations occur within the upper troposphere over the tropics and the lower stratosphere over the high latitudes. ERA-40 reproduces a weaker Antarctic ozone hole, and of less vertical extent, than the independent observations; values in the ozone maximum in the tropical stratosphere are lower for the reanalysis. ERA-40 mixing ratios of water vapor are considerably larger than those for MOZAIC, typically by 20% in the tropical upper troposphere, and they may exceed 60% in the lower stratosphere over high latitudes. The results imply that the Brewer-Dobson circulation in the ECMWF reanalysis system is too fast, as is also evidenced by deficiencies in the way ERA-40 reproduces the water vapor "tape recorder'' signal in the tropical stratosphere. Finally, the paper examines the biases and their temporal variation during the 1990s in the way ERA-40 compares to the independent observations. We also discuss how the evaluation results depend on the instrument used, as well as on the version of the data.

Characteristics and variability of storm tracks in the north Pacific, Bering Sea, and Alaska

The North Pacific and Bering Sea regions represent loci of cyclogenesis and storm track activity. In this paper climatological properties of extratropical storms in the North Pacific/Bering Sea are presented based upon aggregate statistics of individual storm tracks calculated by means of a feature-tracking algorithm run using NCEP–NCAR reanalysis data from 1948/49 to 2008, provided by the NOAA/Earth System Research Laboratory and the Cooperative Institute for Research in Environmental Sciences, Climate Diagnostics Center. Storm identification is based on the 850-hPa relative vorticity field (ζ) instead of the often-used mean sea level pressure; ζ is a prognostic field, a good indicator of synoptic-scale dynamics, and is directly related to the wind speed. Emphasis extends beyond winter to provide detailed consideration of all seasons. Results show that the interseasonal variability is not as large during the spring and autumn seasons. Most of the storm variables—genesis, intensity, track density—exhibited a maxima pattern that was oriented along a zonal axis. From season to season this axis underwent a north–south shift and, in some cases, a rotation to the northeast. This was determined to be a result of zonal heating variations and midtropospheric moisture patterns. Barotropic processes have an influence in shaping the downstream end of storm tracks and, together with the blocking influence of the coastal orography of northwest North America, result in high lysis concentrations, effectively making the Gulf of Alaska the “graveyard” of Pacific storms. Summer storms tended to be longest in duration. Temporal trends tended to be weak over the study area. SST did not emerge as a major cyclogenesis control in the Gulf of Alaska.

Can climate models capture the structure of extratropical cyclones?

Composites of wind speeds, equivalent potential temperature, mean sea level pressure, vertical velocity, and relative humidity have been produced for the 100 most intense extratropical cyclones in the Northern Hemisphere winter for the 40-yr ECMWF Re-Analysis (ERA-40) and the high resolution global environment model (HiGEM). Features of conceptual models of cyclone structure—the warm conveyor belt, cold conveyor belt, and dry intrusion—have been identified in the composites from ERA-40 and compared to HiGEM. Such features can be identified in the composite fields despite the smoothing that occurs in the compositing process. The surface features and the three-dimensional structure of the cyclones in HiGEM compare very well with those from ERA-40. The warm conveyor belt is identified in the temperature and wind fields as a mass of warm air undergoing moist isentropic uplift and is very similar in ERA-40 and HiGEM. The rate of ascent is lower in HiGEM, associated with a shallower slope of the moist isentropes in the warm sector. There are also differences in the relative humidity fields in the warm conveyor belt. In ERA-40, the high values of relative humidity are strongly associated with the moist isentropic uplift, whereas in HiGEM these are not so strongly associated. The cold conveyor belt is identified as rearward flowing air that undercuts the warm conveyor belt and produces a low-level jet, and is very similar in HiGEM and ERA-40. The dry intrusion is identified in the 500-hPa vertical velocity and relative humidity. The structure of the dry intrusion compares well between HiGEM and ERA-40 but the descent is weaker in HiGEM because of weaker along-isentrope flow behind the composite cyclone. HiGEM’s ability to represent the key features of extratropical cyclone structure can give confidence in future predictions from this model.

Cross-hemispheric transport of Central African biomass burning pollutants: Implications for downwind ozone production

Pollutant plumes with enhanced concentrations of trace gases and aerosols were observed over the southern coast of West Africa during August 2006 as part of the AMMA wet season field campaign. Plumes were observed both in the mid and upper troposphere. In this study we examined the origin of these pollutant plumes, and their potential to photochemically produce ozone (O3) downwind over the Atlantic Ocean. Their possible contribution to the Atlantic O3 maximum is also discussed. Runs using the BOLAM mesoscale model including biomass burning carbon monoxide (CO) tracers were used to confirm an origin from central African biomass burning fires. The plumes measured in the mid troposphere (MT) had significantly higher pollutant concentrations over West Africa compared to the upper tropospheric (UT) plume. The mesoscale model reproduces these differences and the two different pathways for the plumes at different altitudes: transport to the north-east of the fire region, moist convective uplift and transport to West Africa for the upper tropospheric plume versus north-west transport over the Gulf of Guinea for the mid-tropospheric plume. Lower concentrations in the upper troposphere are mainly due to enhanced mixing during upward transport. Model simulations suggest that MT and UT plumes are 16 and 14 days old respectively when measured over West Africa. The ratio of tracer concentrations at 600 hPa and 250 hPa was estimated for 14–15 August in the region of the observed plumes and compares well with the same ratio derived from observed carbon dioxide (CO2) enhancements in both plumes. It is estimated that, for the period 1–15 August, the ratio of Biomass Burning (BB) tracer concentration transported in the UT to the ones transported in the MT is 0.6 over West Africa and the equatorial South Atlantic. Runs using a photochemical trajectory model, CiTTyCAT, initialized with the observations, were used to estimate in-situ net photochemical O3 production rates in these plumes during transport downwind of West Africa. The mid-troposphere plume spreads over altitude between 1.5 and 6 km over the Atlantic Ocean. Even though the plume was old, it was still very photochemically active (mean net O3 production rates over 10 days of 2.6 ppbv/day and up to 7 ppbv/day during the first days) above 3 km especially during the first few days of transport westward. It is also shown that the impact of high aerosol loads in the MT plume on photolysis rates serves to delay the peak in modelled O3 concentrations. These results suggest that a significant fraction of enhanced O3 in mid-troposphere over the Atlantic comes from BB sources during the summer monsoon period. According to simulated occurrence of such transport, BB may be the main source for O3 enhancement in the equatorial south Atlantic MT, at least in August 2006. The upper tropospheric plume was also still photochemically active, although mean net O3 production rates were slower (1.3 ppbv/day). The results suggest that, whilst the transport of BB pollutants to the UT is variable (as shown by the mesoscale model simulations), pollution from biomass burning can make an important contribution to additional photochemical production of O3 in addition to other important sources such as nitrogen oxides (NOx) from lightning.

Categorisation of synoptic environments associated with mesoscale convective systems over the UK

Mesoscale convective systems (MCSs) are relatively rare events in the UK but, when they do occur, can be associated with weather that is considered extreme with respect to climatology (as indicated by the number of such events that have been analysed as case studies). These case studies usually associate UK MCSs with a synoptic environment known as the Spanish plume. Here a previously published 17 year climatology of UK MCS events is extended to the present day (from 1998 to 2008) and these events classified according to the synoptic environment in which they form. Three distinct synoptic environments have been identified, here termed the classical Spanish plume, modified Spanish plume, and European easterly plume. Detailed case studies of the two latter, newly defined, environments are presented. Composites produced for each environment further reveal the differences between them. The classical Spanish plume is associated with an eastward propagating baroclinic cyclone that evolves according to idealised life cycle 1. Conditional instability is released from a warm moist plume of air advected northeastwards from Iberia that is capped by warmer, but very dry air, from the Spanish plateau. The modified Spanish plume is associated with a slowly moving mature frontal system associated with a forward tilting trough (and possibly cut-off low) at 500 hPa that evolves according to idealised life cycle 2. As in the classical Spanish plume, conditional instability is released from a warm plume of air advected northwards from Iberia. The less frequent European easterly plume is associated with an omega block centred over Scandinavia at upper levels. Conditional instability is released from a warm plume of air advected westwards across northern continental Europe. Unlike the Spanish plume environments, the European easterly plume is not a warm sector phenomena associated with a baroclinic cyclone. However, in all environments the organisation of convection is associated with the interaction of an upper-level disturbance with a low-level region of warm advection.

The NAO troposphere–stratosphere connection

Using monthly mean data, daily data, and theoretical arguments, relationships between surface pressure variations associated with the North Atlantic Oscillation (NAO), tropopause height, and the strength of the stratospheric vortex are established. An increase in the NAO index leads to a stronger stratospheric vortex, about 4 days later, as a result of increased equatorward refraction of upward-propagating Rossby waves. At tropopause level the effects of the enhanced NAO index and stratospheric polar vortex are opposite, resulting in a lower tropopause over Iceland and a higher tropopause over the Arctic. The raising of the Arctic tropopause leads to a stretching and spinup of the tropospheric column and is therefore associated with a lowering of the surface pressure near the North Pole. For monthly mean data it is found that a standard deviation increase in the NAO index is associated with a 10% increase in the strength of the stratospheric vortex, as measured by potential vorticity at 500 K. A simple theoretical model predicts that this is associated with about 300-m elevation of the Arctic tropopause, as is observed, and a 5-hPa lowering of the surface pressure at the North Pole. The effects of the spinup of the tropospheric column may project on the NAO pattern so that the stratosphere acts as an integrator of the NAO index.