An observing system simulation experiment for climate monitoring with GNSS radio occulation data: setup and testbed study

The long-term stability, high accuracy, all-weather capability, high vertical resolution, and global coverage of Global Navigation Satellite System (GNSS) radio occultation (RO) suggests it as a promising tool for global monitoring of atmospheric temperature change. With the aim to investigate and quantify how well a GNSS RO observing system is able to detect climate trends, we are currently performing an (climate) observing system simulation experiment over the 25-year period 2001 to 2025, which involves quasi-realistic modeling of the neutral atmosphere and the ionosphere. We carried out two climate simulations with the general circulation model MAECHAM5 (Middle Atmosphere European Centre/Hamburg Model Version 5) of the MPI-M Hamburg, covering the period 2001–2025: One control run with natural variability only and one run also including anthropogenic forcings due to greenhouse gases, sulfate aerosols, and tropospheric ozone. On the basis of this, we perform quasi-realistic simulations of RO observables for a small GNSS receiver constellation (six satellites), state-of-the-art data processing for atmospheric profiles retrieval, and a statistical analysis of temperature trends in both the “observed” climatology and the “true” climatology. Here we describe the setup of the experiment and results from a test bed study conducted to obtain a basic set of realistic estimates of observational errors (instrument- and retrieval processing-related errors) and sampling errors (due to spatial-temporal undersampling). The test bed results, obtained for a typical summer season and compared to the climatic 2001–2025 trends from the MAECHAM5 simulation including anthropogenic forcing, were found encouraging for performing the full 25-year experiment. They indicated that observational and sampling errors (both contributing about 0.2 K) are consistent with recent estimates of these errors from real RO data and that they should be sufficiently small for monitoring expected temperature trends in the global atmosphere over the next 10 to 20 years in most regions of the upper troposphere and lower stratosphere (UTLS). Inspection of the MAECHAM5 trends in different RO-accessible atmospheric parameters (microwave refractivity and pressure/geopotential height in addition to temperature) indicates complementary climate change sensitivity in different regions of the UTLS so that optimized climate monitoring shall combine information from all climatic key variables retrievable from GNSS RO data.

Sensitivity of Southern Hemisphere climate to zonal asymmetry in ozone

Climate model simulations of past and future climate invariably contain prescribed zonal mean stratospheric ozone. While the effects of zonal asymmetry in ozone have been examined in the Northern Hemisphere, much greater zonal asymmetry occurs in the Southern Hemisphere during the break up of the Antarctic ozone hole. We prescribe a realistic three-dimensional distribution of ozone in a high vertical resolution atmospheric model and compare results with a simulation containing zonal mean ozone. Prescribing the three dimensional ozone distribution results in a cooling of the stratosphere and upper troposphere comparable to that caused by ozone depletion itself. Our results suggest that changes in the zonal asymmetry of ozone have had important impacts on Southern Hemisphere climate, and will continue to do so in the future.

Linear amplification of marginally neutral baroclinic waves

Baroclinic wave development is investigated for unstable parallel shear flows in the limit of vanishing normal-mode growth rate. This development is described in terms of the propagation and interaction mechanisms of two coherent structures, called counter-propagating Rossby waves (CRWs). It is shown that, in this limit of vanishing normal-mode growth rate, arbitrary initial conditions produce sustained linear amplification of the marginally neutral normal mode (mNM). This linear excitation of the mNM is subsequently interpreted in terms of a resonance phenomenon. Moreover, while the mathematical character of the normal-mode problem changes abruptly as the bifurcation point in the dispersion diagram is encountered and crossed, it is shown that from an initial-value viewpoint, this transition is smooth. Consequently, the resonance interpretation remains relevant (albeit for a finite time) for wavenumbers slightly different from the ones defining cut-off points. The results are further applied to a two-layer version of the classic Eady model in which the upper rigid lid has been replaced by a simple stratosphere.

Analysis of convectively-generated gravity waves in mesoscale model simulations and wind profiler observations

The characteristics of convectively-generated gravity waves during an episode of deep convection near the coast of Wales are examined in both high resolution mesoscale simulations [with the (UK) Met Oce Unified Model] and in observations from a Mesosphere-Stratosphere-Troposphere (MST) wind profiling Doppler radar. Deep convection reached the tropopause and generated vertically propagating, high frequency waves in the lower stratosphere that produced vertical velocity perturbations O(1 m/s). Wavelet analysis is applied in order to determine the characteristic periods and wavelengths of the waves. In both the simulations and observations, the wavelet spectra contain several distinct preferred scales indicated by multiple spectral peaks. The peaks are most pronounced in the horizontal spectra at several wavelengths less than 50 km. Although these peaks are most clear and of largest amplitude in the highest resolution simulations (with 1 km horizontal grid length), they are also evident in coarser simulations (with 4 km horizontal grid length). Peaks also exist in the vertical and temporal spectra (between approximately 2.5 and 4.5 km, and 10 to 30 minutes, respectively) with good agreement between simulation and observation. Two-dimensional (wavenumber-frequency) spectra demonstrate that each of the selected horizontal scales contains peaks at each of preferred temporal scales revealed by the one- dimensional spectra alone.

Temperature trends derived from Stratospheric Sounding Unit radiances: the effect of increasing CO2 on the weighting function

Stratospheric Sounding Units (SSU) on the NOAA operational satellites have been the main source of near global temperature trend data above the lower stratosphere. They have been used extensively for comparison with model-derived trends. The SSU senses in the 15 micron band of CO2 and hence the weighting function is sensitive to changes in CO2 concentrations. The impact of this change in weighting function has been ignored in all recent trend analyses. We show that the apparent trends in global mean brightness temperature due to the change in weighting function vary from about -0.4 K/decade to 0.4 K/decade depending on the altitude sensed by the different SSU channels. For some channels, this apparent trend is of a similar size to the trend deduced from SSU data but ignoring the change in weighting function. In the mid-stratosphere, the revised trends are now significantly more negative and in better agreement with model-calculated trends.

Stratospheric communication of El Nino teleconnections to European winter

The stratospheric role in the European winter surface climate response to El Niño–Southern Oscillation sea surface temperature forcing is investigated using an intermediate general circulation model with a well-resolved stratosphere. Under El Niño conditions, both the modeled tropospheric and stratospheric mean-state circulation changes correspond well to the observed “canonical” responses of a late winter negative North Atlantic Oscillation and a strongly weakened polar vortex, respectively. The variability of the polar vortex is modulated by an increase in frequency of stratospheric sudden warming events throughout all winter months. The potential role of this stratospheric response in the tropical Pacific–European teleconnection is investigated by sensitivity experiments in which the mean state and variability of the stratosphere are degraded. As a result, the observed stratospheric response to El Niño is suppressed and the mean sea level pressure response fails to resemble the temporal and spatial evolution of the observations. The results suggest that the stratosphere plays an active role in the European response to El Niño. A saturation mechanism whereby for the strongest El Niño events tropospheric forcing dominates the European response is suggested. This is examined by means of a sensitivity test and it is shown that under large El Niño forcing the European response is insensitive to stratospheric representation.

Stratospheric Memory and Skill of Extended-Range Weather Forecasts

We use an empirical statistical model to demonstrate significant skill in making extended-range forecasts of the monthly-mean Arctic Oscillation (AO). Forecast skill derives from persistent circulation anomalies in the lowermost stratosphere and is greatest during boreal winter. A comparison to the Southern Hemisphere provides evidence that both the time scale and predictability of the AO depend on the presence of persistent circulation anomalies just above the tropopause. These circulation anomalies most likely affect the troposphere through changes to waves in the upper troposphere, which induce surface pressure changes that correspond to the AO.

Sensitivity of tropospheric forecasts to stratospheric initial conditions

A number of recent papers in the atmospheric science literature have suggested that a dynamical link exists between the stratosphere and troposphere. Numerical modelling studies have shown that the troposphere has a time-mean response to changes to the stratospheric climatological state. In this study the response of the troposphere to an imposed transient stratospheric change is examined. The study uses a high horizontal and vertical resolution numerical weather-prediction model. Experiments compare the tropospheric forecasts of two medium-range forecast ensembles which have identical tropospheric initial conditions and different stratospheric initial conditions. In three case studies described here, stratospheric initial conditions have a statistically significant impact on the tropospheric flow. The mechanism for this change involves, in its most basic step, a change to tropospheric synoptic-scale systems. A consistent change to the tropospheric synoptic-scale systems occurs in response to the stratospheric initial conditions. The aggregated impact of changes to individual synoptic systems maps strongly onto the structure of the Arctic Oscillation, particularly over the North Atlantic storm track. The relationship between the stratosphere and troposphere, while apparent in Arctic Oscillation diagnostics, does not occur on coherent, hemispheric scales.

A Possible Influence of Equatorial Winds on the September 2002 Southern Hemisphere Sudden Warming Event

The stratospheric sudden warming in the Southern Hemisphere (SH) in September 2002 was unexpected for two reasons. First, planetary wave activity in the Southern Hemisphere is very weak, and midwinter warmings have never been observed, at least not since observations of the upper stratosphere became regularly available. Second, the warming occurred in a west phase of the quasi-biennial oscillation (QBO) in the lower stratosphere. This is unexpected because warmings are usually considered to be more likely in the east phase of the QBO, when a zero wind line is present in the winter subtropics and hence confines planetary wave propagation to higher latitudes closer to the polar vortex. At first, this evidence suggests that the sudden warming must therefore be simply a result of anomalously strong planetary wave forcing from the troposphere. However, recent model studies have suggested that the midwinter polar vortex may also be sensitive to the equatorial winds in the upper stratosphere, the region dominated by the semiannual oscillation. In this paper, the time series of equatorial zonal winds from two different data sources, the 40-yr ECMWF Re-Analysis (ERA) and the Met Office assimilated dataset, are reviewed. Both suggest that the equatorial winds in the upper stratosphere above 10 hPa were anomalously easterly in 2002. Idealized model experiments are described in which the modeled equatorial winds were relaxed toward these observations for various years to examine whether the anomalous easterlies in 2002 could influence the timing of a warming event. It is found that the 2002 equatorial winds speed up the evolution of a warming event in the model. Therefore, this study suggests that the anomalous easterlies in the 1–10-hPa region may have been a contributory factor in the development of the observed SH warming. However, it is concluded that it is unlikely that the anomalous equatorial winds alone can explain the 2002 warming event.

The Splitting of the Stratospheric Polar Vortex in the Southern Hemisphere, September 2002: Dynamical Evolution

The polar vortex of the Southern Hemisphere (SH) split dramatically during September 2002. The large-scale dynamical effects were manifest throughout the stratosphere and upper troposphere, corresponding to two distinct cyclonic centers in the upper troposphere–stratosphere system. High-resolution (T511) ECMWF analyses, supplemented by analyses from the Met Office, are used to present a detailed dynamical analysis of the event. First, the anomalous evolution of the SH polar vortex is placed in the context of the evolution that is usually witnessed during spring. Then high-resolution fields of potential vorticity (PV) from ECMWF are used to reveal several dynamical features of the split. Vortex fragments are rapidly sheared out into sheets of high (modulus) PV, which subsequently roll up into distinct synoptic-scale vortices. It is proposed that the stratospheric circulation becomes hydrodynamically unstable through a significant depth of the troposphere–stratosphere system as the polar vortex elongates.

On the influence of stratospheric water vapor changes on the tropospheric circulation

Observations suggest that the mixing ratio of water vapour in the stratosphere has increased by 20–50% between the 1960s and mid-1990s. Here we show that inclusion of such a stratospheric water vapour (SWV) increase in a state-of-the-art climate model modifies the circulation of the extratropical troposphere: the modeled increase in the North Atlantic Oscillation (NAO) index is 40% of the observed increase in NAO index between 1965 and 1995, suggesting that if the SWV trend is real, it explains a significant fraction of the observed NAO trend. Our results imply that SWV changes provide a novel mechanism for communicating the effects of large tropical volcanic eruptions and ENSO events to the extratropical troposphere over timescales of a few years, which provides a mechanism for interannual climate predictability. Finally, we discuss our results in the context of regional climate change associated with changes in methane emissions.

Stratospheric circulation in seasonal forecasting models: implications for seasonal prediction

Accurate seasonal forecasts rely on the presence of low frequency, predictable signals in the climate system which have a sufficiently well understood and significant impact on the atmospheric circulation. In the Northern European region, signals associated with seasonal scale variability such as ENSO, North Atlantic SST anomalies and the North Atlantic Oscillation have not yet proven sufficient to enable satisfactorily skilful dynamical seasonal forecasts. The winter-time circulations of the stratosphere and troposphere are highly coupled. It is therefore possible that additional seasonal forecasting skill may be gained by including a realistic stratosphere in models. In this study we assess the ability of five seasonal forecasting models to simulate the Northern Hemisphere extra-tropical winter-time stratospheric circulation. Our results show that all of the models have a polar night jet which is too weak and displaced southward compared to re-analysis data. It is shown that the models underestimate the number, magnitude and duration of periods of anomalous stratospheric circulation. Despite the poor representation of the general circulation of the stratosphere, the results indicate that there may be a detectable tropospheric response following anomalous circulation events in the stratosphere. However, the models fail to exhibit any predictability in their forecasts. These results highlight some of the deficiencies of current seasonal forecasting models with a poorly resolved stratosphere. The combination of these results with other recent studies which show a tropospheric response to stratospheric variability, demonstrates a real prospect for improving the skill of seasonal forecasts.

On the sensitivity of annular mode dynamics to stratospheric radiative timescales

Long decorrelation timescales of the annular mode are observed in the lower stratosphere. This study uses a simple dynamical model, which has been used extensively to study stratosphere-troposphere coupling, to investigate the origin of the long dynamical timescales. Several long runs of the model are completed, with different imposed thermal damping timescales in the stratosphere. The dynamical timescales of the annular mode are found to be largely insensitive to the input thermal damping timescales, producing similar dynamical timescales in all cases below 50hPa. This result suggests that the hypothesis that long timescales in the lower stratosphere are due to long radiative timescales in this region is false.

Stratospheric temperature and radiative forcing response to 11-year solar cycle changes in irradiance and ozone

The 11-yr solar cycle temperature response to spectrally resolved solar irradiance changes and associated ozone changes is calculated using a fixed dynamical heating (FDH) model. Imposed ozone changes are from satellite observations, in contrast to some earlier studies. A maximum of 1.6 K is found in the equatorial upper stratosphere and a secondary maximum of 0.4 K in the equatorial lower stratosphere, forming a double peak in the vertical. The upper maximum is primarily due to the irradiance changes while the lower maximum is due to the imposed ozone changes. The results compare well with analyses using the 40-yr ECMWF Re-Analysis (ERA-40) and NCEP/NCAR datasets. The equatorial lower stratospheric structure is reproduced even though, by definition, the FDH calculations exclude dynamically driven temperature changes, suggesting an important role for an indirect dynamical effect through ozone redistribution. The results also suggest that differences between the Stratospheric Sounding Unit (SSU)/Microwave Sounding Unit (MSU) and ERA-40 estimates of the solar cycle signal can be explained by the poor vertical resolution of the SSU/MSU measurements. The adjusted radiative forcing of climate change is also investigated. The forcing due to irradiance changes was 0.14 W m−2, which is only 78% of the value obtained by employing the standard method of simple scaling of the total solar irradiance (TSI) change. The difference arises because much of the change in TSI is at wavelengths where ozone absorbs strongly. The forcing due to the ozone change was only 0.004 W m−2 owing to strong compensation between negative shortwave and positive longwave forcings.

Convectively coupled equatorial waves in high resolution Hadley centre climate models

Current global atmospheric models fail to simulate well organised tropical phenomena in which convection interacts with dynamics and physics. A new methodology to identify convectively coupled equatorial waves, developed by NCAS-Climate, has been applied to output from the two latest models of the Met Office/Hadley Centre which have fundamental differences in dynamical formulation. Variability, horizontal and vertical structures, and propagation characteristics of tropical convection and equatorial waves, along with their coupled behaviour in the models are examined and evaluated against a previous comprehensive study of observations. It is shown that, in general, the models perform well for equatorial waves coupled with off-equatorial convection. However they perform poorly for waves coupled with equatorial convection. The vertical structure of the simulated wave is not conducive to energy conversion/growth and does not support the correct physical-dynamical coupling that occurs in the real world. The following figure shows an example of the Kelvin wave coupled with equatorial convection. It shows that the models fail to simulate a key feature of convectively coupled Kelvin wave in observations, namely near surface anomalous equatorial zonal winds together with intensified equatorial convection and westerly winds in phase with the convection. The models are also not able to capture the observed vertical tilt structure and the vertical propagation of the Kelvin wave into the lower stratosphere as well as the secondary peak in the mid-troposphere, particularly in HadAM3. These results can be used to provide a test-bed for experimentation to improve the coupling of physics and dynamics in climate and weather models.