scholarly journals Atmospheric tomography using the Nordic Meteor Radar Cluster and Chilean Observation Network De Meteor Radars: network details and 3D-Var retrieval

2021 ◽  
Vol 14 (10) ◽  
pp. 6509-6532
Author(s):  
Gunter Stober ◽  
Alexander Kozlovsky ◽  
Alan Liu ◽  
Zishun Qiao ◽  
Masaki Tsutsumi ◽  
...  

Abstract. Ground-based remote sensing of atmospheric parameters is often limited to single station observations by vertical profiles at a certain geographic location. This is a limiting factor for investigating gravity wave dynamics as the spatial information is often missing, e.g., horizontal wavelength, propagation direction or intrinsic frequency. In this study, we present a new retrieval algorithm for multistatic meteor radar networks to obtain tomographic 3-D wind fields within a pre-defined domain area. The algorithm is part of the Agile Software for Gravity wAve Regional Dynamics (ASGARD) and called 3D-Var, and based on the optimal estimation technique and Bayesian statistics. The performance of the 3D-Var retrieval is demonstrated using two meteor radar networks: the Nordic Meteor Radar Cluster and the Chilean Observation Network De Meteor Radars (CONDOR). The optimal estimation implementation provide statistically sound solutions and diagnostics from the averaging kernels and measurement response. We present initial scientific results such as body forces of breaking gravity waves leading to two counter-rotating vortices and horizontal wavelength spectra indicating a transition between the rotational k−3 and divergent k-5/3 mode at scales of 80–120 km. In addition, we performed a keogram analysis over extended periods to reflect the latitudinal and temporal impact of a minor sudden stratospheric warming in December 2019. Finally, we demonstrate the applicability of the 3D-Var algorithm to perform large-scale retrievals to derive meteorological wind maps covering a latitude region from Svalbard, north of the European Arctic mainland, to central Norway.

2021 ◽  
Author(s):  
Gunter Stober ◽  
Alexander Kozlovsky ◽  
Alan Liu ◽  
Zishun Qiao ◽  
Masaki Tsutsumi ◽  
...  

Abstract. Ground-based remote sensing of atmospheric parameters is often limited to single station observations of vertical profiles at a certain geographic location. This can be a limiting factor to investigating gravity wave dynamics. In this study we present a new retrieval algorithm for multi-static meteor radar networks to obtain tomographic 3D wind fields within a pre-defined domain area. The algorithm is part of the Agile Software for Gravity wAve Regional Dynamics (ASGARD) called 3DVAR, and based on the optimal estimation technique and Bayesian statistics. The performance of the 3DVAR retrieval is demonstrated using two meteor radar networks, the Nordic Meteor Radar Cluster and the Chilean Observation Network De MeteOr Radars (CONDOR). The optimal estimation implementation provides a statistically sound solution and additional diagnostics from the averaging kernels and measurement response. We present initial scientific results such as body forces of breaking gravity waves leading to two counter-rotating vortices and horizontal wavelength spectra indicating a transition between the vortical κ−3 and divergent κ−5/3 mode at scales of 80–120 km. In addition, we have performed a keogram analysis over extended periods to reflect the latitudinal and temporal impact of a minor sudden stratospheric warming in December 2019. Finally, we demonstrate the applicability of the 3DVAR algorithm to perform large-scale retrievals to derive meteorological wind maps covering a latitude region from Svalbard, north of the European Arctic mainland, to mid-Norway.


2016 ◽  
Vol 9 (3) ◽  
pp. 877-908 ◽  
Author(s):  
Corwin J. Wright ◽  
Neil P. Hindley ◽  
Andrew C. Moss ◽  
Nicholas J. Mitchell

Abstract. Gravity waves in the terrestrial atmosphere are a vital geophysical process, acting to transport energy and momentum on a wide range of scales and to couple the various atmospheric layers. Despite the importance of these waves, the many studies to date have often exhibited very dissimilar results, and it remains unclear whether these differences are primarily instrumental or methodological. Here, we address this problem by comparing observations made by a diverse range of the most widely used gravity-wave-resolving instruments in a common geographic region around the southern Andes and Drake Passage, an area known to exhibit strong wave activity. Specifically, we use data from three limb-sounding radiometers (Microwave Limb Sounder, MLS-Aura; HIgh Resolution Dynamics Limb Sounder, HIRDLS; Sounding of the Atmosphere using Broadband Emission Radiometry, SABER), the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) GPS-RO constellation, a ground-based meteor radar, the Advanced Infrared Sounder (AIRS) infrared nadir sounder and radiosondes to examine the gravity wave potential energy (GWPE) and vertical wavelengths (λz) of individual gravity-wave packets from the lower troposphere to the edge of the lower thermosphere ( ∼  100 km). Our results show important similarities and differences. Limb sounder measurements show high intercorrelation, typically  > 0.80 between any instrument pair. Meteor radar observations agree in form with the limb sounders, despite vast technical differences. AIRS and radiosonde observations tend to be uncorrelated or anticorrelated with the other data sets, suggesting very different behaviour of the wave field in the different spectral regimes accessed by each instrument. Evidence of wave dissipation is seen, and varies strongly with season. Observed GWPE for individual wave packets exhibits a log-normal distribution, with short-timescale intermittency dominating over a well-repeated monthly-median seasonal cycle. GWPE and λz exhibit strong correlations with the stratospheric winds, but not with local surface winds. Our results provide guidance for interpretation and intercomparison of such data sets in their full context.


2011 ◽  
Vol 73 (9) ◽  
pp. 911-920 ◽  
Author(s):  
Manja Placke ◽  
Peter Hoffmann ◽  
Erich Becker ◽  
Christoph Jacobi ◽  
Werner Singer ◽  
...  

2021 ◽  
Author(s):  
Arno Keppens ◽  
Jean-Christopher Lambert ◽  
Daan Hubert ◽  
Steven Compernolle ◽  
Tijl Verhoelst ◽  
...  

<p>Part of the space segment of EU’s Copernicus Earth Observation programme, the Sentinel-5 Precursor (S5P) mission is dedicated to global and European atmospheric composition measurements of air quality, climate and the stratospheric ozone layer. On board of the S5P early afternoon polar satellite, the imaging spectrometer TROPOMI (TROPOspheric Monitoring Instrument) performs nadir measurements of the Earth radiance within the UV-visible and near-infrared spectral ranges, from which atmospheric ozone profile data are retrieved. Developed at the Royal Netherlands Meteorological Institute (KNMI) and based on the optimal estimation method, TROPOMI’s operational ozone profile retrieval algorithm has recently been upgraded. With respect to early retrieval attempts, accuracy is expected to have improved significantly, also thanks to recent updates of the TROPOMI Level-1b data product. This work reports on the initial validation of the improved TROPOMI height-resolved ozone data in the troposphere and stratosphere, as collected both from the operational S5P Mission Performance Centre/Validation Data Analysis Facility (MPC/VDAF) and from the S5PVT scientific project CHEOPS-5p. Based on the same validation best practices as developed for and applied to heritage sensors like GOME-2, OMI and IASI (Keppens et al., 2015, 2018), the validation methodology relies on the analysis of data retrieval diagnostics – like the averaging kernels’ information content – and on comparisons of TROPOMI data with reference ozone profile measurements. The latter are acquired by ozonesonde, stratospheric lidar, and tropospheric lidar stations performing network operation in the context of WMO's Global Atmosphere Watch and its contributing networks NDACC and SHADOZ. The dependence of TROPOMI’s ozone profile uncertainty on several influence quantities like cloud fraction and measurement parameters like sun and scan angles is examined and discussed. This work concludes with a set of quality indicators, enabling users to verify the fitness-for-purpose of the S5P data.</p>


2018 ◽  
Vol 8 (10) ◽  
pp. 1797 ◽  
Author(s):  
Zhuolei Xiao ◽  
Yerong Zhang ◽  
Kaixuan Zhang ◽  
Dongxu Zhao ◽  
Guan Gui

The goal of phase retrieval is to recover an unknown signal from the random measurements consisting of the magnitude of its Fourier transform. Due to the loss of the phase information, phase retrieval is considered as an ill-posed problem. Conventional greedy algorithms, e.g., greedy spare phase retrieval (GESPAR), were developed to solve this problem by using prior knowledge of the unknown signal. However, due to the defect of the Gauss–Newton method in the local convergence problem, especially when the residual is large, it is very difficult to use this method in GESPAR to efficiently solve the non-convex optimization problem. In order to improve the performance of the greedy algorithm, we propose an improved phase retrieval algorithm, which is called the greedy autocorrelation retrieval Levenberg–Marquardt (GARLM) algorithm. Specifically, the proposed GARLM algorithm is a local search iterative algorithm to recover the sparse signal from its Fourier transform magnitude. The proposed algorithm is preferred to existing greedy methods of phase retrieval, since at each iteration the problem of minimizing the objective function over a given support is solved by using the improved Levenberg–Marquardt (ILM) method and matrix transform. A local search procedure such as the 2-opt method is then invoked to get the optimal estimation. Simulation results are given to show that the proposed algorithm performs better than the conventional GESPAR algorithm.


2009 ◽  
Vol 27 (4) ◽  
pp. 1477-1487 ◽  
Author(s):  
H. Takahashi ◽  
M. J. Taylor ◽  
P.-D. Pautet ◽  
A. F. Medeiros ◽  
D. Gobbi ◽  
...  

Abstract. During the Spread F Experiment campaign, under NASA Living with a Star (LWS) program, carried out in the South American Magnetic Equator region from 22 September to 8 November 2005, two airglow CCD imagers, located at Cariri (7.4° S, 36.5° W, geomag. 11° S) and near Brasilia (14.8° S, 47.6° W, geomag. 10° S) were operated simultaneously and measured the equatorial ionospheric bubbles and their time evolution by monitoring the airglow OI 6300 intensity depletions. Simultaneous observation of the mesospheric OH wave structures made it possible to investigate the relationship between the bubble formation in the ionosphere and the gravity wave activity at around 90 km. On the evening of 30 September 2005, comb-like OI 6300 depletions with a distance of ~130 km between the adjacent ones were observed. During the same period, a mesospheric gravity wave with a horizontal wavelength of ~130 km was observed. From the 17 nights of observation during the campaign period, there was a good correlation between the OI 6300 depletion distances and the gravity wave horizontal wavelengths in the mesosphere with a statistically significant level, suggesting a direct contribution of the mesospheric gravity wave to plasma bubble seeding in the equatorial ionosphere.


2019 ◽  
Vol 12 (11) ◽  
pp. 5801-5816 ◽  
Author(s):  
Shayamila Mahagammulla Gamage ◽  
Robert J. Sica ◽  
Giovanni Martucci ◽  
Alexander Haefele

Abstract. We present a new method for retrieving temperature from pure rotational Raman (PRR) lidar measurements. Our optimal estimation method (OEM) used in this study uses the full physics of PRR scattering and does not require any assumption of the form for a calibration function nor does it require fitting of calibration factors over a large range of temperatures. The only calibration required is the estimation of the ratio of the lidar constants of the two PRR channels (coupling constant) that can be evaluated at a single or multiple height bins using a simple analytic expression. The uncertainty budget of our OEM retrieval includes both statistical and systematic uncertainties, including the uncertainty in the determination of the coupling constant on the temperature. We show that the error due to calibration can be reduced significantly using our method, in particular in the upper troposphere when calibration is only possible over a limited temperature range. Some other advantages of our OEM over the traditional Raman lidar temperature retrieval algorithm include not requiring correction or gluing to the raw lidar measurements, providing a cutoff height for the temperature retrievals that specifies the height to which the retrieved profile is independent of the a priori temperature profile, and the retrieval's vertical resolution as a function of height. The new method is tested on PRR temperature measurements from the MeteoSwiss RAman Lidar for Meteorological Observations system in clear and cloudy sky conditions, compared to temperature calculated using the traditional PRR calibration formulas, and validated with coincident radiosonde temperature measurements in clear and cloudy conditions during both daytime and nighttime.


Atmosphere ◽  
2019 ◽  
Vol 10 (2) ◽  
pp. 81 ◽  
Author(s):  
Shaohua Gong ◽  
Guotao Yang ◽  
Jiyao Xu ◽  
Xiao Liu ◽  
Qinzeng Li

A low-frequency inertial atmospheric gravity wave (AGW) event was studied with lidar (40.5° N, 116° E), meteor radar (40.3° N, 116.2° E), and TIMED/SABER at Beijing on 30 May 2012. Lidar measurements showed that the atmospheric temperature structure was persistently perturbed by AGWs propagating upward from the stratosphere into the mesosphere (35–86 km). The dominant contribution was from the waves with vertical wavelengths λ z = 8 − 10   km and wave periods T ob = 6.6 ± 0.7   h . Simultaneous observations from a meteor radar illustrated that MLT horizontal winds were perturbed by waves propagating upward with an azimuth angle of θ = 247 ° , and the vertical wavelength ( λ z = 10   km ) and intrinsic period ( T in = 7.4   h ) of the dominant waves were inferred with the hodograph method. TIMED/SABER measurements illustrated that the vertical temperature profiles were also perturbed by waves with dominant vertical wavelength λ z = 6 − 9   km . Observations from three different instruments were compared, and it was found that signatures in the temperature perturbations and horizontal winds were induced by identical AGWs. According to these coordinated observation results, the horizontal wavelength and intrinsic phase speed were inferred to be ~560 km and ~21 m/s, respectively. Analyses of the Brunt-Väisälä frequency and potential energy illustrated that this persistent wave propagation had good static stability.


2009 ◽  
Vol 2 (2) ◽  
pp. 679-701 ◽  
Author(s):  
G. E. Thomas ◽  
C. A. Poulsen ◽  
A. M. Sayer ◽  
S. H. Marsh ◽  
S. M. Dean ◽  
...  

Abstract. The aerosol component of the Oxford-Rutherford Aerosol and Cloud (ORAC) combined cloud and aerosol retrieval scheme is described and the theoretical performance of the algorithm is analysed. ORAC is an optimal estimation retrieval scheme for deriving cloud and aerosol properties from measurements made by imaging satellite radiometers and, when applied to cloud free radiances, provides estimates of aerosol optical depth at a wavelength of 550 nm, aerosol effective radius and surface reflectance at 550 nm. The aerosol retrieval component of ORAC has several incarnations – this paper addresses the version which operates in conjunction with the cloud retrieval component of ORAC (described by Watts et al., 1998), as applied in producing the Global Retrieval of ATSR Cloud Parameters and Evaluation (GRAPE) data-set. The algorithm is described in detail and its performance examined. This includes a discussion of errors resulting from the formulation of the forward model, sensitivity of the retrieval to the measurements and a priori constraints, and errors resulting from assumptions made about the atmospheric/surface state.


2010 ◽  
Vol 10 (9) ◽  
pp. 4295-4317 ◽  
Author(s):  
D. Wurl ◽  
R. G. Grainger ◽  
A. J. McDonald ◽  
T. Deshler

Abstract. Stratospheric aerosol particles under non-volcanic conditions are typically smaller than 0.1 μm. Due to fundamental limitations of the scattering theory in the Rayleigh limit, these tiny particles are hard to measure by satellite instruments. As a consequence, current estimates of global aerosol properties retrieved from spectral aerosol extinction measurements tend to be strongly biased. Aerosol surface area densities, for instance, are observed to be about 40% smaller than those derived from correlative in situ measurements (Deshler et al., 2003). An accurate knowledge of the global distribution of aerosol properties is, however, essential to better understand and quantify the role they play in atmospheric chemistry, dynamics, radiation and climate. To address this need a new retrieval algorithm was developed, which employs a nonlinear Optimal Estimation (OE) method to iteratively solve for the monomodal size distribution parameters which are statistically most consistent with both the satellite-measured multi-wavelength aerosol extinction data and a priori information. By thus combining spectral extinction measurements (at visible to near infrared wavelengths) with prior knowledge of aerosol properties at background level, even the smallest particles are taken into account which are practically invisible to optical remote sensing instruments. The performance of the OE retrieval algorithm was assessed based on synthetic spectral extinction data generated from both monomodal and small-mode-dominant bimodal sulphuric acid aerosol size distributions. For monomodal background aerosol, the new algorithm was shown to fairly accurately retrieve the particle sizes and associated integrated properties (surface area and volume densities), even in the presence of large extinction uncertainty. The associated retrieved uncertainties are a good estimate of the true errors. In the case of bimodal background aerosol, where the retrieved (monomodal) size distributions naturally differ from the correct bimodal values, the associated surface area (A) and volume densities (V) are, nevertheless, fairly accurately retrieved, except at values larger than 1.0 μm2 cm−3 (A) and 0.05 μm3 cm−3 (V), where they tend to underestimate the true bimodal values. Due to the limited information content in the SAGE II spectral extinction measurements this kind of forward model error cannot be avoided here. Nevertheless, the retrieved uncertainties are a good estimate of the true errors in the retrieved integrated properties, except where the surface area density exceeds the 1.0 μm2 cm−3 threshold. When applied to near-global SAGE II satellite extinction measured in 1999 the retrieved OE surface area and volume densities are observed to be larger by, respectively, 20–50% and 10–40% compared to those estimates obtained by the SAGE~II operational retrieval algorithm. An examination of the OE algorithm biases with in situ data indicates that the new OE aerosol property estimates tend to be more realistic than previous estimates obtained from remotely sensed data through other retrieval techniques. Based on the results of this study we therefore suggest that the new Optimal Estimation retrieval algorithm is able to contribute to an advancement in aerosol research by considerably improving current estimates of aerosol properties in the lower stratosphere under low aerosol loading conditions.


Sign in / Sign up

Export Citation Format

Share Document