An all-sky survey at 230 GHz by MLS on Aura

The Microwave Limb Sounder (MLS) instrument is a small satellite-borne radio telescope. Its purpose is to make limb-scanning measurements of atmospheric composition. One of the gases to which it is sensitive is carbon monoxide (CO), detected via the J = 2 → 1 rotational transition at 230 GHz. CO is present in molecular gas clouds in the Milky Way. Although it was not designed for the purpose, MLS can detect emissions from galactic CO, allowing a map of the 230 GHz radio sky to be constructed. We report the MLS measurements of galactic radio emission and discuss their effect on the atmospheric mission of MLS. The region of the Milky Way with emissions strong enough to significantly affect MLS observations of atmospheric CO is identified. Ground-based radio astronomers have been mapping the sky using CO emission for many years. However, the MLS data are the first such survey to be carried out from space. The MLS survey covers a larger area of the sky than any other 230 GHz survey, but no previously unknown gas clouds are observed.     

A working hypothesis for connections between electrically-induced changes in cloud microphysics and storm vorticity,with possible effects on circulation

This paper outlines, and explores the uncertainties in, hypothesized connections between a series of processes that could explain two long-standing puzzles; those of (1) the observed winter storm vorticity responses to atmospheric energy inputs that change the ionosphere–earth current density, J_z, that appear to involve storm invigoration, and (2) changes in anti-cyclonic blocking and circulation that include the observed colder winters in Great Britain and western Europe at solar minima, and especially at extended solar minima. A working hypothesis for the mechanism responsible for (1) is that the flow of J_z through conductivity gradients, as in stratified cloud layers and fog, especially with sea-salt aerosol haze over the high latitude winter oceans, deposits electric changes on droplets and aerosol particles; most importantly on cloud condensation nuclei (CCN). These electric charges modulate scavenging of the particles in clouds and haze layers, increasing the concentration of small CCN and decreasing the concentration of large CCN. When further cloud formation occurs there is increased concentration of small droplets and decreased concentration of large ones, reducing coalescence and the production of rain. Thus updrafts carry more liquid water above the freezing level, and there the increased production of ice releases more latent heat and invigorates the updraft (the Rosenfeld mechanism), leading to increased vorticity. Here we explore the major uncertainties for the reality of the above chain of physical processes. A consequence of cumulative cyclonic vorticity increases is increases in downstream anti-cyclonic blocking. A further working hypothesis for (2) is that the invigoration may be large enough to contribute to the observed increases in blocking in winters at solar minima (high J_z) in the North Atlantic, that result in colder winters in the UK and northern Europe.     

SCIAMACHY limb measurements of NO2, BrO and OClO. Retrieval of vertical profiles: Algorithm,first results,sensitivity and comparison studies

The Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) measures scattered sun light also in limb viewing mode (i.e. tangential to Earth’s surface and its atmosphere), which allows determining vertical profiles of atmospheric trace gases. First results on the retrieval of NO₂, BrO and OClO profiles from the SCIAMACHY Limb measurements are presented and compared to independent satellite and balloon borne observations.     

Evaluation of the IRI model for the European region

The aim of this paper is to investigate various aspects of the International Reference Ionosphere (IRI) performance in European area and to evaluate its accuracy and efficiency for: long term prediction of the critical frequencies foF2 and the maximum usable frequencies (MUF); using storm-time correction option (ST); the total electron content (TEC) and the maximum observable frequency (MOF) updating. Data of foF2, TEC, MOF are related to 2005. It is obtained that median values of foF2 can be predicted with the mean error σ(med)∼ 0.49 MHz. For median values of MUF absolute σ was 1.39 MHz and relative σ_r was 8.8%. For instanteneous values estimates are increased to 1.58σ(med) MHz for foF2 and could reach 3.84 MHz for MUF. Using correction ST-option and TEC values provided ∼30% improvement but TEC seems to be more preferable. However, from considered parameters of the IRI updating (ST-factor, TEC, MOF) the best results were demonstrated by MOF. Using the IRI2007 to calculate TEC gives 20–50% improvement of TEC correspondence to experimental values but this improvement is not enough to treat TEC without the IRI model adaptation.     

SCIAMACHY solar irradiance observation in the spectral range from 240 to 2380 nm

The SCanning Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY) is part of the payload of ESA’s Environmental Satellite ENVISAT which was launched into a sun-synchronous polar orbit on 2002-03-01. It is the first spaceborne instrument covering a wavelength range from 240 to 2380 nm thus including ultraviolet, visible and near infrared spectral regions.The main purpose of SCIAMACHY is to determine the amount and distribution of a large number of atmospheric trace constituents by measuring the radiance backscattered from the Earth. In addition, several solar observations are performed with daily or orbital frequency.The presented results will cover the following topics: (a) comparison of the solar irradiance measured by SCIAMACHY with data from the instruments SOLSPEC/SOLSTICE/SUSIM and a solar spectrum derived by Kurucz; (b) comparison of the SCIAMACHY solar Mg II index with GOME and NOAA data; (c) correlation of the relative change of solar irradiance measured by SCIAMACHY with the sun spot index.The mean solar irradiance for each of the 8 SCIAMACHY channels agrees with the Kurucz data within ±2–3%. The presented analysis proves that SCIAMACHY is a valuable tool to monitor solar irradiance variations.     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

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°.     

Time dependent charging of layer clouds in the global electric circuit

There is much observational data consistent with the hypothesis that the ionosphere-earth current density (J_z) in the global electric circuit, which is modulated by both solar activity and thunderstorm activity, affects atmospheric dynamics and cloud cover. One candidate mechanism involves J_z causing the accumulation of space charge on droplets and aerosol particles, that affects the rate of scavenging of the latter, notably those of Cloud Condensation Nuclei (CCN) and Ice Forming Nuclei (IFN) (  and ). Space charge is the difference, per unit volume, between total positive and total negative electrical charge that is on droplets, aerosol particles (including the CCN and IFN) and air ions. The cumulative effects of the scavenging in stratiform clouds and aerosol layers in an air mass over the lifetime of the aerosol particles of 1–10 days affects the concentration and size distribution of the CCN, so that in subsequent episodes of cloud formation (including deep convective clouds) there can be effects on droplet size distribution, coagulation, precipitation processes, and even storm dynamics.     

Jovian periodicities (∼10 h, ∼40 min) on Ulysses’ Distant Jupiter Encounter observations around the Halloween CIR events

We analyzed data from four different instruments (HI-SCALE, URAP, SWOOPS, VHM/FGM) onboard Ulysses spacecraft (s/c) and we searched for possible evidence of Jovian emissions when the s/c approached Jupiter during the times of Halloween events (closest time approach/position to Jupiter: February 5, 2004/R = 1683 R_J,θ = ∼49°). In particular, we analyzed extensively the low energy ion measurements obtained by the HI-SCALE experiment in order to examine whether low energy ion/electron emissions show a symmetry, and whether they are observed at north high latitudes upstream from the jovian bow shock, as is known to occur in the region upstream from the south bow shock as well ( Marhavilas et al., 2001). We studied the period from October 2003 to March 2004, as Ulysses moved at distances 0.8–1.2 AU from the planet at north Jovicentric latitudes <75°, and we present here an example of characteristic Jovian periodicities in the measurements around a CIR observed by Ulysses on days ∼348–349/2003 (R = 1894 R_J,θ = 72°). We show that Ulysses observed low energy ion (∼0.055–4.7 MeV) and electron (>∼40 keV) flux and/or spectral modulation with the Jupiter rotation period (∼10 h) as well as variations with the same period in solar wind parameters, radio and magnetic field directional data. In addition, characteristic strong ∼40 min periodic variations were found superimposed on the ∼10 h ion spectral modulation. Both the ∼10 h and ∼40 min ion periodicities in HI-SCALE measurements were present in several cases during the whole period examined (October 2003 to March 2004) and were found to be more evident during some special conditions, for instance during enhanced fluxes around the start (forward shock) and the end (reverse shock) of CIRs. We infer that the Jovian magnetosphere was triggered by the impact of the CIRs, after the Halloween events, and it was (a) a principal source of forward and reverse shock-associated ion flux structures and (b) the cause of generation of ∼10 h quasi-periodic magnetic field and plasma modulation observed by Ulysses at those times.     

Observation of a 27-day solar signature in noctilucent cloud altitude

Previous studies have identified solar 27-day signatures in several parameters in the Mesosphere/Lower thermosphere region, including temperature and Noctilucent cloud (NLC) occurrence frequency. In this study we report on a solar 27-day signature in NLC altitude with peak-to-peak variations of about 400?m. We use SCIAMACHY limb-scatter observations from 2002 to 2012 to detect NLCs. The superposed epoch analysis method is applied to extract solar 27-day signatures. A 27-day signature in NLC altitude can be identified in both hemispheres in the SCIAMACHY dataset, but the signature is more pronounced in the northern hemisphere. The solar signature in NLC altitude is found to be in phase with solar activity and temperature for latitudes $?70°$N. We provide a qualitative explanation for the positive correlation between solar activity and NLC altitude based on published model simulations.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

An investigation of high frequency auroral plasma emissions with a tubular dipole antenna

A band of enhanced amplitudes which follows a local plasma frequency fn in raw high frequency (HF) noise spectra is usually related to plasma emissions in the upper hybrid band (fn, fu). The enhanced band in question occurs permanently in noise spectra recorded on the Intercosmos-19, APEX and CORONAS satellites in the altitude range of 500 km–3000 km. For moderately magnetized plasma with fn > 2fc (fc – electron gyro frequency), the band occurs below fn determined from the topside sounder and impedance data or from electron beam induced spectra. The simulations of an equivalent circuit composed of a dipole antenna in a cold plasma and its preamplifiers, determined the physical origin of the band as the passive circuit resonance, due to inductive character of the antenna in a frequency band (fc, fu). The resonance spectral content is highly structured due to an inflight variability of the circuit impedances. In this report we analyze the noise and impedance spectra which are the most typical in an auroral zone if fn > fc. We focus attention on determination of local electron plasma density, essential for provisional HF mode classification. We found that the natural plasma emission in the upper hybrid band does not manifest itself as the banded natural emission, which may be used for reliable determination of local plasma frequency in the altitude range of 500–3000 km. The fast magnetosonic mode predominates in the auroral emissions. The broadband and multi banded electromagnetic emissions extending from the fast magnetosonic band well above fn > fc are characteristic for the strong wave activity and are much less frequent.     

Observations of non-migrating tides and ionospheric perturbations of O(D) airglow by ISUAL instrument

In this paper, we present the spatial variations of O(¹D) airglow observed by the ISUAL (Imager of Sprites and Upper Atmospheric Lightning) instrument on board the FORMOSAT-2 satellite. With a CCD camera and a 630 nm filter, ISUAL can measure global atmospheric emissions lying between the heights of 80 and 300 km. In days of 3–6 September 2008 and 25–27 February 2009, ISUAL has measured the emissions of O(¹D) airglow with results showing strong longitudinal peak-3 and peak-4 structures. The Lomb-Scargle analyses for these two cases show periods of longitudes of 120° and 90° supporting the DE2 and DE3 non-migrating tides. The 630 nm emissions are enhanced in equatorial regions and are lying along the equator. Over Africa its intensity can sometimes increase up to 80% relative to other longitudes. The perturbation is so strong that non-migrating tides are erased. A case of bimodal distribution with strong emissions at latitudes in equator and mid-latitude in geographic coordinates was observed.     

Development of a regional rain retrieval algorithm for exclusive mesoscale convective systems over peninsular India

The present study emphasize the development of a region specific rain retrieval algorithm by taking into accounts the cloud features. Brightness temperatures (T_bs) from various TRMM Microwave Imager (TMI) channels are calibrated with near surface rain intensity as observed from the TRMM – Precipitation Radar. It shows that T_b–R relations during exclusive-Mesoscale Convective System (MCS) events have greater dynamical range compared to combined events of non-MCS and MCS. Increased dynamical range of T_b–R relations for exclusive-MCS events have led to the development of an Artificial Neural Network (ANN) based regional algorithm for rain intensity estimation. By using the exclusive MCSs algorithm, reasonably good improvement in the accuracy of rain intensity estimation is observed. A case study of a comparison of rain intensity estimation by the exclusive-MCS regional algorithm and the global TRMM 2A12 rain product with a Doppler Weather Radar shows significant improvement in rain intensity estimation by the developed regional algorithm.     

An Artificial Neural Network based approach for estimation of rain intensity from spectral moments of a Doppler Weather Radar

By using a Doppler Weather Radar (DWR) at Shriharikota (13.66°N & 80.23°E), an Artificial Neural Network (ANN) based technique is proposed to improve the accuracy of rain intensity estimation. Three spectral moments of a Doppler spectra are utilized as an input data to an ANN. Rain intensity, as measured by the tipping bucket rain gauges around the DWR station, are considered as a target values for the given inputs. Rain intensity as estimated by the developed ANN model is validated by the rain gauges measurements. With the help of a developed technique, reasonable improvement in the estimation of rain intensity is observed. By using the developed technique, root mean square error and bias are reduced in the range of 34–18% and 17–3% respectively, compared to Z–R approach.     

High resolution limb observations of clouds by the CRISTA-NF experiment during the SCOUT-O3 tropical aircraft campaign

The Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere – New Frontiers (CRISTA-NF) experiment on board the Russian research aircraft Geophysica measures limb emission spectra with an unprecedented vertical and horizontal resolution in the 4–15 μm wavelength region. The IR spectra measured during the SCOUT-O3 Tropical Aircraft Campaign have been analysed with respect of cloud occurrence, cloud vertical and horizontal extent, cloud spatial structures and their utilisation for trace gas retrievals. In addition indicators for ice water content and optical thickness of the clouds have been adopted. These new kinds of measurements in the upper troposphere/lower stratosphere region are especially valuable for the design and development of future space borne high resolution limb sounders.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.     

Comparison of OClO nadir measurements from SCIAMACHY and GOME

The scanning imaging absorption spectrometer for atmospheric chartography was launched successfully onboard ENVISAT on March 1, 2002. It observes the solar radiation transmitted and backscattered from the atmosphere and reflected from the ground in nadir, limb and occultation viewing modes. Chlorine dioxide (OClO), an important indicator for stratospheric chlorine activation, can be measured in the UV spectral range by differential optical absorption spectroscopy (DOAS).

First results of the DOAS retrieval of OClO slant column densities from the SCIAMACHY nadir measurements are presented and compared to measurements of the global ozone monitoring experiment (GOME), which has successfully measured OClO since 1995. While SCIAMACHY operates in the same orbit, it measures ≈30 min earlier than GOME and has an increased spatial resolution (30 × 60 km² compared to 40 × 320 km² for GOME).     

Response of the lower atmosphere to intense geomagnetic storms

In this short paper we examine the possible connection between atmospheric parameters measured at low and middle altitudes and geomagnetic storms occurred in 2000 and 2003. For that, from a chain of stations located near the meridian 60°W we compare the storm time values of temperature and wind speed with their standard deviation 2σ obtained from quiet time values. We observed statistically significant variations at several altitudes during the storm recovery phase and after it, both in neutral wind speed and temperature. The results obtained suggest that atmospheric parameters could be affected by geomagnetic storms.     

Assimilation of Mars Global Surveyor meteorological data

The Mars Global Surveyor (MGS) Mission is more than just a return to Mars. It represents a qualitatively new type of planetary mission. This is true not only because of the capable instrumentation aboard the spacecraft and the choice of orbit and data rate, but also because of a major beneficial change in our understanding of the interplay between observation and modeling. The Thermal Emission Spectrometer (TES), for example, is capable of taking thousands of infrared spectra per day. The 15 micrometer carbon dioxide band in these spectra can be inverted to obtain atmospheric temperature profiles, amounting to some 15,000–25,000 separate measurements each day. In addition, radio occultations produce fewer, but much higher resolution, temperature and pressure profiles. How is such a quantity of data to be handled? Clearly not in the old-fashioned manner in which a single profile is studied at great length and detail. The quantity of data does not allow this; and the quality of the data does not support it. (A single atmospheric profile is bound to contain small-scale spatial and temporal variation—e.g., gravity waves—that cannot be removed unambiguously, as well as possible errors due to noise in the observations or non-uniqueness of the inversions). Furthermore, much of the atmospheric science of interest (winds, for example) depends on quantities like horizontal temperature gradients that are not directly observed. Instead, the data should be examined collectively with the aid of a model that incorporates our knowledge of the governing physics. Described here are the first results of such a data assimilation exercise with TES observations during the aerobraking hiatus period at L_s ≈ 200.