In this study, we discuss the thermospheric density response to the solar activity and QBO. By using the wavelet power spectrum method, the thermospheric density from 1967 to 2013 shows ~11-year period, semiannual and annual variations, while the seasonal variation is usually more significant under high solar activity conditions. Importantly, we investigate the possible link between the density and the QBO, with the aid of GAMDM model and different density residual method. As a result, the difference between the measured density and empirical model seemingly has QBO signal, but the ratio of them revels that the QBO signal does not exist in the thermospheric density. Comprehensively, we thought the stratospheric QBO cannot impact on the thermosphere, and more data and numerical modeling are needed for further validation.

We present in this work a method for estimation of plasma bubble mean zonal drift velocities using keograms generated from images of the OI 6300.0 nm nightglow emission collected from an equatorial station -- Cariri (7.4$^\circ$S, 36.5$^\circ$W), and a mid-latitude station -- Cachoeira Paulista (22.7$^\circ$S, 45$^\circ$W), both in the Brazilian sector. The mean zonal drift velocities were estimated for 239 events recorded from 2000 to 2003 in Cariri, and for 56 events recorded over Cachoeira Paulista from 1998 to 2000. It was found that plasma bubble zonal drift velocities are smaller ($\sim$60 ms$^{-1}$) for events occurring later in the night compared to those occurring earlier ($\sim$150 ms$^{-1}$). The decreasing rate of the zonal drift velocity is of $\sim$10 ms$^{-1}$/h. We have also found that, in general, bubble events appearing first in the west-most region of the keogram are faster than those appearing first in the east-most region of the keograms. Larger zonal drift velocities occur from 19 LT to 23 LT in a longitude range from 37$^\circ$ to 33$^\circ$. The method of velocity estimation using keograms compares favorably against the mosaic method developed by \cite{Arruda:2005}, but the standard deviation of the residuals for the zonal drift velocities from the two methods is $\sim$15 ms$^{-1}$