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.     

CRISTA-NF measurements of water vapor during the SCOUT-O3 Tropical Aircraft Campaign

The new remote sensing experiment CRISTA-NF (Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere – New Frontiers) successfully participated in the SCOUT-O3 Tropical Aircraft Campaign in November and December 2005. CRISTA-NF operated aboard the high-altitude research aircraft M-55 Geophysica. Mid-infrared spectra (4–15 μm) were measured in the limb sounding geometry with high spatial resolution (250 m vertical sampling, 5–15 km along track sampling). Measurements were carried out during transfer flights between Oberpfaffenhofen, Germany, and Darwin, Australia, as well as during several local flights near Darwin. Water vapor volume mixing ratios in the upper troposphere and lower stratosphere were derived from the CRISTA-NF radiance measurements by utilizing a rapid radiative transfer forward model and the optimal estimation retrieval approach. CRISTA-NF water vapor measurements below the hygropause have a total retrieval error of 15–40% (i.e. root mean square of accuracy and precision). The systematic terms are dominating in the retrieval error budget. The contributions of a priori information to the retrieval results are less than 5–10%. The vertical resolution of the observations is about 250–500 m when permitted by instrument sampling. In this paper we present first results for three transfer flights of the campaign. Being generally in good agreement with corresponding ECMWF operational analyzes, the CRISTA-NF measurements show significantly higher variability and local structures in the upper tropospheric water vapor distributions.     

Measurements of the upper troposphere and lower stratosphere during tropical cyclones using the GPS radio occultation technique

Water vapour transport to the upper troposphere and lower stratosphere by deep convective storms affects the radiation balance of the atmosphere and has been proposed as an important component of climate change. The aim of the work presented here is to understand if the GPS radio occultation technique is useful for characterization of this process. Our assessment addresses the question if severe storms leave a significant signature in radio occultation profiles in the upper troposphere/lower stratosphere. Radio occultation data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) were analyzed, focusing on two particular tropical cyclones with completely different characteristics, the hurricane Bertha, which formed in the Atlantic Basin during July 2008 and reached a maximum intensity of Category 3, and the typhoon Hondo, which formed in the south Indian Ocean during 2008 reaching a maximum intensity of Category 4. The result is positive, suggesting that the bending angle of a GPS radio occultation signal contains interesting information on the atmosphere around the tropopause, but not any information regarding the water vapour. The maximum percentage anomaly of bending angle between 14 and 18 km of altitude during tropical cyclones is typically larger than the annual mean by 5–15% and it can reach 20% for extreme cases. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.     

Development of a GNSS water vapour tomography system using algebraic reconstruction techniques

A GNSS water vapour tomography system developed to reconstruct spatially resolved humidity fields in the troposphere is described. The tomography system was designed to process the slant path delays of about 270 German GNSS stations in near real-time with a temporal resolution of 30 min, a horizontal resolution of 40 km and a vertical resolution of 500 m or better. After a short introduction to the GPS slant delay processing the framework of the GNSS tomography is described in detail. Different implementations of the iterative algebraic reconstruction techniques (ART) used to invert the linear inverse problem are discussed. It was found that the multiplicative techniques (MART) provide the best results with least processing time, i.e., a tomographic reconstruction of about 26,000 slant delays on a 8280 cell grid can be obtained in less than 10 min. Different iterative reconstruction techniques are compared with respect to their convergence behaviour and some numerical parameters. The inversion can be considerably stabilized by using additional non-GNSS observations and implementing various constraints. Different strategies for initialising the tomography and utilizing extra information are discussed. At last an example of a reconstructed field of the wet refractivity is presented and compared to the corresponding distribution of the integrated water vapour, an analysis of a numerical weather model (COSMO-DE) and some radiosonde profiles.     

Sensitivity analysis for air temperature profile estimation methods around the tropopause using simulated Aqua/AIRS data

A sensitivity analysis is performed to investigate potential improvements to the accuracy of air temperature profile retrievals near the tropopause. A simple inversion method is employed to identify and remove redundant spectral channels from the retrievals using simulated data for the high-spectral resolution sounder AIRS (Atmospheric Infrared Sounder) on the Aqua satellite. Bayesian optimal theory and inverse technique are applied for the atmospheric temperature profile retrievals, and the 15 μm CO₂ absorption bands (620–750 cm⁻¹) are chosen for this study. Sequentially elimination of redundant channels is directly integrated into the inverse scheme for the temperature profile, in order to accurately retain the valuable channels and remove all the redundant channels, for accurate retrieval of the temperature profile. Also, the tropopause and troposphere are treated differently in the inverse scheme to improve the retrieval accuracy in the tropopause. Results of a sensitivity analysis based on this method, for the Tropical and Middle-Latitude Summer models simulated by MODTRAN4.0, show that the estimated accuracies are improved by 2 K around the tropopause, and are only changed by less than 0.2 K in the troposphere.     

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

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

Validation of space-born Moderate Resolution Imaging Spectroradiometer remote sensors aerosol products using application of ground-based Multi-wavelength Radiometer

The measurements of aerosol optical properties were carried out during April 2006 to March 2011 over Mohal (31.9°N, 77.12°E) in the northwestern Indian Himalaya, using the application of ground-based Multi-wavelength Radiometer (MWR) and space-born Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensors. The average (±standard deviation) values of aerosol optical depth (AOD) at 500 nm, Ångström exponent and turbidity coefficient during the entire measurement period were 0.25 ± 0.09, 1.15 ± 0.42 and 0.12 ± 0.06 respectively. About 86% AOD values retrieved from MODIS remote sensor were found within an uncertainty limit (Δτ_pλ = ±0.05 ± 0.15τ_pλ). In general, the MWR derived AOD values were higher than that of MODIS retrieval with absolute difference ∼0.02. During the entire period of measurement space-born MODIS remote sensor and ground-based MWR observation showed good correspondence with significant correlation coefficient ∼0.78 and root mean square difference ∼0.06. For daily observations the relative difference between these two estimates stood less than 9%. However, satellite-based and ground-based observation showed good correspondence, but further efforts still needed to eliminate systematic errors in the existing MODIS algorithm.     

Atmosphere sounding by GPS radio occultation: First results from TanDEM-X and comparison with TerraSAR-X

On 21 June 2010 the TerraSAR-X satellite was joined by the TanDEM-X satellite. A Global Positioning System (GPS) radio occultation (RO) experiment using the twin satellites has been carried out to estimate the precision of GPS atmospheric soundings. For the Day Of Year (DOY) 330–336, 2011, we analyze phase and amplitude data recorded by GPS receivers separated by a few hundred meters in a low earth orbit and derive collocated atmospheric refractivity profiles. In the altitude range 10–20 km the standard deviation between TerraSAR-X and TanDEM-X refractivity does not exceed 0.15%. The standard deviation is rapidly increasing for lower and higher altitudes; close to the surface and at an altitude of 30 km the standard deviation reaches 0.8% and 0.5%, respectively. Systematic deviations between TerraSAR-X and TanDEM-X refractivity in the considered altitude range (0–30 km) are negligible. The results confirm the anticipated high precision of the GPS RO technique. However, the difference in the retrieved refractivity in the lower troposphere for different Open Loop (OL) signal tracking parameters, altered onboard TanDEM-X for DOY 49–55, 2012, calls for an in depth analysis. At the moment we can not exclude that a potential bias in the OL Doppler model introduces a bias in our retrieved refractivity at altitudes <8$<8$ km.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

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.     

北京上空高分辨率气球探空观测的温度垂直波数谱

使用分辨率为10 m的温度剖面检查2.90～8.01 km和14.65～19.76 km高度范围内归一化温度扰动谱的谱特性,并且将它们与模式谱比较.结果表明,在波数9.8×10-4～2.5×10-2m-1范围内,对流层的垂直波数谱有约-1.9的平均谱斜率;低平流层的垂直波数谱有约-2.2的平均谱斜率.这两个平均谱斜率大大偏离了目前饱和谱模式和普适大气谱模式表示的-3.0和-2.4谱斜率,并且认为是目前气球测量中曾经观测到的最平谱斜率.在波数9.8×10-4～2.5×10-2m-1范围内,对流层和低平流层的平均谱振幅比饱和谱模式的谱振幅分别大24倍和5倍,这也与饱和谱模式以及目前文献中的观测结果有很大不同.这些较大平均谱振幅与最平的平均谱斜率一道共同表明,观测的温度谱并不遵循目前的饱和谱模式和普适大气谱模式.      

Interseismic slip rate of the Garze–Yushu fault belt in the Tibetan Plateau from C-band InSAR observations between 2003 and 2010

InSAR time series techniques can provide high-spatial resolution deformation fields across an active fault belt, even for zones with heavy vegetation coverage. An interseismic deformation map across the Garze–Yushu fault belt in the Tibetan Plateau, ∼300 km by ∼100 km, is derived from C-band Envisat/ASAR imagery collected between 2003 and 2010. Unlike previous research, we obtain a lookup figure which relates the slip rate with the fault locking depth, the dip angle and the rake angle. The estimated slip rate changes significantly with the locking depth and the rake angle, but relatively little with the dip angle. When considering the focal mechanism solutions of historical earthquake along the Garze–Yushu fault, the interseismic slip rate of the Garze–Yushu fault is close to a value of 6.4 mm/yr, which is between the highest (18.2 mm/yr) and the lowest (3.1 mm/yr) slip rate from GPS estimations, but slightly less than the minimum value (∼ 7 mm/yr) from the geological estimations. The earthquake recurrence interval on the Yushu part of Garze–Yushu fault equals 272 yr, and the April 14, 2010 Mw 6.9 earthquake has not completely released the accumulated strain energy between 1738 and 2010.     

Mesospheric doppler wind measurements from Aura Microwave Limb Sounder (MLS)

This paper describes a microwave limb technique for measuring Doppler wind in the Earth’s mesosphere. The research algorithm has been applied to Aura Microwave Limb Sounder (MLS) 118.75 GHz measurements where the O₂ Zeeman lines are resolved by a digital autocorrelation spectrometer. A precision of ∼17 m/s for the line-of-sight (LOS) wind is achieved at 80–92 km, which corresponds to radiometric noise during 1/6 s integration time. The LOS winds from Aura MLS are mostly in the meridional direction at low- and mid-latitudes with vertical resolution of ∼8 km. This microwave Doppler technique has potential to obtain useful winds down to ∼40 km of the Earth’s atmosphere if measurements from other MLS frequencies (near H₂O, O₃, and CO lines) are used. Initial analyses show that the MLS winds from the 118.75 GHz measurements agree well with the TIDI (Thermosphere Ionosphere Mesosphere Energetics and Dynamics Doppler Interferometer) winds for the perturbations induced by a strong quasi 2-day wave (QTDW) in January 2005. Time series of MLS winds reveal many interesting climatological and planetary wave features, including the diurnal, semidiurnal tides, and the QTDW. Interactions between the tides and the QTDW are clearly evident, indicating possible large tidal structural changes after the QTDW events dissipate.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Fault rupture model of the 2008 Dangxiong (Tibet,China) Mw 6.3 earthquake from Envisat and ALOS data

On October 6, 2008, an Mw 6.3 earthquake occurred in Dangxiong county, southern Tibetan Plateau. In this study, Synthetic Aperture Radar (SAR) images from Envisat ASAR C-band descending Track 176 and ALOS PALSAR L-band ascending Track 500 are processed to generate the coseismic deformation caused by the earthquake. To estimate the source model, a downhill simplex non-linear inversion method is used to determine the fault rupture geometry, and an automatic fault discretization technique is employed to divide the fault plane to construct the optimal slip model, in which the uncertainties of the fault parameters are assessed by a Monte Carlo method. The inversion results show that the earthquake strikes almost south–north and has a normal faulting focal mechanism with rake angle and slip of −111.7° and 1.33 m, respectively. Peak slip of 2.15 m is located at a depth of 7.5 km. The estimated geodetic moment is 4.06 × 10¹⁸ N m (Mw 6.37), 71.2% of which is released in the depth range 4.5–11 km. The slip model suggests that coseismic slip also takes place at some fault patches near the earth’s surface and postseismic afterslip occurs below the coseismic rupture area after the earthquake.     

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.     

Improved precision of virtual height measurements with coherent radio pulse sounding based on the maximum likelihood method

The paper describes a new technique that improves precision of the virtual height measurements by a coherent pulse sounding of the ionosphere. Proposed technique is based on the method of maximum likelihood that matches expected and observed spectral domain signatures of the signal intermixed with the noise. Computer simulations show that our technique allows measurements of the echo virtual height with ∼100 m precision even at a much coarser step of the height sampling in the sounder. In experiment, we expect an average 300 m precision of the virtual height measurements for single echoes received during periods of little spread due to ionospheric irregularities.     

Signatures of the Sudden Stratospheric Warming events of January–February 2008 in Seoul,S. Korea

The period January–February 2008 was characterized by four Sudden Stratospheric Warmings (SSWs) in the Northern Hemisphere, of which the last warming, at the end of February 2008, was a major warming. A significant decrease in mesospheric water vapour (H₂O) of more than 2 ppmv (∼40%) was observed by the ground-based microwave (GBMW) radiometer in Seoul, S. Korea [37.3°N, 126.3°E] during the major SSW. A comparison with ground-based mesospheric H₂O observations from the mid-latitude station in Bern [46.9°N, 7°E] revealed an anticorrelation in the mesospheric H₂O data during the major SSW. In addition, prior to the major warming, strong periodic fluctuations were recorded in the Aura MLS vertical temperature distribution between 15 and 0.05 hPa at Seoul. The mesospheric temperature oscillation was found to have a period of ∼10–14 days with a persistency of 3–4 cycles.     

Single photon laser altimeter simulator and statistical signal processing

Spaceborne altimeters are common instruments onboard the deep space rendezvous spacecrafts. They provide range and topographic measurements critical in spacecraft navigation. Simultaneously, the receiver part may be utilized for Earth-to-satellite link, one way time transfer, and precise optical radiometry. The main advantage of single photon counting approach is the ability of processing signals with very low signal-to-noise ratio eliminating the need of large telescopes and high power laser source. Extremely small, rugged and compact microchip lasers can be employed. The major limiting factor, on the other hand, is the acquisition time needed to gather sufficient volume of data in repetitive measurements in order to process and evaluate the data appropriately. Statistical signal processing is adopted to detect signals with average strength much lower than one photon per measurement. A comprehensive simulator design and range signal processing algorithm are presented to identify a mission specific altimeter configuration. Typical mission scenarios (celestial body surface landing and topographical mapping) are simulated and evaluated. The high interest and promising single photon altimeter applications are low-orbit (∼10 km) and low-radial velocity (several m/s) topographical mapping (asteroids, Phobos and Deimos) and landing altimetry (∼10 km) where range evaluation repetition rates of ∼100 Hz and 0.1 m precision may be achieved. Moon landing and asteroid Itokawa topographical mapping scenario simulations are discussed in more detail.