Inversion, error analysis, and validation of GPS/MET occultation data

The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.     

The spatial structure of the dayside ionospheric trough

Tomographic imaging provides a powerful technique for obtaining images of the spatial distribution of ionospheric electron density at polar latitudes. The method, which involves monitoring radio transmissions from the Navy Navigation Satellite System at a meridional chain of ground receivers, has particular potential for complementing temporal measurements by other observing techniques such as the EISCAT incoherent-scatter radar facility. Tomographic reconstructions are presented here from a two-week campaign in November 1995 that show large-scale structuring of the polar ionosphere. Measurements by the EISCAT radar confirm the authenticity of the technique and provide additional information of the plasma electron and ion temperatures. The dayside trough, persistently observed at high latitudes during a geomagnetically quiet period but migrating to lower latitudes with increasing activity, is discussed in relationship to the pattern of the polarcap convection.     

Geophysical phenomena during an ionospheric modification experiment at Tromsø, Norway

We present an analysis of phenomena observed by HF distance-diagnostic tools located in St. Petersburg combined with multi-instrument observation at Tromsø in the HF modified ionosphere during a magnetospheric substorm. The observed phenomena that occurred during the Tromsø heating experiment in the nightside auroral E_s region of the ionosphere depend on the phase of substorm. The heating excited small-scale field-aligned irregularities in the E region responsible for field-aligned scattering of diagnostic HF waves. The equipment used in the experiment was sensitive to electron density irregularities with wavelengths 12–15 m across the geomagnetic field lines. Analysis of the Doppler measurement data shows the appearance of quasiperiodic variations with a Doppler frequency shift, f_d and periods about 100–120 s during the heating cycle coinciding in time with the first substorm activation and initiation of the upward field-aligned currents. A relationship between wave variations in f_d and magnetic pulsations in the Y-component of the geomagnetic field at Tromsø was detected. The analysis of the magnetic field variations from the IMAGE magnetometer stations shows that ULF waves occurred, not only at Tromsø, but in the adjacent area bounded by geographical latitudes from 70.5° to 68° and longitudes from 16° to 27°. It is suggested that the ULF observed can result from superposition of the natural and heater-induced ULF waves. During the substorm expansion a strong stimulated electromagnetic emission (SEE) at the third harmonic of the downshifted maximum frequency was found. It is believed that SEE is accompanied by excitation of the VLF waves penetrating into magneto-sphere and stimulating the precipitation of the energetic electrons (10–40 keV) of about 1-min duration. This is due to a cyclotron resonant interaction of natural precipitating electrons (1–10 keV) with heater-induced whistler waves in the magnetosphere. It is reasonable to suppose that a new substorm activation, exactly above Tromsø, was closely connected with the heater-induced precipitation of energetic electrons.     

Gamma-ray bursts from stellar remnants: probing the Universe at high redshift

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, GRBs should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z  ∼ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied GRB error boxes. The increased distances imply a peak luminosity of 8.3 × 10⁵¹ erg s⁻¹ and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means either that GRBs are caused by a much rarer phenomenon than mergers of binary neutron stars, or that their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.     

Gamma-ray bursts: post-burst evolution of fireballs

The post-burst evolution of fireballs that produce γ-ray bursts (GRBs) is studied, assuming the expansion of fireballs to be adiabatic and relativistic. Numerical results as well as an approximate analytic solution for the evolution are presented. Owing to the adoption of a new relation between t R and γ, our results differ markedly from previous studies. Synchrotron radiation from the shocked interstellar medium is carefully calculated, using a conventional set of equations. The observed X-ray flux of GRB afterglows can be reproduced easily. Although the optical afterglows seem much more complicated, our results can still present a rather satisfactory agreement with observations. We also find that the expansion will no longer be highly relativistic about 4 d after the main GRB. We thus suggest that the marginally relativistic phase of the expansion should be investigated so as to check the afterglows observed a week or more later.     

Oscillations in radiation-driven outflows

In this paper we investigate the dynamical behaviour of radiation-driven winds, specifically winds that arise when Compton scattering transfers momentum from the radiation field to the gas flow. Such winds occur during strong X-ray bursts from slowly accreting neutron stars, and also may be driven from the inner regions of a black hole or neutron star accretion disc when the mass transfer rate is very high. By linearizing the radiation hydrodynamic equations around steady spherical outflow, we evaluate the time-dependent response of these winds to perturbations introduced at their inner boundaries. We find that although radiation-driven winds are generally stable, they act as mechanical filters that should produce quasi-periodic oscillations or peaked noise in their radiation output when perturbations force them stochastically. This behaviour may underlie the photospheric oscillations observed in some strong Type I X-ray bursts.     

Proper motions and variability of the H2 knots in the HH 111 and 121 protostellar outflows

We discuss the formation of spectral features in the decelerating ejecta of gamma-ray bursts, including the possible effect of inhomogeneities. These should lead to blueshifted and broadened absorption edges and resonant features, especially from H and He. An external neutral ISM could produce detectable H and He, as well as Fe X-ray absorption edges and lines. Hypernova scenarios may be diagnosed by Fe Kα and H Lyα emission lines.     

Total flux density radio observations as a tool for understanding AGN behaviour

A wealth of information on the properties of jets in Active Galactic Nuclei (AGNs) can be derived from total flux density observations at high radio frequencies. This includes, for example, the Doppler factor, the Lorentz factor, and the viewing angle of the jet. We have earlier calculated these parameters for a sample of ～80 sources of different AGN types using almost 20 years of 22 and 37 GHz data from Metsähovi Radio Observatory. We have now gathered data for an additional ten years, and studied the long term characteristic variability time scales of a large sample of AGNs using the first order structure function, the discrete autocorrelation function and the Lomb-Scargle periodogram. Some of the results will be presented in this paper. We also stress the importance of long term observations of AGNs, the main reason for this being misinterpretations of source properties due to insufficient time coverage. Only a few observing epochs will too easily lead to incorrect conclusions about variability, continuum spectra, and the general detectability of the source, not to mention the exclusion of interesting objects from further studies. This is particularly important when considering, for example, the Planck satellite for which the quality of the main mission product, the accurate cosmic microwave background anisotropy maps, depends heavily on the elimination of foreground sources such as AGNs.     

Planetary,solar and astrophysical relativistic electrons: Common energization mechanisms?

Direct observations or deduced analysis indicate clearly that formation of intense fluxes of relativistic electrons is an important ingredient in the evolution of numerous active magnetized plasma systems. Examples of relativistic electron energization include the recovery phase of a planetary magnetic storm, post solar flare coronal activity and the afterglow of gamma ray bursts. It is suggested that there exists a universal mechanism, which may explain electron energization at the vastly different magnetized plasma environments. The favorite configuration consists of an inhomogeneous magnetic field anchored at a given magnetic structure and excitation of whistler waves due to external injection of low-energy non-isotropic electrons. The energization proceeds as a bootstrap process due to interaction with the propagating whistler waves along the inhomogeneous magnetic field.     

SWAS observations of water vapor in the Venus mesosphere

We present the first detections of the ground-state H₂¹⁶O (₁₁₀-₁₀₁) rotational transition (at 556.9 GHz) and the ¹³CO (5-4) rotational transition from the atmosphere of Venus, measured with the Submillimeter Wave Astronomy Satellite (SWAS). The observed spectral features of these submillimeter transitions originate primarily from the 70-100 km altitude range, within the Venus mesosphere. Observations were obtained in December 2002, and January, March, and July 2004, coarsely sampling one Venus diurnal period as seen from Earth. The measured water vapor absorption line depth shows large variability among the four observing periods, with strong detections of the line in December 2002 and July 2004, and no detections in January and March 2004. Retrieval of atmospheric parameters was performed using a multi-transition inversion algorithm, combining simultaneous retrievals of temperature, carbon monoxide, and water profiles under imposed constraints. Analysis of the SWAS spectra resulted in measurements or upper limits for the globally averaged mesospheric water vapor abundance for each of the four observation periods, finding variability over at least two orders of magnitude. The results are consistent with both temporal and diurnal variability, but with short-term fluctuations clearly dominating. These results are fully consistent with the long-term study of mesospheric water vapor from millimeter and submillimeter observations of HDO [Sandor, B.J., Clancy, R.T., 2005. Icarus 177, 129-143]. The December 2002 observations detected very rapid change in the mesospheric water abundance. Over five days, a deep water absorption feature consistent with a water vapor abundance of 4.5±1.5 parts per million suddenly gave way to a significantly shallower absorption, implying a decrease in the water vapor abundance by a factor of nearly 50 in less that 48 h. In 2004, similar changes in the water vapor abundance were measured between the March and July SWAS observing periods, but variability on time scales of less than a week was not detected. The mesospheric water vapor is expected to be in equilibrium with aerosol particles, primarily composed of concentrated sulfuric acid, in the upper haze layers of the Venus atmosphere. If true, moderate amplitude (10-15 K) variability in mesospheric temperature, previously noted in millimeter spectroscopy observations of Venus, can explain the rapid water vapor variability detected by SWAS.