Temperature comparison between CHAMP radio occultation and TIMED/SABER measurements in the lower stratosphere

Global Positioning System (GPS) receiver on the CHAllenging Mini-satellite Payload (CHAMP) and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, provide middle atmosphere temperature profiles by Radio Occultation (RO) and limb viewing infrared emission measurements, respectively. These temperature profiles retrieved by two different techniques in the stratosphere are compared with each other using more than 1300 correlative profiles in March, September and December 2005. The over-all mean differences averaged over 15 and 35 km are approximately −2 K and standard deviation is less than 3 K. Below 20 km of altitude, relatively small mean temperature differences ∼1 K are observed in wide latitudinal range except for June (during the SABER nighttime observation). In the middle to low latitudes, between 30°S and 30°N, the temperature difference increases with height from ∼0–1 K at 15 km, to ∼−4 K at 35 km of altitude. Large temperature differences about −4 to −6 K are observed between 60°S and 30°N and 31–35 km of altitude for all months and between 0° and 30°N below 16 km during June (nighttime).     

Vertical wavenumber spectra of atmospheric ozone measured from ozonesonde observations

Vertical profiles of ozone have been measured at balloon altitudes. Our purpose is to examine the character of vertical wavenumber spectra of ozone fluctuations, to assess the possible roles of gravity wave field in ozone fluctuations, and to determine dominant vertical wavelengths of ozone spectra. Vertical wavenumber spectra of 12 ozone fluctuations obtained during June–August 2003 are presented. Results indicate that mean spectral slopes in the wavenumber range from 4.69 × 10⁻⁴ to 2.50 × 10⁻³ cyc/m are about −2.91 in the troposphere and −2.87 in the lower stratosphere, which is close to the slope of −3 predicted by current gravity wave saturation models. The consistency of the observed spectral slopes with the value of −3 predicted by current gravity wave saturation models suggests that the observed ozone fluctuations are due primarily to atmospheric gravity waves. At m = 1/(1000 m) the mean spectral amplitude is over 30 times larger in the lower stratosphere than in the troposphere. Mean vertical wavenumber spectra in area-preserving form reveal dominant vertical wavelengths of ∼2.6 km in the troposphere and ∼2.7 km in the lower stratosphere, which is consistent with the values varying between 1.5 and 3.0 km estimated from the velocity field and temperature field at these heights.     

Lower-mesospheric inversion layers over brazilian equatorial region using TIMED/SABER temperature profiles

Lower-mesospheric inversion layers (MILs) were studied using the temperature profiles observed by TIMED/SABER over Cariri (7.5°S, 36.5°W), Brazil, in 2005. A total 175 MILs were identified with the maximum occurrence in April and October and the minimum in January and July. The lower MIL is located in a height region from 70 to 90 km, with the peak at around 83 ± 4 km with the temperature of 205 ± 5 K, and the thickness of 4–10 km. The results show large amplitudes of MILs during equinoxes and minimum in solstices, with a clear semiannual variation. A general feature of lower MIL in monthly mean profile was observed twice a year, one from February to May, and the other from August to October with a downward shift of the top level. These results suggest that formation and long persistence of MIL is an important factor to investigate propagation of atmospheric gravity waves in the mesosphere-lower thermosphere (MLT) region.     

Generation of vertical–horizontal and horizontal–horizontal gravity gradients using stochastically modified integral estimators

The Earth’s gravity field modelling is an ill-posed problem having a sensitive solution to the error of data. Satellite gravity gradiometry (SGG) is a space technique to measure the second-order derivatives of geopotential for modelling this field, but the measurements should be validated prior to use. The existing terrestrial gravity anomalies and Earth gravity models can be used for this purpose. In this paper, the second-order vertical–horizontal (VH) and horizontal–horizontal (HH) derivatives of the extended Stokes formula in the local north-oriented frame are modified using biased, unbiased and optimum types of least-squares modification. These modified integral estimators are used to generate the VH and HH gradients at 250 km level for validation purpose of the SGG data. It is shown that, unlike the integral estimator for generating the second-order radial derivative of geopotential, the system of equations from which the modification parameters are obtained is unstable for all types of modification, with large cap size and high degree, and regularization is strongly required for solving the system. Numerical studies in Fennoscandia show that the SGG data can be estimated with an accuracy of 1 mE using an integral estimator modified by a biased type least-squares modification. In this case an integration cap size of 2.5° and a degree of modification of 100 for integrating 30′ × 30′ gravity anomalies are required.     

Gravity waves and wind-shear in the MLT at 23°S

We have used the technique suggested by Hocking [Hocking, W. A new approach to momentum flux determinations using SKiYMET meteor radars. Ann. Geophys. 23, 2005.] to derive short period wind variances in the 80–100 km region from meteor radar data. We find that these fluctuating winds, assumed to correspond to gravity waves and turbulence, are closely correlated with the vertical shear of the horizontal tidal winds. This close correlation suggests that in situ wind shear may be a major source of gravity waves and turbulence in the MLT. If this is the case, gravity waves generated in the troposphere and propagating up to the MLT region, generally assumed to constitute an important influence on the climatology of the region, may be a less important source of energy and momentum in the 80–100 km region than has been hitherto believed.     

Kinetic temperature and carbon dioxide from broadband infrared limb emission measurements taken from the TIMED/SABER instrument

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) experiment is one of four instruments on NASA’s Thermosphere–Ionosphere–Energetics and Dynamics (TIMED) satellite. SABER measures broadband infrared limb emission and derives vertical profiles of kinetic temperature (Tk) from the lower stratosphere to approximately 120 km, and vertical profiles of carbon dioxide (CO₂) volume mixing ratio (vmr) from approximately 70 km to 120 km. In this paper we report on SABER Tk/CO₂ data in the mesosphere and lower thermosphere (MLT) region from the version 1.06 dataset. The continuous SABER measurements provide an excellent dataset to understand the evolution and mechanisms responsible for the global two-level structure of the mesopause altitude. SABER MLT Tk comparisons with ground-based sodium lidar and rocket falling sphere Tk measurements are generally in good agreement. However, SABER CO₂ data differs significantly from TIME-GCM model simulations. Indirect CO₂ validation through SABER-lidar MLT Tk comparisons and SABER-radiation transfer comparisons of nighttime 4.3 μm limb emission suggest the SABER-derived CO₂ data is a better representation of the true atmospheric MLT CO₂ abundance compared to model simulations of CO₂ vmr.     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.     

Tidal associated temperature disturbances observed at the middle atmosphere (30–65 km) by a Rayleigh lidar at 23°S

A lidar has been operated in São José dos Campos, Brazil (23.2°S, 45.8°W) since 1972, mainly dedicated to the study of mesospheric sodium at the 589 nm resonant line. The molecular Rayleigh scattering can also be used provided we limit the height to ∼75 km where the sodium scattering begins. Nevertheless, the weak signal obtained only permits the determination of density and temperature profiles by accumulating a large number of shots giving only nocturnal average profiles. Temporal variations in density and temperature on the scale of hours can however, be obtained by performing a superposed epoch analysis for a given time interval and covering a period of several days. In this way we obtained hourly mean profiles grouped by months, seasons and overall, with data acquired from 1993 to 2004. The difference between the hourly temperatures and the nocturnal means shows for some months, with enough data coverage, downward propagating structures that apparently have tidal origin. The seasonal averages show a recurrent feature with high temperatures before and low temperatures after midnight above 50 km. Some similarity is found with the GSWM model, but the observed temperature amplitudes are twice of that for the model.     

Evaluation of uncertainty in gravity wave potential energy calculations through GPS radio occultation measurements

The application of the Global Positioning System (GPS) radio occultation (RO) method to the atmosphere enables the determination of height profiles of temperature, among other variables. From these measurements, gravity wave activity is usually quantified by calculating the potential energy through the integration of the ratio of perturbation and background temperatures between two given altitudes in each profile. The uncertainty in the estimation of wave activity depends on the systematic biases and random errors of the measured temperature, but also on additional factors like the selected vertical integration layer and the separation method between background and perturbation temperatures. In this study, the contributions of different parameters and variables to the uncertainty in the calculation of gravity wave potential energy in the lower stratosphere are investigated and quantified. In particular, a Monte Carlo method is used to evaluate the uncertainty that results from different GPS RO temperature error distributions. In addition, our analysis shows that RO data above 30 km height becomes dubious for gravity waves potential energy calculations.     

A superposed epoch study of the effects of solar wind stream interface events on the upper mesospheric and lower thermospheric temperature

The response of mesosphere and lower thermosphere (MLT) temperature to energetic particle precipitation over the Earth’s polar regions is not uniform due to complex phenomena within the MLT environment. Nevertheless, the modification of MLT temperatures may require an event-based study to be better observed. This work examines the influence of precipitation, triggered by solar wind stream interfaces (SI) event from 2002 to 2007, on polar MLT temperature. We first test the relationship between the ionospheric absorption measured by the SANAE IV (South African National Antarctic Expedition IV) riometer and the layer of energetic particle precipitation from POES (Polar Orbiting Environmental Satellites). The combined particle measurements from POES 15, 16, 17 and 18 were obtained close in time to the pass of the SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) temperature retrieval. Here, a superposed epoch technique is described and implemented to obtain average temperature profiles during SI-triggered particle precipitation. The superposed epoch average shows no significant temperature decrease below 100 km prior to the onset of SI-triggered precipitation, whereas a clear superposed average temperature decrease is observed at 95 km after the SI impact. A case study of SI event also yields similar observations. Results indicate that cooling effects due to the production of mesospheric odd hydrogen might be major contributors to temperature decrease under compressed solar wind stream.     

Semi-stochastic modification of second-order radial derivative of Abel–Poisson’s formula for validating satellite gravity gradiometry data

The geoid can be used to validate the satellite gravity gradiometry data. Validation of such data is important prior to their downward continuation because of amplification of the data errors through this process. In this paper, the second-order radial derivative of Abel–Poisson’s formula is modified stochastically to reduce the effect of the far-zone geoid and generate the second-order radial derivative of geopotential at 250 km level. The numerical studies over Fennoscandia show that this method yields the gradients with an error of 10 mE and when the long wavelength of geoid is removed from the estimator and restored after the computations (remove–compute–restore) the error will be in 1 mE level. We name this method semi-stochastic modification. The best case scenario is found when the degree of modification of the integral formula is 200 and the long wavelength geoid to degree 100 is removed and restored. In this case the geoid should have a resolution of 15′ × 15′ and the integration should be performed over a cap size of 3°.     

On the measurements of the moon’s infrared temperature and its relation to the phase angle

Radiometric measurements of the thermal radiation originating from the moon’s surface were obtained using an infrared detector operating at wavelengths between 8 and 14 μm. The measurements cover a full moon cycle. The variation of the moon’s temperature with the lunar phase angle was established. The lunar temperatures were 391 ± 2.0 K for the full moon, 240 ± 3.5 K for the first quarter, and 236 ± 3 K for the last quarter. For the rest of the phase angles, the lunar temperature varied between 170 and 380 K. Our results are comparable with those obtained previously at these phase angles. For the new moon phase, the obtained temperature was between 120 and 133 K. With the exception of the new moon phase, our measurements at all the phase angles were consistent with those obtained using Earth-based data and those obtained by the Diviner experiment and the Clementine spacecraft. At the new phase, our measurements were comparable with those obtained from the ground but were significantly higher than those obtained by the Diviner and Clementine data. We attribute this inconsistency to either the calibration curve of our detector, which does not perform well at very low temperatures, or to infrared emission from the atmosphere. A simple linear model to predict the lunar temperature as a function of the phase angle was proposed. The experimental errors that affect the measured temperatures are discussed.     

An alternative computation of a gravity field model from GOCE

GOCE is the first satellite with a gravitational gradiometer (SGG). This allows to determine a gravity field model with high spatial resolution and high accuracy. Four of the six independent components of the gravitational gradient tensors (GGT) are measured with high accuracy in the so-called measurement band (MB) from 5 to 100 mHz by the GOCE gradiometer. Based on more than 1 year of GOCE measurements, two gravity field models have been derived. Here, we introduce a strategy for spherical harmonic analysis (SHA) from GOCE measurements, with a bandpass filter applied to the SGG data, combined with orbit analysis based on the integral equation approach, and additional constraints (or stabilization) in the polar areas where no observation is available due to the orbit geometry. In addition, we combined the GOCE SGG part with a set of GRACE normal equations. This improves the accuracy of the gravity field in the long-wavelength parts, due to the complementarity of GOCE and GRACE. Comparison with other models and with external data shows that our results are rather close to the GPS-levelling data in well-selected test regions, with an uncertainty of 4–7 cm, for truncation at degree 200.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

Response of the mesopause temperatures to solar activity over Yakutia in 1999–2013

OH(6-2) rotational temperature trends and solar cycle effects are studied. Observations were carried out at the Maimaga station (63.04°N, 129.51°E) for the period August 1999 to March 2013. Measurements were conducted with an infrared spectrograph. Temperatures were determined from intensity ratios in the P branch of the OH band. The monthly average residuals of temperature after the subtraction of the mean seasonal variation were used for a search for the solar component of temperature response. The dependence of temperatures on solar activity has been investigated using the Ottawa 10.7 cm flux as a proxy. A linear regression fitting on residual temperatures yields a solar cycle coefficient of 4.24 ± 1.39 K/100 solar flux units (SFU). The cross-correlation analyses showed that changes of the residual temperature follow changes of solar activity with a quasi-two year delay (25 months). The temperature response at the delay of 25 months reaches 7 K/100 SFU. The possible reason of the observed delay can be an influence of quasi-biennial oscillations (QBO) of the atmosphere on the relation of temperature and solar activity. The value of the temperature trend after the subtraction of seasonal and solar components is not statistically significant.     

A new approach to solar wind monitoring

The monitoring of solar wind parameters is a key problem of the space weather program. We are presenting a new solution of plasma parameter determination suitable for small and fast solar wind monitors. The first version will be launched during the SPECTR-R project into a highly elongated orbit with apogee ∼350,000 km. The method is based on simultaneous measurements of the total ion flux and ion integral energy spectrum by six identical Faraday cups. Three of them are dedicated to determination of the ion flow direction, the other three (equipped with control grids supplied by a retarding potential) are used for determination of the density, temperature, and speed of the plasma flow. The version under development is primarily designed for the measurements in the solar wind and tail magnetosheath, thus for velocities range from 270 to 750 km/s, temperatures from 1 to 30 eV, and densities up to 200 cm⁻³. However, the instrument design can be simply modified for measurements in other regions with a substantial portion of low-energy plasma as a subsolar magnetosheath, cusp or low-latitude boundary layer. Testing of the engineering model shows that the proposed method can provide reliable plasma parameters with a high time resolution (up to 8 Hz). The paper presents not only the method and its technical realization but it documents all advantages and peculiarities of the suggested approach.     

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

Temperature measurements during the CAMP program

The Cold Arctic Mesopause Program (CAMP) was conducted at ESRANGE, Sweden, in July/August 1982. During the time period of several weeks, the temperature was monitored by ground-based OH emission spectrometers and by stellite radiance measurements. Rocket launchings occurred on the nights of $\text{3}\text{4}$ and $\text{11}\text{12}$ August. On $\text{3}\text{4}$ August, seven rocket payloads were launched during a period of noctilucent cloud sighting over ESRANGE. The presence of the NLC was confirmed by several rocket-borne photometer profiles. The temperature measurements showed that the temperature profiles in the stratosphere and lower mesosphere were near the expected values of high latitude summer models. A large amplitude wave structure with three temperature minima of 139K, 114K and 111K were observed at altitudes between 83 and 94 km. The temperature minimum at 83 km was the location of the observed NLC. The temperature minima caused by the growth of the gravity wave amplitude in the highly stable mesosphere provide the regions for the growth of particles by nucleation to optical scattering size, as well as regions where the nuclei for condensation can be formed through ion chemistry paths.