Estimation of the vertical wavelength of atmospheric gravity waves from airglow imagery

Abstract In the summer of 2010, two imagers were installed in New Mexico with the objective of making stereoscopic observations of atmospheric gravity waves (AGWs). As AGWs propagate vertically, they spatially perturb the airglow emission layers in all three dimensions. Estimates of the vertical wavelength, horizontal wavelength, and the intrinsic frequency are needed to characterize an AGW and quantify its effects on upper atmospheric dynamics. The dispersion relation describes the relationship between vertical and horizontal wavelengths as a function of the intrinsic frequency. Thus, any two of the three aforementioned parameters can be used to determine the third. Mesospheric winds are hard to measure and consequently the intrinsic frequency is difficult to estimate. However, the horizontal wavelength can be directly measured from airglow imagery once the three-dimensional imager field of view is projected onto the two-dimensional image plane. This thesis presents a method to estimate the vertical wavelength using an airglow perturbation model proposed by Anderson et al. (2009). The model is subsequently validated using the observations from ground-based imagers installed in New Mexico. Abstract The perturbed airglow is modeled as a quasi-monochromatic wave and thus, it can be characterized using only a few parameters, one of which is the vertical wavelength. Because the vertical wavelength is embedded in both the phase and the magnitude of this model, two values of the vertical wavelength are estimated by applying two different parameter estimation techniques on the phase and magnitude. The estimation of the vertical wavelength from the phase of the model entails solving an overdetermined system of linear equations by minimizing the sum of the squared residuals. This estimate is then compared to that obtained by iteratively finding the best approximation to the roots of a function, representing the magnitude of the perturbation model. These two techniques are applied on three nights in 2010, and the estimates for the vertical wavelength match to within a few kilometers. Thus, the perturbation model is validated using real data

Estimation of Gravity Wave Momentum Flux with Spectroscopic Imaging

Atmospheric gravity waves play a significant role in the dynamics and thermal balance of the upper atmosphere. In this paper, we present a novel technique for automated and robust calculation of momentum flux of high-frequency quasi-monochromatic wave components from spectroscopic imaging and horizontal radar wind measurements. Our approach uses the two-dimensional (2-D) cross periodogram of two consecutive Doppler-shifted time-differenced (TD) images to identify wave components and estimate intrinsic wave parameters. Besides estimating the average perturbation of dominant waves in the whole field of view, this technique applies 2-D short-space Fourier transform to the TD images to identify localized wave events. With the wave parameters acquired, the momentum flux carried by all vertically propagating wave components is calculated using an analytical model relating the measured intensity perturbation to the wave amplitude. This model is tested by comparing wave perturbation amplitudes inferred from spectroscopic images with those from sodium lidar temperature measurements. The proposed technique enables characterization of the variations in the direction and strength of gravity waves with high temporal resolution for each clear data-taking night. The nightly results provide statistical information for investigating seasonal and geographical variations in momentum flux of gravity waves

Skycorr: A general tool for spectroscopic sky subtraction

Airglow emission lines, which dominate the optical-to-near-IR sky radiation, show strong, line-dependent variability on various time scales. Therefore, the subtraction of the sky background in the affected wavelength regime becomes a problem if plain sky spectra have to be taken at a different time as the astronomical data. A solution of this issue is the physically motivated scaling of the airglow lines in the plain sky data to fit the sky lines in the object spectrum. We have developed a corresponding instrument-independent approach based on one-dimensional spectra. Our code skycorr separates sky lines and sky/object continuum by an iterative approach involving a line finder and airglow line data. The sky lines are grouped according to their expected variability. The line groups in the sky data are then scaled to fit the sky in the science data. Required pixel-specific weights for overlapping groups are taken from a comprehensive airglow model. Deviations in the wavelength calibration are corrected by fitting Chebyshev polynomials and rebinning via asymmetric damped sinc kernels. The scaled sky lines and the sky continuum are subtracted separately. VLT X-Shooter data covering time intervals from two minutes to about one year were selected to illustrate the performance. Except for short time intervals of a few minutes, the sky line residuals were several times weaker than for sky subtraction without fitting. Further tests show that skycorr performs consistently better than the method of Davies (2007) developed for VLT SINFONI data.Comment: 17 pages, 18 figures, accepted for publication in A&

Tomographic Imaging of Airglow from Airborne Spectroscopic Measurements

A description is given of the methodology based on a single, aircraft-mounted spectroscopic imager to tomographically reconstruct airglow perturbations induced by atmospheric gravity waves. In this configuration, the imager passes under the airglow structure to gather multiple-angle views of the wave structure in a relatively short amount of time. Under the assumption that the airglow structure does not change significantly during the acquisition interval, the data can be tomographically inverted to estimate the 2D (horizontal–vertical) airglow structure. We develop an inversion strategy for this image formation task and illustrate its applicability by inverting time-sequential imaging data taken from different vantage points during the ALOHA-93 campaign to reconstruct atmospheric gravity wave structures

Three-dimensional tomographic reconstruction of mesospheric airglow structures using two-station ground-based image measurements

A new methodology is presented to create two-dimensional (2D) and three-dimensional (3D) tomographic reconstructions of mesospheric airglow layer structure using two-station all-sky image measurements. A fanning technique is presented that produces a series of cross-sectional 2D reconstructions, which are combined to create a 3D mapping of the airglow volume. The imaging configuration is discussed and the inherent challenges of using limited-angle data in tomographic reconstructions have been analyzed using artificially generated imaging objects. An iterative reconstruction method, the partially constrained algebraic reconstruction technique (PCART), was used in conjunction with a priori information of the airglow emission profile to constrain the height of the imaged region, thereby reducing the indeterminacy of the inverse problem. Synthetic projection data were acquired from the imaging objects and the forward problem to validate the tomographic method and to demonstrate the ability of this technique to accurately reconstruct information using only two ground-based sites. Reconstructions of the OH airglow layer were created using data recorded by all-sky CCD cameras located at Bear Lake Observatory, Utah, and at Star Valley, Wyoming, with an optimal site separation of ~100 km. The ability to extend powerful 2D and 3D tomographic methods to two-station ground-based measurements offers obvious practical advantages for new measurement programs. The importance and applications of mesospheric tomographic reconstructions in airglow studies, as well as the need for future measurements and continued development of techniques of this type, are discussed. © 2012 Optical Society of America

The MATS satellite mission - gravity wave studies by Mesospheric Airglow/Aerosol Tomography and Spectroscopy

Global three-dimensional data are a key to understanding gravity waves in the mesosphere and lower thermosphere. MATS (Mesospheric Airglow/Aerosol Tomography and Spectroscopy) is a new Swedish satellite mission that addresses this need. It applies space-borne limb imaging in combination with tomographic and spectroscopic analysis to obtain gravity wave data on relevant spatial scales. Primary measurement targets are O-2 atmospheric band dayglow and nightglow in the near infrared, and sunlight scattered from noctilucent clouds in the ultraviolet. While tomography provides horizontally and vertically resolved data, spectroscopy allows analysis in terms of mesospheric temperature, composition, and cloud properties. Based on these dynamical tracers, MATS will produce a climatology on wave spectra during a 2-year mission. Major scientific objectives include a characterization of gravity waves and their interaction with larger-scale waves and mean flow in the mesosphere and lower thermosphere, as well as their relationship to dynamical conditions in the lower and upper atmosphere. MATS is currently being prepared to be ready for a launch in 2020. This paper provides an overview of scientific goals, measurement concepts, instruments, and analysis ideas

衛星搭載分光器による中間圏同心状大気重力波の空間・時間変動

Tohoku University坂野井健課

The Application of Tomographic Reconstruction Techniques to Ill-Conditioned Inverse Problems in Atmospheric Science and Biomedical Imaging

A methodology is presented for creating tomographic reconstructions from various projection data, and the relevance of the results to applications in atmospheric science and biomedical imaging is analyzed. The fundamental differences between transform and iterative methods are described and the properties of the imaging configurations are addressed. The presented results are particularly suited for highly ill-conditioned inverse problems in which the imaging data are restricted as a result of poor angular coverage, limited detector arrays, or insufficient access to an imaging region. The class of reconstruction algorithms commonly used in sparse tomography, the algebraic reconstruction techniques, is presented, analyzed, and compared. These algorithms are iterative in nature and their accuracy depends significantly on the initialization of the algorithm, the so-called initial guess. A considerable amount of research was conducted into novel initialization techniques as a means of improving the accuracy. The main body of this paper is comprised of three smaller papers, which describe the application of the presented methods to atmospheric and medical imaging modalities. The first paper details the measurement of mesospheric airglow emissions at two camera sites operated by Utah State University. Reconstructions of vertical airglow emission profiles are presented, including three-dimensional models of the layer formed using a novel fanning technique. The second paper describes the application of the method to the imaging of polar mesospheric clouds (PMCs) by NASA’s Aeronomy of Ice in the Mesosphere (AIM) satellite. The contrasting elements of straight-line and diffusive tomography are also discussed in the context of ill-conditioned imaging problems. A number of developing modalities in medical tomography use near-infrared light, which interacts strongly with biological tissue and results in significant optical scattering. In order to perform tomography on the diffused signal, simulations must be incorporated into the algorithm, which describe the sporadic photon migration. The third paper presents a novel Monte Carlo technique derived from the optical scattering solution for spheroidal particles designed to mimic mitochondria and deformed cell nuclei. Simulated results of optical diffusion are presented. The potential for improving existing imaging modalities through continual development of sparse tomography and optical scattering methods is discussed

Observation and Modeling of Gravity Wave Propagation through Reflection and Critical Layers above Andes Lidar Observatory at Cerro Pachón, Chile

A complex gravity wave event was observed from 04:30 to 08:10 UTC on 16 January 2015 by a narrow-band sodium lidar and an all-sky airglow imager located at Andes Lidar Observatory (ALO) in Cerro Pachón (30.25∘S, 70.73∘W), Chile. The gravity wave packet had a period of 18–35 min and a horizontal wavelength of about 40–50 km. Strong enhancements of the vertical wind perturbation, exceeding10 m s−1, were found at ∼90 km and ∼103 km, consistent with nearly evanescent wave behavior near a reflection layer. A reduction in vertical wavelength was found as the phase speed approached the background wind speed near ∼93 km. A distinct three-layered structure was observed in the lidar data due to refraction of the wave packet. A fully nonlinear model was used to simulate this event, which successfully reproduced the amplitudes and layered structure seen in observations. The model results provide dynamical insight, suggesting that a double reflection occurring at two separate heights caused the large vertical wind amplitudes, while the three-layered structure in the temperature perturbation was a result of relatively stable regions at those altitudes. The event provides a clear perspective on the filtering processes to which short-period, small-scale gravity waves are subject in mesosphere and lower thermosphere

Gravity wave characteristics in the mesopause region revealed from OH airglow imager observations over Northern Colorado

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Space Physics 119 (2014): 630-645, doi:10.1002/2013JA018955.Using 5 years of all-sky OH airglow imager data over Yucca Ridge Field Station, CO (40.7°N, 104.9°W), from September 2003 to September 2008, we extract and deduce quasi-monochromatic gravity wave (GW) characteristics in the mesopause region. The intrinsic periods are clustered between approximately 4 and 10 min, and many of them are unstable and evanescent. GW occurrence frequency exhibits a clear semiannual variation with equinoctial minima, which is likely related to the seasonal variation of background wind. The anomalous propagation direction in January 2006, with strong southward before major warming starting in 21 January and weak southward propagation afterward, was most likely affected by stratospheric sudden warming. The momentum fluxes show strongly anticorrelated with the tides, with ~180° out of phase in the zonal component. While in the meridional component, the easterly maximum occurred approximately 2–6 h after maximum easterly tidal wind. However, the anticorrelations are both weakest during the summer. The dissipating and breaking of small-scale and high-frequency GW's components could have a potential impact on the general circulation in the mesopause region.This work was carried out at the University of Science and Technology of China, with support from the National Natural Science Foundation of China grants (41025016, 41127901, 41225017, 41074108, and 41121003), the National Basic Research Program of China grant 2012CB825605, the Chinese Academy of Sciences Key Research Program KZZD-EW-01, and the Fundamental Research Funds for the Central Universities.2014-07-3