The influence of the QBO on the propagation of equatorial waves into the stratosphere

The variation of stratospheric equatorial wave characteristics with the phase of the quasi-biennial oscillation (QBO) is investigated using ECMWF Re-Analysis and NOAA outgoing longwave radiation (OLR) data. The impact of the QBO phases on the upward propagation of equatorial waves is found to be consistent and significant. In the easterly phase, there is larger Kelvin wave amplitude but smaller westward-moving mixed Rossby–gravity (WMRG) and n = 1 Rossby (R1) wave amplitude due to reduced propagation from the upper troposphere into the lower stratosphere, compared with the westerly phase. Differences in the wave amplitude exist in a deeper layer in summer than in winter, consistent with the seasonality of ambient zonal winds. There is a strong evidence of Kelvin wave amplitude peaking just below the descending westerly phase, suggesting that Kelvin waves act to bring the westerly phase downward. However, the corresponding evidence for WMRG and R1 waves is less clear. In the lower stratosphere there is zonal variation in equatorial waves. This reflects the zonal asymmetry of wave amplitudes in the upper troposphere, the source for the lower-stratospheric waves. In easterly winters the upper-tropospheric WMRG and R1 waves over the eastern Pacific region appear to be somewhat stronger compared to climatology, perhaps because of the accumulation of waves that are unable to propagate upward into the lower stratosphere. Vertical propagation features of these waves are generally consistent with theory and suggest a mixture of Doppler shifting by ambient flows and filtering. Some lower-stratosphere equatorial waves have a connection with preceding tropical convection, especially for Kelvin and R1 waves in winter.

Using seasonal hindcasts to understand the origin of the equatorial cold tongue bias in CGCMs and its impact on ENSO

The cold equatorial SST bias in the tropical Pacific that is persistent in many coupled OAGCMs severely impacts the fidelity of the simulated climate and variability in this key region, such as the ENSO phenomenon. The classical bias analysis in these models usually concentrates on multi-decadal to centennial time series needed to obtain statistically robust features. Yet, this strategy cannot fully explain how the models errors were generated in the first place. Here, we use seasonal re-forecasts (hindcasts) to track back the origin of this cold bias. As such hindcasts are initialized close to observations, the transient drift leading to the cold bias can be analyzed to distinguish pre-existing errors from errors responding to initial ones. A time sequence of processes involved in the advent of the final mean state errors can then be proposed. We apply this strategy to the ENSEMBLES-FP6 project multi-model hindcasts of the last decades. Four of the five AOGCMs develop a persistent equatorial cold tongue bias within a few months. The associated systematic errors are first assessed separately for the warm and cold ENSO phases. We find that the models are able to reproduce either El Niño or La Niña close to observations, but not both. ENSO composites then show that the spurious equatorial cooling is maximum for El Niño years for the February and August start dates. For these events and at this time of the year, zonal wind errors in the equatorial Pacific are present from the beginning of the simulation and are hypothesized to be at the origin of the equatorial cold bias, generating too strong upwelling conditions. The systematic underestimation of the mixed layer depth in several models can also amplify the growth of the SST bias. The seminal role of these zonal wind errors is further demonstrated by carrying out ocean-only experiments forced by the AOCGCMs daily 10-meter wind. In a case study, we show that for several models, this forcing is sufficient to reproduce the main SST error patterns seen after 1 month in the AOCGCM hindcasts.

Seasonal persistence of ozone and zonal wind anomalies in the equatorial stratosphere

Analysis of the variability of equatorial ozone profiles in the Satellite Aerosol and Gas Experiment‐corrected Solar Backscatter Ultraviolet data set demonstrates a strong seasonal persistence of interannual ozone anomalies, revealing a seasonal dependence to equatorial ozone variability. In the lower stratosphere (40–25 hPa) and in the upper stratosphere (6–4 hPa), ozone anomalies persist from approximately November until June of the following year, while ozone anomalies in the layer between 16 and 10 hPa persist from June to December. Analysis of zonal wind fields in the lower stratosphere and temperature fields in the upper stratosphere reveals a similar seasonal persistence of the zonal wind and temperature anomalies associated with the quasi‐biennial oscillation (QBO). Thus, the persistence of interannual ozone anomalies in the lower and upper equatorial stratosphere, which are mainly associated with the well‐known QBO ozone signal through the QBO-induced meridional circulation, is related to a newly identified seasonal persistence of the QBO itself. The upper stratospheric QBO ozone signal is argued to arise from a combination of QBO‐induced temperature and NOx perturbations, with the former dominating at 5 hPa and the latter at 10 hPa. Ozone anomalies in the transition zone between dynamical and photochemical control of ozone (16–10 hPa) are less influenced by the QBO signal and show a quite different seasonal persistence compared to the regions above and below.

Symmetric stability of compressible zonal flows on a generalized equatorial β plane

Sufficient conditions are derived for the linear stability with respect to zonally symmetric perturbations of a steady zonal solution to the nonhydrostatic compressible Euler equations on an equatorial � plane, including a leading order representation of the Coriolis force terms due to the poleward component of the planetary rotation vector. A version of the energy–Casimir method of stability proof is applied: an invariant functional of the Euler equations linearized about the equilibrium zonal flow is found, and positive definiteness of the functional is shown to imply linear stability of the equilibrium. It is shown that an equilibrium is stable if the potential vorticity has the same sign as latitude and the Rayleigh centrifugal stability condition that absolute angular momentum increase toward the equator on surfaces of constant pressure is satisfied. The result generalizes earlier results for hydrostatic and incompressible systems and for systems that do not account for the nontraditional Coriolis force terms. The stability of particular equilibrium zonal velocity, entropy, and density fields is assessed. A notable case in which the effect of the nontraditional Coriolis force is decisive is the instability of an angular momentum profile that decreases away from the equator but is flatter than quadratic in latitude, despite its satisfying both the centrifugal and convective stability conditions.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part II: combined effects of gravity waves and equatorial planetary waves

This study examines the effect of combining equatorial planetary wave drag and gravity wave drag in a one-dimensional zonal mean model of the quasi-biennial oscillation (QBO). Several different combinations of planetary wave and gravity wave drag schemes are considered in the investigations, with the aim being to assess which aspects of the different schemes affect the nature of the modeled QBO. Results show that it is possible to generate a realistic-looking QBO with various combinations of drag from the two types of waves, but there are some constraints on the wave input spectra and amplitudes. For example, if the phase speeds of the gravity waves in the input spectrum are large relative to those of the equatorial planetary waves, critical level absorption of the equatorial planetary waves may occur. The resulting mean-wind oscillation, in that case, is driven almost exclusively by the gravity wave drag, with only a small contribution from the planetary waves at low levels. With an appropriate choice of wave input parameters, it is possible to obtain a QBO with a realistic period and to which both types of waves contribute. This is the regime in which the terrestrial QBO appears to reside. There may also be constraints on the initial strength of the wind shear, and these are similar to the constraints that apply when gravity wave drag is used without any planetary wave drag. In recent years, it has been observed that, in order to simulate the QBO accurately, general circulation models require parameterized gravity wave drag, in addition to the drag from resolved planetary-scale waves, and that even if the planetary wave amplitudes are incorrect, the gravity wave drag can be adjusted to compensate. This study provides a basis for knowing that such a compensation is possible.

The middle-atmosphere Hadley circulation and equatorial inertial adjustment

In the tropical middle atmosphere the climatological radiative equilibrium temperature is inconsistent with gradient-wind balance and the available angular momentum, especially during solstice seasons. Adjustment toward a balanced state results in a type of Hadley circulation that lies outside the “downward control” view of zonally averaged dynamics. This middle-atmosphere Hadley circulation is reexamined here using a zonally symmetric balance model driven through an annual cycle. It is found that the inclusion of a realistic radiation scheme leads to a concentration of the circulation near the stratopause and to its closing off in the mesosphere, with no need for relaxational damping or a rigid lid. The evolving zonal flow is inertially unstable, leading to a rapid process of inertial adjustment, which becomes significant in the mesosphere. This short-circuits the slower process of angular momentum homogenization by the Hadley circulation itself, thereby weakening the latter. The effect of the meridional circulation associated with extratropical wave drag on the Hadley circulation is considered. It is shown that the two circulations are independent for linear (quasigeostrophic) zonal-mean dynamics, and interact primarily through the advection of temperature and angular momentum. There appears to be no significant coupling in the deep Tropics via temperature advection since the wave-driven circulation is unable to alter meridional temperature gradients in this region. However, the wave-driven circulation can affect the Hadley circulation by advecting angular momentum out of the Tropics. The validity of the zonally symmetric balance model with parameterized inertial adjustment is tested by comparison with a three-dimensional primitive equations model. Fields from a middle-atmosphere GCM are also examined for evidence of these processes. While many aspects of the GCM circulation are indicative of the middle-atmosphere Hadley circulation, particularly in the upper stratosphere, it appears that the circulation is obscured in the mesosphere and lower stratosphere by other processes.

Balance model for equatorial long waves

Geophysical fluid models often support both fast and slow motions. As the dynamics are often dominated by the slow motions, it is desirable to filter out the fast motions by constructing balance models. An example is the quasi geostrophic (QG) model, which is used widely in meteorology and oceanography for theoretical studies, in addition to practical applications such as model initialization and data assimilation. Although the QG model works quite well in the mid-latitudes, its usefulness diminishes as one approaches the equator. Thus far, attempts to derive similar balance models for the tropics have not been entirely successful as the models generally filter out Kelvin waves, which contribute significantly to tropical low-frequency variability. There is much theoretical interest in the dynamics of planetary-scale Kelvin waves, especially for atmospheric and oceanic data assimilation where observations are generally only of the mass field and thus do not constrain the wind field without some kind of diagnostic balance relation. As a result, estimates of Kelvin wave amplitudes can be poor. Our goal is to find a balance model that includes Kelvin waves for planetary-scale motions. Using asymptotic methods, we derive a balance model for the weakly nonlinear equatorial shallow-water equations. Specifically we adopt the ‘slaving’ method proposed by Warn et al. (Q. J. R. Meteorol. Soc., vol. 121, 1995, pp. 723–739), which avoids secular terms in the expansion and thus can in principle be carried out to any order. Different from previous approaches, our expansion is based on a long-wave scaling and the slow dynamics is described using the height field instead of potential vorticity. The leading-order model is equivalent to the truncated long-wave model considered previously (e.g. Heckley & Gill, Q. J. R. Meteorol. Soc., vol. 110, 1984, pp. 203–217), which retains Kelvin waves in addition to equatorial Rossby waves. Our method allows for the derivation of higher-order models which significantly improve the representation of Rossby waves in the isotropic limit. In addition, the ‘slaving’ method is applicable even when the weakly nonlinear assumption is relaxed, and the resulting nonlinear model encompasses the weakly nonlinear model. We also demonstrate that the method can be applied to more realistic stratified models, such as the Boussinesq model.

Observations of diurnal and spatial variability of radiative forcing by equatorial deep convective clouds

From geostationary satellite observations of equatorial Africa and the equatorial east Atlantic during May and June 2000 we explore the radiative forcing by deep convective cloud systems in these regions. Deep convective clouds (DCCs) are associated with a mean radiative forcing relative to non–deep convective areas of −39 W m−2 over the Atlantic Ocean and of +13 W m−2 over equatorial Africa (±10 W m−2 in both cases). We show that over land the timing of the daily cycle of convection relative to the daily cycle in solar illumination and surface temperature significantly affects the mean radiative forcing by DCCs. Displacement of the daily cycle of DCC coverage by 2 hours changes their overall radiative effect by ∼10 W m−2, with implications for the simulation of the radiative balance in this region. The timing of the minimum DCC cover over land, close to noon local time, means that the mean radiative forcing is nearly maximized.

Diabatic balance model for the equatorial atmosphere

Using an asymptotic expansion, a balance model is derived for the shallow-water equations (SWE) on the equatorial beta-plane that is valid for planetary-scale equatorial dynamics and includes Kelvin waves. In contrast to many theories of tropical dynamics, neither a strict balance between diabatic heating and vertical motion nor a small Froude number is required. Instead, the expansion is based on the smallness of the ratio of meridional to zonal length scales, which can also be interpreted as a separation in time scale. The leading-order model is characterized by a semigeostrophic balance between the zonal wind and meridional pressure gradient, while the meridional wind v vanishes; the model is thus asymptotically nondivergent, and the nonzero correction to v can be found at the next order. Importantly for applications, the diagnostic balance relations are linear for winds when inferring the wind field from mass observations and the winds can be diagnosed without direct observations of diabatic heating. The accuracy of the model is investigated through a set of numerical examples. These examples show that the diagnostic balance relations can remain valid even when the dynamics do not, and the balance dynamics can capture the slow behavior of a rapidly varying solution.

The dynamics and relationships of precipitation, temperature and convection boundaries in the dayside auroral ionosphere

A continuous band of high ion temperature, which persisted for about 8 h and zigzagged north-south across more than five degrees in latitude in the dayside (07:00– 15:00MLT) auroral ionosphere, was observed by the EISCAT VHF radar on 23 November 1999. Latitudinal gradients in the temperature of the F-region electron and ion gases (Te and Ti , respectively) have been compared with concurrent observations of particle precipitation and field-perpendicular convection by DMSP satellites, in order to reveal a physical explanation for the persistent band of high Ti , and to test the potential role of Ti and Te gradients as possible markers for the open-closed field line boundary. The north/south movement of the equatorward Ti boundary was found to be consistent with the contraction/expansion of the polar cap due to an unbalanced dayside and nightside reconnection. Sporadic intensifications in Ti , recurring on _10-min time scales, indicate that frictional heating was modulated by time-varying reconnection, and the band of high Ti was located on open flux. However, the equatorward Ti boundary was not found to be a close proxy of the open-closed boundary. The closest definable proxy of the open-closed boundary is the magnetosheath electron edge observed by DMSP. Although Te appears to be sensitive to magnetosheath electron fluxes, it is not found to be a suitable parameter for routine tracking of the open-closed boundary, as it involves case dependent analysis of the thermal balance. Finally, we have documented a region of newly-opened sunward convecting flux. This region is situated between the convection reversal boundary and the magnetosheath electron edge defining the openclosed boundary. This is consistent with a delay of several minutes between the arrival of the first (super-Alfv´enic) magnetosheath electrons and the response in the ionospheric convection, conveyed to the ionosphere by the interior Alfv´en wave. It represents a candidate footprint of the low-latitude boundary mixing layer on sunward convecting open flux

Coordinated Cluster, ground-based instrumentation and low-altitude satellite observations of transient poleward-moving events in the ionosphere and in the tail lobe

During the interval between 8:00-9:30 on 14 January 2001, the four Cluster spacecraft were moving from the central magnetospheric lobe, through the dusk sector mantle, on their way towards intersecting the magnetopause near 15:00 MLT and 15:00 UT. Throughout this interval, the EIS-CAT Svalbard Radar (ESR) at Longyearbyen observed a series of poleward-moving transient events of enhanced F-region plasma concentration ("polar cap patches"), with a repetition period of the order of 10 min. Allowing for the estimated solar wind propagation delay of 75 ( 5) min, the interplanetary magnetic field (IMF) had a southward component during most of the interval. The magnetic footprint of the Cluster spacecraft, mapped to the ionosphere using the Tsyganenko T96 model (with input conditions prevailing during this event), was to the east of the ESR beams. Around 09:05 UT, the DMSP-F12 satellite flew over the ESR and showed a sawtooth cusp ion dispersion signature that also extended into the electrons on the equatorward edge of the cusp, revealing a pulsed magnetopause reconnection. The consequent enhanced ionospheric flow events were imaged by the SuperDARN HF backscatter radars. The average convection patterns (derived using the AMIE technique on data from the magnetometers, the EISCAT and SuperDARN radars, and the DMSP satellites) show that the associated poleward-moving events also convected over the predicted footprint of the Cluster spacecraft. Cluster observed enhancements in the fluxes of both electrons and ions. These events were found to be essentially identical at all four spacecraft, indicating that they had a much larger spatial scale than the satellite separation of the order of 600 km. Some of the events show a correspondence between the lowest energy magnetosheath electrons detected by the PEACE instrument on Cluster (10-20 eV) and the topside ionospheric enhancements seen by the ESR (at 400-700 km). We suggest that a potential barrier at the magnetopause, which prevents the lowest energy electrons from entering the magnetosphere, is reduced when and where the boundary-normal magnetic field is enhanced and that the observed polar cap patches are produced by the consequent enhanced precipitation of the lowest energy electrons, making them and the low energy electron precipitation fossil remnants of the magnetopause reconnection rate pulses.

Long term changes in EUV and X-ray emissions from the solar corona and chromosphere as measured by the response of the Earth’s ionosphere during total solar eclipses from 1932 to 1999

Measurements of the ionospheric E region during total solar eclipses in the period 1932-1999 have been used to investigate the fraction of Extreme Ultra Violet and soft X-ray radiation, phi, that is emitted from the limb corona and chromosphere. The relative apparent sizes of the Moon and the Sun are different for each eclipse, and techniques are presented which correct the measurements and, therefore, allow direct comparisons between different eclipses. The results show that the fraction of ionising radiation emitted by the limb corona has a clear solar cycle variation and that the underlying trend shows this fraction has been increasing since 1932. Data from the SOHO spacecraft are used to study the effects of short-term variability and it is shown that the observed long-term rise in phi has a negligible probability of being a chance occurrence.

Effects of a mid-latitude solar eclipse on the thermosphere and ionosphere: a modelling study

A modelling study is presented which investigates in-situ generated changes of the thermosphere and ionosphere during a solar eclipse. Neutral temperatures are expected to drop by up to 40 degrees K at 240 km height in the totality footprint, with neutral winds of up to 26 m/s responding to the change of pressure. Both temperatures and winds are found to respond with a time lag of 30 min after the passing of the Moon's shadow. A gravity wave is generated in the neutral atmosphere and propagates into the opposite hemisphere at around 300 m/s. The combined effects of thermal cooling and downwelling lead to an overall increase in [O], while [N(2)] initially rises and then for several hours after the eclipse is below the "steady state" level. An enhancement of [NmF2] is found and explained by the atmosphere's contraction during, and the reduced [O]/[N(2)] ratio after the eclipse.

Ionospheric and geomagnetic responses to changes in IMF B_Z: a superposed epoch study

Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices K p, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, f o F2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, δf o F2. For the southward turnings with a change in B z of δB z > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, f o F2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about −60 nT. Four days later, the index has still to fully recover and is at about −25 nT. Both the K p and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like K p, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of δf o Fl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where δB z > 11.5 nT and P ≤ 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in δf o F2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on δB z . At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.

Time-dependent flows in the coupled solar wind-magnetosphere-ionosphere system

Concepts of time-dependent flow in the coupled solar wind-magnetosphere-ionosphere system are discussed and compared with the frequently-adopted steady-state paradigm. Flows are viewed as resulting from departures of the system from equilibrium excited by dayside and nightside reconnection processes, with the flows then taking the system back towards a new equilibrium configuration. The response of the system to reconnection impulses, continuous but unbalanced reconnection and balanced steady-state reconnection are discussed in these terms. It is emphasized that in the time-dependent case the ionospheric and interplanetary electric fields are generally inductively decoupled from each other; a simple mapping of the interplanetary electric field along equipotential field lines into the ionosphere occurs only in the electrostatic steady-state case.