A comparative study of the quasi-16-day wave (QSDW) in the middle atmosphere
using meteor radar observations and reanalysis data from three Brazilian
stations, Sao Joao do Cariri (7.4
The planetary waves (period
A prominent category of the planetary waves consists of normal modes of
periods 2, 5, 10, and 16 days which could be considered as free or resonant
atmospheric modes (Salby, 1981). The 16-day wave is conceived to be a
manifestation of the gravest symmetric westward propagating Rossby normal
mode with zonal wavenumber 1 (Salby, 1981). The wave was first reported
observationally in the MLT by Kingsley et al. (1978) with the help of meteor
wind measurements from Sheffield (53
Seasonal variability of the quasi-16-day wave (QSDW) activity was reported by a handful of past investigations (Mitchell et al., 1999; Luo et al., 2002; Day and Mitchell, 2010). Forbes et al. (1995) simulated the wave behavior from the ground to 140 km and demonstrated direct penetration of the waves from the lower atmosphere as well as a plausible ducting channel through which the wave propagates from the winter to summer hemisphere. Using 4-year meteor radar observations from a mid-latitude station in the Northern Hemisphere, Mitchell et al. (1999) found amplitude maximum around winter–spring and a secondary maximum with comparatively smaller value around summer–fall intervals. Later, using a network of five MF radars distributed from the high to equatorial latitudes in the Northern Hemisphere, Luo et al. (2002) concluded that the wave phenomenon was strong in winter. Using meteor radar measurements over two sites located near both poles, Day and Mitchell (2010) found stronger wave activity in winter as compared to the summer. In contrast to such an occurrence, Williams and Avery (1992) reported stronger wave activity in summer mesopause using observations from a mid-latitude station, although they found significant QSDW amplitude in winter in the stratosphere. With the help of satellite-based Microwave Limb Sounder (MLS) derived temperature observations, McDonald et al. (2011) showed considerable seasonal as well as inter-annual variability of the wave consisting of eastward and westward propagating zonal wavenumber 1 and 2 components. They reported that both eastward and westward components are large in the Northern Hemisphere, whereas eastward components predominate in the Southern Hemisphere.
The phase of the quasi-biennial oscillation (QBO) is believed to cause a wave filtering by supporting the wave propagation in the westerly phase and obstructing the same in easterly phase (Espy et al., 1997; Jacobi et al., 1998). However, a few other investigators concluded that the QBO influence on the QSDW was not a permanent feature rather it was intermittent (Mitchell et al., 1999; Day et al., 2011).
Although a number of investigations carried out observationally and
theoretically over various locations of the globe, our understanding of
these prominent atmospheric normal modes is not complete, especially its
features in the Southern Hemisphere where planetary waves are generally
weaker as compared to the Northern Hemisphere (Volland, 1988). A previous
study by Lima et al. (2006) from Cachoeira Paulista (22.7
The meteor radar systems used for the present study located at Cariri (CA),
Cachoeira Paulista (CP), and Santa Maria (SM) are basically SKiYMET radars.
These identical radars operate at 35.24 MHz with a pulse width of 13
ERA-Interim database is provided by the European Centre for Medium-range
Weather Forecasts (ECMWF). The daily zonal and meridional wind data
available at latitude-longitude grid of 0.25
Standard wavelet technique is utilized to identify the dominant periodicities in the MLT and stratospheric wind data. To calculate the daily variance of the QSDW the data are subjected to a band pass filter of cut off periods of 12 and 20 days. The variance is derived using a 41-day window which is shifted by 1 day progressively. Similarly, the daily amplitude of the wave is estimated by least square fit to the unfiltered data considering a 41-day window and a progressive shift of the same by 1 day. The fit is carried out by constraining the period within the limit of 12–20 days. The daily eastward and westward 16-day wave amplitudes for various wave numbers are estimated in the ERA-Interim horizontal winds by applying two-dimensional fast Fourier transform (FFT) on a window of 32 days with a progressive shift of 1 day.
To find out possible coupling and/or nonlinear interactions of the QSDW with
shorter period planetary waves, bispectral analysis has been carried out. It
can be mentioned that due to nonlinear interaction between two primary
waves, two secondary waves with frequencies equal to sum and difference of
the primary frequencies will be generated. The bispectral analysis is a
useful tool to determine quadratic phase coupling in a triplet formed by any
three of the four wave components (two primary and two secondary waves). A
non-zero value in the bispectrum is expected in case of resonance condition,
i.e., when
The wavelet power spectra at 90 km in the zonal and meridional wind are shown in Fig. 1 to identify the planetary wave components existing during the observational span. The seasons are represented by the abbreviations as LS (late summer), F (fall), W (winter), S (spring), ES (early summer) and divided by bold vertical lines in the present and all the following plots. The solid curves in the plots represent 95 % confidence level. Wave periodicities are found to be more prominent in the zonal wind as compared to the meridional wind for all three stations. Higher periods are more dominant as compared to their low period counterparts. If we consider the period range associated with the QSDW, i.e., 12–20 days, we find it to be strongest in summer in the zonal wind and in winter in the case of meridional wind at Cariri. At Cachoeira Paulista the wave is observed to be strong in winter in the zonal component, whereas it is not evident in the meridional wind. At Santa Maria the wave is intermittently visible over various spans of the year in the zonal wind and it is prominent in the meridional wind in late spring and early summer. Therefore evident from the plots is significant latitudinal variability of the wave prevalence as well as difference in wave activity between the zonal and meridional components implying wave polarization.
Wavelet power spectra estimated at 90 km using meteor radar zonal and meridional winds at CA, CP, and SM. The bold curves represent 95 % confidence level. The seasons, i.e., late summer (LS), fall (F), winter (W), spring (S), early summer (ES) are separated by vertical lines.
Figure 2 shows the variance in the zonal and meridional winds associated with the QSDW in the MLT. The value of the variance is significantly higher in the zonal component as compared to the meridional wind. Another notable feature is that the wave variance is comparable for CA and SM, whereas it is considerably lower in the case of CP. In general, the variance is largest around winter for all the three sites although zonal wind variance at CA shows maximum in late summer. The maximum variance is found to be highly localized. Secondary maxima are also evident during spring and summer at certain altitude bins.
QSDW variance in the MLT estimated from the zonal and meridional winds for the three locations.
Altitudinal mean of daily amplitude of the QSDW is shown in Fig. 3. In general, the amplitude exhibits a number of broad maxima during various times of the year with considerably higher values in the case of zonal wind at CA and SM. On the other hand an additional maximum is visible in fall at CP and in spring at SM in the zonal wind. The highest amplitude among all the three sites is found at CA during late summer. The summer maxima at CA and SM are larger than the winter ones in the zonal wind and winter maxima are found to be larger compared to the summer peaks in the meridional wind. However, it should be mentioned that the mean wave amplitudes are significantly higher in the zonal wind compared to the meridional wind.
Vertical averaged amplitude of the QSDW in the MLT in the
The stratospheric planetary wave activities over the observational sites
illustrated by wavelet profiles of ECMWF zonal and meridional wind at 10 hPa
level (
Wavelet power spectra estimated at 10 hPa pressure level using ECMWF zonal and meridional winds at CA, CP, and SM. The bold curves represent 95 % confidence level.
The variation of the stratospheric winds with respect to longitude at 10 hPa
level filtered by a band pass filter of cut-off periods of 12 and 20 days
corresponding to the latitudes of the present observational stations are
shown in Fig. 5. The horizontal dashed lines denote the longitudes of the
sites 7.5
Longitudinal and temporal variation of the filtered ECMWF zonal and
meridional winds in the period band 12–20 days at 10 hPa pressure level for
the locations 7.5
The amplitudes of the eastward and westward propagating 16-day wave estimated from the ECMWF zonal wind for the zonal wave numbers 1, 2, 3 at 10 hPa level for the three stations are shown in Fig. 6. The amplitude is high in winter–early summer interval and the maximum amplitude peak is observed in winter. In general, wave number 1 component is much higher in amplitude as compared to other components at CA most of the time. At CP wave number 2 amplitude is significantly large and comparable to the wave number 1 most of the times unlike CA. In addition to the wave number 1 and 2 components the eastward propagating wave number 3 component is found to be significant in winter although it is much smaller compared to the other components. Among all the stations the amplitude is found to be maximum at SM for the eastward wave number 1 and 2 components. Therefore from the results it can be noted that the zonal wave number 1 is the strongest component at CA, whereas both wave number 1 and 2 components are dominant at CP and SM. At CA westward components are predominate, at CP both eastward and westward components are significant, and at SM eastward components are found to be dominant which are consistent with the filter profiles.
Amplitude of the eastward and westward propagating 16-day wave at 10 hPa pressure level at three stations corresponding to the zonal wave numbers 1, 2, and 3 derived from the ECMWF zonal wind.
Figure 7 shows the amplitude of the 16-day wave for various zonal wave numbers derived from the ECMWF meridional wind. Evident from the plot is much lower wave amplitude in the meridional wind compared to the zonal wind derived ones which is consistent with the previous analyses results depicted before. At CA the wave amplitudes are comparable among various wave numbers for both eastward and westward components. The wave number 3 component in the meridional wind is found to be comparable with the other two components at CP and SM unlike zonal wind. A common feature between the zonal and meridional components is the largest eastward amplitude at SM.
Same as Fig. 6, but here it is for the meridional wind.
Figure 8 plots the bispectra at 90 km for the three observational sites using
meteor winds. The abscissa is ranged between 12 and 20 days for considering
QSDW and the ordinate is shown between 2 and 12 days to illustrate lower
period planetary waves. It should be mentioned that the absolute values of
the bispectra do not bear any significance, rather the relative value is the one
of importance. The zonal spectrum at CA shows possible interactions of the
QSDW with 4.5-, 6- and 9-day waves. Apart from the longer periods the
meridional spectrum at CA shows coupling with lower periodicities
(< 3.5 days). At CP the wave with period > 8 days is found to be
coupled with the QSDW. In the meridional spectrum wave periods
Bispectra in terms of normalized power calculated from the meteor zonal and meridional winds at 90 km for three stations. Please note that the abscissa is ranged within 12–20 days period to express QSDW and the ordinate is scaled within the range 2–12 days to show shorter periods.
Same as Fig. 8, but it is at 10 hPa pressure level.
Background zonal wind is an important factor to control the propagation of the atmospheric waves. The daily variation of the background zonal wind in the MLT and troposphere–stratosphere (TS) is shown in Fig. 10. Evident from the plot is the increase of magnitude of the wind with the latitude. Also it can be noted that the westerly wind strengthens with the latitude. Common feature from the plot is the strongest westerly wind at low-MLT during winter at all the sites. A secondary peak of westerly wind is found at mid- and high-MLT in early summer at CP and SM. Significantly strong westerly wind during winter at low-MLT may be considered to be responsible for facilitating wave activity in the MLT since westerly wind supports the westward propagating Rossby modes by wave filtering. The daily variability is notably smaller in the TS as compared to the MLT. Strong westerly and easterly winds are found at upper stratosphere (> 40 km) during May–June and November–February intervals, respectively which shows latitudinal increase in magnitude. The altitude region between 2 and 20 km exhibits moderate westerly wind throughout the year. The moderate easterly wind is very consistent throughout the year around 20–30 km at CA.
Temporal variability of the zonal wind in the MLT and TS at three sites.
The present study from three southern hemispheric Brazilian stations depicts the latitudinal variability signature of the QSDW. The features of the wave in the Southern Hemisphere are studied inadequately so far and hence bear significance in the light of the existing literatures and present understanding in this regard.
Charney and Drazin (1961) derived a set of equations which describe the
propagation of the stationary planetary waves. They concluded that the
infinite region of negative refractive index squared caused total reflection
of the waves, while a finite region of negative refractive index squared
between two regions of positive refractive index squared favored partial
reflection of the waves. These inferences could also be applied for
traveling planetary waves by assuming nonzero phase velocity of the waves.
The refractive index is a quantity governed by the zonal wind. According to
the Charney and Drazin (1961) theorem a planetary wave can propagate
vertically upward if it satisfies the condition,
It can be noted that in addition to the westward planetary wave modes, considerable eastward propagating waves are observed at all the stations in the stratosphere which is found to be significant at SM (highest latitude among all stations) especially in winter–spring interval. Therefore strong eastward components of the QSDW at a higher latitude in winter–spring as noted in the present study are not consistent with the theory discussed above. There could be two plausible mechanisms for enhancement of the eastward waves, i.e., (i) propagation through ducting channel from the other hemisphere and (ii) excitation due to barotropic and/or baroclinic instability driven by a westerly jet. If we consider ducting as being the responsible mechanism then similar eastward wave enhancement at lower latitude is also expected which is not found in the present scenario. Therefore we surmise that the ducting mechanism is unlikely to cause eastward propagating wave amplification at higher latitude in the winter–spring interval in the stratosphere, in this case the barotropic and/or baroclinic instability may excite the waves in the presence of the westerly wind. However, ephemeral westward peaks especially in winter–spring period over all the stations could be conceived as a result of ducting favored propagation from the Northern Hemisphere. Nevertheless, with limited observations it is not possible to affirm exact mechanism/s responsible for wave excitation, so we attempted to indicate some probable factors and hence further investigations in this direction are still required to gain more insights into the actual underlying processes.
It is already mentioned that the eastward waves are stronger compared to
the westward counterparts during the winter–spring interval in the stratosphere
at CP and SM. Actually, in the Northern Hemisphere quasi-stationary
planetary waves and westward propagating normal mode Rossby waves are
dominant planetary waves. On the other hand, in the Southern Hemisphere,
apart from these modes, eastward traveling planetary waves are also found to
be significantly strong which are generated due to instability driven by the
strong westerly jet (Shiotani et al., 1990). In the Southern Hemisphere
constructive interference between the quasi-stationary wave of zonal wave
number 1 and eastward propagating wave of wave number 2 could cause the
amplification of the quasi-stationary wave 1 (Shiotani et al., 1990). Such
interaction between the eastward wave number 2 and quasi-stationary wave
number 1 was also verified by Ushimaru and Tanaka (1992). They found
anti-correlation behavior between the aforesaid two waves, indicating a
clear relationship of energy and entropy exchange in the interaction. Since
the westerly jet becomes stronger with latitude in the present study, the
instability also enhances and as a consequence stronger eastward wave
amplitude is observed at higher latitudes. It can be mentioned that wintertime strong quasi-stationary planetary wave amplitudes are very likely to
overestimate the traveling wave amplitudes as derived in the present study.
Therefore observed wintertime amplitude in the stratosphere may not be
purely due to traveling planetary waves and it may be a manifestation of
combined contributions from quasi-stationary and traveling planetary waves.
A previous study from a southern hemispheric mid-latitude station, Grahamstown
(33.3
Using EOS MLS observations, McDonald et al. (2011) reported stronger
eastward QSDW of zonal wave numbers, 1 and 2 as compared to the westward
components in the Southern Hemisphere which is very similar to the present
results. Of late, Lu et al. (2013) studied the characteristics of the
eastward propagating waves in the southern hemispheric winter using
satellite-based data in the middle atmosphere and reported strong wave
activity at higher latitude and sufficiently smaller at lower latitudes due
to negative refractive index squared inhibiting the equatorward propagation
of the excited waves or in other words there was a confinement of the
planetary waves at higher latitude (> 50
Although the QSDW is very strong in summer in the MLT it shows largest amplitude in winter in the case of stratosphere. Taking a glance at Fig. 6 one can note that summer time wave amplitudes are larger for westward waves although they are somewhat smaller compared to the eastward counterparts in winter in the stratosphere. The westerly wind in the upper troposphere is supportive for propagation of the westward waves. But easterly wind in the stratosphere during summer may prevent the further propagation of the waves to the mesosphere. Nevertheless the summer time MLT is found to show considerable wave activity. The summer MLT wave activity can be attributed to propagation of the planetary waves from the winter hemisphere through cross-equatorial ducting. In this context it can be noted that the wind in the MLT is mainly westerly and increases with latitude favoring the filtering of the westward waves from the winter hemisphere through ducting channel. Another possibility of the summer QSDW enhancement in the MLT could be through gravity waves. Gravity waves modulated by the QSDW in the lower atmosphere can propagate to mesosphere and imprint the feature of the QSDW as they dissipate there by depositing energy and momentum in the MLT as explained in the previous literature by Williams and Avery (1992). Unfortunately, contribution of the gravity waves is not looked into in the present case due to a lack of high temporal resolution data. Using numerical simulation, Forbes et al. (1995) inferred that inter-hemispheric propagation through ducting channel is diminished by gravity wave stress to some extent. They also indicated the possibility of in situ excitation of the QSDW in the lower thermosphere. Summer time mesospheric planetary wave generation may be caused by barotropic and/or baroclinic instability due to strong jet stream as discussed before in the case of winter (Plumb, 1983).
The QSDW is considered to have a mighty influence on the background mean flow and hence control the propagation of the other shorter period planetary waves, gravity waves, and tides. Using numerical study, Hagan et al. (1993) demonstrated that existence of a strong lower thermospheric jet in summer could act as a reflecting layer and consequent wave enhancement. Therefore enhanced long period planetary waves, e.g., QSDW etc. have a potential to modulate the propagating short period waves at the timescale of the periods of the longer ones. A recent study by Huang et al. (2013) reported strong interaction between the QSDW and tides (especially the diurnal tide). They also found clear signature of interaction between the newly generated secondary wave and QSDW. With long-term meteor wind observations over CP, Guharay et al. (2015a, b) found predominate coupling of the QSDW with diurnal tide and quasi-2-day wave (QTDW) and subsequent variability in the wave amplitudes. In the present study the results of bispectra clearly reveal appreciable coupling between the QSDW and shorter period planetary waves.
The present work has depicted comparative study of the QSDW utilizing meteor wind observations from three Brazilian stations as well as ECMWF winds. The wave exhibited evident latitudinal variability with increasing amplitude with latitude. Both eastward and westward components of the QSDW are found to be prominent and strong activity in zonal wave number 1 is observed. The wave amplitude exhibits multiple peaks during various times of the year with largest one in summer in the MLT and the largest peak in the stratosphere is found in winter. The dominant wave components vary from the westward to eastward from the tropical to mid-latitude in the stratosphere. The bispectra reveals coupling of the QSDW with other short period planetary waves most probably through nonlinear interactions especially at CA and CP.
The present work is supported by the Department of Space, Government of
India and Conselho Nacional de Desenvolvimento Cientifico e Tecnologio,
CNPq, Government of Brasil. The authors thank the European Center for the
Medium-Range Weather Forecasts (ECMWF) for providing necessary data which
are available at: